Pathology, microbiology, and genetic diversity associated with Erysipelothrix rhusiopathiae and novel Erysipelothrix spp. infections in southern sea otters (Enhydra lutris nereis)

Erysipelothrix spp., including E. rhusiopathiae, are zoonotic bacterial pathogens that can cause morbidity and mortality in mammals, fish, reptiles, birds, and humans. The southern sea otter (SSO; Enhydra lutris nereis) is a federally-listed threatened species for which infectious disease is a major cause of mortality. We estimated the frequency of detection of these opportunistic pathogens in dead SSOs, described pathology associated with Erysipelothrix infections in SSOs, characterized the genetic diversity and antimicrobial susceptibility of SSO isolates, and evaluated the virulence of two novel Erysipelothrix isolates from SSOs using an in vivo fish model. From 1998 to 2021 Erysipelothrix spp. were isolated from six of >500 necropsied SSOs. Erysipelothrix spp. were isolated in pure culture from three cases, while the other three were mixed cultures. Bacterial septicemia was a primary or contributing cause of death in five of the six cases. Other pathology observed included suppurative lymphadenopathy, fibrinosuppurative arteritis with thrombosis and infarction, bilateral uveitis and endophthalmitis, hypopyon, petechia and ecchymoses, mucosal infarction, and suppurative meningoencephalitis and ventriculitis. Short to long slender Gram-positive or Gram-variable bacterial rods were identified within lesions, alone or with other opportunistic bacteria. All six SSO isolates had the spaA genotype-four isolates clustered with spaA E. rhusiopathiae strains from various terrestrial and marine animal hosts. Two isolates did not cluster with any known Erysipelothrix spp.; whole genome sequencing revealed a novel Erysipelothrix species and a novel E. rhusiopathiae subspecies. We propose the names Erysipelothrix enhydrae sp. nov. and Erysipelothrix rhusiopathiae ohloneorum ssp. nov. respectively. The type strains are E. enhydrae UCD-4322-04 and E. rhusiopathiae ohloneorum UCD-4724-06, respectively. Experimental injection of tiger barbs (Puntigrus tetrazona) resulted in infection and mortality from the two novel Erysipelothrix spp. Antimicrobial susceptibility testing of Erysipelothrix isolates from SSOs shows similar susceptibility profiles to isolates from other terrestrial and aquatic animals. This is the first description of the pathology, microbial characteristics, and genetic diversity of Erysipelothrix isolates recovered from diseased SSOs. Methods presented here can facilitate case recognition, aid characterization of Erysipelothrix isolates, and illustrate assessment of virulence using fish models.

Characterization of Type VI secretion system in Edwardsiella ictaluri

Edwardsiella ictaluri is a Gram-negative facultative intracellular fish pathogen causing enteric septicemia of catfish (ESC). While various secretion systems contribute to E. ictaluri virulence, the Type VI secretion system (T6SS) remains poorly understood. In this study, we constructed 13 E. ictaluri T6SS mutants using splicing by overlap extension PCR and characterized them, assessing their uptake and survival in channel catfish (Ictalurus punctatus) peritoneal macrophages, attachment and invasion in channel catfish ovary (CCO) cells, in vitro stress resistance, and virulence and efficacy in channel catfish. Among the mutants, EiΔevpA, EiΔevpH, EiΔevpM, EiΔevpN, and EiΔevpO exhibited reduced replication inside peritoneal macrophages. EiΔevpM, EiΔevpN, and EiΔevpO showed significantly decreased attachment to CCO cells, while EiΔevpN and EiΔevpO also displayed reduced invasion of CCO cells (p < 0.05). Overall, T6SS mutants demonstrated enhanced resistance to oxidative and nitrosative stress in the nutrient-rich medium compared to the minimal medium. However, EiΔevpA, EiΔevpH, EiΔevpM, EiΔevpN, and EiΔevpO were susceptible to oxidative stress in both nutrient-rich and minimal medium. In fish challenges, EiΔevpD, EiΔevpE, EiΔevpG, EiΔevpJ, and EiΔevpK exhibited attenuation and provided effective protection against E. ictaluri wild-type (EiWT) infection in catfish fingerlings. However, their attenuation and protective efficacy were lower in catfish fry. These findings shed light on the role of the T6SS in E. ictaluri pathogenesis, highlighting its significance in intracellular survival, host cell attachment and invasion, stress resistance, and virulence. The attenuated T6SS mutants hold promise as potential candidates for protective immunization strategies in catfish fingerlings.

