Analysis of global Aeromonas veronii genomes provides novel information on source of infection and virulence in human gastrointestinal diseases

Phylogenomics of novel clones of Aeromonas veronii recovered from a freshwater lake reveals unique biosynthetic gene clusters

Aquatic ecosystems serve as crucial reservoirs for pathogens and antimicrobial resistance genes, thus presenting a significant global health risk. Here, we investigated the phylogenomics of Aeromonas veronii from Lake Wilcox in Ontario. Among the 11 bacterial isolates, nine were identified as A. veronii. Notably, 67% of A. veronii isolates were potential human pathogens. Considerable genetic diversity was noted among the A. veronii isolates, suggesting the lake as a reservoir for multiple human pathogenic strains. Comparison of the A. veronii sequenced with global A. veronii genomes highlighted significant genetic diversity and suggests widespread dissemination of strains. All the isolates carried chromosomal genes encoding resistance to β-lactams. Although virulence gene content differed between human and non-human pathogenic strains, type III secretion systems was associated with human pathogenic isolates. The assessment of AMR genes in global isolates showed that β-lactam and tetracycline resistance genes were predominant. Although the machine learning-based pangenome-wide association approach performed did not yield any source-based genes, some genes were enriched in a few isolates from different sources. The mrkABCDF operon that mediates biofilm formation and genes encoding resistance to colistin, chloramphenicol, trimethoprim, and tetracycline were enriched in animal products, whereas macrolide resistance genes and Inc plasmid-types were linked to the aquatic environment. Novel biosynthetic gene clusters were identified, suggesting that A. veronii with varying pathogenic potential could produce unique secondary metabolites. There is a need for continuous tracking of pathogens in aquatic ecosystems to contribute to our understanding of their evolutionary dynamics and the ecological roles of their genetic elements.

IMPORTANCE

Lakes and other aquatic ecosystems can harbor harmful bacteria that can make people sick and resist antibiotics, posing a significant global health risk. In this study, we investigated Aeromonas veronii, a Gram-negative bacteria found in Lake Wilcox in Ontario. We used various techniques, including whole-genome sequencing (WGS), to analyze the bacteria and found that many of the isolates had the potential to cause human disease. We also discovered significant genetic diversity among the isolates, indicating that the lake may be a reservoir for multiple human pathogenic strains. All isolates carried genes that confer resistance to antibiotics, and some virulence genes were associated with human pathogenic isolates. This study highlights the importance of monitoring aquatic ecosystems for harmful bacteria to better understand their evolution, potential for human pathogenicity, and the ecological roles of their genetic elements. This knowledge can inform strategies for preventing the spread of antibiotic-resistant bacteria and protecting public health.

Catching the Big Fish in Big Data: A Meta-Analysis of Zebrafish Kidney scRNA-Seq Datasets Highlights Conserved Molecular Profiles of Macrophages and Neutrophils in Vertebrates

The innate immune system (IIS) is an ancient and essential defense mechanism that protects animals against a wide range of pathogens and diseases. Although extensively studied in mammals, our understanding of the IIS in other taxa remains limited. The zebrafish (Danio rerio) serves as a promising model organism for investigating IIS-related processes, yet the immunogenetics of fish are not fully elucidated. To address this gap, we conducted a meta-analysis of single-cell RNA sequencing (scRNA-seq) datasets from zebrafish kidney marrow, encompassing approximately 250,000 immune cells. Our analysis confirms the presence of key genetic pathways in zebrafish innate immune cells that are similar to those identified in mammals. Zebrafish macrophages specifically express genes encoding cathepsins, major histocompatibility complex class II proteins, integral membrane proteins, and the V-ATPase complex and demonstrate the enrichment of oxidative phosphorylation ferroptosis processes. Neutrophils are characterized by the significant expression of genes encoding actins, cytoskeleton organizing proteins, the Arp2/3 complex, and glycolysis enzymes and have demonstrated their involvement in GnRH and CLR signaling pathways, adherents, and tight junctions. Both macrophages and neutrophils highly express genes of NOD-like receptors, phagosomes, and lysosome pathways and genes involved in apoptosis. Our findings reinforce the idea about the existence of a wide spectrum of immune cell phenotypes in fish since we found only a small number of cells with clear pro- or anti-inflammatory signatures.

Phenotypic and genotypic analysis of antimicrobial resistance and population structure of gastroenteritis-related Aeromonas isolates

BACKGROUND

The population structure and the correlation between antimicrobial resistance (AMR) phenotypes and genotypes in Aeromonas species isolated from patients with gastroenteritis are not well understood. The aims of the study were to: (1) investigate the antimicrobial susceptibility profiles of Aeromonas species isolated from patients with gastroenteritis; (2) explore the relationship between AMR genes and resistance phenotypes; and (3) describe the population structure of these isolates and provide evidence of transmission events among them.

METHODS

This microbiological survey was performed at the Microbiology Laboratory of the Emek Medical Center in Afula, Israel. Cultivation of Aeromonas was attempted from stool samples that tested positive by PCR. Antimicrobial susceptibility testing (AST) was performed using the Sensititre GN3F microdilution panel. Whole genome sequencing (WGS) was done using the Illumina NextSeq500/550 system. Phylogenetic studies involved multi-locus sequence typing (MLST) and core genome (cg) MLST. Resistance mechanisms were identified using the Comprehensive Antibiotic Resistance Database and compared with the AST results.

RESULTS

The study included 67 patient-unique isolates. The species that were identified included A. caviae (n = 58), A. dhakensis (n = 3), A. media (n = 2), A. veronii (n = 2) and A. hydrophila (n = 2). Isolates were almost uniformly susceptible to amikacin, gentamicin, aztreonam, cefepime, ceftazidime, ciprofloxacin and meropenem. All isolates with the exception of 1-2 isolates were resistant to ampicillin, cefazolin and ampicillin-sulbactam which was compatible with the presence of the blaOXA genes. Variable resistance rates were observed to cefuroxime, cefoxitin, ceftriaxone, piperacillin-tazobactam that were not correlated with the presence of other β-lactamase genes. Resistance to tetracycline and trimethoprim-sulfamethoxazole correlated with the presence of tetA and sul1, respectively. The population structure of A. caviae was highly diverse with the minority of the isolates (16/57) clustering into six defined sequence types. A cgMLST-based distance of four genes was found in one pair of isolates, suggesting common source transmission.

CONCLUSIONS

A. caviae is the dominant species related to gastroenteritis and is characterized by a diverse population structure, with almost no evidence for common-source transmission. Resistance rates to most antimicrobial agents were low and partially matched with the presence of resistance genes.

Genomic characterization of Aeromonas spp. isolates from striped catfish with motile Aeromonas septicemia and human bloodstream infections in Vietnam

