Next Article in Journal
Antibiotic Prophylaxis for the Prevention of Urinary Tract Infections in Children: Guideline and Recommendations from the Emilia-Romagna Pediatric Urinary Tract Infections (UTI-Ped-ER) Study Group
Previous Article in Journal
Influence of Multiplex PCR in the Management of Antibiotic Treatment in Patients with Bacteremia
Previous Article in Special Issue
Comparative Genomics of Pseudomonas aeruginosa Strains Isolated from Different Ecological Niches
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genomic Characterization of Aeromonas veronii Provides Insights into Taxonomic Assignment and Reveals Widespread Virulence and Resistance Genes throughout the World

by
José Cleves da Silva Maia
1,2,
Gabriel Amorim de Albuquerque Silva
2,
Letícia Stheffany de Barros Cunha
1,2,
Gisele Veneroni Gouveia
2,
Aristóteles Góes-Neto
3,
Bertram Brenig
4,
Fabrício Almeida Araújo
5,
Flávia Aburjaile
6,
Rommel Thiago Jucá Ramos
5,
Siomar Castro Soares
7,
Vasco Ariston de Carvalho Azevedo
8,
Mateus Matiuzzi da Costa
2,* and
João José de Simoni Gouveia
2,*
1
Graduate Program in Animal Science, Agricultural Sciences Campus, Federal University of Vale of São Francisco (Univasf), Petrolina 56304-917, Pernambuco, Brazil
2
Center for Open Access Genomic Analysis (CALAnGO), Federal University of Vale of São Francisco (Univasf), Petrolina 56304-917, Pernambuco, Brazil
3
Laboratory of Molecular Computational Biology of Fungi (LBMCF), Institute of Biological Sciences, Federal University of Minas Gerais, Belo Horizonte 31270-901, Minas Gerais, Brazil
4
Institute of Veterinary Medicine, University of Göttingen, 37077 Göttingen, Niedersachsen, Germany
5
Biological Engineering Laboratory, Institute of Biological Sciences, Federal University of Pará (UFPA), Belém 66075-110, Pará, Brazil
6
Preventive Veterinary Medicine Department, Veterinary School, Federal University of Minas Gerais, Belo Horizonte 31270-901, Minas Gerais, Brazil
7
Department of Microbiology, Immunology, and Parasitology, Federal University of Triângulo Mineiro, Uberaba 38025-180, Minas Gerais, Brazil
8
Laboratory of Cellular and Molecular Genetics (LGCM), Institute of Biological Sciences, Federal University of Minas Gerais, Belo Horizonte 31270-901, Minas Gerais, Brazil
*
Authors to whom correspondence should be addressed.
Antibiotics 2023, 12(6), 1039; https://doi.org/10.3390/antibiotics12061039
Submission received: 31 March 2023 / Revised: 23 May 2023 / Accepted: 9 June 2023 / Published: 11 June 2023

Abstract

:
Aeromonas veronii is a Gram-negative bacterial species that causes disease in fish and is nowadays increasingly recurrent in enteric infections of humans. This study was performed to characterize newly sequenced isolates by comparing them with complete genomes deposited at the NCBI (National Center for Biotechnology Information). Nine isolates from fish, environments, and humans from the São Francisco Valley (Petrolina, Pernambuco, Brazil) were sequenced and compared with complete genomes available in public databases to gain insight into taxonomic assignment and to better understand virulence and resistance profiles of this species within the One Health context. One local genome and four NCBI genomes were misidentified as A. veronii. A total of 239 virulence genes were identified in the local genomes, with most encoding adhesion, motility, and secretion systems. In total, 60 genes involved with resistance to 22 classes of antibiotics were identified in the genomes, including mcr-7 and cphA. The results suggest that the use of methods such as ANI is essential to avoid misclassification of the genomes. The virulence content of A. veronii from local isolates is similar to those complete genomes deposited at the NCBI. Genes encoding colistin resistance are widespread in the species, requiring greater attention for surveillance systems.

1. Introduction

The genus Aeromonas comprises more than 30 species of Gram-negative bacteria commonly associated with the aquatic environment and is known to infect and cause mortality in fish [1,2]. Among the most important species of the genus is Aeromonas veronii, a mobile mesophilic bacterium that causes furunculosis and motile aeromonas septicemia (MAS) in fish, affecting commercial fish farms and causing economic losses to aquaculture [3,4]. In addition, the role of this species as a pathogen capable of causing intestinal diseases in humans has been increasingly highlighted [5,6].
Several components contribute to the pathogenicity of A. veronii, from lateral and polar flagella, more than one subtype of type IV pili and non-filamentous adhesins, to secretion systems that release toxins and exoenzymes into the environment or deliver directly into target cells [7,8]. These virulence factors help the pathogen during the processes of colonization, adhesion, and destruction of host cells, thus overcoming the immune system [2]. Among them, type I and type II secretion systems (T1SS and T2SS), type IV pili (T4P), and polar flagellum are conserved within the species; whereas type III, type V, and type VI secretion systems (T3SS, T5SS, and T6SS) show variable dispersion among isolates [8]. In addition, important virulent strains may have more than one copy of T3SS, as well as some toxins and iron uptake systems [9].
Species of the genus Aeromonas are already naturally resistant to some beta-lactams [1]. In the A. veronii resistome, antibiotic resistance genes (ARGs) can be grouped according to the mechanism of resistance into at least four categories: antibiotic efflux, antibiotic inactivation, antibiotic target alteration, and antibiotic target replacement [8]. Through horizontal gene transfer (HGT), these bacteria could acquire and disseminate ARGs, especially in the aquatic environment, from where multidrug-resistant strains may emerge [10,11,12]. This situation may involve resistance to antibiotics that are considered the last line of defense against resistant Gram-negative bacteria, such as resistance to colistin encoded by mcr (mobile colistin resistance) genes, since these are increasingly disseminated in the aquatic environment, with probable participation from Aeromonas spp. in the process [10,13,14,15].
Although it is known as one of the main pathogens responsible for cases of MAS in fish [4], A. veronii is amongst the main Aeromonas species associated with infections in humans, being isolated mainly from feces, blood, and wounds [2]. In some cases, this species is the main bacteria isolated from patients with gastroenteritis, and the frequency of infections increases in warmer seasons and in older hosts [5]. In this sense, population aging and the ubiquity of Aeromonas spp. may facilitate a future scenario in which infections by these bacteria become a major health problem [2].
As an emerging pathogen in humans [2,5,6,16], and as likely sources for ARGs [12,13,14], investigating the epidemiology of A. veronii becomes increasingly important; however, care must be taken when identifying this and other species of the genus isolated from animals, humans, and the environment, since classification errors in Aeromonas spp. are common when using phenotypic methods [2,17,18,19]. Regarding this, genome sequencing and the use of in silico tools to classify species are fundamental [17,18], considering that this also enables for the characterization of virulence profile and allows for comparative analyses with more virulent strains isolated from different hosts, leading to the identification of essential components behind the pathogenicity [6,8,9]. Another essential analysis is to annotate and identify ARGs from local isolates and compare them with those present in genomes from other countries to monitor the global spread of resistance to important antibiotics, acting under the context of a One Health approach [20].
Taking into account the importance of genomic sequencing to characterize, classify, and monitor the resistance and virulence of pathogens, this study aims to perform genomic analyses of bacterial isolates from the São Francisco Valley (Petrolina, Pernambuco, Brazil) classified as A. veronii and from complete genomes available in public databases to gain insight into the taxonomic assignment of the species A. veronii and to better understand the virulence and resistance profiles of this species in context of One Health.

2. Results

2.1. MLST, TETRA, ANI, and Phylogeny

