INTRODUCTION

Aeromonas schubertii is a motile, rod-shaped, Gram-negative bacterium characterized by a single polar flagellum and taxonomically assigned to the family Aeromonadaceae. Members of the genus Aeromonas are widely distributed across diverse ecological niches, including soil, freshwater and marine environments, and the intestinal tracts of animals, reflecting their considerable environmental adaptability [1, 2]. A. schubertii has been recognized as an opportunistic pathogen of relevance to human health and has been associated with acute gastroenteritis, with environmental exposure and direct contact representing potential transmission routes [3]. Beyond its public health relevance, A. schubertii has emerged as an important pathogen in aquaculture. Infections have been documented in several economically important aquatic species, including the marine shrimp Penaeus monodon and Litopenaeus vannamei [3, 4], mussels [5], and Nile tilapia (Oreochromis niloticus) [6]. More recently, A. schubertii has been implicated in severe disease outbreaks in snakehead fish, including Channa argus, Channa striata, and Channa maculata, with reported mortality reaching 45% over 40 days [2, 7, 8]. Such outbreaks can cause substantial economic losses through reduced fish survival, increased treatment costs, and decreased farm productivity.

Notably, A. schubertii-associated outbreaks in freshwater fish are characterized by distinctive white nodular lesions distributed across internal organs, an atypical pathological manifestation that is infrequently reported in other bacterial infections of fish [7, 9]. The combination of substantial mortality and characteristic pathological changes highlights the importance of timely and accurate identification of A. schubertii for disease diagnosis, epidemiological surveillance, and outbreak management in aquaculture systems.

Accurate species-level identification of Aeromonas remains a persistent challenge in diagnostic microbiology. Conventional biochemical methods are time-consuming and may produce ambiguous or erroneous results because of phenotypic variability and overlapping characteristics among closely related species [10]. Molecular approaches and whole-genome sequencing provide greater taxonomic resolution; however, their routine application may be constrained by methodological complexity, cost, and reported inconsistencies in species discrimination [11–13]. These limitations can delay reliable species-level identification and complicate epidemiological investigations and disease-management decisions. Routine diagnostic workflows for A. schubertii commonly integrate classical bacteriological techniques with polymerase chain reaction (PCR)-based molecular assays [14–16]. These approaches may involve multiple sequential steps, including bacterial isolation, cultivation, physiological and biochemical characterization, and confirmatory molecular testing, resulting in labor-intensive, time-consuming workflows. Therefore, rapid, reliable, and high-throughput species-level identification remains important for strengthening diagnostic capacity, particularly in Southeast Asian aquaculture settings.

Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) has emerged over the past two decades as a powerful approach for microbial identification in clinical and research laboratories [17–19]. MALDI-TOF MS offers several practical advantages, including rapid analysis, relatively low per-sample operating costs, high-throughput capacity, and the ability to generate identification results within minutes after colony isolation. However, its identification accuracy depends strongly on the quality, diversity, and completeness of the reference spectral database [20]. Consequently, microorganisms absent or poorly represented in commercial databases may receive low-confidence scores, be identified only at the genus level, or be misclassified.

Despite the increasing application of MALDI-TOF MS in aquatic animal diagnostics, reliable species-level identification of A. schubertii remains constrained by its limited species-specific spectral representation in existing reference databases. This limitation is particularly important for isolates associated with white nodular disease in freshwater fish, for which rapid differentiation from closely related Aeromonas species is required for accurate diagnosis and epidemiological surveillance. Although previous work has investigated MALDI-TOF MS-based differentiation of multiple fish-associated Aeromonas species, a dedicated species-specific main spectra profile (MSP) database focused on A. schubertii has not been established. To our knowledge, no previous study has developed and validated a dedicated MSP database specifically for A. schubertii, particularly using isolates associated with white nodular disease in freshwater fish. Therefore, the lack of a curated species-specific spectral reference library represents an important diagnostic and methodological gap that may limit the accuracy and confidence of MALDI-TOF MS-based identification of this emerging aquatic pathogen.

To address this diagnostic gap, the present study aimed to develop and validate a species-specific MSP database for A. schubertii using the MALDI Biotyper Compass platform (Bruker Daltonics). Specifically, the study sought to generate high-quality and reproducible MSPs from representative A. schubertii isolates, assess whether incorporation of the customized MSP database improved species-level identification compared with the existing commercial reference database, evaluate identification repeatability, and determine the discriminatory capacity of the developed database against closely related non-target Aeromonas species using spectral and multivariate analyses. By focusing on A. schubertii, the study was designed to provide a practical framework for improving rapid, reliable identification of this underrepresented pathogen and to support future diagnostic surveillance of A. schubertii-associated disease in freshwater aquaculture.

MATERIALS AND METHODS

Ethical approval

All procedures involving bacterial handling, culture, and MALDI-TOF MS sample preparation were conducted in a Biosafety Level 2 (BSL-2) laboratory in accordance with the institutional biosafety requirements of Chulalongkorn University. The laboratory procedures were reviewed and approved by the Institutional Biosafety Committee of Chulalongkorn University (Approval No. 2031043). Moribund fish used for diagnostic sampling were euthanized at the sampling sites by immersion in clove oil (1 g/L) for approximately 10 min in accordance with the previously described procedure [20]. Internal organs were subsequently collected aseptically for bacteriological examination.

Study period and location

