Comparative Reverse Vaccinology of Piscirickettsia salmonis, Aeromonas salmonicida, Yersinia ruckeri, Vibrio anguillarum and Moritella viscosa, Frequent Pathogens of Atlantic Salmon and Lumpfish Aquaculture

Marine finfish aquaculture is affected by diverse infectious diseases, and they commonly occur as co-infection. Some of the most frequent and prevalent Gram-negative bacterial pathogens of the finfish aquaculture include Piscirickettsia salmonis, Aeromonas salmonicida, Yersinia ruckeri, Vibrio anguillarum and Moritella viscosa. To prevent co-infections in aquaculture, polyvalent or universal vaccines would be ideal. Commercial polyvalent vaccines against some of these pathogens are based on whole inactivated microbes and their efficacy is controversial. Identification of common antigens can contribute to the development of effective universal or polyvalent vaccines. In this study, we identified common and unique antigens of P. salmonis, A. salmonicida, Y. ruckeri, V. anguillarum and M. viscosa based on a reverse vaccinology pipeline. We screened the proteome of several strains using complete available genomes and identified a total of 154 potential antigens, 74 of these identified antigens corresponded to secreted proteins, and 80 corresponded to exposed outer membrane proteins (OMPs). Further analysis revealed the outer membrane antigens TonB-dependent siderophore receptor, OMP assembly factor BamA, the LPS assembly protein LptD and secreted antigens flagellar hook assembly protein FlgD and flagellar basal body rod protein FlgG are present in all pathogens used in this study. Sequence and structural alignment of these antigens showed relatively low percentage sequence identity but good structural homology. Common domains harboring several B-cells and T-cell epitopes binding to major histocompatibility (MHC) class I and II were identified. Selected peptides were evaluated for docking with Atlantic salmon (Salmo salar) and Lumpfish MHC class II. Interaction of common peptide-MHC class II showed good in-silico binding affinities and dissociation constants between −10.3 to −6.5 kcal mol−1 and 5.10 × 10−9 to 9.4 × 10−6 M. This study provided the first list of antigens that can be used for the development of polyvalent or universal vaccines against these Gram-negative bacterial pathogens affecting finfish aquaculture.


Introduction
Bacterial infectious diseases are a common health concern and intractable economic issue affecting global finfish aquaculture [1]. Some of the most frequent marine bacterial pathogens identified in global aquaculture include Aeromonas salmonicida [2][3][4][5][6][7] Moritella viscosa [8,9], Piscirickettsia salmonis [10][11][12], Vibrio anguillarum [13,14] and Yersinia ruckeri [15][16][17][18][19][20]. Some of these pathogens have been reported to affect important fish species like Atlantic salmon and lumpfish [21][22][23][24][25][26]. utilized in the development of universal or polyvalent vaccines. We evaluated the vaccine potentials of these common epitopes by visualizing their interaction with Atlantic salmon and Lumpfish major histocompatibility complex class II (MHC class II). These epitopes provide a library for the development of a polyvalent or universal vaccine against the marine pathogens considered in this study. The approach employed in this study can be used for the development of a universal vaccine against an even wider range of bacterial pathogens affecting the aquaculture industry.