Evaluation of Edwardsiella piscicida basS and basR mutants as vaccine candidates in catfish against edwardsiellosis

Catfish farming is the largest aquaculture industry in the United States and an important economic driver in several southeastern states. Edwardsiella piscicida is a Gram-negative pathogen associated with significant losses in catfish aquaculture. Several Gram-negative bacteria use the BasS/BasR two-component system (TCS) to adapt to environmental changes and the host immune system. Currently, the role of BasS/BasR system in E. piscicida virulence has not been characterized. In the present study, two mutants were constructed by deleting the basS and basR genes in E. piscicida strain C07-087. Both mutant strains were characterized for virulence and immune protection in catfish hosts. The EpΔbasS and EpΔbasR mutants were more sensitive to acidic environments and produced significantly less biofilm than the wild-type. In vivo studies in channel catfish (Ictalurus punctatus) revealed that both EpΔbasS and EpΔbasR were significantly attenuated compared with the parental wild-type (3.57% and 4.17% vs. 49.16% mortalities). Moreover, there was significant protection, 95.2% and 92.3% relative percent survival (RPS), in channel catfish vaccinated with EpΔbasS and EpΔbasR against E. piscicida infection. Protection in channel catfish was associated with a significantly higher level of antibodies and upregulation of immune-related genes (IgM, IL-8 and CD8-α) in channel catfish vaccinated with EpΔbasS and EpΔbasR strains compared with non-vaccinated fish. Hybrid catfish (channel catfish ♀ × blue catfish ♂) challenges demonstrated long-term protection against subsequent challenges with E. piscicida and E. ictaluri. Our findings demonstrate BasS and BasR contribute to acid tolerance and biofilm formation, which may facilitate E. piscicida survival in harsh environments. Further, our results show that EpΔbasS and EpΔbasR mutants were safe and protective in channel catfish fingerlings, although their virulence and efficacy in hybrid catfish warrant further investigation. These data provide information regarding an important mechanism of E. piscicida virulence, and it suggests EpΔbasS and EpΔbasR strains have potential as vaccines against this emergent catfish pathogen.

Comparative genomics of the fish pathogens Edwardsiella ictaluri 93-146 and Edwardsiella piscicida C07-087

Edwardsiella ictaluri and Edwardsiella piscicida are important fish pathogens affecting cultured and wild fish worldwide. To investigate the genome-level differences and similarities between catfish-adapted strains in these two species, the complete E. ictaluri 93-146 and E. piscicida C07-087 genomes were evaluated by applying comparative genomics analysis. All available complete (10) and non-complete (19) genomes from five Edwardsiella species were also included in a systematic analysis. Average nucleotide identity and core-genome phylogenetic tree analyses indicated that the five Edwardsiella species were separated from each other. Pan-/core-genome analyses for the 29 strains from the five species showed that genus Edwardsiella members have 9474 genes in their pan genome, while the core genome consists of 1421 genes. Orthology cluster analysis showed that E. ictaluri and E. piscicida genomes have the greatest number of shared clusters. However, E. ictaluri and E. piscicida also have unique features; for example, the E. ictaluri genome encodes urease enzymes and cytochrome o ubiquinol oxidase subunits, whereas E. piscicida genomes encode tetrathionate reductase operons, capsular polysaccharide synthesis enzymes and vibrioferrin-related genes. Additionally, we report for what is believed to be the first time that E. ictaluri 93-146 and three other E. ictaluri genomes encode a type IV secretion system (T4SS), whereas none of the E. piscicida genomes encode this system. Additionally, the E. piscicida C07-087 genome encodes two different type VI secretion systems. E. ictaluri genomes tend to encode more insertion elements, phage regions and genomic islands than E. piscicida. We speculate that the T4SS could contribute to the increased number of mobilome elements in E. ictaluri compared to E. piscicida. Two of the E. piscicida genomes encode full CRISPR-Cas regions, whereas none of the E. ictaluri genomes encode Cas proteins. Overall, comparison of the E. ictaluri and E. piscicida genomes reveals unique features and provides new insights on pathogenicity that may reflect the host adaptation of the two species.

Identification of Antimicrobial Resistance Determinants in Aeromonas veronii Strain MS-17-88 Recovered From Channel Catfish (Ictalurus punctatus)

