Phylogenetic analysis of the fish pathogen Aeromonas salmonicida underlines the dichotomy between European and Canadian strains for the salmonicida subspecies

Aeromonas salmonicida isolates from Canada demonstrate wide distribution and clustering among mesophilic strains

All the 36 known species to date of the genus Aeromonas are mesophilic except the species Aeromonas salmonicida, which includes both psychrophilic and mesophilic subspecies. For 20 years, more and more mesophilic A. salmonicida strains have been discovered. Only A. salmonicida subsp. pectinolytica has officially been classified as a mesophilic subspecies. Most mesophiles have been isolated in hot countries. We present, for the first time, the characterization of two new mesophilic isolates from Quebec (Canada). Phenotypic and genomic characterizations were carried out on these strains, isolated from dead fish from a fish farm. Isolates 19-K304 and 19-K308 are clearly mesophiles, virulent to the amoeba Dictyostelium discoideum, a surrogate host, and close to strain Y577, isolated in India. To our knowledge, this is the first time that mesophilic strains isolated from different countries are so similar. The major difference between the isolates is the presence of plasmid pY47-3, a cryptic plasmid that sometimes presents in mesophilic strains. More importantly, our extensive phylogenetic analysis reveals two well-defined clades of mesophilic strains with psychrophiles associated with one of these clades. This helps to have a better understanding of the evolution of this species and the apparition of psychrophilic subspecies.

Comparative pan-genomic analysis of 51 Renibacterium salmoninarum indicates heterogeneity in the principal virulence factor, the 57 kDa protein

Renibacterium salmoninarum, a Gram-positive intracellular pathogen, is the causative agent of bacterial kidney disease (BKD), the impacts of which are high mortalities and economic losses for the salmon industry. This study provides novel analyses for the whole-genome sequences of 50 R. salmoninarum isolates and the reference strain ATCC 33209 using a pan-genomic approach to elucidate phylogenomic relationships and identify unique and shared genes associated with pathogenicity and infection mechanisms. Genome size varied from 3,061,638 to 3,155,332 bp; gene count from 3452 to 3580; and predicted coding sequences from 3402 to 3527. Comparative analyses revealed an open, but approaching closed, pan-genome. The pan-genome analysis recovered 4064 genes, with a core genome containing 3306 genes. Phylogenetic analysis of R. salmoninarum showed high genomic homogeneity, apart from one isolate obtained from Salmo trutta in Norway. All genomes presented the 57-kDa protein (p57). Strain ATCC 33209 and the Chilean isolates H-2 and DJ2R presented two copies of the msa gene, while the remaining isolates had one copy. The pan-genome analysis further identified differences in the number of copies and length of the signalling peptide for p57, the principal virulence factor reported for this bacterium. This heterogeneity could be associated with the secretion levels of p57, potentially influencing virulence. Additionally identified were numerous common genes related to iron uptake, the stress response and regulation, and cell signalling-all of which constitute the pathogenic repertoire of R. salmoninarum. This investigation provides information that is applicable in future studies for identifying therapeutic targets and/or for designing new strategies (e.g., vaccines) to prevent BKD infections in salmon farming.

The Pseudomonas aeruginosa Pan-Genome Provides New Insights on Its Population Structure, Horizontal Gene Transfer, and Pathogenicity

The huge increase in the availability of bacterial genomes led us to a point in which we can investigate and query pan-genomes, for example, the full set of genes of a given bacterial species or clade. Here, we used a data set of 1,311 high-quality genomes from the human pathogen Pseudomonas aeruginosa, 619 of which were newly sequenced, to show that a pan-genomic approach can greatly refine the population structure of bacterial species, provide new insights to define species boundaries, and generate hypotheses on the evolution of pathogenicity. The 665-gene P. aeruginosa core genome presented here, which constitutes only 1% of the entire pan-genome, is the first to be in the same order of magnitude as the minimal bacterial genome and represents a conservative estimate of the actual core genome. Moreover, the phylogeny based on this core genome provides strong evidence for a five-group population structure that includes two previously undescribed groups of isolates. Comparative genomics focusing on antimicrobial resistance and virulence genes showed that variation among isolates was partly linked to this population structure. Finally, we hypothesized that horizontal gene transfer had an important role in this respect, and found a total of 3,010 putative complete and fragmented plasmids, 5% and 12% of which contained resistance or virulence genes, respectively. This work provides data and strategies to study the evolutionary trajectories of resistance and virulence in P. aeruginosa.

Study of mesophilic Aeromonas salmonicida A527 strain sheds light on the species' lifestyles and taxonomic dilemma

