Abstract
Despite widespread vaccination, pertussis resurgence persists due to waning immunity and emerging resistant strains in Bordetella pertussis. This study employs a pangenomics-driven reverse vaccinology approach to identify novel vaccine candidates. Analyzing 160 genomes revealed a closed pan-genome, with approximately 60% conserved genes, including 3389 core genes of which 1312 participate in pathogen-specific pathways. Non-homologous proteins were identified by comparison against human and microbiome proteomes, yielding 205 candidates. Essentiality assessment via the Database of Essential Genes (DEG) refined this to 63 non-homologous essential proteins. A multi-criteria selection process evaluated purifiability based on physicochemical properties and transmembrane helices, accessibility (extracellular or secreted localization), and immunogenicity through antigenicity prediction and B-cell epitope mapping. This pipeline culminated in 11 high-potential vaccine targets. The in silico methodology offers rapidity, cost-effectiveness, and reduced side effects compared with conventional vaccinology, though experimental validation is essential for confirmation.
Introduction
Pertussis, commonly known as whooping cough, is a highly contagious respiratory infection caused by the Gram-negative coccobacillus Bordetella pertussis. This pathogen primarily affects the upper and lower respiratory tracts in humans, leading to paroxysmal cough with or without whoops, leukocytosis, alveolitis, and bronchopneumonia in vulnerable populations like infants, which can result in potentially fatal outcomes. 1
Despite the availability of vaccination programs, pertussis remains a major global public health concern. According to the World Health Organization, pertussis exhibits cyclical epidemic peaks every 2 to 5 years, with multiple large-scale outbreaks reported over the past decade. 2 Although reported cases declined markedly during the COVID-19 pandemic, a substantial resurgence has been observed in the post-pandemic period, reaching alarming levels worldwide.3,4
Vaccination has long been the cornerstone of pertussis prevention and has substantially reduced disease incidence worldwide. The shift from DTwP to DTaP vaccines was driven by the latter’s improved safety profile and tolerability due to reduced reactogenicity, making DTaP preferable for pediatric use. 5 Waning immunity, antigenic variations in circulating strains, and diminished natural boosting from subclinical infections have been associated with the resurgence of pertussis in highly vaccinated populations.5,6 Moreover, DTaP vaccines, which target a limited set of antigens, exhibit limitations that may restrict the breadth and durability of the immune response.7,8 Together, these limitations underscore a critical gap, highlighting the urgent need for next-generation vaccines with enhanced and broader protective efficacy.
The persistence and resurgence of pertussis are further illuminated by the genetic evolution of B pertussis, providing critical context for understanding vaccine inefficacy. Comparative genomic studies have reported that genomic variations (polymorphisms and mutations in key virulence genes) affecting virulence-associated factors and surface antigens have increased, underscoring how evolutionary pressures—such as vaccine selection—may drive adaptation.9-11 This genetic evolution limits the effectiveness of vaccines targeting a restricted number of antigens and highlights the limitations of traditional vaccine target discovery and development approaches, which are often protracted and constrained by the intricate biology of pathogens. 12 In contrast, reverse vaccinology harnesses computational tools to systematically screen for antigenic proteins, accelerating discovery and facilitating comprehensive genomic analysis that allows the capture of strain diversity while identifying conserved and functionally relevant targets.13,14
Pangenomics, a subset of this paradigm, examines the complete genomic repertoire across strains to pinpoint conserved elements, such as those involved in essential metabolic processes or signaling pathways, ensuring broad applicability irrespective of strain-specific variations.15,16 By integrating multiple genomes, this approach enables the identification of the core genome and conserved genes shared across all strains, as well as the non-core genome; classified into accessory and singleton genomes; that contributes to phenotypic variability and adaptation. 17
This study addresses these challenges by employing pangenomics reverse vaccinology—a sophisticated, data-driven framework—to identify novel vaccine candidates from B pertussis genomes. By integrating pangenomics with subtractive genomics, we identified and targeted core genome; encoded proteins that are essential for bacterial survival or persistence, while excluding homologous proteins to human and microbiota proteomes to enhance specificity, these proteins are surface-exposed or secreted and immunity-induced for certainty.
Materials and Methods
This in silico study was carried out at the research technology platforms of Mohammed VI University of Sciences and Health (UM6SS) in Casablanca, Morocco. To systematically identify potential vaccine targets, a pangenomics-driven in silico approach was employed. This strategy is based on the comparative analysis of available B pertussis genomes to define conserved, pathogen-specific gene sets, followed by a rational prioritization framework. The overall workflow, illustrated in Figure 1, integrates pangenomics analysis and successive bioinformatic filtering steps to select proteins with high vaccine potential. This approach aims to enhance target specificity while minimizing cross-reactivity with host and microbiota components.
Schematic representation of the pangenomics-driven reverse vaccinology pipeline applied to prioritize targets in B pertussis leading to the identification of 11 putative vaccine and therapeutic targets.
Pangenome Analysis
Prior to initiating the pangenome analysis, a series of preparatory steps are required, starting with the acquisition of genome data. Ten sequencing reads of B pertussis isolates in FASTQ format were retrieved from the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) (Home–SRA–NCBI) (Supplemental Table S1). 18 This retrieval was followed by quality control using FastQC (v0.12.1) (https://github.com/s-andrews/FastQC.git), 19 after which adapters and low-quality reads were trimmed with Fastp (v0.24.1) (https://github.com/OpenGene/fastp.git). 20 The trimmed raw reads, originally downloaded from NCBI-SRA, underwent de novo assembly employing SPAdes (v4.2.0) (https://github.com/ablab/spades.git), with default parameters, to generate a novel complete genome of B pertussis, 21 with the assembly outcomes subsequently assessed via Quast (v5.3.0) (https://github.com/ablab/quast.git) to evaluate standard assembly metrics such as number of contigs, total assembly length, largest contig, N50, and L50. 22 Finally, genome annotation was performed using Prokka (v1.14.5) (https://github.com/tseemann/prokka.git) with the parameters –outdir, –prefix bordpt, and –locustag BP, specifying the output directory, file prefix and locus tag prefix used for the analysis. Prokka used to delineate additional attributes and essential elements (such as coding sequences, rRNA, tRNA, and others) within the de novo assembled genome. 23 Therefore, the pangenome analysis was conducted using Roary (v3.13.0) (https://github.com/sanger-pathogens/Roary.git), analyzing a set of 150 Prokka-generated GFF3 files of complete B pertussis genome (Supplemental Table S2), retrieved from NCBI RefSeq database and the newly annotated genome, all of which were quality-checked prior to analysis using Quast. The purpose of this analysis was to construct the pan-genome and delineate core versus non-core gene sets for comparative genomic interpretation. 24 The pangenome analysis was run with the parameters -e -n -v -p 8 -f with default settings including the default 95 % BLASTP identity threshold and default paralog splitting behavior, which separates paralogous gene copies into distinct clusters and excludes unresolved paralogs from the core alignment.
Identification of Host and Pathogen Metabolic Pathways
Through a systematic workflow (Figure 1), this study represents the first computational and subtractive genomics analysis of various metabolic pathways to identify potential vaccine targets in B pertussis. To derive metabolic pathway-based information, 2 steps were performed; the first step is functional annotation, involved uploading a list of 3389 core genes to an automatic annotation server called BlastKOALA (v3.1) (https://www.kegg.jp/blastkoala/) developed by KEGG (Kyoto Encyclopedia of Genes and Genomes) to correlate each gene with its corresponding KEGG Ortholog (KO) so the metabolic pathway. 25 Data were extracted with reference to the metabolic pathways and identification numbers for the human host. Subsequently, a comparison was conducted between the B pertussis and human pathways, and pathways present in the pathogen but absent in the host were selected as unique to B pertussis.
Identifying the Non-Homologues and Unique Proteins
This step was performed to eliminate proteins from common pathways which serve to minimize potential side effects, the perturbation of the protein function and avoid inducing an immune response against the microbiota. To execute this step, 2-step comparisons were conducted between human and pathogen proteomes to identify non-homologous proteins of B pertussis and between the microbiota and pathogen proteomes to identify unique proteins of B pertussis. Initially, all human protein were retrieved from Refseq database ftp site (v May 23, 2018) (Index of /refseq/H_sapiens/mRNA_Prot) and only proteins from pathogen-specific pathways were subjected to BLASTP (v2.16.0) analysis, with E-value equal to 10−3, against the human proteome. 26 Subsequently, non-homologous proteins were compared using BLASTP, with E-value equal to 10−3, analysis against the microbiota proteome, which yield from NIH Human microbiome project (https://www.ncbi.nlm.nih.gov/datasets/genome/?bioproject=PRJNA43021&assembly_level=3:3). In each case, hits were filtered based on E-value > 10−3, identity < 30%, query coverage < 50%, and others as default. Proteins that did not produce hits meeting these criteria were selected as non-homologous and unique proteins.
Identification of Essential Proteins for Bordetella pertussis
After identifying all non-homologous and unique proteins, we further filtered them based on their essentiality, specifically, targeting proteins critical for B pertussis survival. Essential proteins in cellular organisms are indispensable for viability and replication, making them promising targets for antimicrobial treatments. We compared the list of these proteins to a database of experimentally verified essential genes from approximately 30 bacterial species in the Database of Essential Genes (DEG: v Sep. 1, 2020) (http://origin.tubic.org/deg/public/index.php/download). 27 An E-value threshold of 10−5 was used to screen essential proteins of B pertussis using DEG microbial BLASTP. Significant hits were obtained by filtering based on E-value <= 10−5, identity ⩾ 30% and a minimum query coverage of 50%.
Prioritization of Vaccine Targets
B pertussis essential genes were further evaluated on the basis of their subcellular localization, physicochemical properties, and immunogenicity analysis. This analysis is mandatory to characterize the antigenicity of these essential genes and validate them as vaccine targets.
Subcellular Localization
The subcellular localization prediction and secretion screening was performed on the non-homologous and essential proteins identified in B pertussis. PSORTb (v3.0.3) (https://psort.org/psortb/), was applied first, generating confidence scores from 0 to 10 for each compartment and retaining predictions with the highest score (⩾7.5); proteins scoring below this threshold were further analyzed with CELLO (v2.5) (https://cello.life.nctu.edu.tw/) using normalized probability scores and a reliability index (0-5).28,29 SignalP (v5.0) (https://services.healthtech.dtu.dk/services/SignalP-5.0/) was additionally employed to detect signal peptides and secretion pathways based on probability scores ⩾ 0.7. In cases of disagreement or low confidence, consensus evaluation was applied, prioritizing PSORTb predictions and using CELLO and SignalP for complementary validation. 30 The retained proteins predicted as surface-exposed or secreted (outer membrane, periplasmic, or extracellular), which were advanced to the downstream physicochemical and antigenicity analyses.
Physicochemical Parameters
Initially, information regarding non-homologous, unique, and essential proteins was evaluated to prioritize suitable vaccine targets; however, to minimize the time required for vaccine testing and development, it is preferable to evaluate these vaccine targets by several molecular and structural criteria. 31 This involved calculating the molecular weight, which must be less than 110 kDa, and the GRAVY score, which should be less than or equal to zero. These 2 parameters facilitate the purification of targeted proteins. 32 The stability index, which must be less than or equal to 45, indicates protein stability. For pI, there is no critical value, but an extreme pI value could cause issues; thus, the pI should range between 5 and 9 and differ from human pH. These calculations were performed using the ProtParam (https://web.expasy.org/protparam/) computational web tool. 33 The next criteria involved employing TMHMM (v2.0) (https://services.healthtech.dtu.dk/services/TMHMM-2.0/), which used to predict the number of the transmembrane helices of the proteins. The selection was made if a protein has 1 or 0 transmembrane helices.34,35
Prediction of Antigenicity and Immunogenicity
Proteins classified as PSE or secreted are those that first interact with the immune system 36 ; thus, we predicted their ability to elicit a strong immune response by evaluating their antigenicity using VaxiJen (v2.0) (https://www.ddg-pharmfac.net/vaxijen/VaxiJen/VaxiJen.html). VaxiJen employs to its alignment-independent approach, which evaluates protein antigenicity based on the physicochemical properties of amino acids rather than sequence similarity. Proteins with a score > 0.4 were considered probable antigens with significant immunogenic potential.37,38 Subsequently, we proceeded to predict of potentially immunogenic epitopes in a given protein sequence, that may significantly reduce wet lab effort needed to discover the epitopes required for the design of vaccines and for immunodiagnostics. The goal of B cell epitope prediction was to identify immunogenic epitopes on potential antigens, encompassing both linear (continuous) and discontinuous (conformational) categories, capable of binding to B cell receptors, thereby activating B lymphocytes and eliciting a humoral immune response characterized by the secretion of neutralizing antibodies.39,40 Linear B-cell epitopes were predicted employing BepiPred (v3.0) (https://services.healthtech.dtu.dk/services/BepiPred-3.0/) with a threshold of 0.5 and a minimum residue length of 7 amino acids. 41 For conformational epitopes, which necessitate a 3D protein structure model, the Protein Data Bank (PDB) (https://www.rcsb.org/) was initially queried to ascertain the availability of structures for the target proteins; in cases where no structure existed, a homology model was generated using SWISS-MODEL (https://swissmodel.expasy.org/). 42 The corresponding .pdb files were subsequently imported into the ElliPro tool, applying a protrusion index threshold ⩾0.7 and a minimum residue length of 7 amino acids, which uses a combination of geometric features and the propensity of single amino acids in an antigen for conformational epitope prediction.39,43
Results
Analysis Pangenomics
Using Roary_plots for visualization of pangenomics analysis results give an overview of B pertussis pangenomics type (Figure 2). A total of 5747 genes clusters was found and classified into 3 categories, the core genome (core and soft core) with 3442 genes which present in more than 95% of strains (143 strains) occupied 59.9% of genome which mean a strong conservation of genes and the non-core genome which shell genome with 553 genes present between the 15% and 95% of strains and cloud genome correspond to 1752 of genes which present in less than 15%. The Shell and Cloud genome are occupying, respectively, 9.6% and 30.5% of genome with a total of 60.1%. This distribution, core genome is more than 50% and the ratio non-core genome on core genome strict (without soft genome) equal around 0.75 which is inferior to 1, confirms that the B pertussis pangenomics is closed, which is also confirmed by Parkhill et al. 9
Distribution of 5747 genes into the core (core + soft-core) and the non-core (shell + cloud) genome categories.
Identifying the Non-Homologues and Unique Proteins
Proteins identified as human homologs were excluded in the initial step, as such targets could adversely affect host metabolism. Unique metabolic pathways specific to B pertussis, but absent in the human host, encompass proteins that represent promising candidates for drugs and vaccines. In addition, certain proteins within shared pathways between the pathogen and host may be exclusive to B pertussis; those involved in multiple pathways, particularly if non-homologous, emerge as strong vaccine candidates owing to their essential functions. Submission of the B pertussis core gene list to KEGG yielded annotations for 1993 genes across 48 pathways, including 1312 genes associated with pathogen-specific pathways and 275 genes shared with humans. Proteins within unique pathways, such as lipopolysaccharide biosynthesis, oxidative phosphorylation, and bacterial secretion systems, were deemed robust vaccine targets. Protein sequences from both common and unique pathways were subsequently aligned against the human proteome via NCBI BLASTP to detect non-homologous proteins, resulting in 1010 sequences exhibiting no matches to the human proteome and thus likely to avoid eliciting autoimmune responses. These non-homologous proteins were further compared against the microbiota proteome, with 205 proteins identified as unique to B pertussis and selected for subsequent evaluation.
Identification of Essential Proteins for Bordetella pertussis
Being non-human and involved in unique metabolic pathways are not the sole criteria for selecting favorable vaccine targets. Identifying proteins that regulate key factors, such as nutrient uptake, growth, secretion systems, virulence, and pathogenicity, is of great importance for disrupting pathogen functions and survival.27,44 Such proteins are considered essential for the pathogen. However, not all essential proteins are non-homologous. Therefore, pathogen proteins that are both non-homologous and essential represent more attractive vaccine targets. 45 The non-homologous proteins identified through resistance gene analysis were further screened based on their essentiality. A total of 63 B pertussis proteins were found to be essential, as they are involved in activities significant to the pathogen’s biological pathways (Table 1).
Non-Homologous and Essential B pertussis Proteins Retained After Comparison Against DEG Database.
Prioritization of Vaccine Targets
Sub-Cellular Localization
Subcellular localization is a key functional attribute of a protein, as many cellular processes are compartmentalized; thus, predicting the localization of unknown proteins can offer valuable insights into their roles and assist in selecting candidates for further study. 31 In bacteria, all proteins are synthesized in the cytoplasm, where many remain to perform specific functions, while others contain export signals that direct them to different cellular compartments. In Gram-negative bacteria, compartments include the cytoplasmic membrane, periplasm, inner and outer membranes, and extracellular space. 28 Although most proteins are confined to a single compartment, some span multiple locations.
The subcellular localization of 63 non-homologous essential proteins of B pertussis was evaluated using a combination of prediction tools: PSORTb, followed by cross-validation with CELLO and SignalP. Of these, 27 proteins were predicted to be cytoplasmic, 15 membrane-associated, 18 periplasmic, 9 located in the inner membrane, and 3 extracellular (Figure 3).
A multi-selection screening pipeline for the characterization and prioritization of 63 essential proteins using subcellular localization, signal peptide, transmembrane topology, physicochemical, and antigenicity prediction tools (Made using Napkin).
Subcellular localization is a critical factor in identifying effective drug and vaccine targets; proteins localized in the cytoplasm or membrane are often considered for drug development, while those secreted or exposed on the bacterial surface are ideal candidates for vaccines. 46 The next step involves refining these candidates further through multilevel analysis of physicochemical parameters.
Physicochemical Parameters
We initially screened the secreted or extracellularly exposed proteins of B pertussis for potential vaccine targets, followed by the application of additional prioritization parameters, including physicochemical properties, to refine the selection. These properties were critical in the early stages of vaccine discovery. Most proteins had a molecular weight of less than 110 kDa and a GRAVY score of ⩽ 0, indicating solubility (hydrophilic proteins), which supports their potential for experimental study in vaccine development. In addition, an instability index of ⩽45 signifies protein stability. Although the isoelectric point (pI) lacks a critical threshold, proteins with extreme pI values may pose challenges; thus, a pI range of 5 to 9.6 is preferred, though not strictly codified. Furthermore, estimation of transmembrane helices in the proteins revealed all 11 proteins to have either 0 or 1 transmembrane helices, (except Cytochrome c biogenesis protein) which to its easy ability to be easily purified in experimental analysis (Figure 3 and Table 2).
Analysis of Physicochemical Properties (Molecular Mass, pI, Stability, GRAVY), Subcellular Localization, Transmembrane Helix Number, Signal Peptide Prediction, and Antigenicity.
Prediction of Antigenicity and Immunogenicity
Antigenicity
The proteins classified as secreted or PSE were analyzed using VaxiJen software. This analysis revealed that all 11 proteins have a probability above the threshold of 0.4, indicating that they are immunogenic and capable of generating antibodies, thus making them strong candidates for vaccine targets (Figure 3 and Table 2).
Immunogenicity
Peptide vaccines are more convenient and safer than conventional vaccines, as they include only immunogenic epitopes rather than the full 3-dimensional protein structure. B-cell epitopes are antigenic determinants that are recognized and bound by receptors on the surface of B lymphocytes. A total of 103 linear B-cell epitopes were identified across the 11 examined proteins of B pertussis, employing a threshold score exceeding 0.5 and a minimum epitope length of 7 amino acids (Table 3). For conformational epitope detection, 3-dimensional structures of all proteins were modeled using the SWISS-MODEL homology-modeling server, with the structure exhibiting the highest sequence identity and coverage selected for each protein. The resulting PDB files were subsequently imported into the ElliPro tool, applying a threshold of 0.7 and a minimum epitope length of 7 amino acids. This analysis identified 24 conformational B-cell epitopes from the 11 proteins under investigation (Table 3).
Number of Predicted Linear and Conformational Epitopes for the 11 Candidate Targets.
Discussion
Face to decline of vaccine efficacy and increase of pertussis cases reported around the world, it was suggested the necessary of finding novel therapeutic strategy to improve the vaccines and to combat this resurgence. The purpose of this study was to identifying new vaccine targets against B pertussis, based on pangenomics approach, via a list of 3389 core genes. The efficacy of a vaccine depends on the careful selection of appropriate antigens. Critical criteria for antigen evaluation include (1) non-homology to human proteins, to minimize cross-reactivity and potential adverse effects; (2) accessibility, determined by the antigen’s subcellular location, to ensure it can be targeted by the immune system; and (3) immunogenicity, defined as the ability to stimulate B-cell and/or T-cell responses. These factors collectively enhance the vaccine’s potential to activate protective immunity while avoiding interference with host protein functions.
The pangenomics approach, which is applied in this study, successfully reduced the number of candidate proteins from 1010 to 11. These proteins were selected as vaccine targets for treating B pertussis infections due to their high potential among identified non-homologous proteins and their significant role in the pathogenicity of B pertussis, suggesting that targeting them could disrupt the pathogen’s survival and colonization. Of these 11 proteins, 5 have experimentally confirmed antigenicity and immunogenicity: 2 pertussis toxin (Ptx) subunits S2 (BP_01857) and S5 (BP_01859),6,47 2 LPS assembling proteins (BP_02637 and BP_02906), 48 and the outer membrane protein assembly factors BamB (BP_01207), 48 regarding the outer membrane protein (BP_03126) belonging to the general bacterial porin (GBP) family, needs further characterization to confirm its identity. 49
Two pertussis toxin (Ptx) subunits, S2 (BP_01857) and S5 (BP_01859), are widely recognized as major virulence factors expressed by B pertussis, which contribute significantly to bacterial adhesion and immune modulation, making them key components of B pertussis vaccines (DTaP and DTwP).6,47 Our identification of S2 and S5 through a pangenomics approach confirms their conservation and relevance as vaccine targets. However, vaccine-induced selective pressure has led to observed variability in pertussis toxin genes, which may affect the effectiveness of subunits such as S2 and S5.50,51 Consequently, a vaccination strategy targeting stable genomic regions and continuous monitoring of antigenic variability are required.
Two LPS assembling proteins (BP_02637 and BP_02906) are highly appeared to be the novel vaccine targets, replacing LPS, which are the major component of DTwP and the main responsible for the reactogenicity. 47 Assembling proteins, such as LptD, whose immunogenicity was demonstrated in a study done by Dorji, showed that mice vaccinated with LptD alone exhibited significantly lower bacterial loads compared with those vaccinated solely with DTaP. 48
The outer membrane protein assembly factors BamE (BP_00805) and BamB (BP_01207) are components of the β-barrel assembly machinery (BAM) complex. This complex facilitates the insertion and folding of β-barrel proteins into the outer membrane (OM) of Gram-negative bacteria. 52 Targeting the BAM complex holds potential not only for vaccine development but also for developing potentiators that enhance OM permeability. The immunogenicity of BamB was demonstrated in vivo by Dorji, who showed that mice vaccinated with BamB alone exhibited significantly lower bacterial loads compared with those vaccinated solely with DTaP. 48 However, the immunogenicity of BamE remains to be fully evaluated and warrants further investigation.
Outer membrane protein (BP_03126) belonging to the general bacterial porin (GBP) family are found in the outer membrane of Gram-negative bacteria, which is responsible for antibiotics permeability. 53 Among the proteins of this family in B pertussis, OmpP is the only porin identified as highly immunogenic, making it a promising and essential vaccine candidate with confirmed humoral immune responses in both animals and humans. 49 Further characterization of this outer membrane protein is necessary to confirm its identity and optimize its use in vaccine development.
Hemolysin activation/secretion protein, FhaC (BP_00399),is an outer membrane protein that acts as a transporter in the 2-partner secretion pathway (TPS). Its essentiality has been confirmed in the secretion of filamentous hemagglutinin (FHA). 54 Furthermore, Xu and Hu have noted that secretion systems have changed under vaccine pressure, and they suggested that targeting FhaC is an innovative approach. 55 Therefore, it is suggested that the immunogenicity of FhaC needs to be demonstrated.
The outer membrane usher protein, FimC (BP_00400), is predicted to be an integral outer membrane (OM) β-barrel usher protein involved in fimbriae (Fim2/3) assembly. 56 Despite the critical role of FimC in fimbriae assembly and its potential as a novel vaccine target to disrupt fimbriae assembly and reduce B pertussis colonization, no studies to date have specifically investigated its immunogenicity.
The outer membrane protein (BP_00850), part of the multidrug efflux system, is a relevant element within the microbial repertoire that contributes to multidrug resistance (MDR) and the failure of anti-infectious therapies 57 ; thus, it is considered an attractive target for the development of new vaccines. It has been reported that the multidrug efflux system of B pertussis exhibits low activity, which poses challenges for vaccine manufacturers due to the pathogen’s sensitivity, potentially reducing the yield of biomass for vaccine production. 58
Protein involved in the biogenesis of cytochrome c (BP_00041), an essential component of bacterial energy metabolism, is critical for cellular function. The entire system appears to operate on the outer surface of the cytoplasmic membrane, 59 which limits its potential as a direct vaccine target, but it could be exploited for antimicrobial strategies and used as a potential drug target.
Conclusions
Unlike classical vaccinology methods, which rely on in vitro bacterial cultivation for antigen identification, reverse vaccinology leverages genomic analysis to predict vaccine targets more rapidly and cost-effectively, while enhancing specificity and minimizing risks of autoimmune reactions. These predicted targets subsequently undergo immunogenicity assessments to precisely characterize the elicited immune responses, thereby facilitating advanced testing of potential candidates. In this study, 160 B pertussis genomes, the etiological agent of whooping cough, were analyzed using bioinformatics approaches to identify novel vaccine targets. This methodology yielded 11 proteins as putative vaccine and drug targets, several of which require rigorous in vivo validation to confirm their immunogenicity, protective efficacy, and safety profile before integration into novel vaccine formulations. Collectively, the findings underscore the substantial relevance of this investigation to pertussis vaccine development and highlight the necessity of in vivo testing to optimize the application of these results.
Supplemental Material
sj-docx-1-bbi-10.1177_11779322261442522 – Supplemental material for Pangenomics-Driven Reverse Vaccinology for the Discovery of New Vaccine Candidates Against Bordetella pertussis
Supplemental material, sj-docx-1-bbi-10.1177_11779322261442522 for Pangenomics-Driven Reverse Vaccinology for the Discovery of New Vaccine Candidates Against Bordetella pertussis by Majid Abdelouahab and Diawara Idrissa in Bioinformatics and Biology Insights
Supplemental Material
sj-docx-2-bbi-10.1177_11779322261442522 – Supplemental material for Pangenomics-Driven Reverse Vaccinology for the Discovery of New Vaccine Candidates Against Bordetella pertussis
Supplemental material, sj-docx-2-bbi-10.1177_11779322261442522 for Pangenomics-Driven Reverse Vaccinology for the Discovery of New Vaccine Candidates Against Bordetella pertussis by Majid Abdelouahab and Diawara Idrissa in Bioinformatics and Biology Insights
Supplemental Material
sj-docx-3-bbi-10.1177_11779322261442522 – Supplemental material for Pangenomics-Driven Reverse Vaccinology for the Discovery of New Vaccine Candidates Against Bordetella pertussis
Supplemental material, sj-docx-3-bbi-10.1177_11779322261442522 for Pangenomics-Driven Reverse Vaccinology for the Discovery of New Vaccine Candidates Against Bordetella pertussis by Majid Abdelouahab and Diawara Idrissa in Bioinformatics and Biology Insights
Footnotes
Consent and Ethics
Not applicable (this is an in silico study only).
Author Contributions
Funding
The authors received no financial support for the research, authorship, and/or publication of this article.
Declaration of Conflicting Interests
The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Supplemental Material
Supplemental material for this article is available online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