MLST typing (Multi-Locus Sequence Typing) reveals that all genomes belong to the genus Aeromonas, although only seven NCBI complete genomes have had their sequence type recognized (485, 333, 166, 515, 1000, 160, and 254). All sets of alleles from the complete genomes differ at least in one of the six genes used in typing (gyrB, groL, gltA, metG, ppsA, and recA). Regarding the genomes from the SFV (São Francisco Valley), two sets of alleles repeat twice in each isolate, suggesting they can be clonal: B79 and B80; B44 and B52 (Table S1).
The correlation index of tetranucleotides (TETRAs) clusters together MM_B02 with A. dhakensis and A. hydrophila (TETRA = 1.000 and 0.999, respectively) (Figure S1). All the other SFV (São Francisco Valley) isolates constitute a clear cluster (TETRA ≧ 0.994) encompassing A. veronii and other species (A. fluvialis, A. lacus, A. jandaei, A. australiensis, A. allosaccharophila, and A. finlandensis) (Figure S2). Despite the fact that TETRA analysis seems to not provide a sufficient resolution for unequivocally separating A. veronii from other species (such as A. allosaccharophila and A. filandensis), it can be seen that, considering the threshold of 0.999, the genomes WP2-S18-CRE-03, WP3-W19-ESBL-03, WP8-S18- ESBL-11, WP9-W18-ESBL-04, FC951, HD6448, and MM_B02 (SFV) could be suggested to not belong to the A. veronii species when compared with the reference FDAARGOS-632.
To overcome this, another classification method was performed to estimate the average nucleotide identity (ANI) between paired genomes in four ways: ANIb (BLAST), ANIm (MUMMER), OrthoANIb (BLAST), and OrthoANIu (USEARCH). Considering the threshold of 95%, 33/38 previously classified A. veronii genomes could be suggested to belong to the A. veronii species, and the genomes WP2-S18-CRE-03, WP3-W19-ESBL-03, WP8-S18-ESBL-11, WP9-W18-ESBL-04, and MM_B02 can be suggested to belong to another species. It is interesting to note that using this threshold (95%), all four ANI approaches used are completely concordant (Figure S2). Considering the threshold of 96%, 31/38 previously classified A. veronii genomes could be suggested to belong to A. veronii species. In addition to the aforementioned genomes, 1708-29120 did not reach the 96% threshold in all four ANI analysis, and 71506 genome did not reached the 96% threshold in the OrthoANIU analysis.
The local genome MM_B02 clusters with A. dhakensis KOR1, exhibiting average identity of nucleotides over 0.95 (Figure 1). The other local isolates constitute a cluster which includes the reference A. veronii genome FDAARGOS-632 (ANI > 0.95), and other deposited genomes of the species (Figure 1).
WP3-W19-ESBL-03, WP8-S18-ESBL-11, and WP9-W18-ESBL-04 genomes constitute a cluster with A. allosaccharophila FDAARGOS-933 (ANI > 0.95). The WP2-S18-CRE-03 genome constitutes a cluster with the A. jandaei Riv2 genome (ANI > 0.96) and A. lacus AE122 genome (ANI > 0.95), being related with the first one considering the ANI value. The other isolates used in this study (n = 25) constitute a cluster with the A. veroni reference genome FDAARGOS-632 (ANI > 0.95) (Figure 1). In this way, it can be suggested that the isolates WP3-W19-ESBL-03, WP9-W18-ESBL-04, WP8-S18-ESBL-11, and WP2-S18-CRE-03 do not belong to the A. veronii species, suggesting that they are misclassified in the NCBI database. With A. veronii, 71506 is the most divergent isolate from the others and the only isolate that presents ANI values lower than 96% (0.96) with the four methods being applied.
The isolates B79 and B80 are suggested to be clonal (ANI = 100%), as are the isolates B44 and B52. With the former, the samples were collected in different city points. With the latter, the samples were isolated from different species of hosts in the same location, with B44 isolated from Pseudoplatystoma corruscans and B52 from Phractocephalus hemioliopterus.
A maximum likelihood tree from 841 genes of the core genome was built (Figure 2). It can be observed that the genome MM_B02 from SFV, as with the genomes from NCBI WP2-S18-CRE-03, WP3-W19-ESBL-03, WP8-S18-ESBL-11, and WP9-W18-ESBL-04, is more related to the genomes of A. dhakensis, A. jandaei, and A. allosaccharophila. Inside the branch that includes only the A. veronii genomes, the 71506 isolate is the one that diverges most from the others. From this perspective, the results from phylogenic study are in accord with the ANI results that use 95% as the threshold for species classification. It is important to note that isolates from SFV, different from what was expected, do not constitute an only phylogenetic subgroup, and some of them are even closely related with isolates from China, as, for example, B42 and ZfB1; B48 and 183026; and B45 and GD21SC2322TT (Figure S3). The same can be noted about the isolates from other countries (Figure 2 and Figure S3).

2.2. Virulence

A total of 256 virulence genes were identified, and from these, 239 were present in the local VSF isolates. The number of genes common to all analyzed genomes was 236 (Figure S4), with seventeen exclusive to the complete A. veronii genomes deposited and three exclusive to local SFV genomes, including hlyB, which encodes Hemolysin B.
A total of 190 genes were shared among all isolates regardless of their source of isolation, whether from humans, fish, the environment, or other (Figure S5). In the fish group there were ten exclusive genes, while in the human group there were four exclusive genes. The other groups did not present exclusive genes. Fish isolates are of a bigger quantity in this study and present more exclusive genes than human isolates, including the ones which are involved with secretion systems and exotoxins.
A total of 224 of the 256 virulence genes identified involve functions of adherence, effector delivery systems, and motility, accounting for 59, 78 and 87 genes, respectively (Figure 3). Among the structures behind these functions, those that are present in all isolates are: type IV pili (Tap), polar flagellum, and the type II secretion system (T2SS) (Figure S6). In addition to the elements common to all isolates, present in the genomes from SFV and in some from NCBI are the lateral flagellum, type IV pili (Flp and Msha), type I pili, and type III and type VI secretion systems (T3SS and T6SS). Other genes are related to immune modulation, stress survival, metabolic and nutritional factors, exotoxins and exoenzymes, regulation, biofilm formation, and competitive advantage, among others.
A PCoA (Principal Coordinate Analysis) based on the presence and absence of virulence genes and their matrix reveals that the difference among the genomes is not explained by the type of host or local from which the specimen was collected (Figure S7a), with local SFV isolates showing a similar virulence content to isolates from China (Figure S7b). Interestingly, the single isolate from an asymptomatic host was collected from Hirudo verdana (leech) and presents a similar virulence content to those isolated from diseased fish (Figure S7c).
Plasmids were identified on B44, B52, and B45 local isolates (Table S5). One non-mobilizable plasmid was identified in the three isolates, whereas a conjugative plasmid of 144,051 bp was identified only in isolate B45. In this plasmid, three genes related to the type VI secretion system were identified: vgrG1, AHA_RS09305, and hcp. No genes related to virulence were identified in the 26 plasmids belonging to complete A. veronii genomes downloaded from the NCBI.

2.3. Resistance

There were 60 different genes conferring resistance to 22 classes of antimicrobials in the 33 A. veronii isolates analyzed. This resistance takes place through five mechanisms (Figure S8): antibiotic efflux (15 genes); antibiotic inactivation (27 genes); antibiotic target alteration (10 genes); antibiotic target replacement (6 genes); and antibiotic target protection (2 genes).
All isolates have genes associated with resistance to drugs from the following classes: Peptide, Carbapenem, Phenicol, Fluoroquinolone, Macrolide, Cephalosporin, Diaminopyrimidine, Penam, Penem, Monobactam, Sulfonamide, Cephamycin, Aminocoumarin, Nitroimidazole, and Tetracycline (Figure 4). Resistance to the last four classes is especially due to Efflux Pumps, which contain the largest number of genes in the genomes.
From all the antimicrobial resistance genes identified, 44 were exclusive to the genomes from NCBI, 1 exclusive to the local VSF genomes, and 15 genes were shared by the NCBI and local genomes (Figure S9a). The isolates were assigned to four categories according to the isolation source: human, fish, the environment, and other. A total of 15 genes were shared among them. Fifteen genes were present only in humans, four only in fish, four in the environment, and four in others (Figure S9b).
All analyzed A. veronii genomes harbor genes from the ABC transporters; RND efflux pump; Undercaprinlyl-Pyrophosphate-related proteins; CphA, OXA, and Dfr; and MCR (mcr.3 and mcr.7.1) (Figure S10).
Of the plasmids obtained with the genomes of A. veronii from the NCBI, four were conjugatives, 11 were mobilizable, and 11 were non-mobilizable (Table S6). In total, 29 genes in all were identified in 11 plasmids, with each conferring resistance to 13 classes of antimicrobials. They were grouped into four resistance mechanisms (Figure S11): antibiotic inactivation (27 genes in total), antibiotic target replacement (11 genes in total), antibiotic efflux (4 genes in total), and antibiotic target protection (4 genes in total). Among them, it is worth highlighting three conjugative plasmids: pASP-a58 from the lineage AVNIH1; p158496 of the SW3814 lineage, isolated from an environment in the United States; and pWP8-W19-CRE-03_1, from the WP8-W19-CRE-03 lineage, isolated from an environment in Japan. The first encodes 13 different genes conferring resistance to 10 classes of antimicrobials, the second encodes 8 different genes conferring resistance to 7 classes, and the last one has 5 genes that confer resistance to 5 classes. No gene related to antimicrobial resistance was identified in plasmids from the SFV.

3. Discussion