Aeromonas spp. are commonly found in the aquatic environment and have been responsible for motile Aeromonas septicemia (MAS) in striped catfish, resulting in significant economic loss. These organisms also cause a range of opportunistic infections in humans with compromised immune systems. Here, we conducted a genomic investigation of 87 Aeromonas isolates derived from diseased catfish, healthy catfish and environmental water in catfish farms affected by MAS outbreaks in eight provinces in Mekong Delta (years: 2012-2022), together with 25 isolates from humans with bloodstream infections (years: 2010-2020). Genomics-based typing method precisely delineated Aeromonas species while traditional methods such as aerA PCR and MALDI-TOF were unable identify A. dhakensis. A. dhakensis was found to be more prevalent than A. hydrophila in both diseased catfish and human infections. A. dhakensis sequence type (ST) 656 followed by A. hydrophila ST251 were the predominant virulent species-lineages in diseased catfish (43.7 and 20.7 %, respectively), while diverse STs were found in humans with bloodstream infections. There was evidence of widespread transmission of ST656 and ST251 on striped catfish in the Mekong Delta region. ST656 and ST251 isolates carried a significantly higher number of acquired antimicrobial resistance (AMR) genes and virulence factors in comparison to other STs. They, however, exhibited several distinctions in key virulence factors (i.e. lack of type IV pili and enterotoxin ast in A. dhakensis), AMR genes (i.e. presence of imiH carbapenemase in A. dhakensis), and accessory gene content. To uncover potential conserved proteins of Aeromonas spp. for vaccine development, pangenome analysis has unveiled 2202 core genes between ST656 and ST251, of which 78 proteins were in either outer membrane or extracellular proteins. Our study represents one of the first genomic investigations of the species distribution, genetic landscape, and epidemiology of Aeromonas in diseased catfish and human infections in Vietnam. The emergence of antimicrobial resistant and virulent A. dhakensis strains underscores the needs of enhanced genomic surveillance and strengthening vaccine research and development in preventing Aeromonas diseases in catfish and humans, and the search for potential vaccine candidates could focus on Aeromonas core genes encoded for membrane and secreted proteins.

Analysis of global Aeromonas caviae genomes revealed that strains carrying T6SS are more common in human gastroenteritis than in environmental sources and are often phylogenetically related

Aeromonas caviae is an emerging human enteric pathogen. However, the genomic features and virulence genes of A. caviae strains from human gastroenteritis and other sources have not been fully elucidated. Here, we conducted a genomic analysis of 565 global A. caviae strains isolated from different sources, including 261 strains isolated from faecal samples of gastroenteritis patients, of which 18 genomes were sequenced in this study. The presence of bacterial virulence genes and secretion systems in A. caviae strains from different sources was compared, and the phylogenetic relationship of A. caviae strains was assessed based on the core genome. The complete genome of A. caviae strain A20-9 isolated from a gastroenteritis patient was obtained in this study, from which 300 putative virulence factors and a T4SS-encoding plasmid, pAC, were identified. Genes encoding T4SS were also identified in a novel genomic island, ACI-1, from other T4SS-positive strains. The prevalence of T4SS was significantly lower in A. caviae strains from gastroenteritis patients than in environmental strains (3 %, P<0.0001 vs 14 %, P<0.01). Conversely, the prevalence of T6SS was significantly higher in A. caviae strains isolated from gastroenteritis patients than in environmental strains (25 %, P<0.05 vs 13  %, P<0.01). Four phylogenetic clusters were formed based on the core genome of 565 A. caviae strains, and strains carrying T6SS often showed close phylogenetic relationships. T3SS, aerolysin and thermostable cytotonic enterotoxin were absent in all 565 A. caviae strains. Our findings provide novel information on the genomic features of A. caviae and suggest that T6SS may play a role in A. caviae-induced human gastroenteritis.

Evaluate the potential use of TonB-dependent receptor protein as a subunit vaccine against Aeromonas veronii infection in Nile tilapia (Oreochromis niloticus)

Aeromonas veronii is an emerging bacterial pathogen that causes serious systemic infections in cultured Nile tilapia (Oreochromis niloticus), leading to massive deaths. Therefore, there is an urgent need to identify effective vaccine candidates to control the spread of this emerging disease. TonB-dependent receptor (Tdr) of A. veronii, which plays a role in the virulence factor of the organism, could be useful in terms of protective antigens for vaccine development. This study aims to evaluate the potential use of Tdr protein as a novel subunit vaccine against A. veronii infection in Nile tilapia. The Tdr gene from A. veronii was cloned into the pET28b expression vector, and the recombinant protein was subsequently produced in Escherichia coli strain BL21 (DE3). Tdr was expressed as an insoluble protein and purified by affinity chromatography. Antigenicity test indicated that this protein was recognized by serum from A. veronii infected fish. When Nile tilapia were immunized with the Tdr protein, specific antibody levels increased significantly (p-value <0.05) at 7 days post-immunization (dpi), and peaked at 21 dpi compared to antibody levels at 0 dpi. Furthermore, bacterial agglutination activity was observed in the fish serum immunized with the Tdr protein, indicating that specific antibodies in the serum can detect Tdr on the bacterial cell surface. These results suggest that Tdr protein has potential as a vaccine candidate. However, challenging tests with A.veronii in Nile tilapia needs to be investigated to thoroughly evaluate its protective efficacy for future applications.