The study was conducted over a three-month period, from August to November 2025. The bacterial isolates investigated in this study originated from diseased freshwater fish sampled from July 2019 to July 2020 in Thailand and Vietnam. Thai isolates were obtained from Samut Sakhon and Nakhon Pathom Provinces, whereas Vietnamese isolates originated from Tra Cu District, Tra Vinh Province; Chau Phu District, An Giang Province; and Thanh Binh District, Dong Thap Province. The isolates were submitted to the Center of Excellence in Fish Infectious Diseases, Faculty of Veterinary Science, Chulalongkorn University, Thailand, for diagnostic examination and further characterization. The six isolates selected for MSP database construction were recovered from diseased snakehead fish collected in Nakhon Pathom Province, Thailand, in 2020.

Study design

This laboratory-based diagnostic method development and validation study aimed to construct and evaluate a species-specific MALDI-TOF MS MSP database for A. schubertii. Among 81 bacterial isolates obtained from diseased freshwater fish, 19 were identified as A. schubertii and included in the study. Six representative isolates (CUVETASC03–CUVETASC08) were used for MSP library construction, whereas the complete set of 19 isolates, comprising the six MSP development isolates and 13 additional isolates, was used to assess identification performance. Four closely related non-target Aeromonas species were included as outgroup controls for assessment of spectral discrimination and clustering. Molecular confirmation of the six MSP development isolates was performed using 16S rRNA and gyrB sequence analyses.

Bacterial isolation and identification

Among 81 bacterial isolates obtained from diseased freshwater fish during 2019–2020, 19 were identified as A. schubertii (Table 1). Of these, 17/19 isolates were recovered from striped snakehead (Channa striata), one from Asian seabass (Lates calcarifer), and one from hybrid red tilapia (Oreochromis sp.).

Moribund fish were euthanized at the sampling sites by immersion in clove oil (1 g/L) for approximately 10 min, as previously described [20]. Samples from internal organs, including the liver, spleen, and kidney, were collected aseptically using sterile inoculating loops and streaked onto tryptic soy agar (TSA; Bacto™, Becton, Dickinson and Company [BD], Sparks, MD, USA) supplemented with 5% sheep blood. The inoculated plates were incubated at 28°C for 24 h. Representative colonies with distinct morphologies were selected for preliminary identification using conventional microbiological procedures, including Gram staining and catalase and oxidase tests. The resulting phenotypic characteristics were consistent with classification within the genus Aeromonas. Confirmed isolates were preserved in tryptic soy broth (TSB; Bacto™) supplemented with 20% glycerol and stored at −80°C until further analysis.

PCR assay, DNA sequencing, and phylogenetic analysis based on 16S rRNA and gyrB sequences

Species-level identification of the six isolates selected for MSP development was confirmed by PCR amplification and sequencing of the 16S rRNA and gyrB genes. The 16S rRNA gene was amplified using primers 27F (5′-AGAGTTGATCATGGCTCAG-3′) and 1492R (5′-GGTTCACTTGTTACGACTT-3′) [21], whereas gyrB was amplified using primers gyrB3-F (5′-TCCGGCGGTCTGCACGGCGT-3′) and gyrB14-R (5′-TTGTCCGGGTTGTACTCGTC-3′) [9, 22]. The expected amplicon sizes were approximately 1500 bp for the 16S rRNA gene and 1100 bp for gyrB.

Genomic DNA was extracted from the tested isolates using the boiling method [23]. Amplification of the 16S rRNA gene was performed with an initial denaturation at 94°C for 5 min, followed by 35 cycles of denaturation at 94°C for 1 min, annealing at 56°C for 30 s, and extension at 72°C for 90 s, followed by a final extension at 72°C for 10 min. Amplification of gyrB was performed with an initial denaturation at 94°C for 2 min, followed by 30 cycles of denaturation at 94°C for 30 s, annealing at 56°C for 30 s, and extension at 72°C for 2 min, followed by a final extension at 72°C for 5 min.

PCR products were visualized on 1.5% agarose gels stained with RedSafe™ and purified using the NucleoSpin® Gel and PCR Clean-up Kit (Macherey-Nagel, Germany). Purified amplicons were subjected to Sanger sequencing, and forward and reverse sequences were assembled using BioEdit version 7.1.1. Species identification was confirmed using Basic Local Alignment Search Tool (BLAST) against the National Center for Biotechnology Information (NCBI) database, applying sequence-identity thresholds of ≥99.5% for 16S rRNA and ≥99.0% for gyrB. Sequences were aligned using ClustalW, and maximum-likelihood phylogenetic trees were constructed using MEGA X version 10.2.6 with the Tamura–Nei model and 1,000 bootstrap replicates. Gaps and missing data were handled using the partial-deletion option with a site-coverage cutoff of 95% [24–26].

Preparation of samples for MALDI-TOF MS

All A. schubertii isolates were recovered from frozen stocks and subcultured on TSA (Bacto™) at 28°C for 24 h. Well-isolated colonies were subjected to protein extraction using the standard ethanol–formic acid–acetonitrile protocol for MALDI-TOF MS analysis. Briefly, 1 µL of extracted protein was spotted onto a 96-position polished-steel MALDI target plate, overlaid with α-cyano-4-hydroxycinnamic acid (HCCA) matrix, and air-dried at room temperature. Samples were analyzed using a microflex® LRF MALDI-TOF MS instrument (Bruker Daltonics, Bremen, Germany), and the resulting spectra were compared with the MALDI Biotyper® (Bruker Daltonics) reference library [27].

MALDI-TOF MS extraction and database construction

Six representative A. schubertii isolates (CUVETASC03–CUVETASC08) were selected for MSP library cons-truction to incorporate multiple independent isolates and capture intraspecies spectral variability. These isolates originated from diseased snakehead fish collected in Nakhon Pathom Province, Thailand, in 2020. The isolates were recovered from glycerol stocks, subcultured on TSA, and subjected to ethanol–formic acid–acetonitrile protein extraction.