Nine genomes of SFV isolates were submitted to different tools for taxonomic classification. The tetranucleotide signature (TETRA) revealed that the SFV genome, MM_B02, along with six other NCBI genomes, may not be A. veronii when compared to the reference genome FDAARGOS-632. However, confusing results can be observed with this method, as some genomes are related to more than one species. As pointed out by Richter and Rosselló-Móra [21], TETRA shows irregularities when applied to fragmented genomes, being more inaccurate than other methods, such as ANI, for example.
The average nucleotide identity (ANI) showed greater resolution in genome separation, including for reference sequences from other species (Figure 1). The results of this analysis show that eight isolates from SFV were A. veronii, while one isolate (MM_B02, from human) was misidentified and belongs to the Aeromonas dhakensis species. What stands out is the fact that there is, also, amongst the whole genome data obtained from the NCBI, misidentified genomes were deposited, and they were WP2-S18-CRE-03, WP3-W19-ESBL-03, WP8-S18-ESBL-11, and WP9-W18-ESBL-04. To the best of our knowledge, only the first genome has been explicitly highlighted as erroneously identified in a previously published study [6].
Identification errors in Aeromonas species have already been pointed out in other studies [17,18]. Correct identification of pathogens is critical to understanding their epidemiology and their true impact on human, animal, and environmental health [17]. In this sense, the tools applied to calculate the ANI from genomes, whether complete or drafts, proved to be useful and consistent with the 95% threshold [17] for A. veronii (Figure S1), although 96% is suggested to be more suitable for other species of the genus Aeromonas [22,23].
Interestingly, clonal isolates among A. veronii from the SFV occurred twice (Figure 1), which is noteworthy because they were sampled at different points in the region, as was the case for B79 and B80, and for hosts of different species, as for B44 and B52. This situation reflects what can be observed when analyzing the general phylogeny of the genomes used in this study (Figure 2 and S3), as it is possible to notice that isolates from different countries, or different host species, may be more closely related to each other, for example, isolate B48 and 183026. A similar situation was highlighted by Liu and his collaborators [6], a study which revealed the ease with which A. veronii can be spread locally and globally, in addition to its ability to infect different types of hosts.
Most of the virulence genes are common to SFV and NCBI isolates (Figure 3), and they were mainly involved with adherence, secretion systems, and motility. Most of the genes involved with adherence in this species encode structures such as type IV pili (Tap, Msha, and Flp), which are essential during the host colonization step, especially Msha type IV pili, absent only in isolate FC951, which is necessary for A. veronii adherence and biofilm formation by the species [24]. Interestingly, genes encoding type I pili were also identified in nine isolates, including five from the SFV. This filamentous adhesin is one of the main virulence factors in uropathogenic Escherichia coli and is responsible for mediating attachment to a variety of host surfaces and tissues [25]. In Aeromonas salmonicida, type I pili contributes to the colonization of fish [26]. Most isolates containing this structure are from fish; one of each of the others are from humans, alligators, and the environment.
Most Aeromonas species, including A. veronii, can move via a single polar flagellum [1], thus, genes that encode this component are conserved in the species genomes, according to the results of the analyses. The lateral flagellum, on the other hand, is found with an approximate frequency of 50–60% in mesophilic species associated with diarrheal diseases [27]. Among the whole genomes from the NCBI used in this study, only 44% have a set with more than 30 genes encoding this structure, while all isolates from the SFV have such a set of genes. In addition to motility, the importance of this component for pathogenicity is also due to its participation in the stages of adhesion and colonization of the host’s intestine, invasion, movements on surfaces, and biofilm formation [27,28].
Bacteria use secretion systems to secrete proteins that enhance adherence to host cells, obtain resources from environmental niches, and even directly intoxicate target cells [29]. Aeromonas spp. utilize T1SS, T2SS, T3SS, and T6SS to transport virulence factors, toxins, and effectors during host fish infection [30]. In this work, genes coding for type II, type III, and type VI secretion systems were identified. T2SS is conserved in A. veronii [8,29] and is used to move substrates across the bacterial cell membrane, releasing them directly into the environment, with the substrates then participating in nutrient acquisition, biofilm formation, and overall pathogenicity [31]. The secreted products of this system include proteases, lipases, phosphatases, and toxins, among others [32]. Genes encoding aerolysin (aerA) and metalloprotease (stcE), both products secreted by T2SS [32], were found in all isolates from the SFV and almost all isolates from the NCBI. Aerolysin is a pore-forming toxin [33] and acts by disrupting the fish intestinal barrier, being a key virulence factor in A. veronii [4]. StcE metalloproteases are secreted by Enterohaemorragic-Escherichia coli (EHEC) and are related to reduced mucus levels, facilitating adherence to human hosts [34].
Differently from what happens in T2SS, T3SS and T6SS can both deliver effectors directly into target cells [29]. Gram-negative plants and animal pathogens use T3SS to inject virulence factors directly into the cells of their hosts [35]. T3SS is expressed in important human enteric pathogens, which translocate effectors to the interior of their hosts, which act on the cytoskeleton and cell traffic, in addition to important cellular pathways, aiding in the establishment of infections and blocking phagocytosis [36]. The effectors of this system commonly found in A. veronii are the AexT and AexU toxins [37]; however, the T3SS effectors identified in the present study were phosphatases encoded by the ati2 and aopH genes, with neither of them found in genomes from the SFV, despite the presence of T3SS in these isolates. Genes encoding T6SS are present in six out of eight genomes from the SFV. This secretion system exports effectors which act in both bacterial and eukaryotic cells, leading to the death of competing and/or host cells [38]. Other functions attributed to T6SS involve adaptation to stress conditions caused by temperature or pH changes, as well as the uptake of metal ions in response to oxidative stress [38,39].
Other virulence elements identified in A. veronii genomes include LPS (lipopolysaccharides), LOS (lipooligosaccharides), capsule, genes associated with adherence (Hsp60, IlpA, EF-Tu), genes related to iron uptake (entB, entE and Fur), and luxS. The luxS gene, absent only in the isolate B45 from the SFV, encodes an S-ribosilhomocysteinase that plays a key role in the synthesis of the AI-2 autoinducer [40]. This autoinducer plays an important role in interspecies communication, especially quorum sensing, thus influencing biofilm formation, virulence, stress resistance in different environments, and antibiotic resistance [40,41,42,43].
Genes encoding resistance to 22 classes of antibiotics were identified in the A. veronii genomes used in this study. Part of this resistance can be attributed to the presence of gene coding for beta-lactamases, since many Aeromonas species are naturally resistant to some beta-lactam antibiotics [1]. It is important to draw attention to the Extended Spectrum Beta-lactamases (ESBL), identified only in isolates GD21SC2322TT and 1708-29120, from China, which are responsible for the resistance of some pathogens to a wide spectrum of beta-lactam antibiotics, except for carbapenems [44]. However, this class of antibiotics can be inactivated by Metallo-Beta-Lactamases (MBLs) [45], such as those encoded by the CphA gene family, which is common to all A. veronii genomes investigated in this study (Figure S10).
Even more important is the colistin resistance, which is considered one of the last resources to treat infections caused by multidrug-resistant pathogens [13]. The mcr-3 and mcr-7 genes, which encode phosphoethanolamine transferase, were identified in this work, the last one being present in all isolates. The mcr-3 gene, present only in isolates from the NCBI, was characterized by Yin et al. in 2017, who pointed out its widespread prevalence among Enterobacteriaceae and Aeromonas [13]. The mcr-7 gene was first identified in Klebisiella pneumoniae from chickens and is closely related to mcr-3 from Aeromonas, suggesting a likely origin in this genus [14,15]. Apparently, Aeromonas spp. may be likely reservoirs for mcr genes, disseminating them through horizontal gene transfer in the most diverse aquatic environments in which they are present [13,15].
Several genes identified in the genomes analyzed can confer resistance to antibiotics such as Phenicol, Fluoroquinolone, Macrolide, Tetracycline, Diaminopyrimidine, Sulfonamide, Aminocoumarin, and Peptides, which is due, especially, to the fact that some of them code for antibiotic efflux systems. The ability to rapidly export drugs out of bacterial cells is an important contributing factor to antimicrobial resistance in pathogens [46]. Three families of efflux pumps can be identified in the isolates of this work, namely: MFS, RND, and ABC. The last two are present in all A. veronii genomes and may contribute to multidrug resistance in the species when they are overexpressed, especially the RND efflux pump, and they are considered the most clinically significant [46,47]. MFS transporters, although considered ubiquitous in all domains of life [46], did not have genes identified in all the genomes of this study, being present in only three of the eight isolates from the SFV, and they were B48, B79, and B80. This family tends to be a specific substrate [47], highlighting that genes tet(A) and tet(E), coding for the efflux of tetracyclines, and floR were involved with the efflux of Fenicol, both of which can be transferred horizontally through mobile genetic elements [48,49].
These efflux pumps contribute to prospective resistance, according to CARD information [50], to 18 classes of drugs. The mere presence of efflux pumps may not be enough to confer resistance, since resistance depends on other factors, such as expression [46,47]. However, resistance to some classes of antibiotics may occur due to more than one resistance mechanism encoded in the same genome, as occurs with beta-lactams in isolate GD21SC2322TT, for example.
The spread of antibiotic resistance involves the horizontal transfer of genes through mobile genetic elements, especially plasmids, via conjugation [51]. No genes related to antibiotic resistance were identified in plasmids isolates from the SFV; however, it is necessary to highlight some plasmids deposited with complete genomes from the NCBI. An example of this is the conjugative plasmid pASP-a58, from the environmental isolate AVNIH1, which encodes thirteen genes related to four resistance mechanisms (antibiotic inactivation, target replacement, efflux pump, and target protection) and confers resistance to ten antibiotics classes, including Carbapenem, Phenicol, and Aminoglycoside. Multiple resistance genes on a plasmid such as this help to spread multi-drug resistance more easily [51], especially in the case of conjugative plasmids, since they contain enough elements to carry out their own transfer [52].
Differently from Tekedar et al. [8] and Liu et al. [6], in the present study, T1SS, T4SS, and T5SS were not identified, as was the case in the study by Prediger et al. [7], although in this case a single isolate was characterized. On the other hand, type I pili was only identified in the present study and in the one by Prediger et al. [7]. This difference in the virulence elements found can be caused by several factors, such as tools and methods performed during annotation, database content, coverage percentage, and established identity, among others. This reveals the need to create tools and databases that are personalized with reference genes for some taxa that have gained importance as emerging pathogens, such as Aeromonas [2,4,5,6,16].
The A. veronii species, like other species in the genus, can adapt very well to the selective pressures of antibiotics use, since its pan-genome is open and exhibits high genomic plasticity, being influenced by mobile genetic elements through the horizontal transfer of genes [8,10,16,53]. This genome plasticity also influences, among other aspects, the virulence of Aeromonas spp. [53,54].
Aeromonas veronii is known to infect fish, cause mortality, and negatively impact aquaculture [2,3,4]; however, more attention should be given to the species, as it has increasingly emerged as a pathogen associated with cases of intestinal diseases in humans [5,6,55]. According to the analyses of the present study, it is possible to verify that fish isolates, such as B48, MS-18-37, and X11, are closely related to human isolates (Figure S3). Furthermore, the virulence gene content of isolates from different isolation sources, including diseased humans, are similar (Figure S7). Given its characteristics [1,2], the species has the potential to infect humans from food sources, especially fish [56], and from water used for consumption [57]. A. veronii is a mesophilic species of the genus Aeromonas [2] able to proliferate at higher temperatures, and some studies, such as the one by Yuwono et al. [5], reported a higher frequency of Aeromonas infections in warmer seasons.
In addition to its pathogenic potential, A. veronii may play an important role as a source of antibiotic resistance, maintaining and disseminating antimicrobial resistance genes in the aquatic environment, including mcr genes [10,15,51,58,59]. It is important to note that fish and human isolates practically share resistance to the same classes of antibiotics (Figure 4), and human isolates show an even greater amount of resistance genes, even though they are fewer in number (Figure S9). When thinking about surveillance systems for antimicrobial resistance, which is considered one of the main threats to public health, it is necessary to include data from human and animal pathogens, as well as environmental isolates, following a One Health approach [20]. In this sense, it is also necessary to include emerging species as potential threats to public health; these species already play an important role in the infection of animals or are disseminated in the environment, as with A. veronii.
Finally, isolates from fish native to the São Francisco Valley (SFV) showed a similar virulence profile to fish isolates from other parts of the world, mainly those from China. The antimicrobial resistance content was also similar, although lower, to the content identified for complete genomes of A. veronii deposited at the NCBI, indicating that antimicrobial resistance genes, or resistant isolates, are now widely disseminated throughout the world. As for already deposited genomes, there was also an identification error, which was corrected with genomic sequencing and the application of tools that classify by comparing the entire genome, especially the ANI method, revealing flaws in phenotypic classification systems. In general, the gene content related to virulence and resistance to antibiotics, as well as the ability to infect hosts other than fish, makes A. veronii a species with zoonotic potential that should be on the radar of surveillance systems, especially those based on the One Health initiative.