The Gram-negative bacterium Aeromonas salmonicida contains five subspecies: salmonicida, smithia, achromogenes, masoucida and pectinolytica. Pectinolytica is a mesophilic subspecies with the ability to thrive at a wide range of temperatures, including 37°C, while the four other subspecies are psychrophilic, restricted to lower temperatures. The psychrophilic subspecies are known to infect a wide range of fishes. However, there is no evidence of pathogenicity for the mesophilic subspecies pectinolytica. Study of the differences between the mesophilic and psychrophilic subspecies is hampered by the lack of completely sequenced and closed genomes from the mesophilic subspecies. A previous study reported that insertion sequences, which can induce genomic rearrangements at temperatures around 25°C, could be one of the determinants explaining the differences in lifestyle (mesophilic or psychrophilic) between the subspecies. In this study, the genome of mesophilic strain A527 of A. salmonicida was sequenced, closed and analyzed to investigate the mesophilic-psychrophilic discrepancy. This reference genome supports the hypothesis that insertion sequences are major determinants of the lifestyle differences between the A. salmonicida subspecies. Moreover, the phylogenetic analysis performed to position strain A527 within the taxonomy raises an issue regarding the intraspecies structure of A. salmonicida.

Strong Genomic and Phenotypic Heterogeneity in the Aeromonas sobria Species Complex

Aeromonas sobria is a mesophilic motile aeromonad currently depicted as an opportunistic pathogen, despite increasing evidence of mutualistic interactions in salmonid fish. However, the determinants of its host-microbe associations, either mutualistic or pathogenic, remain less understood than for other aeromonad species. On one side, there is an over-representation of pathogenic interactions in the A. sobria literature, of which only three articles to date report mutualistic interactions; on the other side, genomic characterization of this species is still fairly incomplete as only two draft genomes were published prior to the present work. Consequently, no study specifically investigated the biodiversity of A. sobria. In fact, the investigation of A. sobria as a species complex may have been clouded by: (i) confusion with A. veronii biovar sobria because of their similar biochemical profiles, and (ii) the intrinsic low resolution of previous studies based on 16S rRNA gene sequences and multilocus sequence typing. So far, the only high-resolution, phylogenomic studies of the genus Aeromonas included one A. sobria strain (CECT 4245 / Popoff 208), making it impossible to robustly conclude on the phylogenetic intra-species diversity and the positioning among other Aeromonas species. To further understand the biodiversity and the spectrum of host-microbe interactions in A. sobria as well as its potential genomic diversity, we assessed the genomic and phenotypic heterogeneity among five A. sobria strains: two clinical isolates recovered from infected fish (JF2635 and CECT 4245), one from an infected amphibian (08005) and two recently isolated brook charr probionts (TM12 and TM18) which inhibit in vitro growth of A. salmonicida subsp. salmonicida (a salmonid fish pathogen). A phylogenomic assessment including 2,154 softcore genes corresponding to 946,687 variable sites from 33 Aeromonas genomes confirms the status of A. sobria as a distinct species divided in two subclades, with 100% bootstrap support. The phylogenomic split of A. sobria in two subclades is corroborated by a deep dichotomy between all five A. sobria strains in terms of inhibitory effect against A. salmonicida subsp. salmonicida, gene contents and codon usage. Finally, the antagonistic effect of A. sobria strains TM12 and TM18 suggests novel control methods against A. salmonicida subsp. salmonicida.

Comparative Pan-Genome Analysis of Piscirickettsia salmonis Reveals Genomic Divergences within Genogroups

Piscirickettsia salmonis is the etiological agent of salmonid rickettsial septicemia, a disease that seriously affects the salmonid industry. Despite efforts to genomically characterize P. salmonis, functional information on the life cycle, pathogenesis mechanisms, diagnosis, treatment, and control of this fish pathogen remain lacking. To address this knowledge gap, the present study conducted an in silico pan-genome analysis of 19 P. salmonis strains from distinct geographic locations and genogroups. Results revealed an expected open pan-genome of 3,463 genes and a core-genome of 1,732 genes. Two marked genogroups were identified, as confirmed by phylogenetic and phylogenomic relationships to the LF-89 and EM-90 reference strains, as well as by assessments of genomic structures. Different structural configurations were found for the six identified copies of the ribosomal operon in the P. salmonis genome, indicating translocation throughout the genetic material. Chromosomal divergences in genomic localization and quantity of genetic cassettes were also found for the Dot/Icm type IVB secretion system. To determine divergences between core-genomes, additional pan-genome descriptions were compiled for the so-termed LF and EM genogroups. Open pan-genomes composed of 2,924 and 2,778 genes and core-genomes composed of 2,170 and 2,228 genes were respectively found for the LF and EM genogroups. The core-genomes were functionally annotated using the Gene Ontology, KEGG, and Virulence Factor databases, revealing the presence of several shared groups of genes related to basic function of intracellular survival and bacterial pathogenesis. Additionally, the specific pan-genomes for the LF and EM genogroups were defined, resulting in the identification of 148 and 273 exclusive proteins, respectively. Notably, specific virulence factors linked to adherence, colonization, invasion factors, and endotoxins were established. The obtained data suggest that these genes could be directly associated with inter-genogroup differences in pathogenesis and host-pathogen interactions, information that could be useful in designing novel strategies for diagnosing and controlling P. salmonis infection.