Briefly, bacterial biomass was suspended in 300 µL of high-performance liquid chromatography-grade water and 900 µL of absolute ethanol, centrifuged, air-dried, and subsequently extracted with 30 µL of formic acid and 30 µL of acetonitrile. A 1 µL aliquot of the resulting supernatant was spotted onto a polished-steel MALDI target plate and overlaid with HCCA matrix. Following mass calibration, smoothing, normalization, baseline correction, peak detection, and outlier removal, spectra were analyzed using FlexAnalysis software version 3.0 (Bruker Daltonics) in positive linear mode over a mass range of 2–20 kDa.

For MSP generation, each isolate was analyzed using eight independent spots, with three spectra acquired from each spot, yielding 24 raw spectra per isolate. Spectra showing weak signal intensity, irregular baselines, peak shifting, or poor reproducibility were excluded. A minimum of 20 high-quality spectra per isolate were retained for MSP construction using MBT Compass Explorer® (Bruker Daltonics) [28, 29]. MSPs were accepted when reproducible peak patterns and high intrastrain spectral similarity were consistently observed across technical replicates.

PCA, spectral stacking, and dendrogram analyses were performed using ClinProTools® version 3.0 (Bruker Daltonics), and the resulting MSPs were incorporated into an in-house A. schubertii database [27, 30]. Four closely related non-target species, Aeromonas hydrophila, Aeromonas veronii, Aeromonas jandaei, and Aeromonas caviae, were included as outgroup reference species for discriminatory and clustering analyses. These species were selected because of their phylogenetic relatedness to A. schubertii and their frequent association with freshwater fish diseases [27, 31].

Validation method

The identification performance of the developed MALDI-TOF MS database was evaluated using 19 A. schubertii isolates, with particular emphasis on identification reliability and repeatability. Discriminatory performance was further assessed using the four closely related non-target Aeromonas species as outgroup controls.

MALDI-TOF MS identification scores were interpreted according to the manufacturer's recommendations. Scores <1.70 were classified as unreliable; scores of 1.70–1.99 indicated reliable genus-level identification; scores of 2.00–2.29 indicated reliable species-level identification; and scores ≥2.30 indicated highly reliable species-level identification. Repeatability was assessed by analyzing each isolate in triplicate and calculating the coefficient of variation (CV) of the corresponding log score values. A CV ≤10% was considered indicative of acceptable repeatability. Discriminatory performance was assessed by comparing protein spectral profiles and evaluating clustering patterns between A. schubertii and the four non-target Aeromonas species.

Statistical analysis

Log score values were summarized as mean ± standard deviation from triplicate technical replicates. Repeatability was evaluated using the CV of log score values, with values ≤10% considered indicative of acceptable repeatability. Pearson correlation coefficients were calculated between corresponding processed spectral-intensity profiles of technical replicates to evaluate spectral reproducibility. PCA, spectral stacking, and dendrogram analyses were performed using ClinProTools® version 3.0 (Bruker Daltonics).

RESULTS

MALDI-TOF MS identification performance: commercial versus custom MSP database

The identification performance of MALDI-TOF MS before and after incorporation of the custom MSP database is summarized in Table 2. Using the commercial Bruker database alone, 15/19 (78.95%) Aeromonas schubertii isolates achieved reliable species-level identification, with log scores ranging from 2.01 to 2.22. The remaining 4/19 (21.05%) isolates were identified only at the genus level, with log scores ranging from 1.95 to 1.98. These findings indicated incomplete species-level resolution of A. schubertii using the commercial reference database.

After incorporating the custom MSP database, all 19 isolates (100%) were correctly identified as A. schubertii at the species level, with highly reliable log scores ranging from 2.49 to 2.68. Thus, species-level identification increased from 78.95% with the commercial database to 100% with the custom MSP database, accompanied by higher identification confidence across all tested isolates. These findings demonstrate the diagnostic value of targeted enrichment of reference spectra for improving MALDI-TOF MS identification of A. schubertii.

The six representative A. schubertii isolates used for MSP construction exhibited distinct and reproducible protein mass spectral profiles. Although spectra were acquired across a mass range of 2–20 kDa, the characteristic profiles presented in Figure 1 were visualized predominantly within the 2–16 kDa range. Each isolate exhibited characteristic peak patterns, demonstrating both intraspecies spectral consistency and detectable interstrain variation.

Technical replicate spectra showed high reproducibility, with Pearson correlation coefficients ranging from 0.95 to 0.97. These findings indicated a high degree of similarity among replicate protein fingerprint profiles and supported the reproducibility of the MSPs generated from the six representative isolates.

Spectral discrimination among Aeromonas species

In addition to the six A. schubertii isolates, four closely related non-target species, Aeromonas hydrophila, Aeromonas veronii, Aeromonas jandaei, and Aeromonas caviae, were included as outgroup controls. All isolates generated analyzable MALDI-TOF MS spectra following spectral preprocessing, including baseline correction, smoothing, normalization, peak detection, and outlier removal. Spectral stacking demonstrated distinct protein mass spectral patterns among A. schubertii and the non-target Aeromonas species (Figure 2).

Principal component analysis (PCA)

The A. schubertii isolates CUVETASC03–08 showed a distinct and relatively broad distribution in PCA space, particularly along PC1 and PC2 (Figure 3), indicating considerable intra-species variation among the isolates. Their separation from the reference Aeromonas species suggests that CUVETASC03–08 possess a distinct overall multivariate profile.

Dendrogram analysis