4. Materials and Methods

4.1. Genomes from the São Francisco Valley (SFV)

Nine isolates recovered from fish, humans, and the environment (Table 1) were classified as Aeromonas veronii based on the automated PhoenixTM system (Becton Dickinson, Franklin Lakes, USA). The fish and environmental isolates were collected from 2017–2019, and the human isolate was collected in 2021. All isolates were maintained in a BHI broth added with glycerol in −80C until analysis. DNA was isolated using a salting-out protocol [60]. In summary, 500 µL of lysis solution (10 mM Tris-HCl pH 8; 10 mM EDTA pH 8; 100 mM NaCl; 0.5% SDS; and 2 µg of proteinase K) was added, followed by vortexing and incubation in a dry bath for 4–6 h at 55 °C. Then, 210 µL of TE (pH 7.6) and 240 µL of NaCl 5M were added, with ice incubation for 10 min. Samples were centrifuged for 15 min at (16,000× g), and 750 µL of supernatant was collected and added to new microtubes. In total, 750 µL of isopropyl alcohol (100%) was added, and samples were centrifuged (15min; 16,000× g) and the supernatant then discarded. In total, 500 µL of ethanol (70%) was added to the samples, followed by incubation (5 min) and centrifugation (5 min; 16,000× g). The supernatant was discarded, and 100 µL of H2O + RNAse (10 µg per mL) was added to the samples, which was subsequently incubated at 37 °C for 1 h.
DNA obtained was submitted to sequencing using NEBNext1 Fast DNA Fragmentation and Library Preparation Kit (New England Biolabs Inc., Ipswich, MA, USA) following the manufacturer’s protocol. Whole-genome sequencing was performed using an Illumina HiSeq 2500 platform (Illumina Inc., San Diego, CA, USA) with coverage above 100×.

4.2. NCBI Genomes

All available complete A. veronii genomes (n = 29) were downloaded along with their plasmids from the NCBI (Table S2). Additionally, thirty reference genomes from other Aeromonas species were used for taxonomy, and six randomly selected complete genomes (three A. dhakensis and three A. jandaei) were used for phylogeny analysis (Tables S3 and S4).

4.3. São Francisco Valley Genome Assembly