Aeromonas veronii is a Gram-negative species ubiquitous in different aquatic environments and capable of causing a variety of diseases to a broad host range. Aeromonas species have the capability to carry and acquire antimicrobial resistance (AMR) elements, and currently multi-drug resistant (MDR) Aeromonas isolates are commonly found across the world. A. veronii strain MS-17-88 is a MDR strain isolated from catfish in the southeastern United States. The present study was undertaken to uncover the mechanism of resistance in MDR A. veronii strain MS-17-88 through the detection of genomic features. To achieve this, genomic DNA was extracted, sequenced, and assembled. The A. veronii strain MS-17-88 genome comprised 5,178,226-bp with 58.6% G+C, and it encoded several AMR elements, including imiS, ampS, mcr-7.1, mcr-3, catB2, catB7, catB1, floR, vat(F), tet(34), tet(35), tet(E), dfrA3, and tetR. The phylogeny and resistance profile of a large collection of A. veronii strains, including MS-17-88, were evaluated. Phylogenetic analysis showed a close relationship between MS-17-88 and strain Ae5 isolated from fish in China and ARB3 strain isolated from pond water in Japan, indicating a common ancestor of these strains. Analysis of phage elements revealed 58 intact, 63 incomplete, and 15 questionable phage elements among the 53 A. veronii genomes. The average phage element number is 2.56 per genome, and strain MS-17-88 is one of two strains having the maximum number of identified prophage elements (6 elements each). The profile of resistance against various antibiotics across the 53 A. veronii genomes revealed the presence of tet(34), mcr-7.1, mcr-3, and dfrA3 in all genomes (100%). By comparison, sul1 and sul2 were detected in 7.5% and 1.8% of A. veronii genomes. Nearly 77% of strains carried tet(E), and 7.5% of strains carried floR. This result suggested a low abundance and prevalence of sulfonamide and florfenicol resistance genes compared with tetracycline resistance among A. veronii strains. Overall, the present study provides insights into the resistance patterns among 53 A. veronii genomes, which can inform therapeutic options for fish affected by A. veronii.

Genomic diversity in flavobacterial pathogens of aquatic origin

Flavobacterium species are considered important fish pathogens in wild and cultured fish throughout the world. They can cause acute, subacute, and chronic infections, which are mainly characterized by gill damage, skin lesions, and deep necrotic ulcerations. Primarily, three Flavobacterium species, F. branchiophilum, F. columnare, and F. psychrophilum, have been reported to cause substantial losses to freshwater fish. In this study, we evaluated genomes of 86 Flavobacterium species isolated from aquatic hosts (mainly fish) to identify their unique and shared genome features. Our results showed that F. columnare genomes cluster into four different genetic groups. In silico secretion system analysis identified that all genomes carry type I (T1SS) and type IX (T9SS) secretion systems, but the number of type I secretion system genes shows diversity between species. F. branchiophilum, F. araucananum, F. chilense, F. spartansii, and F. tructae genomes have full type VI secretion system (T6SS). F. columnare, F. hydatis, and F. plurextorum carry partial T6SS with some of the T6SS genes missing. F. columnare, F. araucananum, F. chilense, F. spartansii, F. araucananum, F. tructae, Flavobacterium sp., F. crassostreae, F. succinicans, F. hydatis, and F. plurextorum carry most of the type IV secretion system genes (T4SS). F. columnare genetic groups 1 and 2, Flavobacterium sp., and F. crassostreae encode the least number of antibiotic resistance elements. F. hydatis, F. chilense, and F. plurextorum encode the greatest number of antibiotic resistance genes. Additionally, F. spartansii, F. araucananum, and chilense encode the greatest number of virulence genes while Flavobacterium sp. and F. crassostreae encode the least number of virulence genes. In conclusion, comparative genomics of Flavobacterium species of aquatic origin will help our understanding of Flavobacterium pathogenesis.

Comparative genomics of Aeromonas veronii: Identification of a pathotype impacting aquaculture globally