Dendrogram analysis further demonstrated clear differentiation between A. schubertii and the non-target Aeromonas species. The A. schubertii isolates clustered tightly within a distinct group, whereas A. hydrophila, A. veronii, A. jandaei, and A. caviae were separated from the A. schubertii cluster. The two major clusters merged at a linkage distance of approximately 0.75, whereas most within-cluster distances were <0.40, indicating comparatively high intraspecies spectral similarity and greater interspecies spectral divergence (Figure 4).

The clear separation of the four non-target Aeromonas species from A. schubertii in both PCA and dendrogram analyses supported the discriminatory capability of the custom MSP database.

PCR assay, DNA sequencing, and phylogenetic analysis

PCR amplification of the 16S rRNA and gyrB genes was successfully performed for all six isolates selected for MSP construction (CUVETASC03–CUVETASC08). Clear single amplicons corresponding to the expected sizes of approximately 1500 bp for the 16S rRNA gene and 1100 bp for gyrB were observed, consistent with the positive controls (Figure 5). Subsequent Sanger sequencing followed by BLASTn analysis against the NCBI nucleotide database demonstrated high sequence similarity with reference A. schubertii sequences.

Sequence analysis demonstrated 99.68%–99.80% identity with reference A. schubertii strains across the evaluated 16S rRNA and gyrB sequences (Table 3). These findings supported the species-level taxonomic assignment of the six isolates used for MSP construction. Maximum-likelihood phylogenetic analyses were subsequently performed using approximately 1,500 bp of the 16S rRNA gene and 1100 bp of gyrB.

The maximum-likelihood phylogenetic trees incorporated reference A. schubertii sequences retrieved from GenBank together with selected non-A. schubertii taxa. In both the 16S rRNA- and gyrB-based trees, isolates CUVETASC03–CUVETASC08 clustered with established A. schubertii reference strains and formed a distinct clade, further supporting their taxonomic identification as A. schubertii (Figure 6).

DISCUSSION

The present study developed a species-specific MALDI-TOF MS MSP database for A. schubertii and demonstrated its potential to improve species-level identification of this emerging aquatic pathogen. The custom database addressed limitations observed with the commercial reference database, increasing species-level identification from 78.95% to 100% and producing highly reliable identification scores for all tested isolates. In addition, PCA and dendrogram analyses demonstrated clear spectral differentiation of A. schubertii from the evaluated fish-associated Aeromonas species. Given the increasing importance of A. schubertii in Southeast Asian aquaculture and its reported zoonotic potential, these findings support the value of species-specific spectral enrichment for improving the diagnosis and surveillance of an underrepresented aquatic bacterial pathogen.

MALDI-TOF MS profiling and species differentiation of Aeromonas schubertii

The application of MALDI-TOF MS enabled characterization of A. schubertii isolates based on their protein mass spectral profiles. The generated spectra were highly reproducible among technical replicates, with Pearson correlation coefficients of 0.95–0.97, indicating consistent protein fingerprint profiles. Although gel stack visualization showed similarities among members of the genus Aeromonas, multivariate analyses provided greater discriminatory resolution. PCA demonstrated that A. schubertii isolates formed a distinct cluster separated from the evaluated non-target Aeromonas species, whereas dendrogram analysis similarly showed tight clustering of A. schubertii and separation from Aeromonas hydrophila, Aeromonas veronii, Aeromonas jandaei, and Aeromonas caviae.

The observed differentiation is consistent with previous findings indicating that Aeromonas species can exhibit characteristic protein mass fingerprints that facilitate species-level discrimination by MALDI-TOF MS [29]. The clustering pattern observed for A. schubertii is also consistent with its placement within Aeromonas group 501 reported by Benagli et al. [32], a group characterized by distinct phenotypic and genotypic features [33]. Collectively, the present spectral and clustering results support the capacity of MALDI-TOF MS to discriminate A. schubertii from the closely related species included in this study when appropriate reference spectra are available.

Limitations of commercial MALDI-TOF MS spectral databases

Despite the advantages of MALDI-TOF MS, the commercial Bruker reference database provided incomplete species-level resolution for A. schubertii. Of the 19 isolates examined, 15 (78.95%) were reliably identified at the species level, whereas four (21.05%) achieved only genus-level identification. This finding illustrates an important limitation of reference database-dependent identification: diagnostic performance is influenced by the representation and quality of spectra available for the target organism. Underrepresentation of environmental and aquatic bacterial species in reference libraries may consequently reduce identification confidence for less frequently encountered pathogens such as A. schubertii.

Identification performance may also be affected by methodological differences in sample preparation. Variations in extraction solvents and acid concentrations can influence protein recovery and spectral quality, thereby affecting the resulting identification profiles [34–36]. Organic acid-based extraction procedures preferentially recover abundant low-molecular-weight proteins, including ribosomal and nucleic acid-binding proteins; therefore, differences in extraction chemistry can alter the detectable protein profile [37]. Standardization of both sample preparation procedures and reference spectral libraries is consequently important when MALDI-TOF MS is applied to less well-represented aquatic bacterial species.

Improvement of identification accuracy through a customized MSP database

Development of the custom A. schubertii MSP database substantially improved identification performance. Following incorporation of species-specific reference spectra, all 19 isolates (100%) were correctly identified at the species level, with log scores ranging from 2.49 to 2.68. In comparison, the commercial database yielded species-level identification for only 78.95% of the isolates, with log scores of 2.01–2.22 among those identified at the species level. Thus, targeted addition of high-quality A. schubertii reference spectra improved both the proportion of isolates identified at the species level and their associated identification confidence.

The custom database was evaluated using 19 A. schubertii isolates recovered from multiple freshwater fish hosts, including snakehead fish, Asian seabass, and hybrid red tilapia. This host representation provided an opportunity to evaluate the database against isolates originating from different fish hosts. However, the MSP library itself was constructed from six isolates obtained from diseased snakehead fish in Nakhon Pathom Province, Thailand. Therefore, the present findings should not be interpreted as demonstrating that the MSP library captures the full host-associated or geographic diversity of A. schubertii. Broader validation using isolates from additional hosts and geographic regions remains necessary.

The improvement achieved after database customization is consistent with previous studies demonstrating that MALDI-TOF MS identification of Aeromonas spp. is strongly dependent on reference database completeness [20, 29]. Difficulties in discriminating fish-associated Aeromonas species using commercial databases have previously been reported [29], whereas customized MSP libraries containing additional reference spectra have improved identification performance [20]. The increase in species-level identification from 78.95% to 100% observed in the present study further supports targeted spectral enrichment as a practical strategy for improving MALDI-TOF MS identification of underrepresented aquatic bacterial pathogens.

The capacity to expand MALDI-TOF MS databases also offers opportunities to incorporate spectra representing broader host, geographic, and strain diversity, potentially increasing their utility for epidemiological surveillance and outbreak investigations [38, 39]. Nevertheless, these potential applications require validation using appropriately diverse isolate collections before broader epidemiological utility can be established.

Molecular validation using 16S rRNA and gyrB gene analyses

Molecular analyses provided independent support for the taxonomic identification of the six isolates used for MSP construction. Sequence analysis of the 16S rRNA and gyrB genes demonstrated 99.68%–99.80% identity with reference A. schubertii strains, and maximum-likelihood phylogenetic analyses showed that all six isolates clustered with established A. schubertii reference strains. These findings supported the species-level assignment of the isolates used to generate the custom MSP database.

Although 16S rRNA gene sequencing is widely used for bacterial identification and phylogenetic characterization, its discriminatory capacity may be limited among closely related Aeromonas species, including A. veronii, Aeromonas ichthiosmia, and Aeromonas simiae. In contrast, housekeeping genes such as gyrB can provide greater discriminatory resolution because of their higher sequence variability while retaining conserved functional roles [40]. Furthermore, multilocus sequence analysis based on concatenated housekeeping genes has been reported to provide improved taxonomic and phylogenetic resolution within the genus Aeromonas [41–43]. The combined use of 16S rRNA and gyrB analyses in the present study therefore provided complementary molecular evidence supporting the taxonomic identity of the isolates used for MSP construction.

Implications for aquaculture diagnostics and pathogen surveillance

The findings indicate that MALDI-TOF MS, when supported by an appropriately curated species-specific MSP database, can provide rapid and reliable identification of A. schubertii following bacterial isolation. This capability may be particularly relevant to surveillance of white nodular disease, which has been associated with A. schubertii infection in freshwater fish. Because A. schubertii has also been reported in other aquatic hosts and opportunistic human infections, improved species-level identification may have relevance beyond freshwater fish diagnostics.

Combining MALDI-TOF MS protein fingerprinting with molecular confirmation also provides complementary approaches to bacterial characterization. Molecular methods can support taxonomic confirmation, whereas MALDI-TOF MS can facilitate rapid processing of cultured bacterial isolates once an appropriate reference database has been established. This combination may therefore strengthen laboratory identification and epidemiological surveillance, particularly for aquatic pathogens that are poorly represented in commercial spectral databases.

The customized MALDI-TOF MS approach may also offer practical workflow advantages compared with sequencing-based identification because MALDI-TOF MS enables rapid analysis of cultured colonies and supports high-throughput processing. However, this study did not formally compare turnaround time, workload, and cost. Consequently, the current data cannot quantify the relative economic and operational advantages of the custom MALDI-TOF MS approach over molecular identification methods. Future studies should directly compare diagnostic turnaround time, laboratory workload, and cost-effectiveness between these approaches.

Beyond its immediate application to A. schubertii, the developed MSP database may provide a foundation for expansion of reference spectra for aquatic Aeromonas species. Continued addition of well-characterized isolates could support the development of a broader aquatic Aeromonas MSP reference network. Future technological advances may also permit integration of such databases with portable or field-oriented MALDI-TOF MS platforms; however, the feasibility and diagnostic performance of such applications require dedicated validation.

Limitations and future perspectives

Several limitations of the present study should be considered when interpreting the findings. First, the MSP database was constructed using only six A. schubertii isolates, all originating from diseased snakehead fish collected in Nakhon Pathom Province, Thailand. Although the database was subsequently evaluated using 19 isolates originating from Thailand and Vietnam and from multiple fish hosts, this collection may not adequately represent the broader genetic and proteomic diversity of A. schubertii across different geographic regions, hosts, and production systems. Second, MALDI-TOF MS performance can be influenced by sample preparation procedures, instrument settings, calibration, and laboratory-specific workflows, which may affect the transferability and reproducibility of the database across laboratories.

Third, validation was performed within the study setting and did not include independent external or multicenter laboratory validation. Although the discriminatory analyses included four closely related non-target Aeromonas species, a broader panel of phylogenetically and diagnostically relevant bacterial species would provide a more comprehensive assessment of specificity. Fourth, whole-genome sequencing was not performed; therefore, direct comparison of MALDI-TOF MS identification with genome-based taxonomic approaches was beyond the scope of the study. Finally, formal analyses of diagnostic turnaround time and cost-effectiveness were not undertaken.