Raw data were submitted to an assembly pipeline that combines different approaches (Available online: https://github.com/engbiopct/assembly-hiseq, accessed on 23 April 2022). In summary, FastQC [61] was used in to assess data quality, AdapterRemoval v.2 [62] to remove the adapters, and KmerStream [63] to estimate better k-mer values within the proposed size (7–127, two-by-two). The genome assembly was performed with SPAdes [64], Edena [65], and Unicycler [66], and CD-HIT [67] was used to reduce computational costs and facilitate the assembly of sequences. The best assembly was chosen considering the metrics assessed using QUAST (Table S7) [68].

4.4. Plasmids

Mob-suite [69] was used to identify potential plasmids in the assembled genomes and separate them from chromosome sequences, with the MOB-recon tool used to reconstruct and perform typing of genomes sequenced in this study, and MOB-typer used to perform sequence typing for plasmids downloaded from the NCBI. The Plasmidfinder tool from Abricate version 1.0.1 (Available online: https://github.com/tseemann/abricate, accessed on 12 July 2022) was used to search for replicons.

4.5. Annotation

Prokka version 1.14.6 [70] was used with the genetic code 11 (—gcode 11) and default values for the other under-standard parameters to perform the annotation of possible features and CDSs for genomes.

4.6. Typing and Genomes Classification

The mlst tool (Available online: https://github.com/tseemann/mlst, accessed on 16 May 2022) was used to perform genome typing using multi-locus sequence typing (MLST) following the pubMLST scheme (Available online: https://pubmlst.org/bigsdb?db=pubmlst_aeromonas_seqdef, accessed on 16 May 2022). To calculate the average nucleotide identity (ANI) among the genomes, the pyANI and [71] (Available online: https://github.com/widdowquinn/pyani, accessed on 21 May 2022) and OrthoANI tools were utilized [72] (Available online: https://www.ezbiocloud.net/tools/orthoani, accessed on 2 June 2022 and Available online: https://www.ezbiocloud.net/tools/orthoaniu, accessed on 6 June 2022). The first one, pyANI version 0.2.11, estimated ANI using BLAST (ANIb) and Mummer (ANIm) in the scheme all vs. all. OrthoANIu version 1.2 was used to estimate ANI applying Usearch with the scheme all vs. all, meanwhile OrthoANIb version 1.40 was applied using BLAST with the scheme all vs. all. The reference species was A. veronii genome FDAARGOS-632 (Available online: https://www.ncbi.nlm.nih.gov/assembly/GCF_008693705.1/, accessed on 2 May 2022). The pyANI tool was used to calculate the correlation indexes of tetranucleotides signatures (TETRAs) among the genomes. Genomes were considered clonal in this study if one of the following criteria was met: (1) Present the same alleles for all genes (gyrB, groL, gltA, metG, ppsA, and recA); (2) ANIm = 1; (3) ANIb = 1; (4) OrthoANIb = 100%; (5) OrthoANIu = 100%; (6) TETRA = 1.

4.7. Phylogeny

Prokka output files were submitted to Roary v. 3.13.0 [73], and the core genome was aligned with MAFFT v. 7.505 [74]. IQ-TREE v. 2.2.03 [75] was performed to build the phylogenetic tree from the core genome with the substitution model GTR + F+I + G4 and 1000 replicates from Bootstrap, and iTOL v. 6 (Available online: https://itol.embl.de/, accessed on 15 August 2022) was used to edit and annotate the tree.

4.8. Virulence and Resistance Genes Analysis

Abricate version 1.0.1 (Available online: https://github.com/tseemann/abricate, accessed on 12 July 2022) was used to identify virulence and resistance genes in genomes and plasmids, with a minimum of 60% coverage and identity. The databases used for gene identification were VFDB [76] and CARD [50]. A presence/absence matrix was constructed and used to generate a dissimilarity matrix using the Jaccard index [77,78] and these data were used to perform PCoA (Principal Coordinate Analysis) [79] using the cmdscale function from R (Available online: https://cran.r-project.org/, accessed on 3 January 2022).

5. Conclusions

This study suggests that ANI is the method with the highest resolution for separating species using the genome information once phenotypic methods are not completely reliable for classifying Aeromonas species. Additionally, it can be seen that a widespread distribution of virulence genes identified in isolates originated from different sources and geographic origins. Additionally, the antibiotic resistance profile revealed that A. veronii harbors genes related to several classes of antimicrobial resistance genes, including those that confer resistance to antimicrobials used as the last line of defense against resistant bacteria such as carbapenem and colistin. Taken together, our results pinpoint that A. veronii is an emerging pathogen that deserves more attention from surveillance systems from a One Health perspective, since it can be involved in local and worldwide dispersion of virulence and resistance, thus causing fatal diseases in humans and fishes.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/antibiotics12061039/s1. Figure S1: Heatmap illustrating the relationship between the genomes according to the pairwise calculated tetra nucleotide correlation index. San Francisco Valley (SFV) isolates highlighted in blue. Values closer to red indicate proximity to the maximum index (1,000). The figure was plotted as a standard output of the pyANI tool (Available online: https://github.com/widdowquinn/pyani, accessed on 21 May 2022). Figure S2: Venn diagrams illustrating the number of genomes ranked with the reference for the four methods for estimating average nucleotide identity according to the cutoff points tested: a) cut-off point equal to 95%; b) cut-off point equal to 96%. The figure was plotted with the R language (Avail-able online: https://cran.r-project.org/, accessed on 3 January 2022) using the VennDiagram package (Available online: https://www.rdocumentation.org/packages/VennDiagram/versions/1.7.3/topics/venn.diagram, accessed on 23 August 2022). Figure S3: Phylogeny of SNPs from the core genome including bootstrap values (in red), sample location and isolation source. Misidentified isolates in red; SFV isolates of A. veronii in blue. The tree was edited using iTOL (Available online: https://itol.embl.de/, accessed on 15 August 2022) and host and isolation information was added using inkscape (Available online: https://inkscape.org/pt-br/, accessed on 14 October 2022). Figure S4: Venn Diagram illustrating the amount of virulence genes shared between SFV genomes and NCBI genomes, as well as genes particular to the two groups. The figure was plotted with the R language (Available online: https://cran.r-project.org/, accessed on 3 January 2022) using the VennDi-agram package (Available online: https://www.rdocumentation.org/packages/VennDiagram/versions/1.7.3/topics/venn.diagram, ac-cessed on 23 August 2022). Figure S5: Veen Diagram illustrating the amount of genes shared between isolates from different iso-lation sources, as well as those that are unique to each. The figure was plotted with the R language (Available online: https://cran.r-project.org/, accessed on 3 January 2022) using the VennDiagram package (Available online: https://www.rdocumentation.org/packages/VennDiagram/versions/1.7.3/topics/venn.diagram, ac-cessed on 23 August 2022). Figure S6: Amount of structure genes and virulence-related elements in each Aeromonas veronii isolate. SFV isolates highlighted in blue. The figure was plotted using R (Available online: https://cran.r-project.org/, accessed on 3 January 2022) with the package Pheatmap (Available online: https://www.rdocumentation.org/packages/pheatmap/versions/1.0.12/topics/pheatmap, accessed o 30 August 2022). Figure S7: PCOA from presence/absence of genes plotted as scatter. a) Shape and color according to source of isolation. b) Shapes according to Host and color according to location. c) Shapes according to health status information and colors according to type of isolation source. The figure was plotted using R (Available online: https://cran.r-project.org/, accessed on 3 January 2022) with the package ggplot2 (Available online: https://ggplot2.tidyverse.org/, accessed on 4 January 2022). Figure S8: Amount of gene copies for each resistance mechanism in isolates. SFV isolates highlighted in blue. The figure was plotted using R (Available online: https://cran.r-project.org/, accessed on 3 January 2022) with the package Pheatmap (Available online: https://www.rdocumentation.org/packages/pheatmap/versions/1.0.12/topics/pheatmap, accessed o 30 August 2022). Figure S9: Venn Diagram illustrating the amount of antibiotic resistance-related genes shared between the SFV and NCBI isolates, as well as those that are unique to each group. The figure was plot-ted with the R language (Available online: https://cran.r-project.org/, accessed on 3 January 2022) using the VennDiagram package (Available online: https://www.rdocumentation.org/packages/VennDiagram/versions/1.7.3/topics/venn.diagram, ac-cessed on 23 August 2022). Figure S10: Figure illustrating the amount of genes from each gene family present in each genome. SFV isolates highlighted in blue. The figure was plotted using R (Available online: https://cran.r-project.org/, accessed on 3 January 2022) with the package Pheatmap (Available online: https://www.rdocumentation.org/packages/pheatmap/versions/1.0.12/topics/pheatmap, accessed o 30 August 2022). Figure S11: Number of genes present for each resistance mechanism in plasmids. Conjugative plas-mids in red; mobilizable plasmids in blue; and unmobilizable plasmids in green. The figure was plot-ted using R (Available online: https://cran.r-project.org/, accessed on 3 January 2022) with the pack-age Pheatmap (Available online: https://www.rdocumentation.org/packages/pheatmap/versions/1.0.12/topics/pheatmap, accessed o 30 August 2022). Table S1: Result of mlst typing. Table S2: A. veronii genomes downloaded from the NCBI. Table S3: Reference genomes of species of the genus Aeromonas downloaded from the NCBI. Table S4: Additional reference genomes included in the core genome phylogeny. Table S5: Putative plasmids identified in the SFV (São Francisco Valley) isolates. Table S6: Information on plasmids associated with isolates downloaded from the NCBI. Table S7: Assembly metrics of the genomes sequenced in this study.

Author Contributions

Conceptualization: J.C.d.S.M., J.J.d.S.G., G.V.G. and M.M.d.C.; formal analysis: J.C.d.S.M. and G.A.d.A.S.; investigation: J.C.d.S.M.; methodology: J.C.d.S.M. and J.J.d.S.G.; resources: B.B., F.A.A., F.A., R.T.J.R., S.C.S., V.A.d.C.A., A.G.-N. and M.M.d.C.; software: J.C.d.S.M., R.T.J.R. and S.C.S.; supervision: V.A.d.C.A. and M.M.d.C.; writing—original draft: J.C.d.S.M. and L.S.d.B.C.; writing—review and editing: J.C.d.S.M., G.A.d.A.S., L.S.d.B.C., G.V.G., A.G.-N., B.B., F.A.A., F.A., R.T.J.R., S.C.S., V.A.d.C.A., M.M.d.C. and J.J.d.S.G. All authors have read and agreed to the published version of the manuscript.

Funding

The Coordination for the Improvement of Higher Education Personnel—Brazil (CAPES)—Finance Code 001 has funded this study. The authors also thank the Foundation for the Support of Science and Technology of the State of Pernambuco (FACEPE-APQ-0474-5.05-19) and the National Council for Scientific and Technological Development (CNPq—Process 408695/2022-6) for their financial support.

Institutional Review Board Statement

Bacteria isolation from fish was approved by the Ethics Committee for the use of Animals in Research (no. 0003/170316) from the Federal University of Vale of São Francisco (Univasf), and bacteria isolation from humans was approved by the Human Research Ethics Committee (no. 40360120.4.0000.9547) from the Federal Rural University of Pernambuco (UFRPE).

Informed Consent Statement

Not applicable.

Data Availability Statement

Genome sequence data are available at the National Center for Biotechnology Information (NCBI).

Acknowledgments

The authors would like to thank the Center for Open Access Genomic Analysis (CALAnGO), Brazilian Financier of Studies and Projects (FINEP), and Omics Science Network (RECOM).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Martin-Carnahan, A.; Joseph, S.W. Aeromonas Stanier 1943, 213. In Bergey’s Manual of Systematic Bacteriology; Brenner, D.J., Krieg, N.R., Staley, J.T., Garrity, G.M., Eds.; John Wiley & Sons, Inc.: Philadelphia, PA, USA, 2015; pp. 1–44. ISBN 9781118960608. [Google Scholar]
  2. Fernández-Bravo, A.; Figueras, M.J. An Update on the Genus Aeromonas: Taxonomy, Epidemiology, and Pathogenicity. Microorganisms 2020, 8, 129. [Google Scholar] [CrossRef] [Green Version]
  3. Dong, H.T.; Techatanakitarnan, C.; Jindakittikul, P.; Thaiprayoon, A.; Taengphu, S.; Charoensapsri, W.; Khunrae, P.; Rattanarojpong, T.; Senapin, S. Aeromonas jandaei and Aeromonas veronii Caused Disease and Mortality in Nile Tilapia, Oreochromis niloticus (L.). J. Fish Dis. 2017, 40, 1395–1403. [Google Scholar] [CrossRef] [PubMed]
  4. Ran, C.; Qin, C.; Xie, M.; Zhang, J.; Li, J.; Xie, Y.; Wang, Y.; Li, S.; Liu, L.; Fu, X.; et al. Aeromonas veronii and Aerolysin are Important for the Pathogenesis of Motile Aeromonad Septicemia in Cyprinid Fish. Environ. Microbiol. 2018, 20, 3442–3456. [Google Scholar] [CrossRef]
  5. Yuwono, C.; Wehrhahn, M.C.; Liu, F.; Riordan, S.M.; Zhang, L. The Isolation of Aeromonas Species and Other Common Enteric Bacterial Pathogens from Patients with Gastroenteritis in an Australian Population. Microorganisms 2021, 9, 1440. [Google Scholar] [CrossRef]
  6. Liu, F.; Yuwono, C.; Tay, A.C.Y.; Wehrhahn, M.C.; Riordan, S.M.; Zhang, L. Analysis of Global Aeromonas veronii Genomes Provides Novel Information on Source of Infection and Virulence in Human Gastrointestinal Diseases. BMC Genom. 2022, 23, 166. [Google Scholar] [CrossRef]
  7. Prediger, K.D.C.; Dallagassa, C.B.; Moriel, B.; Vizzotto, B.S.; Volanski, W.; Souza, E.M.; Pedrosa, F.O.; Weiss, V.; Alberton, D.; Guizelini, D.; et al. Virulence Characteristics and Antimicrobial Resistance of Aeromonas veronii Biovar Sobria 312M, a Clinical Isolate. Braz. J. Microbiol. 2020, 51, 511–518. [Google Scholar] [CrossRef]
  8. Tekedar, H.C.; Kumru, S.; Blom, J.; Perkins, A.D.; Griffin, M.J.; Abdelhamed, H.; Karsi, A.; Lawrence, M.L. Comparative Genomics of Aeromonas veronii: Identification of a Pathotype Impacting Aquaculture Globally. PLoS ONE 2019, 14, e0221018. [Google Scholar] [CrossRef] [Green Version]
  9. Yang, B.; Chen, C.; Sun, Y.; Cao, L.; Zhang, D.; Sun, W.; Zhang, L.; Wang, G.; Shan, X.; Kang, Y.; et al. Comparative Genomic Analysis of Different Virulence Strains Reveals Reasons for the Increased Virulence of Aeromonas veronii. J. Fish Dis. 2021, 44, 11–24. [Google Scholar] [CrossRef] [PubMed]
  10. Sheng, T.; Song, G.; Yue, T.; Zhang, J.; Wang, W.; Yang, Z.; Lu, Q. Whole-Genome Sequencing and Antimicrobial Resistance Analysis of Multidrug-Resistant Aeromonas veronii strain JC529 from a Common Carp. J. Glob. Antimicrob. Resist. 2021, 27, 118–122. [Google Scholar] [CrossRef]
  11. Yang, Q.; Zhao, M.; Wang, K.Y.; Wang, J.; He, Y.; Wang, E.L.; Liu, T.; Chen, D.F.; Lai, W. Multidrug-Resistant Aeromonas veronii Recovered from Channel Catfish (Ictalurus punctatus) in China: Prevalence and Mechanisms of Fluoroquinolone Resistance. Microb. Drug Resist. 2017, 23, 473–479. [Google Scholar] [CrossRef] [PubMed]
  12. Bello-López, J.M.; Cabrero-Martínez, O.A.; Ibáñez-Cervantes, G.; Hernández-Cortez, C.; Pelcastre-Rodríguez, L.I.; Gonzalez-Avila, L.U.; Castro-Escarpulli, G. Horizontal Gene Transfer and Its Association with Antibiotic Resistance in the Genus Aeromonas spp. Microorganisms 2019, 7, 363. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Yin, W.; Li, H.; Shen, Y.; Liu, Z.; Wang, S.; Shen, Z.; Zhang, R.; Walsh, T.R.; Shen, J.; Wang, Y. Beijing Novel Plasmid-Mediated Colistin Resistance Gene mcr-3 in Escherichia coli Wenjuan. MBio 2017, 8, e00543-17. [Google Scholar] [PubMed] [Green Version]
  14. Yang, Y.Q.; Li, Y.X.; Lei, C.W.; Zhang, A.Y.; Wang, H.N. Novel Plasmid-Mediated Colistin Resistance Gene mcr-7.1 in Klebsiella pneumoniae. J. Antimicrob. Chemother. 2018, 73, 1791–1795. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Cherak, Z.; Loucif, L.; Moussi, A.; Rolain, J.M. Epidemiology of Mobile Colistin Resistance (mcr) Genes in Aquatic Environments. J. Glob. Antimicrob. Resist. 2021, 27, 51–62. [Google Scholar] [CrossRef]
  16. Pessoa, R.B.G.; de Oliveira, W.F.; dos Santos Correia, M.T.; Fontes, A.; Coelho, L.C.B.B. Aeromonas and Human Health Disorders: Clinical Approaches. Front. Microbiol. 2022, 13, 1266. [Google Scholar] [CrossRef]
  17. Beaz-Hidalgo, R.; Hossain, M.J.; Liles, M.R.; Figueras, M.J. Strategies to Avoid Wrongly Labelled Genomes Using as Example the Detected Wrong Taxonomic Affiliation for Aeromonas Genomes in the Genbank Database. PLoS ONE 2015, 10, e0115813. [Google Scholar] [CrossRef]
  18. Awan, F.; Dong, Y.; Liu, J.; Wang, N.; Mushtaq, M.H.; Lu, C.; Liu, Y. Comparative Genome Analysis Provides Deep Insights into Aeromonas hydrophila Taxonomy and Virulence-Related Factors. BMC Genomics 2018, 19, 712. [Google Scholar] [CrossRef] [Green Version]
  19. Lamy, B.; Laurent, F.; Verdier, I.; Decousser, J.W.; Lecaillon, E.; Marchandin, H.; Roger, F.; Tigaud, S.; de Montclos, H.; Kodjo, A. Accuracy of 6 Commercial Systems for Identifying Clinical Aeromonas Isolates. Diagn. Microbiol. Infect. Dis. 2010, 67, 9–14. [Google Scholar] [CrossRef]
  20. Ceric, O.; Tyson, G.H.; Goodman, L.B.; Mitchell, P.K.; Zhang, Y.; Prarat, M.; Cui, J.; Peak, L.; Scaria, J.; Antony, L.; et al. Enhancing the One Health Initiative by Using Whole Genome Sequencing to Monitor Antimicrobial Resistance of Animal Pathogens: Vet-LIRN Collaborative Project with Veterinary Diagnostic Laboratories in United States and Canada. BMC Vet. Res. 2019, 15, 130. [Google Scholar] [CrossRef] [Green Version]
  21. Richter, M.; Rosselló-Móra, R. Shifting the Genomic Gold Standard for the Prokaryotic Species Definition. Proc. Natl. Acad. Sci. USA 2009, 106, 19126–19131. [Google Scholar] [CrossRef] [Green Version]
  22. Colston, S.M.; Fullmer, M.S.; Beka, L.; Lamy, B.; Gogarten, J.P. Bioinformatic Genome Comparisons for Taxonomic and Phylogenetic Assignments Using Aeromonas as a Test Case. MBio 2014, 5, e02136-14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Beaz-Hidalgo, R.; Latif-Eugenín, F.; Hossain, M.J.; Berg, K.; Niemi, R.M.; Rapala, J.; Lyra, C.; Liles, M.R.; Figueras, M.J. Aeromonas aquatica sp. nov., Aeromonas finlandiensis sp. nov. and Aeromonas lacus sp. nov. Isolated from Finnish Waters Associated with Cyanobacterial Blooms. Syst. Appl. Microbiol. 2015, 38, 161–168. [Google Scholar] [CrossRef] [PubMed]
  24. Hadi, N.; Yang, Q.; Barnett, T.C.; Tabei, S.M.B.; Kirov, S.M.; Shaw, J.G. Bundle-Forming Pilus Locus of Aeromonas veronii bv. sobria. Infect. Immun. 2012, 80, 1351–1360. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Werneburg, G.T.; Thanassi, D.G. Pili Assembled by the Chaperone/Usher Pathway in Escherichia coli and Salmonella. EcoSal Plus 2018, 8, 631–632. [Google Scholar] [CrossRef] [Green Version]
  26. Dacanay, A.; Boyd, J.M.; Fast, M.D.; Knickle, L.C.; Reith, M.E. Aeromonas salmonicida Type I Pilus System Contributes to Host Colonization but Not Invasion. Dis. Aquat. Organ. 2010, 88, 199–206. [Google Scholar] [CrossRef]
  27. Kirov, S.M.; Tassell, B.C.; Semmler, A.B.T.; O’Donovan, L.A.; Rabaan, A.A.; Shaw, J.G. Lateral Flagella and Swarming Motility in Aeromonas Species. J. Bacteriol. 2002, 184, 547–555. [Google Scholar] [CrossRef] [Green Version]
  28. Gavín, R.; Merino, S.; Altarriba, M.; Canals, R.; Shaw, J.G.; Tomás, J.M. Lateral Flagella Are Required for Increased Cell Adherence, Invasion and Biofilm Formation by Aeromonas spp. FEMS Microbiol. Lett. 2003, 224, 77–83. [Google Scholar] [CrossRef] [Green Version]
  29. Green, E.R.; Mecsas, J. Bacterial Secretion Systems: An Overview. Microbiol. Spectr. 2016, 4, 213–239. [Google Scholar] [CrossRef] [Green Version]
  30. Mekasha, S.; Linke, D. Secretion Systems in Gram-Negative Bacterial Fish Pathogens. Front. Microbiol. 2021, 12, 782673. [Google Scholar] [CrossRef]
  31. Naskar, S.; Hohl, M.; Tassinari, M.; Low, H.H. The Structure and Mechanism of the Bacterial Type II Secretion System. Mol. Microbiol. 2021, 115, 412–424. [Google Scholar] [CrossRef]
  32. Korotkov, K.V.; Sandkvist, M. Architecture, Function, and Substrates of the Type II Secretion System. EcoSal Plus. 2019, 8. [Google Scholar] [CrossRef] [PubMed]
  33. Chakraborty, T.; Schmid, A.; Notermans, S.; Benz, R. Aerolysin of Aeromonas sobria: Evidence for Formation of Ion-Permeable Channels and Comparison with Alpha-Toxin of Staphylococcus aureus. Infect. Immun. 1990, 58, 2127–2132. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Hews, C.L.; Tran, S.L.; Wegmann, U.; Brett, B.; Walsham, A.D.S.; Kavanaugh, D.; Ward, N.J.; Juge, N.; Schüller, S. The StcE Metalloprotease of Enterohaemorrhagic Escherichia coli Reduces the Inner Mucus Layer and Promotes Adherence to Human Colonic Epithelium ex vivo. Cell. Microbiol. 2017, 19, e12717. [Google Scholar] [CrossRef] [Green Version]
  35. Dey, S.; Chakravarty, A.; Guha Biswas, P.; De Guzman, R.N. The Type III Secretion System Needle, Tip, and Translocon. Protein Sci. 2019, 28, 1582–1593. [Google Scholar] [CrossRef] [PubMed]
  36. Pinaud, L.; Sansonetti, P.J.; Phalipon, A. Host Cell Targeting by Enteropathogenic Bacteria T3SS Effectors. Trends Microbiol. 2018, 26, 266–283. [Google Scholar] [CrossRef] [PubMed]
  37. Silver, A.C.; Graf, J. Prevalence of Genes Encoding the Type Three Secretion System and the Effectors AexT and AexU in the Aeromonas veronii Group. DNA Cell Biol. 2009, 28, 383–388. [Google Scholar] [CrossRef]
  38. Yu, K.W.; Xue, P.; Fu, Y.; Yang, L. T6SS Mediated Stress Responses for Bacterial Environmental Survival and Host Adaptation. Int. J. Mol. Sci. 2021, 22, 478. [Google Scholar] [CrossRef]
  39. Si, M.; Wang, Y.; Zhang, B.; Zhao, C.; Kang, Y.; Bai, H.; Wei, D.; Zhu, L.; Zhang, L.; Dong, T.G.; et al. The Type VI Secretion System Engages a Redox-Regulated Dual-Functional Heme Transporter for Zinc Acquisition. Cell Rep. 2017, 20, 949–959. [Google Scholar] [CrossRef] [Green Version]
  40. Schauder, S.; Shokat, K.; Surette, M.G.; Bassler, B.L. The LuxS Family of Bacterial Autoinducers: Biosynthesis of a Novel Quorum-Sensing Signal Molecule. Mol. Microbiol. 2001, 41, 463–476. [Google Scholar] [CrossRef]
  41. Liu, L.; Wu, R.; Zhang, J.; Li, P. Overexpression of LuxS Promotes Stress Resistance and Biofilm Formation of Lactobacillus paraplantarumL-ZS9 by Regulating the Expression of Multiple Genes. Front. Microbiol. 2018, 9, 2628. [Google Scholar] [CrossRef]
  42. Wang, Y.; Liu, B.; Grenier, D.; Yi, L. Regulatory Mechanisms of the LuxS/AI-2 System and Bacterial Resistance. Am. Soc. Microbiol. 2019, 63, e01186-19. [Google Scholar] [CrossRef] [Green Version]
  43. Ma, Y.; Hao, L.; Ke, H.; Liang, Z.; Ma, J.; Liu, Z.; Li, Y. LuxS/AI-2 in Streptococcus agalactiae Reveals a Key Role in Acid Tolerance and Virulence. Res. Vet. Sci. 2017, 115, 501–507. [Google Scholar] [CrossRef] [PubMed]
  44. Castanheira, M.; Simner, P.J.; Bradford, P.A. Extended-Spectrum β-Lactamases: An Update on Their Characteristics, Epidemiology and Detection. JAC Antimicrob. Resist. 2021, 3, dlab092. [Google Scholar] [CrossRef]
  45. Bahr, G.; González, L.J.; Vila, A.J. Metallo-β-Lactamases in the Age of Multidrug Resistance: From Structure and Mechanism to Evolution, Dissemination, and Inhibitor Design. Chem. Rev. 2021, 121, 7957–8094. [Google Scholar] [CrossRef] [PubMed]
  46. Du, D.; Wang-Kan, X.; Neuberger, A.; van Veen, H.W.; Pos, K.M.; Piddock, L.J.V.; Luisi, B.F. Multidrug Efflux Pumps: Structure, Function and Regulation. Nat. Rev. Microbiol. 2018, 16, 523–539. [Google Scholar] [CrossRef] [PubMed]
  47. Blair, J.M.A.; Richmond, G.E.; Piddock, L.J.V. Multidrug Efflux Pumps in Gram-Negative Bacteria and Their Role in Antibiotic Resistance. Future Microbiol. 2014, 9, 1165–1177. [Google Scholar] [CrossRef] [PubMed]
  48. Cloeckaert, A.; Baucheron, S.; Chaslus-Dancla, E. Nonenzymatic Chloramphenicol Resistance Mediated by IncC Plasmid R55 Is Encoded by a floR Gene Variant. Antimicrob. Agents Chemother. 2001, 45, 2381–2382. [Google Scholar] [CrossRef] [Green Version]
  49. Shi, Y.; Shi, Y.; Tian, Z.; Tian, Z.; Gillings, M.R.; Zhang, Y.; Zhang, Y.; Zhang, H.; Huyan, J.; Huyan, J.; et al. Novel Transposon Tn 6433 Variants Accelerate the Dissemination of tet(E) in Aeromonas in an Aerobic Biofilm Reactor under Oxytetracycline Stresses. Environ. Sci. Technol. 2020, 54, 6781–6791. [Google Scholar] [CrossRef]
  50. McArthur, A.G.; Waglechner, N.; Nizam, F.; Yan, A.; Azad, M.A.; Baylay, A.J.; Bhullar, K.; Canova, M.J.; De Pascale, G.; Ejim, L.; et al. The Comprehensive Antibiotic Resistance Database. Antimicrob. Agents Chemother. 2013, 57, 3348–3357. [Google Scholar] [CrossRef] [Green Version]
  51. Von Wintersdorff, C.J.H.; Penders, J.; Van Niekerk, J.M.; Mills, N.D.; Majumder, S.; Van Alphen, L.B.; Savelkoul, P.H.M.; Wolffs, P.F.G. Dissemination of Antimicrobial Resistance in Microbial Ecosystems through Horizontal Gene Transfer. Front. Microbiol. 2016, 7, 173. [Google Scholar] [CrossRef] [Green Version]
  52. Shintani, M.; Sanchez, Z.K.; Kimbara, K. Genomics of Microbial Plasmids: Classification and Identification Based on Replication and Transfer Systems and Host Taxonomy. Front. Microbiol. 2015, 6, 242. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Zhong, C.; Han, M.; Yang, P.; Chen, C.; Yu, H.; Wang, L.; Ning, K. Comprehensive Analysis Reveals the Evolution and Pathogenicity of Aeromonas, Viewed from Both Single Isolated Species and Microbial Communities. mSystems 2019, 4, e00252-19. [Google Scholar] [CrossRef] [Green Version]
  54. Talagrand-Reboul, E.; Colston, S.M.; Graf, J.; Lamy, B.; Jumas-Bilak, E. Comparative and Evolutionary Genomics of Isolates Provide Insight into the Pathoadaptation of Aeromonas. Genome Biol. Evol. 2020, 12, 535–552. [Google Scholar] [CrossRef] [Green Version]
  55. Pessoa, R.B.G.; Marques, D.S.C.; Lima, R.O.H.A.; Oliveira, M.B.M.; Lima, G.M.S.; Maciel de Carvalho, E.V.M.; Coelho, L.C.B.B. Molecular Characterization and Evaluation of Virulence Traits of Aeromonas spp. Isolated from the Tambaqui Fish (Colossoma macropomum). Microb. Pathog. 2020, 147, 104273. [Google Scholar] [CrossRef] [PubMed]
  56. Li, T.; Raza, S.H.A.; Yang, B.; Sun, Y.; Wang, G.; Sun, W.; Qian, A.; Wang, C.; Kang, Y.; Shan, X. Aeromonas veronii Infection in Commercial Freshwater Fish: A Potential Threat to Public Health. Animals 2020, 10, 608. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Didugu, H.; Thirtham, M.; Nelapati, K.; Reddy, K.K.; Kumbhar, B.S.; Poluru, A.; Pothanaboyina, G. A Study on the Prevalence of Aeromonas spp. and Its Enterotoxin Genes in Samples of Well Water, Tap Water, and Bottled Water. Vet. World 2015, 8, 1237–1242. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Snyman, Y.; Whitelaw, A.C.; Barnes, J.M.; Maloba, M.R.B.; Newton-Foot, M. Characterisation of Mobile Colistin Resistance Genes (mcr-3 and mcr-5) in River and Storm Water in Regions of the Western Cape of South Africa. Antimicrob. Resist. Infect. Control 2021, 10, 1–9. [Google Scholar] [CrossRef]
  59. Shao, S.; Hu, Y.; Cheng, J.; Chen, Y. Research Progress on Distribution, Migration, Transformation of Antibiotics and Antibiotic Resistance Genes (ARGs) in Aquatic Environment. Crit. Rev. Biotechnol. 2018, 38, 1195–1208. [Google Scholar] [CrossRef]
  60. Gouveia, J.J.S.; Regitano, L.C.A. Extração de DNA. In Protocolos em Biologia Molecular Aplicada à Produção Animal; Regitano, L.C.A., Niciura, S.C.M., Ibelli, A.M.G., Gouveia, J.J., Eds.; Embrapa: São Carlos, Brazil; pp. 3–8. ISBN 9788586763127.
  61. Trivedi, U.H.; Cézard, T.; Bridgett, S.; Montazam, A.; Nichols, J.; Blaxter, M.; Gharbi, K. Quality Control of Next-Generation Sequencing Data without a Reference. Front. Genet. 2014, 5, 111. [Google Scholar] [CrossRef] [Green Version]
  62. Lindgreen, S. AdapterRemoval: Easy Cleaning of next-Generation Sequencing Reads. BMC Res. Notes 2012, 5, 337. [Google Scholar] [CrossRef] [Green Version]
  63. Melsted, P.; Halldórsson, B.V. KmerStream: Streaming Algorithms for k-Mer Abundance Estimation. Bioinformatics 2014, 30, 3541–3547. [Google Scholar] [CrossRef] [Green Version]
  64. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D.; et al. SPAdes: A New Genome Assembly Algorithm and Its Applications to Single-Cell Sequencing. J. Comput. Biol. 2012, 19, 455–477. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Hernandez, D.; François, P.; Farinelli, L.; Østerås, M.; Schrenzel, J. De Novo Bacterial Genome Sequencing: Millions of Very Short Reads Assembled on a Desktop Computer. Genome Res. 2008, 18, 802–809. [Google Scholar] [CrossRef] [Green Version]
  66. Wick, R.R.; Judd, L.M.; Gorrie, C.L.; Holt, K.E. Unicycler: Resolving Bacterial Genome Assemblies from Short and Long Sequencing Reads. PLoS Comput. Biol. 2017, 13, e1005595. [Google Scholar] [CrossRef] [Green Version]
  67. Li, W.; Godzik, A. Cd-Hit: A Fast Program for Clustering and Comparing Large Sets of Protein or Nucleotide Sequences. Bioinformatics 2006, 22, 1658–1659. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Gurevich, A.; Saveliev, V.; Vyahhi, N.; Tesler, G. QUAST: Quality Assessment Tool for Genome Assemblies. Bioinformatics 2013, 29, 1072–1075. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Robertson, J.; Nash, J.H.E. MOB-Suite: Software Tools for Clustering, Reconstruction and Typing of Plasmids from Draft Assemblies. Microb. Genom. 2018, 4, e000206. [Google Scholar] [CrossRef]
  70. Seemann, T. Prokka: Rapid Prokaryotic Genome Annotation. Bioinformatics 2014, 30, 2068–2069. [Google Scholar] [CrossRef] [Green Version]
  71. Pritchard, L.; Glover, R.H.; Humphris, S.; Elphinstone, J.G.; Toth, I.K. Genomics and Taxonomy in Diagnostics for Food Security: Soft-Rotting Enterobacterial Plant Pathogens. Anal. Methods 2016, 8, 12–24. [Google Scholar] [CrossRef]
  72. Lee, I.; Kim, Y.O.; Park, S.C.; Chun, J. OrthoANI: An Improved Algorithm and Software for Calculating Average Nucleotide Identity. Int. J. Syst. Evol. Microbiol. 2016, 66, 1100–1103. [Google Scholar] [CrossRef]
  73. Page, A.J.; Cummins, C.A.; Hunt, M.; Wong, V.K.; Reuter, S.; Holden, M.T.G.; Fookes, M.; Falush, D.; Keane, J.A.; Parkhill, J. Roary: Rapid Large-Scale Prokaryote Pan Genome Analysis. Bioinformatics 2015, 31, 3691–3693. [Google Scholar] [CrossRef] [Green Version]
  74. Katoh, K.; Misawa, K.; Kuma, K.I.; Miyata, T. MAFFT: A Novel Method for Rapid Multiple Sequence Alignment Based on Fast Fourier Transform. Nucleic Acids Res. 2002, 30, 3059–3066. [Google Scholar] [CrossRef] [Green Version]
  75. Minh, B.Q.; Schmidt, H.A.; Chernomor, O.; Schrempf, D.; Woodhams, M.D.; Von Haeseler, A.; Lanfear, R.; Teeling, E. IQ-TREE 2: New Models and Efficient Methods for Phylogenetic Inference in the Genomic Era. Mol. Biol. Evol. 2020, 37, 1530–1534. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Chen, L.; Yang, J.; Yu, J.; Yao, Z.; Sun, L.; Shen, Y.; Jin, Q. VFDB: A Reference Database for Bacterial Virulence Factors. Nucleic Acids Res. 2005, 33, 325–328. [Google Scholar] [CrossRef] [Green Version]
  77. Jaccard, P. Étude Comparative de La Distribution Florale Dans Une Portion Des Alpes et Du Jura. Bull. Soc. Vaud. Sci. Nat. 1901, 37, 547–579. [Google Scholar] [CrossRef]
  78. Gower, J.C.; Legendre, P. Metric and Euclidean Properties of Dissimilarity Coefficients. J. Classif. 1986, 3, 5–48. [Google Scholar] [CrossRef]
  79. Gower, J.C. Some Distance Properties of Latent Root and Vector Methods Used in Multivariate Analysis. Biometrika 1966, 53, 325. [Google Scholar] [CrossRef]
Figure 1. Heatmap illustrating the genomes relationship according to the estimated pairwise nucleotide average identity. Warmer colors represent values closer to 1. A. veronii isolates from SFV are highlighted in blue; misidentified isolates are highlighted in red.
Figure 1. Heatmap illustrating the genomes relationship according to the estimated pairwise nucleotide average identity. Warmer colors represent values closer to 1. A. veronii isolates from SFV are highlighted in blue; misidentified isolates are highlighted in red.
Antibiotics 12 01039 g001
Figure 2. Tree built from core genome SNPs (841 genes). Misidentified isolates are highlighted in red; A. veronii isolates from SFV are highlighted in blue.
Figure 2. Tree built from core genome SNPs (841 genes). Misidentified isolates are highlighted in red; A. veronii isolates from SFV are highlighted in blue.
Antibiotics 12 01039 g002
Figure 3. Gene numbers belonging to each category of function related to the virulence in each genome. A.veronii isolates from SFV highlighted in blue.
Figure 3. Gene numbers belonging to each category of function related to the virulence in each genome. A.veronii isolates from SFV highlighted in blue.
Antibiotics 12 01039 g003
Figure 4. Genes number conferring drug class resistance in each isolate. Isolates from SFV are highlighted in blue.
Figure 4. Genes number conferring drug class resistance in each isolate. Isolates from SFV are highlighted in blue.
Antibiotics 12 01039 g004
Table 1. Information related to the Aeromonas veronii isolates from fishes, humans, and the environment in the São Francisco Valley.
Table 1. Information related to the Aeromonas veronii isolates from fishes, humans, and the environment in the São Francisco Valley.
IDIsolatesIsolatesHostLocationPatientPlasmids
B42PRJNA590960B42_Aer_Pi6_2MCPhractocephalus
hemioliopterus
CODEVASF *kidney-
B44PRJNA594379B44_Aer_Su2_1ASPseudoplatystoma corruscansCODEVASF *kidneyYes
B45PRJNA590964B45_Pi_23_1MCPhractocephalus
hemioliopterus
CODEVASF *kidneyYes
B47PRJNA590505B47_Aer_Pi4_2HPASPhractocephalus
hemioliopterus
CODEVASF *kidney-
B48PRJNA594423B48_Aer_Pi5_1HPASPhractocephalus
hemioliopterus
CODEVASF *kidney-
B52PRJNA595563B52_Aer_Pi17_Phractocephalus
hemioliopterus
CODEVASF *kidneyYes
1HPAS
B79PRJNA590958B79_Aer_onGT1Oreochromis
niloticus
Instituto Federalkidney-
B80PRJNA594352B80_Aer_Environmental (water)CODEVASF *Not applicable-
WatCTCBM21
MM_B02-MM_B02Homo sapiensTeaching HospitalSoft tissues infection-
* Development Company of the São Francisco and Parnaiba Valleys.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Maia, J.C.d.S.; Silva, G.A.d.A.; Cunha, L.S.d.B.; Gouveia, G.V.; Góes-Neto, A.; Brenig, B.; Araújo, F.A.; Aburjaile, F.; Ramos, R.T.J.; Soares, S.C.; et al. Genomic Characterization of Aeromonas veronii Provides Insights into Taxonomic Assignment and Reveals Widespread Virulence and Resistance Genes throughout the World. Antibiotics 2023, 12, 1039. https://doi.org/10.3390/antibiotics12061039

AMA Style

Maia JCdS, Silva GAdA, Cunha LSdB, Gouveia GV, Góes-Neto A, Brenig B, Araújo FA, Aburjaile F, Ramos RTJ, Soares SC, et al. Genomic Characterization of Aeromonas veronii Provides Insights into Taxonomic Assignment and Reveals Widespread Virulence and Resistance Genes throughout the World. Antibiotics. 2023; 12(6):1039. https://doi.org/10.3390/antibiotics12061039

Chicago/Turabian Style

Maia, José Cleves da Silva, Gabriel Amorim de Albuquerque Silva, Letícia Stheffany de Barros Cunha, Gisele Veneroni Gouveia, Aristóteles Góes-Neto, Bertram Brenig, Fabrício Almeida Araújo, Flávia Aburjaile, Rommel Thiago Jucá Ramos, Siomar Castro Soares, and et al. 2023. "Genomic Characterization of Aeromonas veronii Provides Insights into Taxonomic Assignment and Reveals Widespread Virulence and Resistance Genes throughout the World" Antibiotics 12, no. 6: 1039. https://doi.org/10.3390/antibiotics12061039

APA Style

Maia, J. C. d. S., Silva, G. A. d. A., Cunha, L. S. d. B., Gouveia, G. V., Góes-Neto, A., Brenig, B., Araújo, F. A., Aburjaile, F., Ramos, R. T. J., Soares, S. C., Azevedo, V. A. d. C., Costa, M. M. d., & Gouveia, J. J. d. S. (2023). Genomic Characterization of Aeromonas veronii Provides Insights into Taxonomic Assignment and Reveals Widespread Virulence and Resistance Genes throughout the World. Antibiotics, 12(6), 1039. https://doi.org/10.3390/antibiotics12061039

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop