Abstract
Streptococcus agalactiae is a causative agent of streptococcosis disease in various fish species, including Nile tilapia (Oreochromis niloticus Linn.). Vaccination is an effective disease prevention and control method, but limitations remain for protecting against catastrophic mortality of fish infected with different strains of streptococci. Immunoproteomics analysis of S. agalactiae was used to identify antigenic proteins and construct a chimeric multiepitope vaccine. Epitopes from five antigenic proteins were shuffled in five helices of a flavodoxin backbone, and in silico analysis predicted a suitable RNA and protein structure for protein expression. 45F2 and 42E2 were identified as the best candidates for a chimeric multiepitope vaccine. Recombinant plasmids were constructed to produce a recombinant protein vaccine and DNA vaccine system. Overexpressed proteins were determined to be 30 kDa and 25 kDa in the E. coli and TK1 systems, respectively. The efficacy of the chimeric multiepitope construct as a recombinant protein vaccine and DNA vaccine was evaluated in Nile tilapia, followed by S. agalactiae challenge at 1 × 107 CFU/mL. Relative percentage survival (RPS) and cumulative mortality were recorded at approximately 57–76% and 17–30%, respectively. These chimeric multiepitope vaccines should be applied in streptococcosis disease control and developed into a multivalent vaccine to control multiple diseases.
Similar content being viewed by others
Introduction
Tilapia is a globally economically important aquaculture fish species, particularly in tropical and subtropical countries, such as China and Thailand1. Highly intensive farming systems can increase stress and disease outbreaks in cultured fish2. Bacterial infectious diseases, such as streptococcosis disease caused by the gram-positive pathogen Streptococcus agalactiae, have severe and devastating effects in tilapia aquaculture3. Infected fish can develop various symptoms, including an anorexia phase, hyperemic gills, dermal hemorrhages, dark skin pigment, eye lesions, spinal curvature, erratic swimming, and diffuse epithelial tissue proliferation symptoms3,4.
S. agalactiae is categorized into 11 serotypes: serotypes I-XI5. Serotypes Ia, Ib, and III are the most commonly found serotypes in infected fish6,7. Serotypes Ia and III have been commonly isolated in Thailand7. Vaccines have been used to control disease emergence, including whole-cell inactivated vaccines8,9, live attenuated vaccines10, recombinant vaccines11,12, and DNA vaccines13,14. Although whole-cell vaccines exhibit excellent protection in tilapia, they have limitations in controlling heterologous S. agalactiae serotypes. Subunit vaccines can surpass this limitation using common antigens present in all serotypes, such as ornithine carbamoyl transferase (OCT), pilus island (PI)-1 ancillary protein 115, CAMP factor, R5 protein, enolase, hemolysin (cyLE)16, fibrinogen-binding protein A (FbsA)17, and surface immunogenic protein (Sip)14.
Multivalent and multiepitope vaccines combining at least three segments or epitopes conjugated by linkers have been presented as alternative disease prevention and control strategies18. Various bioinformatics approaches, such as immunoinformatics, molecular dynamics simulation, and protein-protein interaction studies, have been applied to design appropriate and effective multivalent and multiepitope subunit vaccines19. Each individual epitope in a chimeric polypeptide vaccine may provide a high efficacy vaccine by inducing and enhancing robust and specific humoral responses in addition to other cellular responses, particularly opsonization activity20. Moreover, proper linkers have been considered to minimize steric hindrance of each chimeric epitope and enhance epitope presentation to the host immune system21.
Chimeric multiepitope vaccines were generated by combining five different segments of antigenic genes of S. agalactiae on the surface of flavodoxin. Two best chimeric multiepitope vaccines were created by molecular modeling analysis and were produced as a recombinant protein vaccine and a DNA vaccine that were shown to effectively protect against streptococcosis disease in tilapia with different immune response patterns. This platform will elucidate the development of vaccines that combine multiple epitopes from different pathogens to create multivalent vaccines that effectively control fish diseases by single vaccination.
Results
Immunogenic protein characterization
Proteins bound to a S. agalactiae antibody were eluted from protein A agarose and divided into two fractions. The first fraction was subjected to 4–20% gradient SDS-PAGE to observe the protein features and compare the protein profile from serotypes Ia and III. The second fraction was subjected to LC-MS/MS mass spectrometry to identify the immunogenic proteins. The protein profile from the immunoprecipitation on 1D-SDS-PAGE demonstrated that the major protein (approximately 55 kDa) corresponded to rabbit immunoglobulin. However, several bacterial proteins could not be bound to rabbit immunoglobulin and were removed through the flow-through fraction (FT), whereas the protein that specifically bound to the anti-S. agalactiae antibody could be detected in the eluted fraction (Fig. 1).
Analysis of the immunogenic protein profile of S. agalactiae with a 4–20% gradient NuPAGE gel. The flow-through (FT) fraction represents the unbound S. agalactiae proteins, and the eluted fraction represents the immunogenic proteins of S. agalactiae. Markedly different protein patterns between the flow-through fraction and eluted fraction were revealed.
Comparative immunoproteomics analysis of S. agalactiae serotypes Ia and III was determined by LC-MS/MS and assessed by a Venn diagram (Supplementary Fig. 1). One hundred proteins were matched and identified between serotype Ia and serotype III via in-house protein databases, resulting in 79 shared proteins between serotype Ia and serotype III. The protein expression levels of the 79 common proteins were determined by hierarchical clustering (HCL). Two groups of immunogenic proteins were demonstrated based on their abundance, and 37 proteins were overexpressed in serotype III, whereas there was a lower abundance of 39 immunogenic proteins in serotype III than in serotype Ia (Fig. 2). Regarding specific antigen-antibody interactions, 10 and 11 proteins were uniquely identified in serotypes Ia and III, respectively (Supplementary Figs. 1, 2).
Heat map with hierarchical clustering (HCL) of normalized protein abundance reveals the 79 differentially expressed immunogenic proteins. The expression value showed in the relative intensities ranges from the highest protein abundance (red) to the lowest protein abundance (green) expression value.
Linear β-cell epitope prediction and chimeric vaccine design
The epitopes of immunogenic proteins were predicted by the BCPREDS server based on B cell epitopes to be used in chimeric multiepitope vaccine construction. In this study, not only immunogenic proteins from the immunoproteomics analysis were used but also other subunit vaccine candidates were subjected to epitope prediction and combined to produce a chimeric multiepitope vaccine. The amino acid sequences of the C-β protein (BAC), surface protein rib (Rib), LPXTG cell wall anchor domain-containing protein (SPB1), surface immunogenic protein (Sip), and cell surface protein (CSF) showed a different number of epitopes with a BCPredScore of 1.0 (Table 1). However, to construct a chimeric multiepitope vaccine, all epitopes were randomly shuffled on 6 linkers. A total of 7,200 chimeric multiepitope vaccine models were designed.
Molecular modeling of a chimeric multiepitope vaccine
Protein conformation is important for chimeric multiepitope design regarding whether suitable folding can display proper epitopes and maintain high stability. Therefore, a structural domain with an α/β conformation that produces α-helices and parallel β-strands alternatively throughout the backbone likely provides potential bioactivity22. Considering α/β fold structure, flavodoxin from Escherichia coli [PDB accession code: 3CHY]23 was utilized as a linker to combine the epitope fragments from five antigenic proteins.
Predicted epitopes were randomly displayed on the α-helix structure of flavodoxin, generating 1,440 designed models due to the variance of 6 epitopes of BAC and of 2 epitopes of Sip protein. After joining, protein conformation was examined by molecular modeling with 7,200 constructs. I-TASSER and stereochemical qualitative allowance manifested from 45F2 and 42E2 showed appropriate potential tertiary structure with optimal C-scores between −5 and 2. 45F2 and 42E2 also demonstrated the highest score of the amino acid allowance region in the Ramachandran plot. The 45F2 multiepitope model represented 90.2%, 8.5%, and 0.7% of residues located in the most favored, allowed, and disallowed regions, respectively. Meanwhile, the Ramachandran plot regions for the 42E2 designed model comprised 83.0%, 11.6%, and 0.7%, respectively, of the residues (Supplementary Fig. 3). The epitope arrangements in 45F2 and 42E2 were represented in a 3D structure of chimeric proteins, showing that all chosen epitopes were exposed to the protein surface. The five epitopes were displayed as α-helical layers surrounding 3CHY linkers, which appeared as five-stranded parallel β-sheets at the structure’s center, with the order 21345 (Fig. 3).
Schematic diagram and predicted 3D model of designed chimeric multiepitope vaccines. (A) The schematic diagram depicting the chimeric multiepitope vaccine of the 45F2 and 42E2 constructs containing 5 different B cell epitopes in α-helical structures (Rip1, BAC6, BAC5, SPB1_1, Sip2, and CSF1 illustrated in violet, blue, brown, green, yellow and orange, respectively) linked with 6 fragments of a β-pleated sheet of the flavodoxin backbone (gray color). (B) The predicted 3D protein structures of the 45F2 and 42E2 models, which were the two best vaccine candidates, are shown as α/β proteins with a flavodoxin fold. Their colors are represented as colors in the schematic diagram, and their N terminus and C terminus are indicated by arrows. These vaccines were designed based on the desirable construction of the TIM-barrel structure, as shown in Fig. 4B.
Codon optimization of chimeric multiepitope vaccines and plasmid construction
The ectopic expression of bacterial protein in the fish cells may not be achieved due to different codon utilization in the bacterial system. Subsequently, codon optimization of the chimeric multiepitope vaccine was analyzed by GeneArt™’s gene optimization according to ISO 9001 standards (registration no. 1210024212) to apply the codon bias of Oreochromis niloticus. The region of an ideal GC content range—between 30% to 70%—was well optimized. Moreover, negative cis-acting sites included internal TATA-boxes, chi-sites and ribosomal sites; AT-rich or GC-rich sequence stretches; RNA instability motifs; repeat sequences; RNA secondary structures; and splice donor and acceptor sites in higher eukaryotes, which were successfully removed from these chimeric multiepitope DNA vaccine sequences. The best two predicted chimeric multiepitope vaccines were designated 45F2 and 42E2.
Codon adaptation index (CAI) presented 45F2 and 42E2 scores that matched in codon utilization with that of Nile tilapia of 0.92 and 0.93, respectively. The codon quality distribution index of 45F2 and 42E2 demonstrated that the codons within the DNA sequence were distributed frequently in 90–100 positions at 77% and 78% (Supplementary Fig. 4A–D). The average GC content of both chimeric multiepitope vaccines was 56% (Supplementary Fig. 4E,F). Single-stranded RNA-folding prediction revealed the minimum free energy (MFE) secondary structure of 45F2 and 42E2 (Supplementary Fig. 5), together with the free energy of the thermodynamic ensemble at -194.44 kcal/mol and −181.16 kcal/mol, respectively.
N- and O-linked glycosylation sites of 45F2 showed N-glycosylation at N161 (in the CSF protein fragment), whereas 42E2 showed glycosylation at N129 (in the CSF protein fragment). Meanwhile, there were eight (O122, O130, O132, O133, O135, O141, O169, and O171) and four (O92, O93, O95 and O150) putative O-glycosylation sites in 45F2 and 42E2, respectively.
The ProtParam server demonstrated a theoretical pI of 4.1 and a molecular mass of 20 kDa for 45F2 and 42E2. The total number of negatively (Asp and Glu) and positively (Arg and Lys) charged amino acid residues of 45F2 was 33 and 13 residues, while for 42E2, there were 31 and 12 residues, respectively.
The estimated half-life of both chimeric multiepitope constructs was approximately 30 h in mammalian reticulocytes (in vitro), more than 20 h in yeast (in vivo), and over 10 h in E. coli (in vitro). 45F2 showed aliphatic index and grand average of hydropathicity values of 65.32 and −0.401, respectively, whereas 42E2 showed values of 67.93 and −0.296, respectively. The 45F2 and 42E2 proteins were indicated to be stable proteins, as represented by instability indexes of 31.31 and 25.53, respectively.
Antigenicity of the 45F2 and 42E2 chimeric multiepitope vaccines was predicted as 0.7538% and 0.7424% at a 0.4% threshold for the bacterial model, consistent with ANTIGENpro server prediction by representing 0.936 and 0.923, respectively. These results indicate that both vaccine candidates have high potential antigenic properties. Conformational B cell epitopes from the 3D protein structure computed by the DiscoTope server demonstrated 11 B cell epitope residues in both 45F2 and 42E2 at a −3.7 threshold (Table 2)24. Interestingly, the number of epitopes was reduced when computed at the −2.5 and −1.0 thresholds, with 45F2 showing 6 and 3 B cell epitope residue regions, respectively, while 42E2 contained only 2 and 1 B cell epitope residue regions, respectively (Table 2). Recombinant plasmids harboring 42E2 and 45F2 were constructed, namely, pET28a (+)_42E2 or _45F2 and pcDNA3.1 (+)_42E2 or _45F2, which were used to determine the recombinant chimeric multiepitope vaccine expression (Fig. 4).
Four plasmid constructs of chimeric multiepitope vaccines. (A,C) are the expression plasmids of 42E2 and 45F2 in pET28a, a prokaryotic bacterial system; (B,D) represent expression plasmids based on pcDNA3.1, a eukaryotic tilapia cell culture system.
In vitro determination of chimeric multiepitope vaccine expression
Determination of recombinant chimeric multiepitope protein expression was tested in a bacterial expression system and a fish cell (TK-1) culture expression system. These results demonstrated that both chimeric multiepitope proteins could be expressed in both systems, with the expression detectable within 3 h in E. coli (30 kDa) and within 7 days post-transfection in TK-1 cells (25 kDa) (Fig. 5). Larger-sized chimeric multiepitope proteins in the E. coli expression system resulted from an additional tag at the N-terminus, which was contained in the pET28 expression vector.
Western blot analysis of chimeric multiepitope vaccine proteins from the bacterial expression system and fish cell (TK1) culture system. (A) 42E2 and (B) 45F2 proteins expressed in TK1 cells and the E. coli Rosetta-gami B (DE3) pLysS strain at approximately 25 kDa and 30 kDa, respectively (M = prestained protein marker; N = noninduced).
Vaccine efficacy
After vaccination, fish were challenged with S. agalactiae, and infected fish showed clinical signs of streptococcosis disease, such as swirling swimming, opaque eye, exophthalmia and abscess. These moribund fish were collected, and bacteria were re-isolated, showing that they were infected with S. agalactiae serotype III (Supplementary Fig. S7).
DNA vaccine efficacy testing showed that fish immunized with either 45F2 or 42E2 had cumulative mortality rates of 16.67 ± 5.77% and 16.67 ± 15.27%, respectively, which were not significantly different from those of the FKC-vaccinated fish (P > 0.05). However, in the control group [empty vector; pcDNA3.1(+)], 70.00 ± 10.00% mortality was observed at 21 days post-challenge (Fig. 6A).
Analysis of chimeric multiepitope vaccine efficacy. (A) Cumulative mortality for chimeric multiepitope DNA vaccines and recombinant protein vaccines compared with that for the FKC (positive control) and pcDNA3.1(+) (negative control) groups. (B) Relative percentage survival (RPS) of chimeric multiepitope vaccine groups based on that of the control pcDNA3.1 group at 21 days post challenge. Data are represented as the means ± SDs (n = 3). Statistical analysis was performed via a one way ANOVA, compared with the control group (pcDNA3.1). The different letters above the bars indicate significant differences (P < 0.05).
The recombinant chimeric multiepitope protein vaccination showed that 45F2 and 42E2 produced cumulative mortality rates of 30.00 ± 10.00% and 26.67 ± 5.77%, respectively, which were significantly lower than those of the negative control group, at 70% (P < 0.05) (Fig. 6A).
The 45F2 and 42E2 DNA vaccines demonstrated similar patterns of RPS, with 76.19 ± 8.24% and 76.19 ± 21.82%, respectively, which were not significantly different from those of the FKC-immunized fish (76.19 ± 21.82%). However, they were significantly higher than those of the recombinant protein vaccines, which showed 61.90 ± 8.24% and 57.14 ± 14.28% for 42E2 and 45F2, P < 0.05, respectively (Fig. 6B).
Immune response
To determine the immune response, dot blot analysis of serum prepared from 42E2- or 45F2-vaccinated fish was used. It was demonstrated that the DNA vaccine could gradually activate the production of fish antibodies from the 1st to the 4th week. The pattern of antibody response differed from that for the recombinant protein vaccine, with the highest activation of antibody production being significantly produced in the 2nd - 3rd week and suddenly dropping in the 4th week. The highest induction was observed in FKC-immunized fish (Fig. 7). Dot blot analysis of vaccinate fish sera against whole cell lysate of S. agalactiae serotype Ia and III demonstrated that fish vaccinated with recombinant protein vaccine 42E2 and 45F2 showed cross-reactivity to whole cell lysate of S. agalactiae serotype Ia and III (Supplementary Fig. S8).
Dot blot analysis of serum from vaccinated fish. (A) The detectable antibody from vaccinated fish serum is shown on the membrane compared with that for the mock, FKC (formalin killed cell) and pcDNA3.1 group (negative control) against recombinant protein of multiepitope vaccine themselves of FKC vaccine at 1 month. (B) The integrated density from all the dot blot results were converted to the values using the ImageJ 1.x server. Data are represented as the means ± SDs (n = 2). The different letters above the bars indicate significant differences within treatment (P < 0.05).
Discussion
For the reverse vaccinology approach, computational analysis using a variety of bioinformatics tools is robust and beneficial when identifying appropriate vaccine candidates18. Bacterial genomics and proteomics analysis indeed help researchers analyze proteins, short domains, and pathogenic epitopes that provide high immunogenicity and high antigenicity for multimeric vaccine development25. Therefore, immunoproteomics should be applied as a preliminary process to screen antigenic proteins and minimize potential candidates for vaccine development21.
Several immunogenic proteins in this study were described previously, such as C5a peptidase and laminin-binding surface protein (lmb), which are cell surface proteins that have an important function in chemoattractant activities and are proteins promoting invasion of group B streptococcus (GBS)26,27. However, the current immunoproteomics analysis from this study identified new immunoreactive proteins, such as bacteriocin transport accessory protein, dihydrofolate reductase, SSU ribosomal protein S8p, transposase TnpA, 1,4-alpha-glucan, cell wall surface anchor family protein, and the GTP-binding protein Era. As expected, most of these are cell surface proteins, which are suggested to be associated with bacterial virulence27,28. Subsequently, the identified immunoreactive proteins may be used in further vaccine development.
Multiepitope vaccines are an interesting issue since constructed vaccines designed by in silico analysis may elicit cellular immunity and provide effective responses25,29. It is known that immunodominant B cells could strongly induce both cellular and humoral immunity; thus, evaluation of B cell epitopes was performed to identify potent epitopes before integrating them to produce a multiepitope vaccine. Moreover, this vaccine type is more efficient than whole antigens for controlling Staphylococcus spp. infections20,30. From the present study, linear B cell epitope prediction was assessed and identified 11 potent epitopes from 5 common immunogenic and virulence proteins that were present in serotypes Ia and III. The previous studies supported that one of the chosen proteins, Sip, represented a highly conserved protein among GBS isolates and showed cross-protective immunity against GBS infections11,13,14,31. Prediction of candidate antigenic proteins can be used to select the bacterial strains that carry antigenic genes, as well as to determine high expression levels in the target host and the accessibility of particular antigens in host organisms16. Therefore, these selected immunogenic proteins might be suitable for consideration in a rational vaccine design.
Rational chimeric multiepitope vaccine design was achieved by randomly combining epitopes from 5 immunogenic proteins and conjugating with core structures of flavodoxin (β-1–5–3CHY) to produce a secondary structure with α/β folding. In addition to the α/β-type folding of flavodoxin, it was also useful to construct our chimeric multiepitope vaccine by forcing the 5 chosen epitope segments to fit within 5 α-helix loops and protrude out of the 3D-folded structure since that configuration benefits protein solubility by exposure to water molecules23. Additionally, this linker may promote the solubility of the constructed vaccine and help enhance the recognition of the vaccine by the host’s immune system, which contributes to vaccine efficacy. 45F2 and 42E2 presented the most favored region of protein folding, with the stereochemical quality representing the disallowed region at only 0.7%, which is acceptable since the minimum quality should be less than 2%32. It is suggested that in silico analysis could design a chimeric multiepitope vaccine that could probably manifest effective properties33,34.
To achieve a high level of protein expression in Nile tilapia, codon optimization was conducted to improve the transcription and translation capability by removing all possible cis-acting sequence motifs, which may have a negative impact on protein expression. Both proteins had a CAI > 0.8 and a codon with frequent distribution (CFD) > 30%, which are acceptable for high expression in the target organism18,21. The GC content of 45F2 and 42E2 was optimized between 30–70% and had a suitable thermodynamic ensemble free energy, which allowed RNA folding and thermodynamic stability35,36. The overall points suggested that the modeled 45F2 chimeric multiepitope vaccine was clearly the best candidate vaccine.
Numerous effective single-serotype GBS vaccines have been reported, including vaccines for controlling streptococcosis in tilapia10,11,13,14,17,37,38. However, it is known that single-serotype whole-cell inactivated vaccines have limitations during cross-prevention against different serotypes. For instance, a S. iniae vaccine (serotype I) could not protect Atlantic salmon from infection by S. iniae (serotype II)39. Meanwhile, mixed-serotype vaccines (serotypes IV and VII) could promote antiserum levels and enhance the survival rate of newborn pups against streptococcal infection40. Although formalin-killed vaccines generally provide highly protective effects compared with those of subunit vaccines and DNA vaccines, the subunit and DNA vaccines may replace the original formalin-killed vaccines or inactivated vaccines due to their promising efficacy, which are similar to those of inactivated vaccines, and longer shelf life18,33. Evidence suggests that DNA and subunit vaccines can efficiently trigger the immune system and promote protective efficacy, with an RPS value greater than 50%11,13,14,17. Nevertheless, these vaccines have limitations, such as their mass production costs, and they may require various optimizations to obtain the highest stable storage conditions18,41.
Regarding this idea, a chimeric multiepitope vaccine composed of different epitopes from different proteins common in both serotypes Ia and III was generated to achieve broad protection against different serotypes and increase their stability. Interestingly, the designed chimeric multiepitope DNA vaccine and protein vaccine exhibited effective prevention in Nile tilapia against S. agalactiae, with efficacy similar to that of the whole-cell inactivated vaccine. This evidence supports the strategy of rational vaccine design through B cell recognition using in silico analysis. Importantly, immunoproteomics analysis could assist the preliminary determination of suitable immunogenic proteins for vaccine development due to the distinct antigenic determinants that can mediate dissimilar immune responses. The criterion in immunogenic protein selection for vaccine development has focused on the ability of a particular protein to induce an immune response. Among 79 identified proteins shared in both serotypes, in addition to providing the highest BCPRED scores (Table 1), the 5 proteins chosen were also reported as virulence proteins and used as vaccine candidates for streptococcosis disease prevention13. For example, C-β protein (BAC) can lead to antibody production through Fc region binding of human IgA42. Sip protein has been shown to mediate protection against streptococcal infection11,13,14. Additionally, the chosen immunogenic protein should be conserved among Streptococcus spp., so it would be suitable for application in cross-reactive prevention among S. agalactiae serotypes11,43. Moreover, it should be mentioned that peptide vaccines or epitopes with only 30 amino acid residues may trigger immune responses through binding directly to MHC-I or MHC-II molecules. These molecules localize to nonprofessional antigen-presenting cells. Vaccines containing proteins with longer amino acid sequences can enhance the presentation of epitopes to dendritic cells due to T cell induction25,44,45. Herein, the comparative efficacy of both the 45F2 and 42E2 DNA and recombinant protein vaccines indicated that the DNA vaccine provided a higher efficacy than the recombinant protein vaccine. This result suggests that the DNA vaccine can prolong the activation of the immune response by triggering both humoral and cellular immune responses46,47. Moreover, the clearance rate of the recombinant protein vaccine in the host system may be faster than that of the DNA vaccine. This difference implies that the DNA vaccine can enter the host cell to produce chimeric multiepitope protein, with that protein existing in the host system for longer than the recombinant protein vaccine, thus enhancing its bioavailability.
Taken together, these data indicate that the antigen combination has shown promise for streptococcosis disease control in Nile tilapia. This research demonstrated a novel platform for rational vaccine design based on chimeric vaccine development that used flavodoxin with a Tim-barrel structure as a template. Our chimeric protein backbone is suitable for presenting epitopes to be recognized by the host immune system. With 5 epitopes, it could activate antibody production and demonstrated promising protection against bacterial disease similar to that of a whole-cell inactivated vaccine. This platform will promote the production of multivalent vaccines to control multiple diseases and for other applications in the future.
Materials and Methods
Experimental fish, bacterial strain and antibody
All male S. agalactiae-free Nile tilapia (Oreochromis niloticus Linn.) were obtained from a commercial GAP farm in Thailand. The experiments were conducted in accordance with guidelines approved by the National Research Council of Thailand. The experimental fish were anesthetized with clove oil to minimize stress during vaccination and challenge testing. S. agalactiae serotypes Ia and III were cultured as described previously7. S. agalactiae serotype III was used for polyclonal antibody (pAb) production, which was kindly provided by Prof. Ikuo Hirono, TUMSAT, Japan. Antibody against IgM of Nile tilapia was kindly provided by Assist. Prof. Eakapol Wangkahart. Mahasarakham University, Thailand.
Immunoproteomics analysis
S. agalactiae was grown in BHI broth at 30 °C with agitation until reaching exponential phase. Bacterial cells were collected by centrifugation, lysed in 100 µL of lysis buffer [Tris-buffered saline (TBS) with 1% Tween-20 and 0.01% lysozyme] and incubated at 50 °C for 20 min following sonication on ice. Protein A agarose beads (Cell Signaling, USA) were added to the bacterial protein lysate, and nonspecific proteins were removed by 10 min of centrifugation at 10,000 × g at 4 °C. Clarified supernatant was supplemented with 5% glycerol and then with a pAb specific to S. agalactiae serotype III (1:500 dilution). Then, 30 µL of protein A agarose beads were added to separate bound immunogenic proteins, and the bound proteins were separated by acetone precipitation [1:5 (v/v)]. Precipitated proteins were solubilized in 20 mM Tris-HCl with 0.5% SDS, and a Lowry assay was used to measure the protein concentration. The protein profile was assessed by fractionating 25 µg of protein on a NuPAGE 4–12% Bis-Tris Protein gel (ThermoFisher, USA).
3 µg of immunogenic protein was mixed with a lysis buffer (0.1% RapidGest SF in 20 mM ammonium bicarbonate) and 5 mM DTT in 10 mM ammonium bicarbonate at 60 °C for 3 h. This step was followed by incubation with 15 mM iodoacetamide (IAA) in 10 mM ammonium bicarbonate at room temperature for 45 min in the dark. The protein solution was cleaned up by a Zeba Spin Desalting Column before digestion with 50 ng of sequencing-grade trypsin (Promega, Germany) at 37 °C for 6 h. Tryptic peptides were dried at 44 °C under a vacuum and then protonated with 0.1% formic acid in LC water before injection into an LC-MS/MS.
The tryptic peptides’ immunoproteomics profiles were analyzed using an Ultimate™ 3000 Nano/Capillary LC System (Dionex, UK) and Hybrid quadrupole Q-Tof impact II™ (Bruker Daltonics GmbH, Germany) equipped with a Nano-CaptiveSpray ion source. First, 500 nL of extracted peptide was subjected to a trapping column (Thermo Scientific, PepMap100, C18, 300 μm i.d. × 5 mm) through a full loop injection before being resolved in an analytical column (PepSwift C18 Nano Column, 100 μm × 15 cm, i.d.) at 60 °C. The linear gradient method was used to elute peptides with mobile phase A (0.1% formic acid in water) and mobile phase B (0.1% formic acid in 80% acetonitrile) at a 0.35 µL/min constant flow rate into the mass spectrometer. Electrospray ionization was conducted at 1.6 kV using CaptiveSpray. Mass spectra (MS) and MS/MS spectra were fully acquired in positive ion mode (Compass 1.9 for otofSeries software, Bruker Daltonics). Mass accuracy was assessed using positive detection mode after internal calibration with sodium trifluoroacetate (Na-TFA) within 1.6 ppm. Raw LC-MS/MS spectra were collected using CompassXport Version 3.0.9.2 (Bruker Daltonics GmbH, Germany) to convert all spectra into the mzXML data format. The mzXML files of the LC-MS/MS datasets for label-free quantification of peptides were evaluated based on the MS profile by Maxquant software.
Chimeric multiepitope vaccine design
The linear B cell epitope was predicted by BCPRED48. The SCOP and CATH databases were used to design an appropriate chimeric multiepitope vaccine structure49 using BAC (accession no. [AN]: BAE45252), Rip (AN: EAO72273), Sip (AN: AUP09114), CFS (AN: AIK73093), and SPB1 (AN: WP_000913277). A 3D structure was rendered by I-TASSER (Iterative Threading ASSEmbly Refinement) using the qualifying C-score value as a confidence score32. To refine the tertiary structure, the derived I-TASSER results in the PDB files were prepared using the GalaxyRefine server, which performed a repeated structure perturbation, and the best structural relaxation candidates were chosen19. Moreover, to obtain the best chimeric multiepitope vaccine candidates, the residues were determined according to residue stereochemical quality for all the refined chimeric multiepitope models and validated by the PROCHECK program v.3.5.4 to generate Ramachandran plots50.
Codon optimization
Amino acid sequences were reverse-translated to nucleotide sequences using Nile tilapia codon usage (Oreochromis niloticus [gbvrt]: 113). The codon adaptation index (CAI) of the designed vaccine candidates’ nucleotides was analyzed by an optimizer program (http://genomes.urv.es/OPTIMIZER/) and combined with GeneArtTM’s gene optimization process (Thermo Fisher Scientific, USA). The secondary structure of the single-stranded RNA folding and free energy of the thermodynamic ensemble were calculated by the RNAfold web server51. The optimized DNA sequence was synthesized by GeneArt® Gene Synthesis (Thermo Fisher Scientific, USA)52.
Chimeric multiepitope vaccine characterization
N-linked and O-linked glycosylation sites were evaluated using NetNGlyc 1.0 Server53 and NetOGlyc 4.0 Server54. The theoretical pI (isoelectric point), MW (molecular weight), composition of positive and negative residues, estimated half-life, extinction coefficient, aliphatic index, and grand average of hydropathicity (GRAVY) were revealed using the ProtParam server of ExPASy55. Antigenicity was analyzed with the VaxiJen v2.0 server56 and ANTIGENpro software57. The DiscoTope 2.0 server was employed to define discontinuous B cell epitopes at the default threshold of −1.0 to −3.724.
Chimeric multiepitope vaccine preparation
The DNA vaccine and recombinant vaccine expression vector were constructed by inserting the synthesized nucleotides of the 45F2 and 42E2 genes into the pET28a (+) and pcDNA3.1 (+) vectors at the BamHI and XhoI sites. For prokaryotic expression, the pET28a vectors harboring chimeric multiepitope vaccines—45F2 or 42E2—were transformed into E. coli Rosetta-gami (DE3) pLysS strains (Novagen), and protein expression was induced at 30 °C for 3 h with 0.1 mM IPTG.
To verify the ectopic expression of the chimeric multiepitope DNA vaccine, pcDNA3.1(+)_42E2 or _45F2 was transfected into TK1 (Tilapia Kidney 1) tilapia cells using Effectene Transfection Reagent (QIAGEN, Germany). The transfected fish cell cultures were maintained with Leibovitz’s L-15 media containing 10% FBS and penicillin-streptomycin at 25 °C, and DNA vaccine expression was determined after 1 week.
Recombinant chimeric multiepitope protein was purified by Ni-NTA agarose beads (Qiagen) with a gradient concentration buffer of imidazole ranging from 5 mM to 500 mM. Subsequently, the gel filtration chromatography method was performed by fast protein liquid chromatography (FPLC) incorporated with a HiPrep 16/60&26/60 Sephacryl S-300 High-Resolution column (GE Healthcare, USA) using a 1 × PBS buffer with a 1 mL/min flow rate. Recombinant protein detection was confirmed by SDS-PAGE analysis and Western blot analysis using an anti-His tag antibody (recombinant protein vaccine) or an anti-flag (rabbit IgG) (DNA vaccine) and anti-rabbit antibody conjugated to AP (alkaline phosphatase).
Vaccine efficacy analysis
To evaluate vaccine performance, Nile tilapia (O. niloticus) were immunized with chimeric multiepitope vaccines (recombinant protein and DNA vaccines), followed by bacterial challenge. A total of 6 experimental groups, namely, 1) the 45F2 recombinant protein vaccine, 2) 42E2 recombinant protein vaccine, 3) 45F2 DNA vaccine, 4) 42E2 DNA vaccine, 5) formalin-killed (FKC) S. agalactiae vaccine58, and 6) pcDNA3.1(+) [empty vector], were conducted in triplicate. Before vaccination, 25 streptococcosis-free Nile tilapia (60 ± 5 g) were transferred into 18 glass aquarium tanks containing 30 L of water for one week. After a week of acclimatization, fish were vaccinated according to above mentioned groups. All fish were maintained under running and aerated water at 30 ± 3 °C and fed with commercial pellet feed twice a day.
For the chimeric multiepitope protein vaccination, purified 45F2 and 42E2 proteins were mixed with Montanide ISA 763 (Seppic, France) in a 7:3 ratio prior to intraperitoneal injection with 200 µg of protein per fish. For the chimeric multiepitope DNA vaccine, plasmid DNA of 45F2 and 42E2 were purified by ultracentrifugation using a CsCl gradient59 and dissolved in TE buffer (pH 8.0) to obtain a concentration of 0.1 µg/µL. The DNA vaccine was applied to the fish with 10 µg of DNA through intramuscular injection. FKC and pcDNA3.1(+) were used as positive and negative controls, respectively. The schedule of vaccine efficacy analysis and immune response analysis was demonstrated in Supplementary Fig. 9.
For the immune response analysis, blood was drawn from the caudal vein to separate serum for the immunoblotting assay, and those fish were transferred to another separate tank. The analysis was performed every week, using 3 fish in each treatment from the 1st week to the 4th week.
After one month of vaccination, 10 vaccinated fish in each treatment group were taken from among the remaining fish for serum collection and anesthetized with eugenol before challenge with S. agalactiae (serotype III) at 1 × 107 CFU/mL through IP administration. Mortality and clinical signs of infected tilapia were recorded daily for 3 weeks. The brain, head kidney, and liver were collected from moribund fish for bacterial isolation and identification7. Cumulative mortality and relative percentage survival (RPS) were calculated60. A one way analysis of variance (ANOVA) was used for statistical analysis and P < 0.05 was considered significant.
Dot-blot immunoassay
To detect the antibody response after immunization, antibody production was evaluated through dot blot analysis using the Minifold® I dot blot system (GE Healthcare, Germany). Briefly, 20 µL of purified 42E2, 45F2 proteins, or a whole-cell lysate of S. agalactiae (10 µg/mL) were spotted on a nitrocellulose membrane and blocked with blocking solution (0.1% BSA in TBST) before adding 10 µL of serum of the different treatment groups as above mentioned. Then, the membrane was probed with a primary antibody (anti-IgM at 1:5,000) for 1.5 h, followed by washing 3 times with TBST buffer and 45 min of incubation with an anti-mouse IgG HRP-linked Ab (1:10,000). Subsequently, the signal was detected with a ChemiDoc™ Imaging System (Bio-Rad) after adding a substrate reagent (PerkinElmer, USA). The integrated density of the dot blot was analyzed by ImageJ (version 1.x)61.
References
Peng, K. C. et al. Five different piscidins from Nile tilapia, Oreochromis niloticus: analysis of their expression and biological functions. PLOS ONE. 7, e50263, https://doi.org/10.1371/journal.pone.0050263 (2012).
Munang’andu, H. M., Mutoloki, S. & Evensen, O. Prevention and control of viral disease in aquaculture. In Aquaculture Virology. Academic Press: San Diego, CA, USA. 79–93 (2016).
Su, Y. L., Feng, J., Li, Y. W., Bai, J. S. & Li, A. X. Development of a quantitative PCR assay for monitoring Streptococcus agalactiae colonization and tissue tropism in experimentally infected tilapia. J. Fish Dis. 39, 229–238 (2016).
Shoemaker, C. A., Xu, D. H., Evans, J. J. & Klesius, P. H. Parasites and Diseases. (ed. Lim, C. & Webster, C.) 561-582 (Haworth, 2006).
Sabine, B. et al. Streptococcus agalactiae serotype distribution and antimicrobial susceptibility in pregnant women in Gabon, Central Africa. Sci Rep. 5, 17281, https://doi.org/10.1038/srep17281 (2015).
Kawasaki, M. et al. Microevolution of Streptococcus agalactiae ST-261 from Australia indicates dissemination via imported tilapia and ongoing adaption to marine hosts or environment. Appl Environ Microbiol. 84, e00859–18, https://doi.org/10.1128/AEM.00859-18 (2018).
Kannika, K. et al. Molecular serotyping, virulence gene profiling and pathogenicity of Streptococcus agalactiae isolated from tilapia farms in Thailand by multiplex PCR. J Appl Microbiol. 122, 1497–1507 (2017).
Pasnik, D. J., Evans, J. J. & Klesius, P. H. A microwave-irradiated S. agalactiae vaccine provides partial protection against experimental challenge in Nile tilapia Oreochromis niloticus. World J Vaccines. 4, 184–189 (2014).
Pretto-Giordano, L. G., Muller, E. E., Klesius, P. & Da Silva, V. G. Efficacy of an experimentally inactivated S. agalactiae vaccine in Nile tilapia (Oreochromis niloticus) reared in Brazil. Aquacult Res. 41, 1539–1544 (2010).
Li, L. P., Wang, R., Liang, W. W. & Huang, T. Development of live attenuated Streptococcus agalactiae vaccine for tilapia via continuous passage in vitro. Fish Shellfish Immunol. 45, 955–963 (2015).
He, Y., Wang, K. Y., Xiao, D. & Chen, D. F. A recombinant truncated surface immunogenic protein (tSip) plus adjuvant FIA confers active protection against Group B streptococcus infection in tilapia. Vaccine. 32, 7025–7032 (2014).
Nur-Nazifah, M., Sabri, M. Y. & Siti-Zahrah, A. Development and efficacy of feed- based recombinant vaccine encoding the cell wall surface anchor family protein of S. agalactiae against Streptococcosis in Oreochromis sp. Fish Shellfish Immunol. 37, 193–200 (2014).
Huang, L. Y., Wang, K. Y., Xiao, D. & Chen, D. F. Safety and immunogenicity of an oral DNA vaccine encoding Sip of Streptococcus agalactiae from Nile tilapia Oreochromis niloticus delivered by live attenuated Salmonella typhimurium. Fish Shellfish Immunol. 38, 34–41 (2014).
Ma, Y. H. Protective efficacy of cationic-PlGA microspheres loaded with DNA vaccine encoding the sip gene of S. agalactiae in tilapia. Fish Shellfish Immunol. 66, 345–353 (2017).
Dramsi, S. et al. Assembly and role of pili in group B streptococci. Mol Microbiol. 60, 1401–1413 (2006).
Maione, D., Margarit, I., Rinaudi, C. D. & Masignani, V. dentification of universal Group B streptococcus vaccine by multiple genome screen. Science. 309, 148–150 (2005).
Yi, T., Li, Y. W., Liu, L., Xiao, X. X. & Li, A. X. Protection of Nile tilapia (Oreochromis niloticus L.) against Streptococcus agalactiae following immunization with recombinant FbsA and alpha-enolase. Aquaculture. 428, 35–40 (2014).
Nefasat, N. et al. Designing an efficient multi-epitope peptide vaccine against Vibrio cholera via combined immunoinformatics and protein interaction based approaches. Comput Biol Chem. 62, 82–95 (2016).
Hajighahramani, N. et al. Immunoinformatics analysis and in silico designing of a novel multi-epitope peptide vaccine against Staphylococcus aureus. Infect. Genet. Evol. 48, 83–94 (2017).
Zhao, Z. et al. Multiple B-cell epitope vaccine induces a Staphylococcus enterotoxin B- specific lgG1 protective response against MRSA infection. Sci Rep. 5, 12371, https://doi.org/10.1038/srep12371 (2015).
Farhadi, T. et al. Designing of complex multi-epitope peptide vaccine based on Omps of Klebsiella pneumonia: an in silico approach. Int J Pept Res Ther. 21, 325–341 (2015).
Mortenson, D. E. et al. Evaluation of β‐amino acid replacements in protein loops: effects on conformational stability and structure. Chem Bio Chem. 19, 604–612 (2018).
Volz, K. & Matsumura, P. Crystal structure of Escherichia coli CheY refined at 1.7-A resolution. J Biol Chem. 266, 15511–15519 (1991).
Kringelum, J. V., Lundegaard, C., Lund, O. & Nielsen, M. Reliable B cell epitope predictions: impacts of method development and improved benchmarking. PLoS Comput Biol. 8, e1002829, https://doi.org/10.1371/journal.pcbi.1002829 (2012).
Nefasat, N., Ghasemi, Y., Javadi, G., Khoshnoud, M. J. & Omidinia, E. A novel multi- epitope peptide vaccine against cancer: An in silico approach. Journal of Theoretical Biology. 349, 121–134 (2014).
Brown, C. K. et al. Structure of the streptococcal cell wall C5a peptidase. PNAS. 102, 18391–18396 (2005).
Safadi, R. A. et al. Enhanced expression of lmb gene encoding laminin-binding protein in Streptococcus agalactiae strains harboring IS 1548 in scpB-lmb intergenic region. PLoS One. 5, e10794, https://doi.org/10.1371/journal.pone.0010794 (2010).
Foster, T. J., Geoghegan, J. A., Ganesh, V. K. & Hook, M. Adhesion, invasion and evasion: the many functions of the surface proteins of Staphylococcus aureus. Nat Rev Microbiol. 12, 49–62 (2017).
Zhang, L. Multi-epitope vaccines: a promising strategy against tumors and viral infections. Cell Mol Immunol. 15, 182–184 (2018).
Otto, M. Novel targeted immunotherapy approaches for staphylococcal infection. Expert Opin Biol Ther. 10, 1049–1059 (2010).
Brodeur, R. B. et al. Identification of group B streptococcal sip protein, which elicits cross- protective immunity. Infect Immun. 68, 5610–5618 (2000).
Yang, J. et al. The I-TASSER Suite: protein structure and function prediction. Nat Methods. 12, 7–8 (2015).
Farhani, I., Nezafat, N. & Mahmoodi, S. Designing a novel multi-epitope peptide vaccine against pathogenic Shigella spp. based immunoinformatics approaches. International Journal of Peptide Research and Therpeutics, https://doi.org/10.1007/s10989-018-9698-5 (2018).
Pandey, R. K., Bhatt, T. K. & Prajapati, V. K. Novel immunoinformatics approaches to design multi-epitope subunit vaccine for malaria by investigating Anopheles salivary protein. Sci. Rep. 8, 1125, https://doi.org/10.1038/s41598-018-19456-1 (2018).
Ali, M., Pandey, R. K., Khatoon, N., Narula, A. & Mishra, A. Exploring dengue genome to construct a multi-epitope based subunit vaccine by utilizing immunoinformatics approach to battle against dengue infection. Sci. Rep. 7, 9232, https://doi.org/10.1038/s41598-017-09199-w (2017).
Trotta, E. On the normalization of the minimum free energy of RNAs by sequence length. PLoS ONE. 9, e113380, https://doi.org/10.1371/journal.pone.0113380 (2014).
Nawi, M. F., Mohd Yusoff, S., Yusof, H., Abdullah, S. Z. & Saad, M. Z. Efficacy of feed- based adjuvant vaccine against Streptococcus agalactiae in Oreochromis spp. In Malaysia. Aquacult Res. 45, 87–96 (2012).
Suwannasang, A. et al. Growth, immune responses and protection of Nile tilapia Oreochromis niloticus immunized with formalin-killed Streptococcus agalactiae serotype Ia and III. vaccines. Songklanakarin J. Sci. Technol. 39, 429–437 (2017).
Bachrach, G., Zlotkin, A., Hurvitz, A., Evans, D. L. & Eldar, A. Recovery of Streptococcus iniae from diseased fish previously vaccinated with a Streptococcus vaccine. Applied and Environmental Microbiology. 67, 3756–3758 (2001).
Paoletti, L. C. & Kasper, D. L. Conjugate vaccines against group B Streptococcus types IV and VII. The Journal of Infectious Diseases. 186, 123–126 (2002).
Liu, G. et al. Development of Streptococcus agalactiae vaccines for tilapia. Dis Aquat Org. 122, 163–170 (2016).
Areschoug, T., Carlemalm, M. S., Karlsson, I. & Lindahl, G. Streptococcal β protein has separate binding sites for human factor h and IgA-Fc. The Journal of Biological Chemistry. 277, 12642–12648 (2002).
Brodeur, B. R. et al. Identification of group B streptococcal Sip protein, which elicits cross-protective immunity. Infect Immun. 68, 5610–5618 (2000).
Melief, C. J. & van der Burg, S. H. Immunotherapy of established (pre)malignant disease by synthetic long peptide vaccines. Nat Rev Cancer. 8, 351–360 (2008).
Nagpal, G., Usmani, S. S. & Raghava, G. P. S. A web resource for designing subunit vaccine against major pathogenic species of bacteria. Frontiers in Immunology. 9, 2280, https://doi.org/10.3389/fimmu.2018.02280 (2018).
Moss, R. B. Prospects for control of emergimg infectious dyseases with plasmid DNA vaccines. Journal of Immune Based Therapies and Vaccines. 7, 13620, https://doi.org/10.1186/1476-8518-7-3 (2009).
Kutzler, M. A. & Weiner, D. B. DNA vaccines: ready for prime time? Nature Reviews Genetics. 9, 776–788 (2008).
EL-Manzalawy, Y., Dobbs, D. & Honavar, V. Predicting linear B-cell epitopes using string kernels. J. Mol. Recogn. 21, 243–255 (2008).
Fox, N. K., Brenner, S. E. & Chandonia, J. M. SCOPe: Structural Classification of Proteins extended, integrating SCOP and ASTRAL data and classification of new structures. Nucleic Acids Res. 42, 304–309 (2014).
Laskowski, R. A., Rullmannn, J. A., MacArthur, M. W., Kaptein, R. & Thornton, J. M. AQUA and PROCHECK-NMR: programs for checking the quality of protein structures solved by NMR. J Biomol NMR. 8, 477–486 (1996).
Lorenz, R. et al. ViennaRNA Package 2.0. Algorithms for Molecular Biology. 6, 26 (2011).
Buegl et al. DNA synthesis and biological security. Nature Biotechnology. 25, 627–629 (2007).
Gupta, R., Jung, E. & Brunak, S. Prediction of N-glycosylation sites in human proteins. in preparation. Available online, http://www.cbs.dtu.dk/services/NetNGlyc/ (2004).
Steentoft, C. et al. Precision mapping of the human O-GalNAc glycoproteome through SimpleCell technology. EMBO J. 32, 1478–1488 (2013).
Gasteiger, E. et al. Protein identification and analysis tools on the ExPASy Server. (ed. Walker, J. M.) 571–607 (Humana, 2005).
Doytchinova, I. A. & Flower, D. R. Vaxijen: a server for prediction of protective antigens, tumor antigens and subunit vaccines. BMC Bioinf. 8, 1–7 (2007).
Magnan, C. N. et al. High-throughput prediction of protein antigenicity using protein microarray data. Bioinformatics. 26, 2936–2943 (2010).
Asenciosa, Y. O. et al. First report of Streptococcus agalactiae isolated from Oreochromis niloticus in Piura, Peru: Molecular identification and histopathological lesions. Aquaculture Reports. 4, 74–79 (2016).
Wright, J. J., Lee, S., Zaikova, E., Walsh, D. A. & Hallam, S. J. DNA extraction from 0.22 µM sterivex filters and Cesium chloride density gradient centrifugation. J Vis Exp. 31, e1352 (2009).
Sholichah, L., Yuhana, M., Lusiastuti, A. M. & Prihadi, T. H. Potency and efficacy test of a vaccine in addition with adjuvant against Koi Herpesvirus in Koi (Cypronus carpio). Indonesian Aquaculture Journal. 11, 41–47 (2016).
Schneider, C. A., Rasband, W. S. & Eliceiri, K. W. NIH Image to ImageJ: 25 years of image analysis. Nature methods. 9, 671–675 (2012).
Acknowledgements
A.P. is supported by Science Achievement Scholarship of Thailand, SAST, and The Capacity Building of K.U. Students on Internationalization Program: KUCSI, Kasetsart University. This work was partially supported by the Center for Advanced Studies for Agriculture and Food, Institute for Advanced Studies, Kasetsart University Under the Higher Education Research Promotion and National Research University Project of Thailand, Office of the Higher Education Commission, Ministry of Education, Thailand, KURDI (KSD-1-59) and Kasetsart University for research funding*.
Author information
Authors and Affiliations
Contributions
S.U., N.A., S.R., H.K. and I.H. conceived and designed the experiments. S.K. performed immunoproteomics experiment. O.S. conducted bioinformatics analysis. A.P. carried out molecular analysis and animal testing.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
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 license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license 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 license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Pumchan, A., Krobthong, S., Roytrakul, S. et al. Novel Chimeric Multiepitope Vaccine for Streptococcosis Disease in Nile Tilapia (Oreochromis niloticus Linn.). Sci Rep 10, 603 (2020). https://doi.org/10.1038/s41598-019-57283-0
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-019-57283-0
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