Fishing Innate Immune System Properties through the Transcriptomic Single-Cell Data of Teleostei

The innate immune system is the first line of defense in multicellular organisms. Danio rerio is widely considered a promising model for IIS-related research, with the most amount of scRNAseq data available among Teleostei. We summarized the scRNAseq and spatial transcriptomics experiments related to the IIS for zebrafish and other Teleostei from the GEO NCBI and the Single-Cell Expression Atlas. We found a considerable number of scRNAseq experiments at different stages of zebrafish development in organs such as the kidney, liver, stomach, heart, and brain. These datasets could be further used to conduct large-scale meta-analyses and to compare the IIS of zebrafish with the mammalian one. However, only a small number of scRNAseq datasets are available for other fish (turbot, salmon, cavefish, and dark sleeper). Since fish biology is very diverse, it would be a major mistake to use zebrafish alone in fish immunology studies. In particular, there is a special need for new scRNAseq experiments involving nonmodel Teleostei, e.g., long-lived species, cancer-resistant fish, and various fish ecotypes.

Emerging Aeromonas enteric infections: their association with inflammatory bowel disease and novel pathogenic mechanisms

Aeromonas species are emerging human enteric pathogens. This study examines the isolation of Aeromonas and other enteric bacterial pathogens from patients with and without inflammatory bowel disease (IBD). This study also investigates the intestinal epithelial pathogenic mechanisms of Aeromonas veronii. The isolation rates of seven enteric bacterial pathogens from 2,279 patients with IBD and 373,276 non-IBD patients were compared. An A. veronii strain (AS1) isolated from intestinal biopsies of a patient with IBD was used for pathogenic mechanism investigation, and Escherichia coli K12 was used as a bacterial control. HT-29 cells were used as a model of human intestinal epithelium. A significantly higher isolation of Aeromonas species was found in patients with IBD as compared to non-IBD patients (P = 0.0001, odds ratio = 2.11). A. veronii upregulated 177 inflammatory genes and downregulated 52 protein-coding genes affecting chromatin assembly, multiple small nuclear RNAs, multiple nucleolar RNAs, and 55 cytoplasmic tRNAs in HT-29 cells. These downregulation effects were unique to A. veronii and not observed in HT-29 cells infected with E. coli K12. A. veronii induced intestinal epithelial apoptosis involving the intrinsic pathway. A. veronii caused epithelial microvilli shortening and damage and epithelial production of IL-8. In conclusion, this study for the first time reports the association between IBD and Aeromonas enteric infection detected by bacterial cultivation. This study also reports that A. veronii damages intestinal epithelial cells via multiple mechanisms, of which the downregulating cytoplasmic tRNA, small nuclear RNA, and small nucleolar RNA are novel bacterial pathogenic mechanisms. IMPORTANCE This study for the first time reports the association between inflammatory bowel disease (IBD) and Aeromonas enteric infection detected by bacterial pathogen cultivation, highlighting the need of clinical and public health attention. The finding that patients with IBD are more susceptible to Aeromonas enteric infection suggests that detection of Aeromonas enteric infection should be routinely performed for the diagnosis and treatment of IBD. This study also reports novel bacterial pathogenic mechanisms employed by Aeromonas veronii. Through comparative transcriptomic analysis and other techniques, this study revealed the pathogenic mechanisms by which A. veronii causes damage to intestinal epithelial cells. Among the various pathogenic mechanisms identified, the downregulating tRNA, small nuclear and nucleolar RNAs in human intestinal epithelial cells are novel bacterial pathogenic mechanisms.