Future studies should expand the MSP database using a larger collection of molecularly and genomically characterized A. schubertii isolates representing diverse hosts, geographic regions, and epidemiological settings. Multicenter validation should assess interlaboratory reproducibility and database transferability across different MALDI-TOF MS workflows. Inclusion of a broader range of non-target Aeromonas species and other phenotypically similar aquatic bacteria would further strengthen assessment of diagnostic specificity. Integration with whole-genome sequencing and multilocus sequence approaches may also help characterize strain diversity and determine whether spectral variation corresponds to genetically distinct lineages. In addition, prospective diagnostic studies comparing MALDI-TOF MS with conventional and sequencing-based methods in terms of accuracy, turnaround time, workload, and cost would help establish its practical value for routine aquaculture diagnostics.

CONCLUSION

This study developed and validated a species-specific MALDI-TOF MS MSP database for A. schubertii, addressing the limited representation of this emerging aquatic pathogen in the commercial reference database. The custom MSP database increased species-level identification from 78.95% with the commercial database to 100%, with all 19 tested isolates achieving highly reliable log scores of 2.49–2.68. High reproducibility of technical replicate spectra, with Pearson correlation coefficients of 0.95–0.97, further supported the consistency of the generated protein fingerprints. PCA and dendrogram analyses clearly differentiated A. schubertii from the four evaluated non-target Aeromonas species. Molecular characterization based on 16S rRNA and gyrB sequences, together with phylogenetic analyses, further supported the taxonomic identity of the isolates used to construct the MSP.

The developed MSP database provides a practical approach for improving MALDI-TOF MS-based identification of A. schubertii following bacterial isolation and may strengthen laboratory diagnosis and surveillance of white nodular disease in freshwater aquaculture. A major strength of the study was the integration of proteomic profiling with molecular confirmation and multivariate spectral analyses. In addition, evaluation of the database using 19 A. schubertii isolates from multiple fish hosts, and inclusion of four closely related non-target Aeromonas species provided evidence of identification performance and spectral discrimination within the evaluated isolate collection. Expanding MALDI-TOF MS reference libraries also provides a framework for improving identification of other aquatic bacterial pathogens underrepresented in commercial databases.

The species-specific MSP database substantially improved MALDI-TOF MS identification of A. schubertii within the evaluated isolate collection and provided clear spectral differentiation from the tested non-target Aeromonas species. These findings support targeted enrichment of MALDI-TOF MS reference databases as a promising strategy for improving the identification of underrepresented aquatic bacterial pathogens. With broader isolate representation and independent multicenter validation, the developed database could contribute to more rapid and reliable diagnosis and surveillance of A. schubertii-associated disease in freshwater aquaculture.

DATA AVAILABILITY

The custom MSP database and associated MALDI-TOF MS spectral data are available from the corresponding authors upon reasonable request for research and validation purposes.

GENERATIVE AI DECLARATION

During the preparation of this manuscript, the authors used generative artificial intelligence tools solely for language editing, grammar correction, and improvement of clarity and readability. The authors critically reviewed and revised all AI-assisted content and take full responsibility for the accuracy, integrity, and originality of the final manuscript. No generative AI tools were used to generate, analyze, or interpret the scientific data or results.

AUTHORS’ CONTRIBUTIONS

MFM: Conceptualization, investigation, methodology, formal analysis, and writing – original draft. HB: Investigation. MM: Investigation and writing – review and editing. SFA: Investigation and writing – review and editing. RHS: Investigation and writing – review and editing. CR: Conceptualization, supervision, methodology validation, funding acquisition, and writing – review and editing. All authors have read and approved the final manuscript.

COMPETING INTERESTS

The authors declare that they have no competing interests.

PUBLISHER’S NOTE

Veterinary World remains neutral with regard to jurisdictional claims in the published institutional affiliations.

ACKNOWLEDGMENTS

The authors sincerely thank Assoc. Prof. Dr. Nuvee Prapasarakul, Diagnosis and Monitoring of Animal Pathogen Research Unit, Faculty of Veterinary Science, Chulalongkorn University, Bangkok, Thailand, and the staff of the Department of Veterinary Microbiology, Faculty of Veterinary Science, Chulalongkorn University, for their valuable technical assistance and support throughout the study. The Researchers would like to thank the Deanship of Graduate Studies and Scientific Research at Qassim University for financial support. This study was funded by the Thailand Science Research and Innovation Fund and Chulalongkorn University (FOOD_FF_68_201_3100_017), the 90th Anniversary of Chulalongkorn University Fund (Ratchadaphiseksomphot Endowment Fund), and a postdoctoral fellowship awarded to M.M. by the Ratchadapisek Somphot Fund, Graduate School, Chulalongkorn University, Thailand.

REFERENCES

1. Li T, Raza SHA, Yang B, Sun Y, Wang G, Sun W. Aeromonas veronii infection in commercial freshwater fish: a potential threat to public health. Animals 2020;10(4):608.. [Google Scholar]
2. Budiyansah H, Mursalim MF, Chokmangmeepisarn P, Rung-ruangkijkrai T, Mabrok M, Khan MIR. Unveiling the threat: Aeromonas schubertii, an emerging bacterial pathogen inducing white nodule lesions in snakehead fish (Channa striata) at a nursery farm in Thailand. Aquac Int 2025;33(6):502.. [Google Scholar]
3. Trullàs C, Kitiyodom S, Pinpimai K, Chokmangmeepisarn P, Kamble MT, Rodkhum C. Pathogenicity and virulence of Aeromonas schubertii in farmed Asian seabass (Lates calcarifer) in Thailand. BMC Vet Res 2025;21(1):493.. [Google Scholar]
4. Sangpo P, Thitamadee S, Dong HT, Senapin S. Aeromonas schubertii, a novel bacterium recovered from an AHPND-affected farm, is lethal to whiteleg shrimp, Penaeus vannamei. Microb Pathog 2020;149:1-8.. [Google Scholar]
5. Latif-Eugenín F, Beaz-Hidalgo R, Figueras MJ. First record of the rare species Aeromonas schubertii from mussels: phenotypic and genetic reevaluation of the species and a review of the literature. Arch Microbiol 2016;198(4):333-346.. [Google Scholar]
6. Ren Z, Cai Y, Wang S, Liu S, Li A, Xiong Y. First case of Aeromonas schubertii infection in brackish water wild Nile tilapia, Oreochromis niloticus, in China. Aquaculture 2019;501:247-254.. [Google Scholar]
7. Tahir NDM, Matori MF, Gan HM. First report of Aeromonas schubertii infection in striped snakehead Channa striata (Bloch, 1793) fingerlings in Malaysia. J Med Vet 2024;7(1):33-40.. [Google Scholar]
8. Pham Thi HH, Kim DH, Quach Van CT, Nguyen PT, Nguyen TL. Prevalence and antibiotic resistance of Aeromonas schubertii causing internal white spot disease on snakehead fish, Channa striata, in the Mekong Delta, Vietnam. J World Aquac Soc 2023;54(5):1260-1276.. [Google Scholar]
9. Mursalim MF, Budiyansah H, Raharjo HM, Debnath PP, Sakulworakan R, Chokmangmeepisarn P. Diversity and antimicrobial susceptibility profiles of Aeromonas spp. isolated from diseased freshwater fishes in Thailand. J Fish Dis 2022;45(8):1149-1163.. [Google Scholar]
10. Tao LT, Gao DP, Liu YL, Sun WW, Shan XF. The genus Aeromonas in aquaculture: a comprehensive review of prevalence, virulence, and antibiotic resistance with an emphasis on key pathogenic species. Rev Aquac 2026;18(1):e70097.. [Google Scholar]
11. Vieira T, Costa GL, Yamada AY, Bertani AMJ, Medeiros MIC, Moreno LZ. Genomic and proteomic identification of Aeromonas isolates from diverse sources: comparative performance of WGS and MALDI-TOF MS in species-level resolution. Front Bacteriol 2025;4:1708292.. [Google Scholar]
12. Fono-Tamo EUK, Kamika I, Dewar JB, Lekota KE. Comparative genomics revealed a potential threat of Aeromonas rivipollensis G87 strain and its antibiotic resistance. 2023. [Google Scholar]
13. Sakulworakan R, Chokmangmeepisarn P, Dinh-Hung N, Sivaramasamy E, Hirono I, Chuanchuen R. Insight into the whole genome of Aeromonas veronii isolated from freshwater fish by resistome analysis reveals antibiotic-resistant traits. Front Microbiol 2021;12:733668.. [Google Scholar]
14. Luo X, Liao G, Fu X, Liang H, Niu Y, Lin Q. A novel and effective therapeutic method for treating Aeromonas schubertii infection in Channa maculata. Animals 2024;14(6):957.. [Google Scholar]
15. Badenas-Alzugaray D, Tristancho-Baró A, García-Lechuz JM, Burillo-Navarrete N, Sanz-Sanz S, Milagro-Beamonte AM. Identification and genomic characterization of Aeromonas dhakensis from a human sample. Acta Microbiol Hell 2025;70(3):34.. [Google Scholar]
16. Liu C, Guo YM, Cao JZ, Zhang DF, Chang OQ, Li K. Detection and quantification of Aeromonas schubertii in Channa maculata by TaqMan MGB probe fluorescence real-time quantitative PCR. J Fish Dis 2019;42(1):109-117.. [Google Scholar]
17. Recio G, Fernández-Bravo A, Latif-Eugenin F, Tena D, Rezusta A, Figueras M. Identification of 273 clinical Aeromonas strains with a gold standard method and MALDI-TOF: a review on the limitations of the method. In: Aeromonas—a ubiquitous bacterium that infects humans and animals. London: IntechOpen; 2024. [Google Scholar]
18. Oviaño M, Rodríguez-Sánchez B. MALDI-TOF mass spectrometry in the 21st-century clinical microbiology laboratory. Enferm Infecc Microbiol Clin 2021;39(4):192-200.. [Google Scholar]
19. Kitagawa D, Suzuki Y, Abe N, Ui K, Suzuki K, Yamashita T. Comparison of MALDI-TOF mass spectrometry and rpoB gene sequencing for the identification of clinical isolates of Aeromonas spp. Heliyon 2022;8(11):e11585.. [Google Scholar]
20. Mursalim MF, Raharjo HM, Budiyansah H, Chokmangmeepisarn P, Mabrok M, Yindee J. Development and validation of a custom matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) peptide database for the rapid and accurate identification of Aeromonas spp. isolated from diseased fish. BMC Vet Res 2026;22:157.. [Google Scholar]
21. Shin HB, Yoon J, Lee Y, Kim MS, Lee K. Comparison of MALDI-TOF MS, housekeeping gene sequencing, and 16S rRNA gene sequencing for identification of Aeromonas clinical isolates. Yonsei Med J 2015;56(2):550-555.. [Google Scholar]
22. Dhanapala PM, Kalupahana RS, Kalupahana AW, Wijesekera DPH, Kottawatta SA, Jayasekera NK. Characterization and antimicrobial resistance of environmental and clinical Aeromonas species isolated from freshwater ornamental fish and associated farming environment in Sri Lanka. 2021. [Google Scholar]
23. Shin SK, Lee Y, Kwon H, Rhee JS, Kim JK. Validation of the direct boiling method for simple and efficient genomic DNA extraction and PCR-based macroalgal species determination. J Phycol 2021;57(4):1368-1372.. [Google Scholar]
24. Hoel S, Vadstein O, Jakobsen AN. Species distribution and prevalence of putative virulence factors in mesophilic Aeromonas spp. isolated from fresh retail sushi. Front Microbiol 2017;8:931.. [Google Scholar]
25. Janda JM, Abbott SL. 16S rRNA gene sequencing for bacterial identification in the diagnostic laboratory: pluses, perils, and pitfalls. J Clin Microbiol 2007;45(9):2761-2764.. [Google Scholar]
26. Beaz-Hidalgo R, Hossain MJ, Liles MR, Figueras MJ. Strategies to avoid wrongly labelled genomes, using as an example the detected wrong taxonomic affiliation for Aeromonas genomes in the GenBank database. PLoS One 2015;10(1):e0115813.. [Google Scholar]
27. Pérez-Sancho M, Cerda I, Fernández-Bravo A, Domínguez L, Figueras MJ, Fernández-Garayzabal JF. Limited performance of MALDI-TOF for identification of fish Aeromonas isolates at species level. J Fish Dis 2018;41(10):1485-1493.. [Google Scholar]
28. Nellessen CM, Nehl DB. An easy adjustment of instrument settings (‘Peak MALDI’) improves identification of organisms by MALDI-TOF mass spectrometry. Sci Rep 2023;13(1):15018.. [Google Scholar]
29. Marczyk M, Polanska J. Methods for quality control of low-resolution MALDI-ToF spectra. Volume 3. Angers: SCITEPRESS; 2014. p. 172-177. [Google Scholar]
30. Murugaiyan J, Lewin A, Kamal E, Bakuła Z, van Ingen J, Ulmann V. MALDI spectra database for rapid discrimination and subtyping of Mycobacterium kansasii. Front Microbiol 2018;9:587.. [Google Scholar]
31. Piamsomboon P, Jaresitthikunchai J, Hung TQ, Roytrakul S, Wongtavatchai J. Identification of bacterial pathogens in cultured fish with a custom peptide database constructed by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS). BMC Vet Res 2020;16(1):52.. [Google Scholar]
32. Benagli C, Demarta A, Caminada A, Ziegler D, Petrini O, Tonolla M. A rapid MALDI-TOF MS identification database at genospecies level for clinical and environmental Aeromonas strains. PLoS One 2012;7(10):e48441.. [Google Scholar]
33. Fernández-Bravo A, Figueras JM. An update on the genus Aeromonas: taxonomy, epidemiology, and pathogenicity. Microorganisms 2020;8:1-39.. [Google Scholar]
34. Alatoom AA, Cunningham SA, Ihde SM, Mandrekar J, Patel R. Comparison of direct colony method versus extraction method for identification of Gram-positive cocci by use of Bruker Biotyper matrix-assisted laser desorption ionization-time of flight mass spectrometry. J Clin Microbiol 2011;49(8):2868-2873.. [Google Scholar]
35. Matsuda N, Matsuda M, Notake S, Yokokawa H, Kawamura Y, Hiramatsu K. Evaluation of a simple protein extraction method for species identification of clinically relevant staphylococci by matrix-assisted laser desorption ionization-time of flight mass spectrometry. J Clin Microbiol 2012;50(12):3862-3866.. [Google Scholar]
36. Kim SW, Nho SW, Im SP, Lee JS, Jung JW, Lazarte JM. Rapid MALDI biotyper-based identification and cluster analysis of Streptococcus iniae. J Microbiol 2017;55(4):260-266.. [Google Scholar]
37. Niyompanich S, Srisanga K, Jaresitthikunchai J, Roytrakul S, Tungpradabkul S. Utilization of whole-cell MALDI-TOF mass spectrometry to differentiate Burkholderia pseudomallei wild-type and constructed mutants. PLoS One 2015;10(12):e0144128.. [Google Scholar]
38. McLean K, Palarea-Albaladejo J, Currie CG, Imrie LHJ, Manson EDT, Fraser-Pitt D. Rapid and robust analytical protocol for E. coli STEC bacteria subspecies differentiation using whole-cell MALDI mass spectrometry. Talanta 2018;182:162-170.. [Google Scholar]
39. Rothen J, Pothier JF, Foucault F, Blom J, Nanayakkara D, Li C. Subspecies typing of Streptococcus agalactiae based on ribosomal subunit protein mass variation by MALDI-TOF MS. Front Microbiol 2019;10:471.. [Google Scholar]
40. Yáñez MA, Catalán V, Apráiz D, Figueras MJ, Martínez-Murcia AJ. Phylogenetic analysis of members of the genus Aeromonas based on gyrB gene sequences. Int J Syst Evol Microbiol 2003;53(3):875-883.. [Google Scholar]
41. Beaz-Hidalgo R, Alperi A, Bujan N, Romalde JL, Figueras M. Comparison of phenotypical and genetic identification of Aeromonas strains isolated from diseased fish. Syst Appl Microbiol 2010;33:149-153.. [Google Scholar]
42. Martínez-Murcia AJ, Monera A, Saavedra MJ, Oncina R, López-Álvarez M, Lara E. Multilocus phylogenetic analysis of the genus Aeromonas. Syst Appl Microbiol 2011;34:189-199.. [Google Scholar]
43. Soler L, Yáñez MA, Chacón MR, Aguilera-Arreola MG, Catalán V, Figueras MJ. Phylogenetic analysis of the genus Aeromonas based on two housekeeping genes. Int J Syst Evol Microbiol 2004;54(5):1511-1519.. [Google Scholar]