Acquisition of Proteomes
Protein sequences from complete and annotated genomes utilized are listed in Table  S1 (See Supplementary File S5). P. salmonis, A. salmonicida, M. viscosa, V. anguillarum and Y. ruckeri proteomes were obtained in FASTA format from National Centre for Biotechnology Information (NBCI) public database (https://www.ncbi.nlm.nih.gov/) (accessed on 10 January 2022).

Subcellular Localization
Potential protective antigens must be able to interact with the host molecular receptor that will trigger a signal cascade and ultimately induce a protective adaptive immune response. Therefore, outer membrane proteins and secreted proteins were selected for further analyses based on the pipeline designed for this study (Figure 1). To identify outer membrane proteins and secreted proteins, the protein sequences retrieved from NCBI were analyzed for their subcellular localization using the dynamic Vaxign analysis (http://www.violinet.org/vaxign/) (accessed on 28 December 2021) [60][61][62] and PSORTb V.3.0.2 server (https://www.psort.org/psortb/) (accessed on 28 December 2021) [63]. The Vaxign web server utilizes the vaxign pipeline to dynamically calculate the possibilities of using the proteins as vaccine candidates [52]. Parameters were set to analyze subcellular localization, transmembrane helices, and adhesion probability. Self-antigen screening was performed on vaxign by matching the bacteria antigen to human antigens. Vaxign is not specifically trained to identify antigens for fish, therefore, fish MHC alleles are not included in the vaxign database. The best alternative for self-antigen screening was against human MHC alleles since they have the closest sequence identity to fish MHCs [49]. The subcellular locations assessed were outer membrane, cytoplasmic membrane, cell wall, periplasmic and extracellular regions. The subcellular localization probability scores range from 0 to 1 on Vaxign and 0 to 10 on PSORTb for each cellular location. Only proteins localized on the outer membrane and extracellular regions that have transmembrane helices score of ≤ 1 and no similarity to eukaryotes (human) MHC were selected for further analysis.

Sequence and Structural Alignment
The multiple sequence alignment for the common antigens of the bacterial pathogens was conducted using Clustal 2.1 (clustal omega) with default parameters set to generate a percent identity matrix [66]. The tridimensional protein structure of selected antigens was modeled using HHpred [67]. Structural Alignment of Multiple Proteins (STAMP) on Visual molecular dynamics (VMD) was further used to compare the tridimensional protein structures of the antigens (http://www.ks.uiuc.edu/Research/vmd/) (accessed on 10 January 2022). STAMP aligns protein sequences based on a three-dimensional structure.
Its algorithm minimizes the alpha carbon (Cα) distance between aligned residues of each molecule by applying globally optimal rigid-body rotations and translations [68,69]. This structural alignment method requires that multiple alignments should be used for a single domain evolutionary analysis and not for the alignment of larger multi-domain proteins unless the multi-domain proteins display homology over all domains. Any attempt to align unrelated multiple domain structures or proteins that have no similarities with STAMP may result in alignment failure [68]. The structural homology (QH) values, root mean square deviation (RMSD) and percent identity were generated using Visual Molecular Dynamic (VMD) tools. Structural conservation between proteins (Q) was determined by calculating Structural homology between protein structures [70]. Structures that showed a Q = 1 were considered identical. Structural alignment with Q values <0.3 was considered not well-aligned [70].

Identification of Putative Domains
To further characterize the selected antigens, the presence of putative domains for all common antigens identified in this study was generated by querying the sequence against NCBI conserved domain database (CDD) using the position-specific iterated (PSI) BLAST.

Detection of B and T Cell Epitopes
B-cell epitopes were predicted for selected antigens using tools from the Immune Epitope Database (IEDB) analysis resource BepiPred linear epitope prediction [71]. The BepiPred linear epitope predicts the location of linear B-cell epitopes using a combination of a hidden Markov model and a propensity scale method [71]. The residues with scores > 0.35, sensitivity > 0.49 and specificity > 0.75 were considered to be part of a B cell epitope [71]. Vaxitop, a tool in the Vaxign server [60,61] was used to identify T cell epitopes binding to human MHCs (Table S1). A p-value of~0.05 was used as a cutoff.

Atlantic Salmon and Lumpfish Major Histocompatibility Homology Modelling
Homology modeling of Atlantic salmon and Lumpfish Major Histocompatibility complex class II (MHC class II) was carried out on SWISS-MODEL webserver (https: //swissmodel.expasy.org/) (accessed on 10 January 2022) [72]. Atlantic salmon MHC II α (Q5ZQM4_SALSA) and β (Q95IS0_SALSA) chains were retrieved from the SWISS-MODEL repository. FASTA files of the protein sequence of Lumpfish MHC II α (Accession number: XP_034409868.1) and β (Accession number: XP_034409871.1) chains were uploaded to the webserver program and the tridimensional models were generated. Subsequently, the obtained MHC class II chains were docked together to produce the MHC II complex using the ClusPro V2.0 webserver (https://cluspro.bu.edu/login.php) (accessed on 10 January 2022) [73]. All parameters were set to default settings for the protein-protein docking. The quality of the model was checked using the Procheck webserver (https://www.ebi.ac.uk/ thornton-srv/software/PROCHECK/) (accessed on 10 January 2022) [74].

Epitope Docking of Common Epitopes with Atlantic salmon and Lumpfish MHC Complex
To date, all the available programs to predict epitopes based on MHC complexes do not contain fish MHC alleles. Therefore, in this study, human MHC alleles were used to predict the epitopes. However, it has been reported that human MHC alleles have low sequence identities with fish MHC alleles of~35 to 38% to human MHC sequences [49]. To validate that the predicted common epitopes identified in this study can bind to the Atlantic salmon MHC class II, protein-peptide docking was performed using the online program CABSdock (http://biocomp.chem.uw.edu.pl/CABSdock) (accessed on 10 January 2022). This program performs protein-peptide docking using simulation to search for the binding site on the receptors, which allows for full flexibility of the peptide as well as small fluctuations of the receptor backbone [75]. The online webserver PRODIGY (https://wenmr.science.uu. nl/prodigy/) (accessed on 10 January 2022) was then used to predict binding affinities and dissociation constants of the protein-peptide complexes [76].

Pipeline Development
The reverse vaccinology pipeline used in this study was categorized into selection and characterization phases, respectively ( Figure 1). The selection phase consisted of bacterial proteome acquisition, screening, and selection of common antigens. A total of 5 in-silico tools and one manual filtration were used to achieve this. The characterization phase involved a further description of the common antigens, epitopes and investigating the degree of their similarity across the pathogens used in this study (Table 1). A total of 9 in-silico tools were employed for the characterization phase ( Figure 1).

Subcellular Localization
The Vaxign subcellular localization analysis ( Figure 2) revealed the number of proteins localized in various bacterial cellular regions. As expected, the cytoplasmic region has the highest number of proteins, while outer membrane proteins (OMPs) and secreted proteins were less abundant. On average, we detected approximately 3146 cytoplasmic proteins, 76 OMPs, and 57 secreted proteins for all bacterial proteomes analyzed in this study. M. viscosa showed the highest number of OMPs and secreted proteins, in contrast to P. salmonis which showed the lowest OMPs and V. anguillarum, which showed the lowest number of secreted proteins.  (Table 1). Moritella viscosa has no standard deviation (SD)bar because only one proteome was available for this study.

Adhesin and Antigenicity
The selected OMPs and secreted proteins were evaluated for their ability to function as adhesins and for their antigenicity index, which makes them potential vaccine candidates [49]. The Vaxign vaccine design and Vaxijen informatic packages generated adhesin and antigenicity scores, respectively, for the OMPs and secreted protein. A total 161 immunogenic proteins were identified in this study, 80 of them are secreted proteins, and 81 are exposed OMPs. Hypothetical proteins, which are proteins not yet characterized on NCBI, were not included in this analysis. Based on the description of the immunogenic proteins identified, some proteins are common to the bacterial pathogens used in this study, while others are unique to each bacteria pathogen ( Figure 3). The outer membrane protein assembly factor (BamA), the LPS assembly protein (LptD), and dependent siderophore receptor (TonB) are the 3 common immunogenic proteins identified from the outer membrane proteins across all bacteria proteomes, while flagellar hook assembly protein FlgD and flagellar basal-body rod protein were the two common antigens identified from the secreted immunogenic proteins identified ( Figure 3 and Table 2). Outer membrane transport protein was a common antigen identified in P. salmonis, A. salmonicida, V. anguillarum and M. viscosa but absent in Y. ruckeri. Porin, hemoglobin/transferrin/lactoferrin family receptor and flagellar basal body L-ring protein (FlgH) were common antigens identified in A. salmonicida, V. anguillarum, Y. ruckeri and M. viscosa but absent in P. salmonis. Flagellar filament capping protein FliD was a common antigen identified in P. salmonis, A. salmonicida, Y. ruckeri and M. viscosa but absent in V. anguillarum. Other proteins include the OmpA family protein, outer membrane protein assembly factor BamE, TolC family outer membrane protein, domain-containing proteins, TraV, maltoporin, TonB-dependent receptor, DUF2860 family protein, membrane protein, M23 family metallopeptidase, carbohydrate porin, conjugal transfer protein TraF, murein hydrolase activator NlpD, outer membrane beta-barrel protein, MipA/OmpV family protein, flagellar hook-associated protein FlgK and FlgL, flagellin, phage tail, and others were common to two or three bacteria pathogens (Tables S2 and S3).

Comparative Topology
Multiple sequence alignments of the common antigens identified in this study showed low percentage identity but high structural homology across the various bacterial pathogens (Tables 3 and 4). The 3-dimensional secondary structures of the identified common antigens (OMP assembly factor BamA, TonB dependent siderophore receptor, LPS assembly protein LptD, flagellar basal-body rod protein FlgG and flagellar hook assembly protein (FlgD) showed similar structures (Figures 4 and 5). The structural homology (Q H ) values of the common antigens range from 0.24 to 0.88 with corresponding RMSD values ranging from 0.90 to 6.27, indicating a good structural homology (Tables 3 and 4). Furthermore, the PSI-BLAST analysis of the protein sequence for these common antigens revealed that they have the same conserved domains across the bacterial pathogens (Table S4).   Table S5 for a complete list of the conserved regions highlighted).  Table S5 for a complete list of the conserved regions highlighted).

Epitope Prediction and Docking to Fish MHC Complexes
A total of 64 epitopes was predicted from the common proteins identified using P. salmonis LF 89 as the model for this epitope prediction (Table S4 and Figure 6). Sequence alignment revealed 13 epitopes are conserved across the bacterial pathogens studied (Figure 7 and Table 5). These common epitopes were selected for protein-peptide docking with Atlantic salmon and Lumpfish MHC II. The model for alpha and beta chains for Atlantic salmon MHC II (Figure 8) showed 91.3% (336/368) of all residues in favored (98%) regions and 98.9% (364/368) of all residues in allowed (>99.8%) regions. Protein-peptide docking of the common epitopes with Atlantic salmon MHC II (Figure 8) showed acceptable predicted binding affinities ranging from −10.3 to −6.5 kcal mol −1 and dissociation constants ranging from 5.10 × 10 −9 to 9.4 × 10 −6 M ( Table 5). Protein-peptide docking of the common epitopes with lumpfish MHC II ( Figure S1) showed acceptable predicted binding affinities ranging from −10.7 to −6.5 kcal mol −1 and dissociation constants ranging from 5.0 × 10 −9 to 9.4 × 10 −6 M ( Table 5).

Discussion
Bacterial pathogens affecting finfish aquaculture are relentless, and their impact is enhanced by co-infections [37] and the increasing diversity of pathogen species and serotypes [36]. Commercial polyvalent vaccines available in the aquaculture industry are not providing optimal protection in all finfish species. This comparative reverse vaccinology study aimed to identify common antigens and epitopes across various marine Gram-negative bacteria affecting finfish aquaculture and provide crucial information for the development of effective polyvalent or universal vaccines to prevent bacterial infectious diseases outbreaks in aquaculture.
Gram-negative pathogens synthesize several extracellular and secreted proteins including virulence factors like fimbria, secretion systems, and toxins essential for pathogenicity that can serve as effective antigens [77,78]. In this study, we identified antigens that are exposed onto the bacteria outer membrane or secreted to the extracellular milieu (Table S2). Coincidently, some of these proteins such as BamA, LptD, TonB, FlgD, FlgG, aerolysin, hemolysin have been shown to confer immune protection in fish (e.g., Olive flounder, Zebrafish, Catfish, Grouper fish, Turbot,) against Edwardsiella tarda, Y. ruckeri, V. anguillarum, V. alginolyticus and A. salmonicida [58,[79][80][81]. Adhesin and antigenic characteristics are also important features of a potential vaccine candidate [49]. Secreted and outer membrane proteins that scored adhesin and antigenic probability values of ≥ 0.51 and ≥ 0.4, respectively, were considered for the development of potential vaccine candidates in this study ( Table 2). Five common antigens, including the LPS assembly protein LptD, the outer membrane protein assembly factor BamA, the TonB-dependent siderophore receptor, the flagellar hook assembly protein FlgD, and the flagellar basal-body rod protein FlgG were identified to be conserved across all the strains and species used in this study (Figure 1). These antigens are identical proteins on NCBI across strains of the same species, but a low percentage identity of the antigens across strains of the same species but a low percentage identity of the antigens across different bacterial species used in this study (Tables 3  and 4). However, we hypothesized that these antigens might have important structural or functional similarities because they have the same protein description ( Table 2). Since structural alignments are usually more specific to protein and are more evolutionarily conserved than the amino acid sequence [82], further analyses were conducted to investigate the extent of their structural homology using STAMP on VMD. STAMP has been used to identify the structurally conserved epitope "HAFYLQYKNVKVDFA" among superoxide dismutase homolog structures of Mycobacterium tuberculosis (superoxide dismutase 1IDS), Aspergillus fumigatus (Mn superoxide dismutase 1KKC), and human (Mn superoxide dismutase 2QKC) associated with autoimmune atopic dermatitis [83]. The common proteins identified across all bacterial pathogens in this study were successfully aligned by STAMP on VMD (Figures 4 and 5), which is an indication of homology in their domains. This finding is further supported by the results of the (PSI)-BLAST of the common proteins identified in this study, which showed the presence of similar domains (Table S5). Further analysis of the common antigens revealed a total of 64 novel B cells and T cells epitopes binding to several MHC class I and II alleles (Tables S1 and S4). Sequence alignment of these antigens across P. salmonis A. salmonicida, M. viscosa, Y. ruckeri and V. anguillarum revealed external domains harboring common epitopes (Figure 7 and Table 5) across these bacterial pathogens. These epitopes were predicted using human MHC alleles because the database does not contain fish MHC alleles for teleost. Thus, to evaluate if these epitopes can bind to Atlantic salmon MHC class II, protein-peptide docking was performed to calculate the binding affinities and dissociation constants ( Table 5). This analysis showed that these epitopes can successfully bind to MHC class II from Atlantic salmon and could be used for the polyvalent vaccine design against furunculosis, vibriosis, winter ulcerative disease, piscirickettsiosis, and enteric red mouth disease.
The common antigens in this study play important roles in bacteria cellular activities and are highly conserved in Gram-negative bacteria. Outer membrane protein assembly factor BamA is a key component of the beta-barrel assembly machinery (BAM) complex, which is a key player in OMP assembly [84]. PSI-blast of BamA revealed 5 flexible polypeptidetransport-associated (PORTA) domains (Table S4), which guide substrates into the barrel scaffold and facilitate β-strand formation and assembly. BamA has been identified in a previous study as a potential vaccine candidate against Acinetobacter baumannii [85]. Mice intraperitoneally immunized with 20 µg BamA formulated with 2% Al (OH) 3 adjuvant and intranasally challenged with 10 9 CFU dose −1 of A. baumannii 45 days post-immunization, showed 80% survival. In addition, serum from BamA immunized animals used in passive immunization assays conferred 60% protection to mice after intranasal challenge with a lethal dose of A. baumannii. A. baumannii was not detected in lungs of BamA immunized mice after 7 days post-challenge. Reduction of pro-inflammatory cytokines (TNF-α, IL-6 and IL-1β) and increase in levels of anti-inflammatory cytokine (IL-10) were also observed [85]. These results suggest that BamA might be a good vaccine candidate and could protect the fish against the bath challenge.
LptD has been characterized in bacteria as an essential and well conserved outer membrane protein that mediates the final transport of lipopolysaccharide (LPS) to outer leaflet [86]. Outer membrane protein assembly factor BamA and LPS assembly protein LptD has been identified in silico as a promising vaccine candidate against Moraxella catarrhalis responsible for respiratory tract infections and middle ear infections in children and adults [87]. Numerous potential gonorrhea vaccine targets including BamA and LptD have been identified using proteomics-driven discovery [88]. LptD has also been identified as an immunogenic protein in Vibrio parahaemolyticus and Salmonella enteritca Serovar Typhi using immunoproteomics [89,90]. Subsequently, LptD has been shown to confer 100% protection against a lethal intraperitoneal challenge with 10 7 CFU dose −1 of V. parahaemolyticus in mice and clearance of the bacteria in vaccinated mice [86].
TonB-dependent siderophore receptor mediates the transport of siderophores into the periplasm in Gram-negative bacteria [91]. Other iron-related antigens identified though not common across all strains studied here, include the siderophore amonabactin TonB-dependent receptor and TonB-dependent hemoglobin/transferrin/lactoferrin. TonBdependent siderophore receptor has been described as an immunogenic protein because of the role it plays in iron transport during pathogenesis which is essential in the virulence of bacteria [92,93]. The potential of the TonB receptor as a vaccine candidate has been explored against some Gram-negative bacteria pathogens. For instance, flounder immunized with 100 µL purified recombinant TonB receptor showed strong protective immunity with a relative percent of survival (RPS) of 80.6% and production of serum-specific antibodies against a lethal dose of P. fluorescens challenge [92]. The protective properties of Burkholderia mallei attenuated strain which is a TonB mutant deficient in iron acquisition, have been evaluated in acute inhalational infection models of murine glanders and melioidosis, and it demonstrates great potential as a backbone strain for future vaccine development against both infections [94]. HgbA, a TonB-dependent hemoglobin receptor in Haemophilus ducreyi has been implicated in the virulence of this pathogen in humans and has the potential to be used in vaccine development [95]. In another study, HgbA showed immune protection against chancroid infection in swine model challenged, with 5 × 10 4 CFU of H. ducreyi 3 weeks post-vaccination [96]. These studies support our findings and indicate that TonB-dependent siderophore receptors are a potential candidate for fish vaccine development.
Common secreted antigens identified are the Flagellar basal-body rod protein FlgG and flagellar hook assembly protein FlgD. Flagellins have been shown to bind to the Toll-like receptor 5 (TLR5), which induces innate immune system, triggering a cascade of signaling pathways response [97]. FlgG and FlgD have been explored as potential vaccine candidates in various experiments. Flagellar basal-body rod protein FlgG used in a recombinant vaccine with chaperones Hsp60 and Hsp70 showed high protection, with a relative percent survival (RPS) of 95% against piscirickettsiosis in Atlantic salmon [98].
Olive flounder (Paralichthys olivaceus) challenged with 5 × 10 6 CFU dose −1 of Edwardsiella tarda after vaccination with oral vaccines expressing OmpA-FlgD, FlgD, and OmpA were protected from edwardsiellosis, with survival rates of 82.5%, 55% and 50%, respectively. In this experiment, flounder fed the FlgD-expressing vaccine, showed a significantly increased expression of TLR5M, IL-1β, and IL-12p40, suggesting that the FlgD may be a ligand of olive flounder TLR5M receptor or closely related to the TLR5M pathway [99]. Zebrafish and turbot showed high cellular mediated immune response with high expressions of MHC class I, MHC class II, CD4 and CD8ά genes, which are indications of T cells activation to flagellar hook assembly protein FlgD against E. tarda with RPS of 70% [100]. These studies support our findings and indicate that secreted antigens Flagellar basal-body rod protein FlgG and flagellar hook assembly protein FlgD are potential vaccine candidates that should be explored in fish vaccinology.
An ideal vaccine should induce both innate immunity and long-term stimulation of humoral and cell-mediated adaptive immunity by producing effector and memory cells [101]. This immune response could be achieved by using antigens that possess high epitope density [102]. In this study, 13 common epitopes were identified in the common antigens across the five bacteria species analyzed that could be used as part of polyvalent vaccines (Table 5 and Figure 9). Most in-silico tools for reverse vaccinology currently available are not adapted for fish vaccinology. For example, Vaxitop software used to predict T cell epitopes is not specifically trained to predict peptides based on binding to the fish MHC alleles, but based on binding to human and other mammalian MHC alleles, like mice [49]. In the absence of fish MHC alleles to predict epitopes, epitope predictions were conducted using human alleles, which have an identity of~35 to 38% to fish MHC sequences [49]. Atlantic salmon and lumpfish MHC class II was modeled ( Figure 8 and Figure S1), and a protein-peptide docking was performed with the identified common epitopes to confirm binding to MHC class II (Figure 9). Binding affinity was considered high when Kd < 0.1 nM (1.0 × 10 −10 ), medium when Kd is between 0.1 nM to 1 µM (1.0 × 10 −10 -1.0× 10 −6 ), and low when Kd >1 µM [103]. Furthermore, according to thermodynamics, the binding of a ligand/epitope to a protein only occurs when the system's change in the Gibbs free energy (∆G) is negative and the magnitude indicates the stability and binding affinity of the protein-ligand complex [104]. Our results showed that all the predicted epitopes could bind to lumpfish MHC class II complex with high negative ∆G values with medium to high binding affinities. Hence, the predicted epitopes have a high potential of being recognized by Atlantic salmon and Lumpfish MHC class II to elicit a humoral response in lumpfish. Aside from the common antigens identified, several other antigens were common to two, three or four bacteria species used in this study. Some of these antigens have been reported in other studies to be potential vaccine candidates for fish or other animals. For instance, OmpA was identified to be common in P. salmonis, M. viscosa, A. salmonicida (as porin OmpAI and OmpAII), and V. anguillarum (as porin OmpA). OmpA has been reported to provide 73.3% survival to Labeo fimbriatus that were orally vaccinated with recombinant OmpA encapsulated in chitosan nanoparticles and challenged against E. tarda [105]. A new hybrid OmpA developed by DNA shuffling was demonstrated as an ideal polyvalent vaccine against infections caused V. alginolyticus and E. tarda [106]. Olive flounder fed feeds expressing OmpA, flagellar FlgD, or a fusion antigen of the two for one week showed 55%, 50% and 82.5%, respectively [99]. OmpA protein has been characterized and observed to be highly immunogenic in fish. Common carp vaccinated with recombinant OmpA protein showed a relative percentage survival of 54.3% against E. tarda challenge and elicited a high antibody production in immunized fish [107]. Another antigen identified in this study to be common in up to 3 bacteria species (P. salmonis, V. anguillarum and M. viscosa) and has been reported in other studies to be potential vaccines candidate is the TolC family outer membrane protein. TolC has been identified to be good vaccine candidates against E. tarda and F. columnare in a computer aided vaccine design [108]. Purified recombinant TolC has been shown to confer significant protection with a RPS of 73.68% in hybrid grouper fish against V. harveyi challenge and competitively inhibit the invasion of V. harveyi to grouper embryonic cells in vitro [109]. The M23 family metallopeptidase identified in V. anguillarum and A. salmonicida was evaluated as potential candidate for a broad-spectrum vaccine against Neisseria meningitidis, the causative agent of meningococcal disease in humans [110]. Murein hydrolase activator NlpD identified in Y. ruckeri and V. anguillarum has been reported to be essential for Y. pestis virulence and the nlpD mutant can serve as a basis for developing a potent live vaccine against bubonic and pneumonic plague [111]. OmpV (MipA) family protein identified in Y. rucker and V. anguillarum has been reported to induce high IgG production and protection against Salmonella enterica serovar Typhimurium systemic and gastrointestinal disease in mice. Additionally, it has been shown that OmpV activates innate immune cells, such as monocytes, macrophages, and intestinal epithelial cells, in a Toll-like receptor 2-dependent manner [112]. OmpV can be recognized by both TLR1/2 and TLR2/6 heterodimers, this increases its potential to act as a good adjuvant in other vaccine formulations [112].
Several bacterial toxins and virulent factors identified in this study are unique to specific bacteria species and can be explored in the development of a monovalent or polyvalent multiepitope vaccine in fish vaccinology. Some of these antigens have been reported to be potential vaccine candidates. For instance, hemolysin and aerolysin are extracellular virulent factors identified in V. anguillarum and M. viscosa, respectively. Catfish immunized with recombinant hemolysin (rHly) and aerolysin induced humoral immune response as evidenced by immunoblotting and cell agglutination; immunized catfish showed 71-78% survival when challenged with virulent A. hydrophila and monitored for 2-5 weeks [113]. In another experiment, catfish immunized with the T6SS effector, hemolysin co-regulated protein (Hcp) resulted in an increased survival (46.67%) in common carp during a 10-day challenge time compared to non-vaccinated fish (7.14%). The vaccinated fish also showed significantly increased levels of IgM antibody in serum and cytokines (interleukin-1β and tumor necrosis factor-α) in the kidney, spleen and gills [114]. Mass spectrometric analysis of extracellular products (ECP) from A. hydrophila that conferred 100% protection to catfish vaccinated with ECP when challenged with the pathogen two weeks post-vaccination showed several putative proteins that may serve as important immunogens, including chitinase, chitodextrinase, outer membrane protein 85, putative metalloprotease, extracellular lipase, hemolysin, and elastase [115]. Other enzymes like collagenase, lytic polysaccharide monooxygenase, sphingomyelin phosphodiesterase, methyltransferase, triacylglycerol lipase, exo-alpha-sialidase, deoxyribonuclease I, nucleotidyltransferase, class C betalactamase and SGNH/GDSL hydrolase family protein identified in this study as potential vaccine candidates have been implicated in the virulence of bacterial pathogen [116][117][118][119][120][121][122][123][124][125], and could be considered for a design of an effective polyvalent vaccine for finfish.

Conclusions
There is a need for effective polyvalent vaccines in aquaculture that prevent the emergence of novel virulent variants and reduce co-infections. The increasing number of bacterial genomes available and tools for reverse vaccinology offers a platform for the identification of novel common and uncommon antigens across bacterial pathogens that allows the development of a potential universal anti-bacterial vaccine. We have provided a reverse vaccinology pipeline to identify common and unique antigens across strains of frequent marine bacterial pathogens of finfish, including P. salmonis, A. salmonicida, Y. ruckeri, V. anguillarum and M. viscosa. Potential B cell and T cell epitopes from these common antigens were also identified and docked to Atlantic salmon and Lumpfish MHC class II. Based on these analyses, we identified five common antigens, which contain 13 exposed epitopes that could interact with CD4 + T and B cells. These epitopes could be considered as a base to develop an anti-bacterial polyvalent vaccine for finfish, but additional epitopes from individual antigens need to be incorporated to complete a possible fully functional universal vaccine. In addition, docking with other fish species would be needed to ensure broad immunogenicity. This study contributes to the antigen selection and design of a universal vaccine across a wider range of bacteria affecting the aquaculture industry.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/vaccines10030473/s1, Files S1-S5; Table S1: Categories of MHC HLA-A (MHC class I) and HLA-DP and HLA-DQ (MHC class II) sequences present in the Vaxitop program used in this study. Table S2: Comparative reverse vaccinology outlining outer membrane antigens identified in this study. Table S3: Comparative reverse vaccinology outlining secreted antigens identified in this study. Table S4: Selected B and T (MHC I and II) cell epitopes identified for the common antigens. Table S5. Conserved putative domains of qualifying vaccine candidates.