Aeromonas veronii is a gram-negative species abundant in aquatic environments that causes disease in humans as well as terrestrial and aquatic animals. In the current study, 41 publicly available A. veronii genomes were compared to investigate distribution of putative virulence genes, global dissemination of pathotypes, and potential mechanisms of virulence. The complete genome of A. veronii strain ML09-123 from an outbreak of motile aeromonas septicemia in farm-raised catfish in the southeastern United States was included. Dissemination of A. veronii strain types was discovered in dispersed geographical locations. Isolate ML09-123 is highly similar to Chinese isolate TH0426, suggesting the two strains have a common origin and may represent a pathotype impacting aquaculture in both countries. Virulence of strain ML09-123 in catfish in a dose-dependent manner was confirmed experimentally. Subsystem category disposition showed the majority of genomes exhibit similar distribution of genomic elements. The type I secretion system (T1SS), type II secretion system (T2SS), type 4 pilus (T4P), and flagellum core elements are conserved in all A. veronii genomes, whereas the type III secretion system (T3SS), type V secretion system (T5SS), type VI secretion system (T6SS), and tight adherence (TAD) system demonstrate variable dispersal. Distribution of mobile elements is dependent on host and geographic origin, suggesting this species has undergone considerable genetic exchange. The data presented here lends insight into the genomic variation of A. veronii and identifies a pathotype impacting aquaculture globally.

Transposon mutagenesis and identification of mutated genes in growth-delayed Edwardsiella ictaluri

Background

Edwardsiella ictaluri is a Gram-negative facultative intracellular anaerobe and the etiologic agent of enteric septicemia of channel catfish (ESC). To the catfish industry, ESC is a devastating disease due to production losses and treatment costs. Identification of virulence mechanisms of E. ictaluri is critical to developing novel therapeutic approaches for the disease. Here, we report construction of a transposon insertion library and identification of mutated genes in growth-delayed E. ictaluri colonies. We also provide safety and efficacy of transposon insertion mutants in catfish.

Results

An E. ictaluri transposon insertion library with 45,000 transposants and saturating 30.92% of the TA locations present in the E. ictaluri genome was constructed. Transposon end mapping of 250 growth-delayed E. ictaluri colonies and bioinformatic analysis of sequences revealed 56 unique E. ictaluri genes interrupted by the MAR2xT7 transposon, which are involved in metabolic and cellular processes and mostly localized in the cytoplasm or cytoplasmic membrane. Of the 56 genes, 30 were associated with bacterial virulence. Safety and vaccine efficacy testing of 19 mutants showed that mutants containing transposon insertions in hypothetical protein (Eis::004), and Fe-S cluster assembly protein (IscX, Eis::039), sulfurtransferase (TusA, Eis::158), and universal stress protein A (UspA, Eis::194) were safe and provided significant protection (p < 0.05) against wild-type E. ictaluri.

Conclusions

The results indicate that random transposon mutagenesis causing growth-delayed phenotype results in identification bacterial virulence genes, and attenuated strains with transposon interrupted virulence genes could be used as vaccine to activate fish immune system.

Comparative Genomics of Aeromonas hydrophila Secretion Systems and Mutational Analysis of hcp1 and vgrG1 Genes From T6SS

Virulent Aeromonas hydrophila causes severe motile Aeromonas septicemia in warmwater fishes. In recent years, channel catfish farming in the U.S.A. and carp farming in China have been affected by virulent A. hydrophila, and genome comparisons revealed that these virulent A. hydrophila strains belong to the same clonal group. Bacterial secretion systems are often important virulence factors; in the current study, we investigated whether secretion systems contribute to the virulent phenotype of these strains. Thus, we conducted comparative secretion system analysis using 55 A. hydrophila genomes, including virulent A. hydrophila strains from U.S.A. and China. Interestingly, tight adherence (TaD) system is consistently encoded in all the vAh strains. The majority of U.S.A. isolates do not possess a complete type VI secretion system, but three core elements [tssD (hcp), tssH, and tssI (vgrG)] are encoded. On the other hand, Chinese isolates have a complete type VI secretion system operon. None of the virulent A. hydrophila isolates have a type III secretion system. Deletion of two genes encoding type VI secretion system proteins (hcp1 and vgrG1) from virulent A. hydrophila isolate ML09-119 reduced virulence 2.24-fold in catfish fingerlings compared to the parent strain ML09-119. By determining the distribution of genes encoding secretion systems in A. hydrophila strains, our study clarifies which systems may contribute to core A. hydrophila functions and which may contribute to more specialized adaptations such as virulence. Our study also clarifies the role of type VI secretion system in A. hydrophila virulence.

Universal Stress Proteins Contribute Edwardsiella ictaluri Virulence in Catfish

