- Open Access
In silico prediction of Gallibacterium anatis pan-immunogens
Veterinary Research volume 45, Article number: 80 (2014)
The Gram-negative bacterium Gallibacterium anatis is a major cause of salpingitis and peritonitis in commercial egg-layers, leading to reduced egg production and increased mortality. Unfortunately, widespread multidrug resistance and antigenic diversity makes it difficult to control infections and novel prevention strategies are urgently needed. In this study, a pan-genomic reverse vaccinology (RV) approach was used to identify potential vaccine candidates. Firstly, the genomes of 10 selected Gallibacterium strains were analyzed and proteins selected on the following criteria; predicted surface-exposure or secretion, none or one transmembrane helix (TMH), and presence in six or more of the 10 genomes. In total, 42 proteins were selected. The genes encoding 27 of these proteins were successfully cloned in Escherichia coli and the proteins expressed and purified. To reduce the number of vaccine candidates for in vivo testing, each of the purified recombinant proteins was screened by ELISA for their ability to elicit a significant serological response with serum from chickens that had been infected with G. anatis. Additionally, an in silico prediction of the protective potential was carried out based on a protein property prediction method. Of the 27 proteins, two novel putative immunogens were identified; Gab_1309 and Gab_2312. Moreover, three previously characterized virulence factors; GtxA, FlfA and Gab_2156, were identified. Thus, by combining the pan-genomic RV approach with subsequent in vitro and in silico screening, we have narrowed down the pan-proteome of G. anatis to five vaccine candidates. Importantly, preliminary immunization trials indicated an in vivo protective potential of GtxA-N, FlfA and Gab_1309.
Poultry meat and eggs are considered very important and sustainable sources of animal protein worldwide , making efficient strategies to prevent and control the spread of poultry diseases highly important . Gallibacterium anatis is a Gram-negative, non-motile, encapsulated coccobacillus of the Pasteurellaceae family , and commonly associated with poultry . Besides constituting a part of the normal microflora of the upper respiratory tract and lower genital tract in chickens , it is also considered a major cause of salpingitis and peritonitis in egg-laying chickens -. Thus, G. anatis infections lead to a drop in egg production and increased mortality in commercial layers . Unfortunately, widespread multiple-drug resistance  hinders treatment with traditional antimicrobial agents, while substantial antigenic diversity  among disease-causing field isolates hampers disease prevention by classical vaccination with inactivated whole cell bacterins. Hence, novel prevention strategies are urgently needed.
The sequencing of the first bacterial genome in 1995  initiated the genomic era and catalyzed a shift from conventional culture-based approaches to genome-based vaccinology . This gave rise to the Reverse Vaccinology (RV) approach , in which bioinformatics tools are used to analyze genome sequences to identify genes encoding likely protective antigens. The concept of RV was initially applied to Neisseria meningitidis serogroup B (MenB) , for which conventional vaccine development approaches had failed in producing an efficacious vaccine. Based on the genomic sequence of MenB strain MC58 , five universal vaccine candidates were identified , and the resulting 4CMenB vaccine (Bexsero®) is now approved in the EU for active immunization of individuals aged over two months against disease caused by MenB . Since this pioneering MenB project, the RV approach has been applied to a variety of other important pathogens . However, the increased availability of multiple genomes for the same bacterial species has shown that genomic variability in bacteria is much more extensive than initially anticipated. Thus, analysis of the genome of a single strain often fails to address intra-species genetic variability and limits the effectiveness of genome-wide screens for vaccine candidates. To overcome this, a pan-genomic RV model utilizing the global gene repertoire for a species was proposed by Tettelin et al. . Pan-genomic RV was first applied to vaccine development in Group B Streptococcus, and this study demonstrated the importance of sequencing multiple strains of a single pathogen for the identification of vaccine antigens . The application of in silico and in vitro predictions has not only enabled a much more rational selection of vaccine candidates, but has also shown promise at reducing the number of experimental animals needed to verify the effectiveness of vaccine candidates.
We report here the use of a pan-genomic RV approach for identification of novel and conserved immunogens of G. anatis. By implementing different in silico approaches and in vitro assays, we screened the Gallibacterium pan-proteome, resulting in a final selection of five proteins with a high predicted potential as vaccine candidates. Importantly, preliminary in vitro immunization results indicate protective potential of at least three of these candidates including FlfA, which has previously been tested and confirmed highly protective against homologous challenge in chickens . Together, these results provide an important step in the development of a new and broadly protective vaccine against G. anatis.
Materials and methods
Animal ethics statement
All work on experimental animals was carried out with the approval of the Danish National Animal Ethics Committee (Approval no. 2012-15-2934-00339 and 2012-15-2934-00923).
Gallibacterium strains and growth conditions
The 10 Gallibacterium strains included in the study are listed in Table 1. The strains were selected based on their pathogenic potential, prevalence in the field and genetic diversity, in order to provide as much diversity as possible within the species. Bacteria were incubated at 37 °C on brain heart infusion (BHI) agar supplemented with 5% citrated bovine blood in a closed plastic bag or in BHI broth with aeration.
RV in silico analysis
Sequencing and assembly of G. anatis strains 12656–12 and F149T was performed as described previously . The remaining genomes were sequenced using the Illumina Genome Analyzer IIx (CD genomics, New York, USA). Reads were assembled using VelvetOptimiser 2.0 . All 10 genomes were annotated using Prokka v.1.0 . The subcellular localization of the 31 564 annotated proteins from the 10 genomes was predicted using standalone PSORTb v.3.0 . The presence of N-terminal signal peptides was predicted using SignalP v.3.0  and the number of transmembrane helices (TMHs) was predicted by TMHMM v.2.0 . The protein conservation among the strains was analyzed using BLASTp v.2.2.22  with default parameters. The theoretical molecular masses and isoelectric points were calculated using the pepstats tool in EMBOSS . In total, 42 proteins were selected (Table 2).
Cloning and small-scale protein expression
Each of the selected genes was amplified from the G. anatis 12656–12 genome by PCR and cloned into the Gateway entry vector pENTR™/SD/D-TOPO (Invitrogen). Primers were designed using Oligo Explorer 1.2 (Gene Link™, Hawthorne, NY, USA) as described previously . Areas with high predicted hydrophobicity in the N and/or C terminus were removed, as were predicted signal peptides. In addition, GtxA was cloned as two parts (N- and C-terminal) due to its size. The E. coli strains and plasmids used in this study are listed in Additional file 1. Genes were cloned and small-scale expressed as described in Additional file 2 using Gateway cloning and ligation-independent cloning (LIC) systems. Altogether, 37 expression clones were constructed for 36 of the 42 selected proteins (two clones were made for GtxA). The primer sequences used for gene amplification, and the final expression vector chosen for protein expression from each gene, are listed in Additional file 3.
Large-scale protein expression and purification
All proteins were expressed and purified in large-scale from E. coli Rosetta 2 (DE3) cells (Novagen, Madison, WI, USA). Large-scale expression was performed in a custom-made large-scale expression system (LEX) (Harbinger Biotech, Toronto, Canada) as described previously  and in Additional file 2. Of the 37 expression clones, 27 recombinant proteins were successfully purified; the majority (17) of these proteins had a purity > 90%.
Production of antiserum against G. anatis 12656–12 in chickens
Two Lohmann brown chickens (21 weeks old) were purchased from a commercial breeder with high biosecurity standards. The chickens were kept under free indoor housing conditions and were provided with fresh water and feed ad libitum. The chickens were swabbed for the presence of G. anatis by a cloacal swab. After two weeks of acclimatization the chickens were challenged with 105 colony forming units (CFU) of G. anatis 12656–12 by injection into the peritoneal cavity as previously described  and re-infected 2 weeks after the first infection. Blood for serum purification was collected from the brachial vein prior to the first infection (pre-immune antiserum) and one week after the second infection (hyper-immune antiserum).
Enzyme-Linked Immunosorbent Assay (ELISA)
The putative immunogenicity of each of the purified recombinant proteins was assessed by indirect ELISA as described previously , using pooled anti-G. anatis pre-immune and hyper-immune antiserum. Briefly, Nunc-Immuno™ MicroWell™ 96-Well Plates (Thermo Scientific, Waltham, MA, USA) were coated overnight at 4 °C with 0.5 μg recombinant protein (48 wells per protein) diluted in carbonate-bicarbonate buffer (pH 9.6) (Sigma-Aldrich, St. Louis, MO, USA). Each well was then washed; this and all subsequent washing steps consisted of three washes in 350 μL wash buffer (PBS + 0.05% Tween 20). The wells were blocked for 2 h at room temperature in 200 μL blocking solution (PBS containing 0.05% Tween 20 and 2% bovine serum albumin (BSA)) and washed. The antibody titers were assayed by serial 3-fold dilutions of chicken serum ranging from 1:200 to 1:48600. All dilutions were prepared in triplicate in dilution buffer (PBS containing 0.05% Tween 20 and 0.1% BSA), 100 μL were added to each well and plates were incubated for 1 h at 37 °C. For each assay, 12 control wells were included, which contained pure dilution buffer; secondary antibody was added to 6 of these wells as a measure of background, and the other 6 wells remained blank as a negative control for the ELISA. Following incubation, the wells were washed and 100 μL polyclonal goat anti-chicken IgG (Fc):HRP (AbD Serotec, Puchheim, Germany), diluted 1:4000 in diluting buffer, were added to each well and the plates incubated for a further 1 h at 37 °C and then washed. To detect the binding, 100 μL of 3,3′,5,5′-Tetramethylbenzidine (TMB) liquid substrate (Sigma) were added to each well. The plates were incubated for 2 min and then the reaction was stopped by addition of 100 μL 1 M HCl. The absorbance was read immediately at 450 nm in a PowerWave XS spectrophotometry (BioTek Instruments, Winooski, VT, USA).
The antibody titers were calculated for the measured absorbances at 450 nm , using the “Antibody Titers” online data analysis tool . To compare and rank the ELISA results, a P/N ratio (P = hyper-immune serum, N = pre-immune serum) of mean antibody titers was calculated . All statistical analysis was performed using SAS version 9.3 (SAS Institute, Cary, NC, USA), and differences between groups assessed using a one-way t-test. The recognition of recombinant protein by hyper-immune serum was deemed significant at P < 0.05, indicating that the protein was expressed in vivo during G. anatis infection and elicited a specific immune response.
VacFinder® in silico protein analysis
To further predict the protective potential of each of the expressed proteins, each of the proteins was analyzed using the proprietary VacFinder® in silico technology platform (Evaxion Biotech, LLC, USA). VacFinder® is a data-driven machine learning method trained by protein property pattern recognition on known and protective B-cell protein antigens (except for exotoxins), aiming at identifying novel and protective B-cell protein antigens with a neutralizing opsonizing profile. The machine-learned prediction is based on specific protein property features of protein sequences rather than sequence similarity, allowing antigen classification based solely on protein properties . The output is a list of proteins from the proteome ranked by their ability to elicit a highly protective antibody response.
Immunization of layer chickens with recombinant proteins
24 Isa Brown layer chickens (16 weeks old) were purchased from a commercial breeder with high biosecurity standards. The chickens were swabbed for the presence of G. anatis by a cloacal swab. The chickens were randomly divided in eight groups of four each and allowed to acclimatize for one week after arrival. The chickens were kept under free indoor housing conditions and provided with fresh water and feed ad libitum. Each group was immunized subcutaneously with 100 μg of one of the five selected recombinant proteins (GtxA-N, FlfA, Gab_1309, Gab_2156 or Gab_2312) mixed in 0.5 mL of SEC buffer (50 mM NaP, 150 mM NaCl, 0.5 mM TCEP, 10% glycerol; pH 7.5) and 0.5 mL of Freund’s incomplete adjuvant (Sigma-Aldrich). As a control (non-immunized), a group of four chickens was immunized with a placebo (SEC buffer and Freund’s incomplete adjuvant). Two weeks after the immunization all chickens were infected intraperitoneally with 1.5 × 106 CFU of G. anatis 12656–12 as described previously . Forty-eight hours after infection the chickens were euthanized and a post mortem examination was conducted. To assess the protective effect of the immunization, the lesions found in peritoneum of each bird were scored according to three parameters: (i) inflammatory reaction, (ii) amount of exudate, and (iii) clarity of the peritoneal surfaces. All parameters were scored on a scale from 0–3, thus giving a maximum score of 9. Furthermore, the presence of Gallibacterium was detected by swabbing the peritoneum with a sterile cotton swab and streaking it onto BHI blood agar. The scorings of the lesions in the peritoneum were analyzed by a Mann Whitney U test and P < 0.05 were deemed significant.
Multiple sequence alignments
Multiple amino acid sequence alignments of the Gab_1309 and Gab_2312 proteins and their orthologs were prepared using MAFFT v7.130b  and formatted using Jalview 2.8.0b1 .
Genbank accession numbers
The genome sequence of G. anatis 12656–12 has recently been made available  and was submitted to Genbank (BioProject ID: 213810, accession number AVOX00000000). The nucleotide sequence accession numbers for the genes included in this study are listed in Table 2. Genome sequence reads from the nine Gallibacterium strains used for Gab 1309 and Gab_2312 multiple sequence alignment were submitted to the NCBI Sequence Read Archive (SRA)  and can be retrieved using the study accession number SRP029613.
RV in silico prediction of candidate vaccine antigens
For the identification and selection of putative immunogens, the genomes of 10 Gallibacterium strains were analyzed (Table 1). Based on the central premise that protective antigens should be accessible to the host immune system, proteins predicted to be surface-exposed or secreted were selected from the G. anatis pan-proteome. Moreover, proteins with more than one TMH were discarded, based on the premise that they are unlikely to be transported beyond the inner membrane. In addition, these proteins have the highest rate of expression failure during subsequent procedures  or are less likely to be over-expressed in E. coli. Finally, proteins present in six or more of the 10 G. anatis genomes were identified. A protein was considered present if a significant full length match (E-value < 10−8) was obtained. From a total of 31564 proteins, 162 proteins were predicted as extracellular proteins and 482 proteins as outer membrane proteins. Of these, 42 proteins were present in six or more of the 10 genomes and these proteins were selected for further studies (Table 2).
Cloning, expression and purification
Each of the 42 selected genes was amplified from the G. anatis strain 12656–12 genome by PCR. This strain was chosen for gene cloning as it is a well-characterized and highly pathogenic strain originally isolated from the liver of a chicken with septicaemia ,. Using different high-throughput cloning strategies, 37 proteins were successfully purified at a small-scale from E. coli and of these, 27 were successfully purified in the large-scale trials (Table 3). The 10 remaining proteins were lost during the large-scale purification process due to low yield, impurities or lack of expression in E. coli in large-scale cultures.
Screening of immunogenic potential by ELISA
In order to reduce the number of antigens to be tested in vaccine trials, the immunogenic potential of the 27 purified proteins was assessed in vitro by an indirect ELISA approach using pools of anti-G. anatis 12656–12 pre-immune and hyper-immune antiserum from two chickens. The chickens were, as expected, found to be positive for G. anatis on the cloacal mucosa prior to the immunizations. The mean antibody titers were calculated for each protein using both pre-immune and hyper-immune antiserum. To rank the proteins, P/N values between the two antibody titers were calculated for each protein (Table 3). A significant and specific reaction with hyper-immune antiserum was identified for 12 of the 27 recombinant proteins tested (Figure 1), namely GtxA (both halves of the protein), FlfA, Gab_0091, Gab_0523, Gab_1309, Gab_1575, Gab_1654, Gab_1755, Gab_2156, Gab_2312 and Gab_2348. These results indicate (i) that these 12 proteins are expressed in vivo by G. anatis 12656–12, (ii) that the proteins are recognized by the chicken immune system, and (iii) that the proteins elicit a specific antibody response. The remaining proteins were not significantly recognized by the hyper-immune antiserum.
Screening of protective potential by VacFinder®
To provide an additional prediction of the protective potential of the 27 proteins available in recombinant form, the protein sequences were analyzed in silico with the VacFinder® algorithm. The entire G. anatis strain 12656–12 proteome (approximately 2500 proteins) was analyzed and ranked for protective potential. Proteins ranked in the Top20 are those predicted to have the highest protective potential (Table 3). Of the 27 proteins cloned and expressed in E. coli, only Gab_1309 and Gab_2312 were ranked within the Top20 VacFinder® hits.
Conservation of Gab_1309 and Gab_2312
Two of the 27 proteins, Gab_1309 and Gab_2312, elicited a significant ELISA response and were ranked in the Top20 proteins with protective potential in the G. anatis 12656– 12 proteome as determined by VacFinder®, suggesting that these proteins should be prioritized in further analyses. Both proteins were expressed in E. coli with an N-terminus His6-tag and purified in soluble form. Due to the high degree of predicted hydrophobicity within the C-terminal region of Gab_2312, as well as the size of the whole protein (325.7 kDa), only a small part (338 residues) was included in the final recombinant protein. Homologs of Gab_1309 and Gab_2312 were identified within the genomes of all the 10 strains included in this study, as well as in the recently sequenced G. anatis strain UMN179 (UMN_1211 and UMN_1565) , further supporting their potential as broadly protective immunogens. To determine the level of protein sequence conservation of Gab_1309 and Gab_2312, multiple alignments were conducted against the protein homologs from other G. anatis strains. The Gab_1309 protein sequences were almost identical across all strains (Additional file 4). On the other hand, the Gab_2312 protein sequences varied considerably in length from 2157 residues in G. anatis strain F149 to 5202 residues in G. genomospecies 2 strain CCM5976. Despite this difference, several well-conserved sections could be identified within the sequences, indicating a degree of conservation within the structure and the presence of common epitopes. Additional file 5 shows the alignment between the 338 residues included in the recombinant protein and the corresponding parts in the protein homologs.
Potential protective capacity
To evaluate the protective potential in vivo of the five most promising vaccine candidates (GtxA-N, FlfA, Gab_1309, Gab_2156, or Gab_2312), groups of four chickens were each immunized with one of the recombinant proteins, followed by intraperitoneal challenge with G. anatis 12656–12. G. anatis was recovered in pure culture from chickens within all the immunized groups, as well as from the non-immunized group. However, two chickens from the group immunized with GtxA-N and one chicken from the group immunized with FlfA were culture-negative. A significant lower lesion score was found in the group immunized with GtxA-N (P = 0.02), FlfA (P = 0.04) and Gab_1309 (P = 0.02) when compared to the non-immunized group (Figure 2), indicating a protective potential of these three proteins. The two groups immunized with Gab_2156 or Gab_2312 did not show a significant difference when compared with the non-immunized group.
Novel prevention and treatment strategies are urgently needed to prevent G. anatis infections in the reproductive tract of chickens. In this study, a pan-genomic RV approach , was applied to identify novel and potentially broadly protective immunogens from G. anatis. The screening procedures and the main results are summarized in Figure 3. Of the 42 in silico predicted immunogens (Table 2), 27 proteins were successfully expressed in E. coli (Table 1), and of these, two novel proteins, Gab_1309 and Gab_2312, elicited a significant ELISA response and were also ranked in the Top20 of the G. anatis 12656– 12 proteome by VacFinder®. Furthermore, these two proteins were present in all 10 strains included in this study as well as in the recently sequenced genome from G. anatis strain UMN179 .
Gab_1309 is a predicted lipoprotein, and this annotation is supported by the DOLOP database . The protein sequence of Gab_1309 shows sequence similarity to the NlpD lipoprotein from Yersinia pestis. By use of the Conserved Domain Database (CDD) feature in BLASTp , an N-terminal LysM domain and a C-terminal M23 peptidase domain can be identified in Gab_1309. The LysM domain is widespread among several bacterial species and is involved in peptidoglycan-binding and bacterial cell wall degradation ,, while the M23 family of endopeptidases is thought to be involved in bacterial cell wall separation . Further studies are warranted to determine the specific function and role of Gab_1309 in G. anatis pathogenesis.
The C-terminal section of the second promising novel vaccine candidate identified in this study, Gab_2312, demonstrated 34% sequence similarity (conserved or identical residues) to the 200 kDa extracellular matrix protein adhesin A (EmaA) from Aggregatibacter actinomycetemcomitans. EmaA is widespread among A. actinomycetemcomitans strains  and belongs to the family of trimeric autotransporter adhesins (TAAs) . TAAs are a group of homotrimeric virulence-related proteins in Gram-negative bacteria that primarily act as adhesins . Interestingly, the Neisseria adhesin A (NadA) protein from MenB, which is one of the proteins in the recently approved 4CMenB vaccine (Bexsero®), also belongs to this family of TAAs , further suggesting a protective potential of Gab_2312. The prediction of Gab_2312 as a TAA protein was further confirmed using the recently developed domain annotation workflow of TAAs (daTAAs) .
Although Gab_1309 and Gab_2312 were the only proteins identified by both the ELISA screening and the in silico VacFinder® analysis, a number of proteins elicited a significant response as measured by ELISA, indicating that these proteins are expressed and immunogenic during G. anatis infection. It is currently unclear whether proteins that do not stimulate a natural immune response should be included as vaccine candidates. As our intraperitoneal immunization procedure bypassed the mucosal surface of e.g. the salpinx, this may not have induced a broad based natural immunological response, which in turn could explain why some proteins failed to induce a significant antibody response. However, studies have shown that high antigenicity of a protein does not necessarily lead to protection , questioning the validity of the antigen-specific titer as a marker of the overall ability to stimulate a protective immune response in vivo . On the other hand, serological titers have been demonstrated to provide useful information about the reactivity of certain antigens to antibodies and thereby give indications about their potential as vaccine candidates . In support of the ELISA antigen screening approach, three of the 10 proteins identified by the ELISA screening; GtxA, FlfA and Gab_2156, have recently been demonstrated as important, well-conserved G. anatis virulence-factors. GtxA is a large, cytolytic RTX toxin responsible for both the haemolytic and leukotoxic activity of G. anatis,, while FlfA and Gab_2156 are both subunits of F17-like fimbriae . To this, FlfA has been suggested to play a role in the tissue tropism of G. anatis by promoting adhesion to the epithelial lining in the reproductive tract and beyond during the typical course of disease . The identification of GtxA, FlfA and Gab_2156 by ELISA further supports the predictive potential and importance of this screening approach, where particularly the N-terminal part of the GtxA-N elicited an extremely high antibody response in comparison to the remaining 26 proteins investigated (Figure 1 and Table 3), thereby suggesting that this part of the protein is particularly promising as a vaccine candidate. This suggestion is supported by a previous study by , demonstrating the protective potential of the N-terminal portion of the RTX protein, ApxI, against infections caused by Actinobacillus pleuropneumoniae.
The use of hyper-immune antiserum to screen for antigenic potential has previously been used successfully to identify putative immunogens from Bacillus anthracis, and in that study, the control sera (pre-immune antiserum) consisted of a mixture of naïve animal sera . In our study, the pre-immune antiserum and hyper-immune serum was obtained from the same chickens before and after G. anatis 12656–12 infection. However, the chickens were found to be positive for G. anatis in the trachea and cloaca prior to infection. Since G. anatis is a typical part of the microflora in these organ systems , we find that this reflects the natural in vivo condition, which allows us to identify those proteins significantly expressed and exposed to the chicken immune system during infection. Thus, the ELISA screening applied in this study can be considered a semi-in vivo-evaluation and an important screening assay to evaluate antigenic potentials of recombinant proteins.
Of the 27 proteins tested by ELISA, 15 were not recognized by hyper-immune antiserum. However, these proteins might still possess immunogenic potential, as it is well-established that in vitro assessments may miss some proteins . Moreover, a recent study focused interest on non-immunodominant protein regions for vaccine development and suggested the importance of analyzing major surface-exposed proteins for the presence of sub-regions that elicit protective immunity as a complement to the RV approach . Additionally, the insolubility of many of the proteins tested by ELISA might also lead to false-negative titer-values, as it has been shown that linear epitopes are less likely to be immunodominant and elicit a good B-cell response if the protein possesses a higher-order structure . On the other hand, the folding of a recombinant protein is not necessarily equivalent to the native structure of the protein, and moreover, folded recombinant proteins might expose epitopes normally hidden , and thereby lead to a false-positive ELISA response.
To provide an additional prediction of the protective potential, the sequences of 27 proteins were ranked within the G. anatis 12656–12 proteome by use of the in silico VacFinder® platform, in order to add an extra selective criterion. The ranking by VacFinder® is based on the predicted ability of a protein to elicit a high protective antibody response. Thus, the identified proteins are not affected by possible solubility problems or competing immunodominant epitopes. Moreover, VacFinder® does not consider general practicalities of cloning, expression and purification. The combination of VacFinder® prediction and ELISA screening identified two proteins, namely the above-mentioned Gab_1309 and Gab_2312. However, the previously identified and well-conserved virulence factors, GtxA, FlfA and Gab_2156, were not predicted as highly protective by VacFinder®. Leaving out GtxA, which is an exotoxin and not possible for VacFinder® to predict, FlfA and Gab_2156 were respectively ranked #58 and #41 out of the whole proteome. Hence, the Top20 of the proteome by VacFinder® is predicted to elicit the highest protective antibody response, yet this study demonstrates that proteins with a lower predicted protectiveness can be highly immunogenic.
Finally, the protective potential of the five most promising vaccine candidates; GtxA-N, FlfA, Gab_1309, Gab_2156 and Gab_2312, was evaluated in vivo. The results from these preliminary in vivo immunization trials indicated that at least GtxA-N, FlfA and Gab_1309 are promising vaccine candidates with a good protective potential. These results correspond well with previous results demonstrating the in vivo protective potential of FlfA . Further studies, including larger group sizes, different doses and repeated immunization, are needed to confirm the protective potential of the recombinant proteins. Moreover, the ability of the recombinant proteins to elicit a broadly cross-protective immune response against heterologous strains should also be investigated. Still, the results presented in this paper provide an important step towards the development of a new and broadly protective G. anatis vaccine.
- G. anatis 12656–12:
Gallibacterium anatis biovar haemolytica strain 12656–12
Neisseria meningitides seropgroup B
Enzyme-linked immunosorbent assay
Trimeric autotransporter adhesin
AVEC: Annual report 2011. Assoc Poult Processors Poult Trade EU Ctries. 2011, 41: 1-52.
Marangon S, Busani L: The use of vaccination in poultry production. Rev Sci Tech. 2007, 26: 265-274.
Christensen H, Bisgaard M, Bojesen AM, Mutters R, Olsen JE: Genetic relationships among avian isolates classified asPasteurella haemolytica, “Actinobacillus salpingitidis” orPasteurella anatiswith proposal ofGallibacterium anatisgen. nov., comb. nov. and description of additional genomospecies withinGallibacteriumgen. nov.Int J Syst Evol Microbiol 2003, 53:275–287.,
Bisgaard M, Korczak BM, Busse H-J, Kuhnert P, Bojesen AM, Christensen H: Classification of the taxon 2 and taxon 3 complex of Bisgaard withinGallibacteriumand description ofGallibacterium melopsittacisp. nov.,Gallibacterium trehalosifermentanssp. nov. andGallibacterium salpingitidissp. nov.Int J Syst Evol Microbiol 2009, 59:735–744.,
Bisgaard M: Incidence ofPasteurella haemolyticain the respiratory tract of apparently healthy chickens and chickens with infectious bronchitis: characterisation of 213 strains.Avian Pathol 1977, 6:285–292.,
Bojesen AM, Nielsen SS, Bisgaard M: Prevalence and transmission of haemolyticGallibacteriumspecies in chicken production systems with different biosecurity levels.Avian Pathol 2003, 5:503–510.,
Mirle C, Schöngarth M, Meinhart H, Olm U: Untersuchungen zu auftreten and bedeutung vonPasteurella haemolyticainfektionen bei hennen unter besonderer berücksichtigung von erkrankungen des legeapparates.Mh Vet Med 1991, 46:545–549 (in German).,
Jordan FTW, Williams NJ, Wattret A, Jones T: Observations on salpingitis, peritonitis and salpingoperitonitis in a layer breeder flock. Vet Rec. 2005, 157: 573-577.
Neubauer C, De Souza-Pilz M, Bojesen AM, Bisgaard M, Hess M: Tissue distribution of haemolyticGallibacterium anatisisolates in laying birds with reproductive disorders.Avian Pathol 2009, 38:1–7.,
Bojesen AM, Christensen JP, Bisgaard M:Gallibacterium infections and other avian Pasteurellaceae. Poultry diseases. Edited by: Pattison M, Mcmullin PF, Bradbury JM, Alexander DJ. 2008, Saunders Elsevier, Edinburg, 160-163. 6,
Bojesen AM, Vazquez ME, Bager RJ, Ifrah D, Gonzalez C, Aarestrup FM: Antimicrobial susceptibility and tetracycline resistance determinant genotyping ofGallibacterium anatis.Vet Microbiol 2011, 148:105–110.,
Schumacher T, Vazquez ME, Leon-Kempis MR: Analysis of strains of Gallibacterium anatis by agar gel precipitation and determination of cross reactivity [Poster]. 2008, International Pasteurellaceae Society Meeting, Sorrento, Italy
Fleischmann RD, Adams MD, White O, Clayton RA, Kirkness EF, Kerlavage AR, Bult CJ, Tomb JF, Dougherty BA, Merrick JM: Whole-genome random sequencing and assembly ofHaemophilus influenzaeRd.Science 1995, 269:496–512.,
Rinaudo CD, Telford JL, Rappuoli R, Seib KL: Vaccinology in the genome era. J Clin Invest. 2009, 119: 2515-2525.
Rappuoli R: Reverse vaccinology, a genome-based approach to vaccine development. Vaccine. 2001, 19: 2688-2691.
Pizza M, Scarlato V, Masignani V, Giuliani MM, Aricò B, Comanducci M, Jennings GT, Baldi L, Bartolini E, Capecchi B, Galeotti CL, Luzzi E, Manetti R, Marchetti E, Mora M, Nuti S, Ratti G, Santini L, Savino S, Scarselli M, Storni E, Zuo P, Broeker M, Hundt E, Knapp B, Blair E, Mason T, Tettelin H, Hood DW, Jeffries AC: Identification of vaccine candidates against serogroup B meningococcus by whole-genome sequencing. Science. 2000, 287: 1816-1820.
Tettelin H, Saunders NJ, Heidelberg J, Jeffries AC, Nelson KE, Eisen JA, Ketchum KA, Hood DW, Peden JF, Dodson RJ, Nelson WC, Gwinn ML, DeBoy R, Peterson JD, Hickey EK, Haft DH, Salzberg SL, White O, Fleischmann RD, Dougherty BA, Mason T, Ciecko A, Parksey DS, Blair E, Cittone H, Clark EB, Cotton MD, Utterback TR, Khouri H, Qin H, et al: Complete genome sequence ofNeisseria meningitidisserogroup B strain MC58.Science 2000, 287:1809–1815.,
Giuliani MM, Adu-Bobie J, Comanducci M, Aricò B, Savino S, Santini L, Brunelli B, Bambini S, Biolchi A, Capecchi B, Cartocci E, Ciucchi L, Di Marcello F, Ferlicca F, Galli B, Luzzi E, Masignani V, Serruto D, Veggi D, Contorni M, Morandi M, Bartalesi A, Cinotti V, Mannucci D, Titta F, Ovidi E, Welsch JA, Granoff D, Rappuoli R, Pizza M: A universal vaccine for serogroup B meningococcus. Proc Natl Acad Sci U S A. 2006, 103: 10834-10839.
Carter N: Multicomponent meningococcal serogroup B vaccine (4CMenB; Bexsero®): a review of its use in primary and booster vaccination. BioDrugs. 2013, 3: 263-274.
Donati C, Rappuoli R: Reverse vaccinology in the 21st century: improvements over the original design. Ann N Y Acad Sci. 2013, 1285: 115-132.
Tettelin H, Masignani V, Cieslewicz MJ, Donati C, Medini D, Ward NL, Angiuoli SV, Crabtree J, Jones AL, Durkin AS, Deboy RT, Davidsen TM, Mora M, Scarselli M, Margarit y Ros I, Peterson JD, Hauser CR, Sundaram JP, Nelson WC, Madupu R, Brinkac LM, Dodson RJ, Rosovitz MJ, Sullivan SA, Daugherty SC, Haft DH, Selengut J, Gwinn ML, Zhou L, Zafar N, et al: Genome analysis of multiple pathogenic isolates ofStreptococcus agalactiae: implications for the microbial “pan-genome”.Proc Natl Acad Sci U S A 2005, 102:13950–13955.,
Maione D, Margarit I, Rinaudo CD, Masignani V, Mora M, Scarselli M, Tettelin H, Brettoni C, Iacobini ET, Rosini R, D’Agostino N, Miorin L, Buccato S, Mariani M, Galli G, Nogarotto R, Nardi-Dei V, Vegni F, Fraser C, Mancuso G, Teti G, Madoff LC, Paoletti LC, Rappuoli R, Kasper DL, Telford JL, Grandi G: Identification of a universal Group Bstreptococcusvaccine by multiple genome screen.Science 2005, 309:148–150.,
Medini D, Donati C, Tettelin H, Masignani V, Rappuoli R: The microbial pan-genome. Curr Opin Genet Dev. 2005, 15: 589-594.
Bager RJ, Nesta B, Pors SE, Soriani M, Serino L, Boyce JD, Adler B, Bojesen AM: The fimbrial protein FlfA fromGallibacterium anatisis a virulence factor and vaccine candidate.Infect Immun 2013, 81:1964–1973.,
Bojesen AM, Torpdahl M, Christensen H, Olsen JE, Bisgaard M: Genetic diversity ofGallibacterium anatisisolates from different chicken flocks.J Clin Microbiol 2003, 41:2737–2740.,
Kristensen BM, Frees D, Bojesen AM: Expression and secretion of the RTX-toxin GtxA among members of the genusGallibacterium.Vet Microbiol 2011, 153:116–123.,
Johnson TJ, Danzeisen JL, Trampel D, Nolan LK, Seemann T, Bager RJ, Bojesen AM: Genome analysis and phylogenetic relatedness ofGallibacterium anatisstrains from poultry.PLoS One 2013, 1:e54844.,
The VelvetOptimiser 2.0de novosequence assembler. In [http://www.vicbioinformatics.com/software.velvetoptimiser.shtml]
Seemann T: Prokka: rapid prokaryotic genome annotation. Bioinformatics. 2014, 30: 2068-2069.
Yu NY, Wagner JR, Laird MR, Melli G, Rey S, Lo R, Dao P, Sahinalp SC, Ester M, Foster LJ, Brinkman FS: PSORTb 3.0: improved protein subcellular localization prediction with refined localization subcategories and predictive capabilities for all prokaryotes. Bioinformatics. 2010, 13: 1608-1615.
Bendtsen JD, Nielsen H, von Heijne G, Brunak S: Improved prediction of signal peptides: SignalP 3.0. J Mol Biol. 2004, 4: 783-795.
Krogh A, Larsson B, von Heijne G, Sonnhammer EL: Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J Mol Biol. 2001, 3: 567-580.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410.
Rice P, Longden I, Bleasby A: EMBOSS: the European molecular biology open software suite. Trends Genet. 2000, 16: 276-277.
Kristensen BM, Frees D, Bojesen AM: GtxA fromGallibacterium anatis, a cytolytic RTX-toxin with a novel domain organisation.Vet Res 2010, 41:25.,
Bager RJ, Persson G, Nesta B, Soriani M, Serino L, Jeppsson M, Nielsen TK, Bojesen AM: Outer membrane vesicles reflect environmental cues inGallibacterium anatis.Vet Microbiol 2013, 167:565–572.,
Al-Hasani K, Boyce J, McCarl VP, Bottomley S, Wilkie I, Adler B: Identification of novel immunogens inPasteurella multocida.Microb Cell Fact 2007, 6:3.,
Gräslund S, Sagemark J, Berglund H, Dahlgren L-G, Flores A, Hammarström M, Johansson I, Kotenyova T, Nilsson M, Nordlund P, Weigelt J: The use of systematic N- and C-terminal deletions to promote production and structural studies of recombinant proteins. Protein Expr Purif. 2008, 2: 210-221.
Bojesen AM, Nielsen OL, Christensen JP, Bisgaard M: In vivostudies ofGallibacterium anatisinfection in chickens.Avian Pathol 2004, 33:145–152.,
Raida MK, Nylén J, Holten-Andersen L, Buchmann K: Association between plasma antibody response and protection in rainbow troutOncorhynchus mykissimmersion vaccinated againstYersinia ruckeri.PLoS One 2011, 6:e18832.,
Villumsen KR, Dalsgaard I, Holten-Andersen L, Raida MK: Potential role of specific antibodies as important vaccine induced protective mechanism againstAeromonas salmonicidain rainbow trout.PLoS One 2012, 7:e46733.,
The MyAssays® Antibody Titer Calculation Software.In [http://www.myassays.com/antibody-titers.assay]
Case JT, Ardans AA, Bolton DC, Reynolds BJ: Optimization of parameters for detecting antibodies against infectious bronchitis virus using an enzyme-linked immunosorbent assay: temporal response to vaccination and challenge with live virus. Avian Dis. 1983, 27: 196-210.
Doytchinova IA, Flower DR: VaxiJen: a server for prediction of protective antigens, tumor antigens and subunit vaccines. BMC Bioinformatics. 2007, 8: 4-
Katoh K, Standley DM: MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013, 30: 772-780.
Waterhouse AM, Procter JB, Martin DMA, Clamp M, Barton GJ: Jalview Version 2 - a multiple sequence alignment editor and analysis workbench. Bioinformatics. 2009, 25: 1189-1191.
Kudirkiene E, Christensen H, Bojesen AM: Draft genome sequence ofGallibacterium anatisbiovarhaemolytica12656–12 Liver, an isolate obtained from the liver of a septicaemic chicken.Genome Announc 2013, 1:e00810-13.,
The NCBI Sequence Read Archive (SRA). In [https://www.ncbi.nlm.nih.gov/sra]
Korepanova A, Gao FP, Hua Y, Qin H, Nakamoto RK, Cross TA: Cloning and expression of multiple integral membrane proteins fromMycobacterium tuberculosisinEscherichia coli.Protein Sci 2005, 14:148–158.,
Babu MM, Sankaran K: DOLOP–database of bacterial lipoproteins. Bioinformatics. 2002, 18: 641-643.
Tidhar A, Flashner Y, Cohen S, Levi Y, Zauberman A, Gur D, Aftalion M, Elhanany E, Zvi A, Shafferman A, Mamroud E: The NlpD lipoprotein is a novelYersinia pestisvirulence factor essential for the development of plague.PLoS One 2009, 4:e7023.,
Marchler-Bauer A, Lu S, Anderson JB, Chitsaz F, Derbyshire MK, DeWeese-Scott C, Fong JH, Geer LY, Geer RC, Gonzales NR, Gwadz M, Hurwitz DI, Jackson JD, Ke Z, Lanczycki CJ, Lu F, Marchler GH, Mullokandov M, Omelchenko MV, Robertson CL, Song JS, Thanki N, Yamashita RA, Zhang D, Zhang N, Zheng C, Bryant SH: CDD: a Conserved Domain Database for the functional annotation of proteins. Nucleic Acids Res. 2011, 39: D225-D229.
Bateman A, Bycroft M: The structure of a LysM domain fromE. colimembrane-bound lytic murein transglycosylase D (MltD).J Mol Biol 2000, 299:1113–1139.,
Buist G, Steen A, Kok J, Kuipers OP: LysM, a widely distributed protein motif for binding to (peptido)glycans. Mol Microbiol. 2008, 68: 838-847.
Bernhardt TG, de Boer PA: Screening for synthetic lethal mutants inEscherichia coliand identification of EnvC (YibP) as a periplasmic septal ring factor with murein hydrolase activity.Mol Microbiol 2004, 52:1255–1269.,
Mintz KP: Identification of an extracellular matrix protein adhesin, EmaA, which mediates the adhesion ofActinobacillus actinomycetemcomitansto collagen.Microbiology 2004, 150:2677–2688.,
Linke D, Riess T, Autenrieth IB, Lupas A, Kempf VAJ: Trimeric autotransporter adhesins: variable structure, common function. Trends Microbiol. 2006, 14: 264-270.
Łyskowski A, Leo JC, Goldman A: Structure and biology of trimeric autotransporter adhesins. Adv Exp Med Biol. 2011, 715: 143-158.
Serruto D, Bottomley MJ, Ram S, Giuliani MM, Rappuoli R: The new multicomponent vaccine against meningococcal serogroup B, 4CMenB: immunological, functional and structural characterization of the antigens. Vaccines. 2012, 30 (Suppl 2): B87-B97.
Szczesny P, Lupas A: Domain annotation of trimeric autotransporter adhesins–daTAA. Bioinformatics. 2008, 24: 1251-1256.
Ooka H, Terkawi MA, Goo Y-K, Luo Y, Li Y: Babesia microti: molecular and antigenic characterizations of a novel 94-kDa protein (BmP94).Exp Parasitol 2011, 127:287–293.,
Chao C-C, Huber ES, Porter TB, Zhang Z, Ching W-M: Analysis of the cross-reactivity of various 56 kDa recombinant protein antigens with serum samples collected afterOrientia tsutsugamushiinfection by ELISA.Am J Trop Med Hyg 2011, 84:967–972.,
Seah JN, Frey J, Kwang J: The N-terminal domain of RTX toxin ApxI ofActinobacillus pleuropneumoniaeelicits protective immunity in mice.Infect Immun 2002, 70:6464–6467.,
Ariel N, Zvi A, Grosfeld H, Gat O, Inbar Y, Velan B, Cohen S, Shafferman A: Search for potential vaccine candidate open reading frames in theBacillus anthracisvirulence plasmid pXO1:in silicoandin vitroscreening.Infect Immun 2002, 70:6817–6827.,
Stålhammar-Carlemalm M, Waldemarsson J, Johnsson E, Areschoug T, Lindahl G: Nonimmunodominant regions are effective as building blocks in a streptococcal fusion protein vaccine. Cell Host Microbe. 2007, 2: 427-434.
Ito H, Nakashima T, So T, Hirata M, Inoue M: Immunodominance of conformation-dependent B-cell epitopes of protein antigens. Biochem Biophys Res Commun. 2003, 308: 770-776.
Keyburn AL, Boyce JD, Vaz P, Bannam TL, Ford ME, Parker D, Di Rubbo A, Rood JI, Moore RJ: NetB, a new toxin that is associated with avian necrotic enteritis caused by Clostridium perfringens. PLoS Pathog. 2008, 2: e26-
Gileadi O, Burgess-Brown NA, Colebrook SM, Berridge G, Savitsky P, Smee CE, Loppnau P, Johansson C, Salah E, Pantic NH: High throughput production of recombinant human proteins for crystallography. Methods Mol Biol. 2008, 426: 221-246.
Vernet E, Kotzsch A, Voldborg B, Sundström M: Screening of genetic parameters for soluble protein expression inEscherichia coli.Protein Expr Purif 2011, 77:104–111.,
Gasteiger E, Hoogland C, Gattiker A, Duvaud S, Wilkins MR, Appel RD, Bairoch A: Protein identification and analysis tools in the ExPASy server. The Proteomics Protocols Handbook. Edited by: Walker JM. 2005, Humana Press, Totowa, NJ, 571-607.
This work was supported by the Danish Research Council for Technology and Production, grant 09–065909, and a PhD stipend from the Faculty of Life Sciences, Copenhagen University. The authors acknowledge the excellent technical assistance of Kate Rainczuk, for her help with the Gateway cloning, and Sara Bjørn, for her help with the LIC cloning.
The authors declare that they have no competing interests.
RJB contributed to the design of the study, carried out the experiments and drafted the manuscript. EK and TS carried out the genome assembly, annotation and submission to Genbank. IP prepared the pan-proteome and carried out the in silico selection of the putative immunogens. JB and BA participated in the small-scale cloning and expression of the recombinant proteins. TKN participated in the large-scale expression and purification of the recombinant proteins. SEP carried out the in vivo experiments. AHM carried out the VacFinder predictions. AMB conceived the study, participated in its design and its coordination, and helped to draft the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: E. coli strains and plasmids. E. coli strains and plasmid used in this study to clone and express recombinant proteins . (DOCX 12 KB)
Additional file 3: Primer sequences, PCR products and expression vectors. A table summarizing the primer sequences used for gene amplification, the size of the resulting PCR product and the final expression vector chosen for protein expression from each of the G. anatis 12656– 12 genes that were successfully small-scaled cloned and expressed ,,. (DOCX 19 KB)
Additional file 5: Multiple sequence alignment of recombinant Gab_2312. A multiple alignment between the 338-residue recombinant Gab_2312 protein sequence and Gab_2312 protein homologs was conducted using MAFFT (v7.130b)  and formatted using Jalview 2.8.0b1 . Amino acids were colored in blue based on their conservation (dark blue = fully conserved). The two head domains included in the recombinant protein, which were predicted by the domain annotation workflow of TAAs (daTAAs) , are marked with red and green boxes, respectively. (PDF 134 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.
The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.
To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.
The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Bager, R.J., Kudirkiene, E., da Piedade, I. et al. In silico prediction of Gallibacterium anatis pan-immunogens. Vet Res 45, 80 (2014). https://doi.org/10.1186/s13567-014-0080-0
- Vaccine Candidate
- Protective Potential
- Cloacal Swab
- LysM Domain
- Reverse Vaccinology