Immunological effects of recombinant Lactobacillus casei expressing pilin MshB fused with cholera toxin B subunit adjuvant as an oral vaccine against Aeromonas veronii infection in crucian carp

Aeromonas veronii is a zoonotic agent capable of infecting fish and mammals, including humans, posing a serious threat to the development of aquaculture and public health safety. Currently, few effective vaccines are available through convenient routes against A. veronii infection. Herein, we developed vaccine candidates by inserting MSH type VI pili B (MshB) from A. veronii as an antigen and cholera toxin B subunit (CTB) as a molecular adjuvant into Lactobacillus casei and evaluated their immunological effect as vaccines in a crucian carp (Carassius auratus) model. The results suggested that recombinant L. casei Lc-pPG-MshB and Lc-pPG-MshB-CTB can be stably inherited for more than 50 generations. Oral administration of recombinant L. casei vaccine candidates stimulated the production of high levels of serum-specific immunoglobulin M (IgM) and increased the activity of acid phosphatase (ACP), alkaline phosphatase (AKP) superoxide dismutase (SOD), lysozyme (LZM), complement 3 (C3) and C4 in crucian carp (carassius auratus) compared to the control group (Lc-pPG612 group and PBS group) without significant changes. Moreover, the expression levels of interleukin-10 (IL-10), interleukin-1β (IL-1β), tumour necrosis factor-α (TNF-α) and transforming growth factor-β (TGF-β) genes in the gills, liver, spleen, kidney and gut of crucian carp orally immunized with recombinant L. casei were significantly upregulated compared to the control groups, indicating that recombinant L. casei induced a significant cellular immune response. In addition, viable recombinant L. casei can be detected and stably colonized in the intestine tract of crucian carp. Particularly, crucian carp immunized orally with Lc-pPG-MshB and Lc-pPG-MshB-CTB exhibited higher survival rates (48% for Lc-pPG-MshB and 60% for Lc-pPG-MshB-CTB) and significantly reduced loads of A. veronii in the major immune organs after A. veronii challenge. Our findings indicated that both recombinant L. casei strains provide favorable immune protection, with Lc-pPG-MshB-CTB in particular being more effective and promising as an ideal candidate for oral vaccination.

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

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.

Bacterial diversity of middle ear cholesteatoma by 16S rRNA gene sequencing in China

In this study, the bacterial diversity of acquired middle ear cholesteatoma (MEC) was evaluated to reveal its pathogenesis and provides a guide for the use of antibiotics. Twenty-nine cases of acquired MEC and eight cases of healthy middle ears undergoing cochlear implantation (CI) were evaluated. Full-length 16S rRNA gene sequencing was performed to profile the bacterial communities in lesions and healthy tissues of the middle ear. ACE (P = 0.043) and Chao1 (P = 0.039) indices showed significant differences in alpha diversity (P < 0.05). Analysis of PERMANOVA/Anosim using the Bray-Curtis distance matrix results suggested that the between-group differences were greater than the within-group differences (R = 0.238, P < 0.05, R² = 0.066, P < 0.05). Bacterial community analysis revealed that Alphaproteobacteria at the class level and Caulobacterales and Sphingomonadales at the order level were significantly different (P < 0.05). In the LefSe (Linear discriminant analysis effect size) analysis, Porphyromonas bennonis was elevated, and Bryum argenteum and unclassified Cyanobacteriales were reduced at the species level in MEC (P < 0.05). Fifteen metabolic pathways were found to be significantly different between the two groups by analysing the abundance of metabolic pathways in level 2 of the Kyoto Encyclopaedia of Genes and Genomes (KEGG). Seven and eight metabolic pathways were significantly elevated in the MEC and control groups, respectively (P < 0.05). The role of bacteria in the pathogenesis of acquired MEC was further refined through analysis of metabolic pathways. These findings indicate that the acquired MEC and healthy middle ear contain more diverse microbial communities than previously thought.