Edwardsiella ictaluri is an intracellular Gram-negative facultative pathogen causing enteric septicemia of catfish (ESC), a common disease resulting in substantial economic losses in the U.S. catfish industry. Previously, we demonstrated that several universal stress proteins (USPs) are highly expressed under in vitro and in vivo stress conditions, indicating their importance for E. ictaluri survival. However, the roles of these USPs in E. ictaluri virulence is not known yet. In this work, 10 usp genes of E. ictaluri were in-frame deleted and characterized in vitro and in vivo. Results show that all USP mutants were sensitive to acidic condition (pH 5.5), and EiΔusp05 and EiΔusp08 were very sensitive to oxidative stress (0.1% H₂O₂). Virulence studies indicated that EiΔusp05, EiΔusp07, EiΔusp08, EiΔusp09, EiΔusp10, and EiΔusp13 were attenuated significantly compared to E. ictaluri wild-type (EiWT; 20, 45, 20, 20, 55, and 10% vs. 74.1% mortality, respectively). Efficacy experiments showed that vaccination of catfish fingerlings with EiΔusp05, EiΔusp07, EiΔusp08, EiΔusp09, EiΔusp10, and EiΔusp13 provided complete protection against EiWT compared to sham-vaccinated fish (0% vs. 58.33% mortality). Our results support that USPs contribute E. ictaluri virulence in catfish.

Comparative Analysis of the Flavobacterium columnare Genomovar I and II Genomes

Columnaris disease caused by Gram-negative rod Flavobacterium columnare is one of the most common diseases of catfish. F. columnare is also a common problem in other cultured fish species worldwide. F. columnare has three major genomovars; we have sequenced a representative strain from genomovar I (ATCC 49512, which is avirulent in catfish) and genomovar II (94-081, which is highly pathogenic in catfish). Here, we present a comparative analysis of the two genomes. Interestingly, F. columnare ATCC 49512 and 94-081 meet criteria to be considered different species based on the Average Nucleotide Identity (90.71% similar) and DNA–DNA Hybridization (42.6% similar). Genome alignment indicated the two genomes have a large number of rearrangements. However, function-based comparative genomics analysis indicated that the two strains have similar functional capabilities with 2,263 conserved orthologous clusters; strain ATCC 49512 has 290 unique orthologous clusters while strain 94-081 has 391. Both strains carry type I secretion system, type VI secretion system, and type IX secretion system. The two genomes also have similar CRISPR capacities. The F. columnare strain ATCC 49512 genome contains a higher number of insertion sequence families and phage regions, while the F. columnare strain 94-081 genome has more genomic islands and more regulatory gene capacity. Transposon mutagenesis using Tn4351 in pathogenic strain 94-081 yielded six mutants, and experimental infections of fish showed hemolysin and glycine cleavage protein mutants had 15 and 10% mortalities, respectively, while the wild-type strain caused 100% mortalities. Our comparative and mutational analysis yielded important information on classification of genomovars I and II F. columnare as well as potential virulence genes in F. columnare strain 94-081.

Comparative Genomics and Transcriptional Analysis of Flavobacterium columnare Strain ATCC 49512

Flavobacterium columnare is a Gram-negative fish pathogen causing columnaris disease in wild and cultured fish species. Although the pathogen is widespread in aquatic environments and fish worldwide, little is known about biology of F. columnare and mechanisms of columnaris disease pathogenesis. Previously we presented the complete genome sequence of F. columnare strain ATCC 49512. Here we present a comparison of the strain ATCC 49512 genome to four other Flavobacterium genomes. In this analysis, we identified predicted proteins whose functions indicate F. columnare is capable of denitrification, which would enable anaerobic growth in aquatic pond sediments. Anaerobic growth of F. columnare ATCC 49512 with nitrate supplementation was detected experimentally. F. columnare ATCC 49512 had a relatively high number of insertion sequences and genomic islands compared to the other Flavobacterium species, suggesting a larger degree of horizontal gene exchange and genome plasticity. A type VI subtype III secretion system was encoded in F. columnare along with F. johnsoniae and F. branchiophilum. RNA sequencing proved to be a valuable technique to improve annotation quality; 41 novel protein coding regions were identified, 16 of which had a non-traditional start site (TTG, GTG, and CTT). Candidate small noncoding RNAs were also identified. Our results improve our understanding of F. columnare ATCC 49512 biology, and our results support the use of RNA sequencing to improve annotation of bacterial genomes, particularly for type strains.