Comparative Genomics and Transcriptional Analysis of Flavobacterium columnare Strain ATCC 49512

Gliding motility proteins GldJ and SprB contribute to Flavobacterium columnare virulence

Flavobacterium columnare causes columnaris disease in fish. Columnaris disease is incompletely understood, and adequate control measures are lacking. The type IX secretion system (T9SS) is required for F. columnare gliding motility and virulence. The T9SS and gliding motility machineries share some, but not all, components. GldN (required for gliding and for secretion) and PorV (involved in secretion but not required for gliding) are both needed for virulence, implicating T9SS-mediated secretion in virulence. The role of motility in virulence is uncertain. We constructed and analyzed sprB, sprF, and gldJ mutants that were defective for motility but that maintained T9SS function to understand the role of motility in virulence. Wild-type cells moved rapidly and formed spreading colonies. In contrast, sprB and sprF deletion mutants were partially defective in gliding and formed nonspreading colonies. Both mutants exhibited reduced virulence in rainbow trout fry. A gldJ deletion mutant was nonmotile, secretion deficient, and avirulent in rainbow trout fry. To separate the roles of GldJ in secretion and in motility, we generated gldJ truncation mutants that produce nearly full-length GldJ. Mutant gldJ563, which produces GldJ truncated at amino acid 563, was defective for gliding but was competent for secretion as measured by extracellular proteolytic activity. This mutant displayed reduced virulence in rainbow trout fry, suggesting that motility contributes to virulence. Fish that survived exposure to the sprB deletion mutant or the gldJ563 mutant exhibited partial resistance to later challenge with wild-type cells. The results aid our understanding of columnaris disease and may suggest control strategies.IMPORTANCEFlavobacterium columnare causes columnaris disease in many species of freshwater fish in the wild and in aquaculture systems. Fish mortalities resulting from columnaris disease are a major problem for aquaculture. F. columnare virulence is incompletely understood, and control measures are inadequate. Gliding motility and protein secretion have been suggested to contribute to columnaris disease, but evidence directly linking motility to disease was lacking. We isolated and analyzed mutants that were competent for secretion but defective for motility. Some of these mutants exhibited decreased virulence. Fish that had been exposed to these mutants were partially protected from later exposure to the wild type. The results contribute to our understanding of columnaris disease and may aid development of control strategies.

Complete genome sequence and potential pathogenic assessment of Flavobacterium plurextorum RSG-18 isolated from the gut of Schlegel's black rockfish, Sebastes schlegelii

Flavobacterium plurextorum is a potential fish pathogen of interest, previously isolated from diseased rainbow trout (Oncorhynchus mykiss) and oomycete-infected chum salmon (Oncorhynchus keta) eggs. We report here the first complete genome sequence of F. plurextorum RSG-18 isolated from the gut of Schlegel's black rockfish (Sebastes schlegelii). The genome of RSG-18 consists of a circular chromosome of 5,610,911 bp with a 33.57% GC content, containing 4858 protein-coding genes, 18 rRNAs, 63 tRNAs and 1 tmRNA. A comparative analysis was conducted on 11 Flavobacterium species previously reported as pathogens or isolated from diseased fish to confirm the potential pathogenicity of RSG-18. In the SEED classification, RSG-18 was found to have 36 genes categorized in 'Virulence, Disease and Defense'. Across all Flavobacterium species, a total of 16 antibiotic resistance genes and 61 putative virulence factors were identified. All species had at least one phage region and type I, III and IX secretion systems. In pan-genomic analysis, core genes consist of genes linked to phages, integrases and matrix-tolerated elements associated with pathology. The complete genome sequence of F. plurextorum RSG-18 will serve as a foundation for future research, enhancing our understanding of Flavobacterium pathogenicity in fish and contributing to the development of effective prevention strategies.

Genomic evidence of genetic diversity and functional evolution in Flavobacterium columnare

Flavobacterium columnare is the causative agent of columnaris disease in freshwater fish. Columnaris disease can cause heavy economic losses in aquaculture. In this study, whole-genome sequencing was used to characterize this pathogen. F. columnare isolate AH-01 had a circular chromosome and plasmid that encoded a total of 3,022 genes. Isolate GX-01 only had a circular chromosome and encoded 2,965 genes. Genomic islands, prophage regions, and CRISPR/Cas systems were identified in both genomes. Both genomes presented evidence of gene variation and horizontal transfer, both of which are the essential components of genetic diversity, genome plasticity, and functional evolution. Single-gene phylogeny and comparative genome analyses were performed to investigate the variation and evolution of this pathogen. Genetic analysis of 16S rRNA and housekeeping gene sequences significantly clustered 55 F. columnare isolates into four clades. The intragroup identity of the 16S rRNA gene exceeded 99%, while the intergroup identity was below the species delineation threshold. We discovered significant translocation, inversion, and rearrangement events that influenced local synteny within each group. Notably, the observed alignments varied considerably among all the studied groups. The core genomes of all strains with available sequences comprised 747 genes, corresponding to approximately 25% of the genome. Core genome multilocus sequence typing, genome-wide orthology and phylogenetic analyses, and average nucleotide identity suggested that the currently existing F. columnare was an assemblage of several distinct species, with levels of divergence at least equivalent to those between recognized bacterial species. The present investigation provided genomic evidence of gene variation and horizontal transfer, which were the basis of genetic diversity, genome plasticity, and functional evolution. The findings supported a proposed new taxonomic perspective on F. columnare.

Pangenome analysis provides insights into the genetic diversity, metabolic versatility, and evolution of the genus Flavobacterium

Members of the genus Flavobacterium are widely distributed and produce various polysaccharide-degrading enzymes. Many species in the genus have been isolated and characterized. However, few studies have focused on marine isolates or fish pathogens, and in-depth genomic analyses, particularly comparative analyses of isolates from different habitat types, are lacking. Here, we isolated 20 strains of the genus from various environments in South Korea and sequenced their full-length genomes. Combined with published sequence data, we examined genomic traits, evolution, environmental adaptation, and putative metabolic functions in total 187 genomes of isolated species in Flavobacterium categorized as marine, host-associated, and terrestrial including freshwater. A pangenome analysis revealed a correlation between genome size and coding or noncoding density. Flavobacterium spp. had high levels of diversity, allowing for novel gene repertories via recombination events. Defense-related genes only accounted for approximately 3% of predicted genes in all Flavobacterium genomes. While genes involved in metabolic pathways did not differ with respect to isolation source, there was substantial variation in genomic traits; in particular, the abundances of tRNAs and rRNAs were higher in the host-associdated group than in other groups. One genome in the host-associated group contained a Microviridae prophage closely related to an enterobacteria phage. The proteorhodopsin gene was only identified in four terrestrial strains isolated for this study. Furthermore, recombination events clearly influenced genomic diversity and may contribute to the response to environmental stress. These findings shed light on the high genetic variation in Flavobacterium and functional roles in diverse ecosystems as a result of their metabolic versatility. IMPORTANCE The genus Flavobacterium is a diverse group of bacteria that are found in a variety of environments. While most species of this genus are harmless and utilize organic substrates such as proteins and polysaccharides, some members may play a significant role in the cycling for organic substances within their environments. Nevertheless, little is known about the genomic dynamics and/or metabolic capacity of Flavobacterium. Here, we found that Flavobacterium species may have an open pangenome, containing a variety of diverse and novel gene repertoires. Intriguingly, we discovered that one genome (classified into host-associated group) contained a Microviridae prophage closely related to that of enterobacteria. Proteorhodopsin may be expressed under conditions of light or oxygen pressure in some strains isolated for this study. Our findings significantly contribute to the understanding of the members of the genus Flavobacterium diversity exploration and will provide a framework for the way for future ecological characterizations.

Investigating the Ability of Edwardsiella ictaluri and Flavobacterium covae to Persist within Commercial Catfish Pond Sediments under Laboratory Conditions

Two prevalent bacterial diseases in catfish aquaculture are enteric septicemia of catfish and columnaris disease caused by Edwardsiella ictaluri and Flavobacterium covae, respectively. Chronic and recurring outbreaks of these bacterial pathogens result in significant economic losses for producers annually. Determining if these pathogens can persist within sediments of commercial ponds is paramount. Experimental persistence trials (PT) were conducted to evaluate the persistence of E. ictaluri and F. covae in pond sediments. Twelve test chambers containing 120 g of sterilized sediment from four commercial catfish ponds were inoculated with either E. ictaluri (S97-773) or F. covae (ALG-00-530) and filled with 8 L of disinfected water. At 1, 2, 4-, 6-, 8-, and 15-days post-inoculation, 1 g of sediment was removed, and colony-forming units (CFU) were enumerated on selective media using 6 × 6 drop plate methods. E. ictaluri population peaked on Day 3 at 6.4 ± 0.5 log₁₀ CFU g^-1. Correlation analysis revealed no correlation between the sediment physicochemical parameters and E. ictaluri log₁₀ CFU g^-1. However, no viable F. covae colonies were recovered after two PT attempts. Future studies to improve understanding of E. ictaluri pathogenesis and persistence, and potential F. covae persistence in pond bottom sediments are needed.

Variations in antibiotic resistance genes and microbial community in sludges passing through biological nutrient removal and anaerobic digestion processes in municipal wastewater treatment plants

Antimicrobial resistance (AMR) represents a relentless, silent pandemic. Contributing to this are wastewater treatment plants (WWTPs), a potential source of antibiotic resistance genes' (ARGs) transmission to the environment, threatening public health. The presence of ARGs in pathogenic bacteria and their release into the environment by WWTPs threatens the public health. The current study investigated changes in ARGs' abundance in biological nutrient removal (BNR) processes and anaerobic digestion (AD) reactors of two WWTPs. Also, microbial community structure, which is known to shape the distribution and abundance of ARGs, was also analyzed. The relative abundance of eight ARGs (tetX, tetA, tetM, TEM, sul1, sul2, ermB and qnrD) was quantified as ARGs' copies/16 S rRNA gene copies using quantitative polymerase chain reaction (qPCR). Microbial community composition was assessed by 16 S rRNA microbiome sequencing analysis. TetX was prevalent among the eight ARGs, followed by TEM and sul1. However, its abundance was decreased in the AD sludges compared to BNR sludges. Proteobacteria was the major bacterial phylum found in all the sludge samples, while Arcobacter, 12up and Acidovorax were the predominant genera. Acinetobacter and Flavobacterium were significantly more abundant in the BNR sludges, while 12up and Aeromonas were predominant in AD sludges. Principal component analysis (PCA) revealed a clear difference in dominant ARGs and bacteria between the sludges in the processes of BNR and AD of the two WWTPs. Clinically relevant bacterial genera, Klebsiella and Enterococcus, found in both the BNR and AD sludges, were significantly correlated with the tetX gene. Throughout this study, the relationship between microbial communities and specific ARGs was revealed, illustrating that the composition of the microbial community could play a vital role in the abundance of ARGs. These results will better inform future studies aimed at controlling the spread of ARGs and their potential hosts from WWTPs.

Whole genome analysis of Flavobacterium aziz-sancarii sp. nov., isolated from Ardley Island (Antarctica), revealed a rich resistome and bioremediation potential

Despite being one of the most isolated regions in the world, Antarctica is at risk of increased contamination with potentially toxic elements and other toxic chemicals through anthropogenic interventions. In this study, a psychrotolerant bacterium was isolated using the lake water collected from Ardley Island (Antarctica), which can grow at temperatures between 4 and 30 °C and pH values between 6.0 and 9.0. The isolate, named AC, had protease, amylase, and lipase activities with no NaCl tolerance and could degrade 1-5% diesel fuel. Multilocus sequence analysis (MLSA) using 16S rRNA, gyrB, tuf, and rpoD genes resulted in 92.91-98.6% sequence similarities between the isolate AC and other Flavobacterium spp. Whole genome analysis indicated that the genome length of Flavobacterium sp. AC is 5.8 Mbp with a GC content of 34.04% and 1274 genes predicted. The strain AC branched independently from other Flavobacterium spp. in the phylogenetic and phylogenomic trees and ranked a new species named Flavobacterium aziz-sancarii. Genome mining identified several cold-inducible genes, including stress-associated genes such as cold-shock proteins, chaperones, carotenoid biosynthetic genes, or oxidative-stress response genes. In addition, virulence, gliding motility, and biofilm-related genes were determined. Its genome contains 35 and 88 open-reading frames related to potentially toxic element and antibiotic resistance, respectively. F. aziz-sancarii showed a remarkable tolerance of Cr and Ni, with minimal inhibitory concentration values of 2.88 and 2.81 mM, respectively. Pb, Cu, and Zn exposure resulted in moderate toxicity (2.14-2.41 mM), while Cd showed the highest inhibitory effect in bacterial growth (0.74 mM). Antibiotic susceptibility testing indicated multidrug-resistant phenotype in correlation to in silico prediction of antibiotic resistance genes. Overall, our results contribute to biodiversity of Antarctica and provide new insights into resistome profile of Antarctic microorganisms. Additionally, the diesel degradation feature of F. aziz-sancarii offers potential use for the bioremediation of hydrocarbon-contaminated polar ecosystems.

Flavobacterium columnare ferric iron uptake systems are required for virulence

Flavobacterium columnare, which causes columnaris disease, is one of the costliest pathogens in the freshwater fish-farming industry. The virulence mechanisms of F. columnare are not well understood and current methods to control columnaris outbreaks are inadequate. Iron is an essential nutrient needed for metabolic processes and is often required for bacterial virulence. F. columnare produces siderophores that bind ferric iron for transport into the cell. The genes needed for siderophore production have been identified, but other components involved in F. columnare iron uptake have not been studied in detail. We identified the genes encoding the predicted secreted heme-binding protein HmuY, the outer membrane iron receptors FhuA, FhuE, and FecA, and components of an ATP binding cassette (ABC) transporter predicted to transport ferric iron across the cytoplasmic membrane. Deletion mutants were constructed and examined for growth defects under iron-limited conditions and for virulence against zebrafish and rainbow trout. Mutants with deletions in genes encoding outer membrane receptors, and ABC transporter components exhibited growth defects under iron-limited conditions. Mutants lacking multiple outer membrane receptors, the ABC transporter, or HmuY retained virulence against zebrafish and rainbow trout mirroring that exhibited by the wild type. Some mutants predicted to be deficient in multiple steps of iron uptake exhibited decreased virulence. Survivors of exposure to such mutants were partially protected against later infection by wild-type F. columnare.

Siderophores Produced by the Fish Pathogen Flavobacterium columnare Strain MS-FC-4 Are Not Essential for Its Virulence

Flavobacterium columnare causes columnaris disease in wild and aquaculture-reared freshwater fish. F. columnare virulence mechanisms are not well understood, and current methods to control columnaris disease are inadequate. Iron acquisition from the host is important for the pathogenicity and virulence of many bacterial pathogens. F. columnare iron acquisition has not been studied in detail. We identified genes predicted to function in siderophore production for ferric iron uptake. Genes predicted to encode the proteins needed for siderophore synthesis, export, uptake, and regulation were deleted from F. columnare strain MS-FC-4. The mutants were examined for defects in siderophore production, for growth defects in iron-limited conditions, and for virulence against zebrafish and rainbow trout. Mutants lacking all siderophore activity were obtained. These mutants displayed growth defects when cultured under iron-limited conditions, but they retained virulence against zebrafish and rainbow trout similar to that exhibited by the wild type, indicating that the F. columnare MS-FC-4 siderophores are not required for virulence under the conditions tested. IMPORTANCE Columnaris disease, which is caused by Flavobacterium columnare, is a major problem for freshwater aquaculture. Little is known regarding F. columnare virulence factors, and control measures are limited. Iron acquisition mechanisms such as siderophores are important for virulence of other pathogens. We identified F. columnare siderophore biosynthesis, export, and uptake genes. Deletion of these genes eliminated siderophore production and resulted in growth defects under iron-limited conditions but did not alter virulence in rainbow trout or zebrafish. The results indicate that the F. columnare strain MS-FC-4 siderophores are not critical virulence factors under the conditions tested but may be important for survival under iron-limited conditions in natural aquatic environments or aquaculture systems.

Ubiquity and Diversity of Cold Adapted Denitrifying Bacteria Isolated From Diverse Antarctic Ecosystems

Nitrogen cycle has been poorly investigated in Antarctic ecosystems. In particular, how extreme conditions of low temperature, dryness, and high radiation select the microorganisms involved in the cycle is not yet understood. Denitrification is an important step in the nitrogen cycle in which nitrate is reduced stepwise to the gases NO, N2O, and N2. Denitrification is carried out by a wide group of microorganisms spread in the phylogenetic tree. The aim of this work was to isolate and characterize denitrifying bacteria present in different cold environments from Antarctica. Bacterial isolates were obtained from lake, meltwater, sea, glacier ice, ornithogenic soil, and penguin feces samples from King George Island, Fildes peninsula in the Antarctic. Samples were taken during the deicing season in five sampling campaigns. From all the samples we were able to isolate denitrifying strains. A total of 199 bacterial isolates with the capacity to grow in anaerobic mineral media reducing nitrate at 4°C were obtained. The characterization of the isolates by 16S rRNA gene sequence analysis showed a high predominance of the genus Pseudomonas, followed by Janthinobacterium, Flavobacterium, Psychrobacter, and Yersinia. Other minor genera detected were Cryobacterium, Iodobacter, Kaistella, and Carnobacterium. The capacity to denitrify was not previously described for most of the bacteria related to our isolates and in many of them denitrifying genes were not present suggesting the presence of new genes in this extreme environment. Our work demonstrates the ubiquity of denitrification in the Maritime Antarctica and gives important information linking denitrification at cold temperature with taxa in an unequivocal way.

Streambed microbial communities in the transition zone between groundwater and a first-order stream as impacted by bidirectional water exchange

The input of nitrate and other agricultural pollutants in higher-order streams largely derives from first-order streams. The streambed as the transition zone between groundwater and stream water has a decisive impact on the attenuation of such pollutants. This reactivity is not yet well understood for lower-order agricultural streams, which are often anthropogenically altered and lack the streambed complexity allowing for extensive hyporheic exchange. Reactive hot spots in such streambeds have been hypothesized as a function of hydrology, which controls the local gaining (groundwater exfiltration) or losing (infiltration) of stream water. However, streambed microbial communities and activities associated with such reactive zones remain mostly uncharted. In this study, sediments of a first-order agriculturally impacted stream in southern Germany were investigated. Along with a hydraulic dissection of distinct gaining and losing reaches of the stream, community composition and the abundance of bacterial communities in the streambed were investigated using PacBio long-read sequencing of bacterial 16S rRNA gene amplicons, and qPCR of bacterial 16S rRNA and denitrification genes (nirK and nirS). We show that bidirectional water exchange between groundwater and the stream represents an important control for sediment microbiota, especially for nitrate-reducing populations. Typical heterotrophic denitrifiers were most abundant in a midstream net losing section, while up- and downstream net gaining sections were associated with an enrichment of sulfur-oxidizing potential nitrate reducers affiliated with Sulfuricurvum and Thiobacillus spp. Dispersal-based community assembly was found to dominate such spots of groundwater exfiltration. Our results indicate a coupling of N- and S-cycling processes in the streambed of an agricultural first-order stream, and a prominent control of microbiology by hydrology and hydrochemistry in situ. Such detailed local heterogeneities in exchange fluxes and streambed microbiomes have not been reported to date, but seem relevant for understanding the reactivity of lower-order streams.

Detecting Flavobacterial Fish Pathogens in the Environment via High-Throughput Community Analysis

Diseases caused by the fish pathogens Flavobacterium columnare and Flavobacterium psychrophilum are major contributors of preventable losses in the aquaculture industry. The persistent and difficult-to-control infections caused by these bacteria make timely intervention and prophylactic elimination of pathogen reservoirs important measures to combat these disease-causing agents. In this study, we present two independent assays for detecting these pathogens in a range of environmental samples. Natural water samples were inoculated with F. columnare and F. psychrophilum over 5 orders of magnitude, and pathogen levels were detected using Illumina MiSeq sequencing and droplet digital PCR. Both detection methods accurately identified pathogen-positive samples and showed good agreement in quantifying each pathogen. Additionally, the real-world application of these approaches was demonstrated using environmental samples collected at a rainbow trout (Oncorhynchus mykiss) aquaculture facility. These results show that both methods can serve as useful tools for surveillance efforts in aquaculture facilities, where the early detection of these flavobacterial pathogens may direct preventative measures to reduce disease occurrence. IMPORTANCE Early detection of a deadly disease outbreak in a population can be the difference between mass mortality or mitigated effects. In the present study, we evaluated and compared two molecular techniques for detecting economically impactful aquaculture pathogens. We demonstrate that one of these techniques, 16S rRNA gene sequencing using Illumina MiSeq technology, provides the ability to accurately detect two freshwater fish pathogens, F. columnare and F. psychrophilum, while simultaneously profiling the native microbial community. The second technique, droplet digital PCR, is commonly used for pathogen detection, and the results obtained using the assays we designed with this method served to validate those obtained using the MiSeq method. These two methods offer distinct advantages. The MiSeq method pairs pathogen detection and microbial community profiling to answer immediate and long-term fish health concerns, while the droplet digital PCR method provides fast and highly sensitive detection that is useful for surveillance and rapid clinical responses.

Deciphering the Molecular Basis for Attenuation of Flavobacterium columnare Strain Fc1723 Used as Modified Live Vaccine against Columnaris Disease

Vaccines are widely employed in aquaculture to prevent bacterial infections, but their use by the U.S. catfish industry is very limited. One of the main diseases affecting catfish aquaculture is columnaris disease, caused by the bacterial pathogen Flavobacterium columnare. In 2011, a modified-live vaccine against columnaris disease was developed by selecting mutants that were resistant to rifampin. The previous study has suggested that this vaccine is stable, safe, and effective, but the mechanisms that resulted in attenuation remained uncharacterized. To understand the molecular basis for attenuation, a comparative genomic analysis was conducted to identify specific point mutations. The PacBio RS long-read sequencing platform was used to obtain draft genomes of the mutant attenuated strain (Fc1723) and the parent virulent strain (FcB27). Sequence-based genome comparison identified 16 single nucleotide polymorphisms (SNP) unique to the mutant. Genes that contained mutations were involved in rifampin resistance, gliding motility, DNA transcription, toxin secretion, and extracellular protease synthesis. The results also found that the vaccine strain formed biofilm at a significantly lower rate than the parent strain. These observations suggested that the rifampin-resistant phenotype and the associated attenuation of the vaccine strain result from the altered activity of RNA polymerase (RpoB) and possible disrupted protein secretion systems.

Flavobacterium erciyesense sp. nov., a putative non-pathogenic fish symbiont

During a study aiming at isolation and genomic characterization of non-pathogenic fish symbionts, strain F-328^T was isolated from skin mucus of healthy-seeming rainbow trout (Oncorhynchus mykiss) and subjected to polyphasic characterization, including a comprehensive genome analysis. Strain F-328^T was Gram stain negative, none-gliding, strictly aerobic, and rod shaped. The 16S rRNA gene sequence of strain F-328^T exhibited the highest level of identity (98.9%) with F. turcicum F-339^T. The major fatty acids were iso-C_15:0, C_15:1 ω6c, and summed feature 3 (C_16:1 ω7c, C_16:1 ω6c). The polar lipid profile consisted mainly of phosphatidylethanolamine and aminolipids, while the predominant menaquinone was MK-6. The DNA G + C content and genome size of the strain were 35.3% and 3.4 Mb, respectively. Strain F-328^T shows 93.8% average nucleotide identity and 53% digital DNA-DNA hybridization identity with the closest type strain F. turcicum F-339^T. In addition, strain F-328^T was found to carry antimicrobial resistance genes, which confer resistance to several antimicrobials, including aminoglycoside, macrolides, and streptogramin. Differential phenotypic properties, together with the phylogenetic inference, demonstrate that strain F-328^T should be classified as a novel species of the genus Flavobacterium for which the name Flavobacterium erciyesense sp. nov. is proposed, with F-328^T as the type strain (= JCM 34201^T = KCTC 82261^T).

Spatial Variation of Cladophora Epiphytes in the Nan River, Thailand

Cladophora is an algal genus known to be ecologically important. It provides habitats for microorganisms known to provide ecological services such as biosynthesis of cobalamin (vitamin B12) and nutrient cycling. Most knowledge of microbiomes was obtained from studies of lacustrine Cladophora species. However, whether lotic freshwater Cladophora microbiomes are as complex as the lentic ones or provide similar ecological services is not known. To illuminate these issues, we used amplicons of 16S rDNA, 18S rDNA, and ITS to investigate the taxonomy and diversity of the microorganisms associated with replicate Cladophora samples from three sites along the Nan River, Thailand. Results showed that the diversity of prokaryotic and eukaryotic members of Cladophora microbiomes collected from different sampling sites was statistically different. Fifty percent of the identifiable taxa were shared across sampling sites: these included organisms belonging to different trophic levels, decomposers, and heterotrophic bacteria. These heterogeneous assemblages of bacteria, by functional inference, have the potential to perform various ecological functions, i.e., cellulose degradation, cobalamin biosynthesis, fermentative hydrogen production, ammonium oxidation, amino acid fermentation, dissimilatory reduction of nitrate to ammonium, nitrite reduction, nitrate reduction, sulfur reduction, polyphosphate accumulation, denitrifying phosphorus-accumulation, and degradation of aromatic compounds. Results suggested that river populations of Cladophora provide ecologically important habitat for microorganisms that are key to nutrient cycling in lotic ecosystems.

Transcriptome architecture and regulation at environmental transitions in flavobacteria: the case of an important fish pathogen

The family Flavobacteriaceae (phylum Bacteroidetes) is a major component of soil, marine and freshwater ecosystems. In this understudied family, Flavobacterium psychrophilum is a freshwater pathogen that infects salmonid fish worldwide, with critical environmental and economic impact. Here, we report an extensive transcriptome analysis that established the genome map of transcription start sites and transcribed regions, predicted alternative sigma factor regulons and regulatory RNAs, and documented gene expression profiles across 32 biological conditions mimicking the pathogen life cycle. The results link genes to environmental conditions and phenotypic traits and provide insights into gene regulation, highlighting similarities with better known bacteria and original characteristics linked to the phylogenetic position and the ecological niche of the bacterium. In particular, osmolarity appears as a signal for transition between free-living and within-host programs and expression patterns of secreted proteins shed light on probable virulence factors. Further investigations showed that a newly discovered sRNA widely conserved in the genus, Rfp18, is required for precise expression of proteases. By pointing proteins and regulatory elements probably involved in host-pathogen interactions, metabolic pathways, and molecular machineries, the results suggest many directions for future research; a website is made available to facilitate their use to fill knowledge gaps on flavobacteria.

Beyond denitrification: The role of microbial diversity in controlling nitrous oxide reduction and soil nitrous oxide emissions

Many biotic and abiotic processes contribute to nitrous oxide (N₂ O) production in the biosphere, but N₂ O consumption in the environment has heretofore been attributed primarily to canonical denitrifying microorganisms. The nosZ genes encoding the N₂ O reductase enzyme, NosZ, responsible for N₂ O reduction to dinitrogen are now known to include two distinct groups: the well-studied Clade I which denitrifiers typically possess, and the novel Clade II possessed by diverse groups of microorganisms, most of which are non-denitrifiers. Clade II N₂ O reducers could play an important, previously unrecognized role in controlling N₂ O emissions for several reasons, including: (1) the consumption of N₂ O produced by processes other than denitrification, (2) hypothesized non-respiratory functions of NosZ as an electron sink or for N₂ O detoxification, (3) possible differing enzyme kinetics of Clade II NosZ compared to Clade I NosZ, and (4) greater nosZ gene abundance for Clade II compared to Clade I in soils of many ecosystems. Despite the potential ecological significance of Clade II NosZ, a census of 800 peer-reviewed original research articles discussing nosZ and published from 2013 to 2019 showed that the percentage of articles evaluating or mentioning Clade II nosZ increased from 5% in 2013 to only 22% in 2019. The census revealed that the slowly spreading awareness of Clade II nosZ may result in part from disciplinary silos, with the percentage of nosZ articles mentioning Clade II nosZ ranging from 0% in Agriculture and Agronomy journals to 32% in Multidisciplinary Sciences journals. In addition, inconsistent nomenclature for Clade I nosZ and Clade II nosZ, with 17 different terminologies used in the literature, may have created confusion about the two distinct groups of N₂ O reducers. We provide recommendations to accelerate advances in understanding the role of the diversity of N₂ O reducers in regulating soil N₂ O emissions.

Aromatic compound-degrading taxa in an anoxic coal seam microbiome from the Surat Basin, Australia

Methane is an important energy resource internationally, and a large proportion of this methane is produced by microbial communities living in coal seams. Despite the value of this resource for human energy security, our understanding of the metabolic roles played by specific taxa during the biodegradation of coal to methane in situ is quite limited. In order to develop a greater understanding of microbial catabolism on coal, a community from a coal seam in the Surat Basin, Australia, was incubated on 10 different aromatic organic compounds: coronene, benzo[a]pyrene, pyrene, phenanthrene, naphthalene, ethylbenzene, phenol, benzoate, vanillate and syringate. Each of these aromatic compounds either occurs in coal or is a possible product of the coal biodegradation process. 16S rRNA sequencing revealed substantial changes to each community in response to each aromatic carbon substrate provided. Abundant taxa from these substrate-specific communities were identified and their probable catabolic roles proposed based on literature searches of related taxa. This study is the first to link specific coal seam taxa to aromatic substrates available in coal seam environments. Two conceptual models of the putative degradation pathways and key taxa responsible are proposed.

Comparative genomics of Flavobacterium columnare unveils novel insights in virulence and antimicrobial resistance mechanisms

This study reports the comparative analyses of four Flavobacterium columnare isolates that have different virulence and antimicrobial resistance patterns. The main research goal was to reveal new insights into possible virulence genes by comparing the genomes of bacterial isolates that could induce tissue damage and mortality versus the genome of a non-virulent isolate. The results indicated that only the genomes of the virulent isolates possessed unique genes encoding amongst others a methyl-accepting chemotaxis protein possibly involved in the initial colonization of tissue, and several VgrG proteins engaged in interbacterial competition. Furthermore, comparisons of genes unique for the genomes of the highly virulent (HV) carp and trout isolates versus the, respectively, low and non-virulent carp and trout isolates were performed. An important part of the identified unique virulence genes of the HV-trout isolate was located in one particular gene region identified as a genomic island. This region contained araC and nodT genes, both linked to pathogenic and multidrug-resistance, and a luxR-gene, functional in bacterial cell-to-cell communication. Furthermore, the genome of the HV-trout isolate possessed unique sugar-transferases possibly important in bacterial adhesion. The second research goal was to obtain insights into the genetic basis of acquired antimicrobial resistance. Several point-mutations were discovered in gyrase-genes of an isolate showing phenotypic resistance towards first and second-generation quinolones, which were absent in isolates susceptible to quinolones. Tetracycline-resistance gene tetA was found in an isolate displaying acquired phenotypic resistance towards oxytetracycline. Although not localized on a prophage, several flanking genes were indicative of the gene’s mobile character.

Prophages and Past Prophage-Host Interactions Revealed by CRISPR Spacer Content in a Fish Pathogen

The role of prophages in the evolution, diversification, or virulence of the fish pathogen Flavobacterium columnare has not been studied thus far. Here, we describe a functional spontaneously inducing prophage fF4 from the F. columnare type strain ATCC 23463, which is not detectable with commonly used prophage search methods. We show that this prophage type has a global distribution and is present in strains isolated from Finland, Thailand, Japan, and North America. The virions of fF4 are myoviruses with contractile tails and infect only bacterial strains originating from Northern Finland. The fF4 resembles transposable phages by similar genome organization and several gene orthologs. Additional bioinformatic analyses reveal several species in the phylum Bacteroidetes that host a similar type of putative prophage, including bacteria that are important animal and human pathogens. Furthermore, a survey of F. columnare Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) spacers indicate a shared evolutionary history between F. columnare strains and the fF4 phage, and another putative prophage in the F. columnare strain ATCC 49512, named p49512. First, CRISPR spacer content from the two CRISPR loci (types II-C and VI-B) of the fF4 lysogen F. columnare ATCC 23463 revealed a phage terminase protein-matching spacer in the VI-B locus. This spacer is also present in two Chinese F. columnare strains. Second, CRISPR analysis revealed four F. columnare strains that contain unique spacers targeting different regions of the putative prophage p49512 in the F. columnare strain ATCC 49512, despite the geographical distance or genomovar of the different strains. This suggests a common ancestry for the F. columnare prophages and different host strains.

The role of denitrification genes in anaerobic growth and virulence of Flavobacterium columnare

AIMS

Comparative genomics analyses indicated that the Flavobacterium columnare genome has unique denitrification genes relative to Flavobacterium psychrophilum and Flavobacterium johnsoniae, including nasA (nitrate reductase), nirS (nitrite reductase), norB (nitric oxide reductase) and nosZ (nitrous oxide reductase). The current study determines the roles of nasA, nirS, norB and nosZ in anaerobic growth, nitrate reduction, biofilm formation and virulence.

METHODS AND RESULTS

Four in-frame deletion mutants in virulent F. columnare strain 94-081 were constructed by allelic exchange using pCP29 plasmid. Compared with parent strain 94-081, FcΔnasA,FcΔnirS and FcΔnosZ mutants did not grow as well anaerobically, whereas the growth of FcΔnorB strain was similar to the parent strain (FcWT). Exogenous nitrate was not significantly consumed under anaerobic conditions in FcΔnasA, FcΔnirS and FcΔnosZ compared to parent strain 94-081. Under anaerobic conditions, Fc∆nasA, Fc∆norB and Fc∆nosZ formed significantly less biofilm than the wild type strain at 24 and 96 h, but FcΔnirS was not significantly affected. The nitrite reductase mutant FcΔnirS was highly attenuated in catfish, whereas FcΔnasA, FcΔnorB and FcΔnosZ had similar virulence to FcWT.

CONCLUSIONS

These results show, for the first time, that denitrification genes enable F. columnare to grow anaerobically using nitrate as an electron acceptor, and nitrite reductase contributes to F. columnare virulence.

SIGNIFICANCE AND IMPACT OF THE STUDY

These findings indicate potential for F. columnare to grow in nitrate-rich anaerobic zones in catfish production ponds, and they suggest that a Fc∆nirS strain could be useful as a safe live vaccine if it protects catfish against columnaris disease.

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.

Symbiosis, Selection, and Novelty: Freshwater Adaptation in the Unique Sponges of Lake Baikal

Freshwater sponges (Spongillida) are a unique lineage of demosponges that secondarily colonized lakes and rivers and are now found ubiquitously in these ecosystems. They developed specific adaptations to freshwater systems, including the ability to survive extreme thermal ranges, long-lasting dessication, anoxia, and resistance to a variety of pollutants. Although spongillids have colonized all freshwater systems, the family Lubomirskiidae is endemic to Lake Baikal and plays a range of key roles in this ecosystem. Our work compares the genomic content and microbiome of individuals of three species of the Lubomirskiidae, providing hypotheses for how molecular evolution has allowed them to adapt to their unique environments. We have sequenced deep (>92% of the metazoan "Benchmarking Universal Single-Copy Orthologs" [BUSCO] set) transcriptomes from three species of Lubomirskiidae and a draft genome resource for Lubomirskia baikalensis. We note Baikal sponges contain unicellular algal and bacterial symbionts, as well as the dinoflagellate Gyrodinium. We investigated molecular evolution, gene duplication, and novelty in freshwater sponges compared with marine lineages. Sixty one orthogroups have consilient evidence of positive selection. Transporters (e.g., zinc transporter-2), transcription factors (aristaless-related homeobox), and structural proteins (e.g. actin-3), alongside other genes, are under strong evolutionary pressure in freshwater, with duplication driving novelty across the Spongillida, but especially in the Lubomirskiidae. This addition to knowledge of freshwater sponge genetics provides a range of tools for understanding the molecular biology and, in the future, the ecology (e.g., colonization and migration patterns) of these key species.

Comparative genomic analyses reveal diverse virulence factors and antimicrobial resistance mechanisms in clinical Elizabethkingia meningoseptica strains

Three human clinical isolates of bacteria (designated strains Em1, Em2 and Em3) had high average nucleotide identity (ANI) to Elizabethkingia meningoseptica. Their genome sizes (3.89, 4.04 and 4.04 Mb) were comparable to those of other Elizabethkingia species and strains, and exhibited open pan-genome characteristics, with two strains being nearly identical and the third divergent. These strains were susceptible only to trimethoprim/sulfamethoxazole and ciprofloxacin amongst 16 antibiotics in minimum inhibitory tests. The resistome exhibited a high diversity of resistance genes, including 5 different lactamase- and 18 efflux protein- encoding genes. Forty-four genes encoding virulence factors were conserved among the strains. Sialic acid transporters and curli synthesis genes were well conserved in E. meningoseptica but absent in E. anophelis and E. miricola. E. meningoseptica carried several genes contributing to biofilm formation. 58 glycoside hydrolases (GH) and 25 putative polysaccharide utilization loci (PULs) were found. The strains carried numerous genes encoding two-component system proteins (56), transcription factor proteins (187~191), and DNA-binding proteins (6~7). Several prophages and CRISPR/Cas elements were uniquely present in the genomes.

Gliding Motility and Expression of Motility-Related Genes in Spreading and Non-spreading Colonies of Flavobacterium columnare

Gliding motility facilitates the movement of bacteria along surfaces in many Bacteroidetes species and results in spreading colonies. The adhesins required for the gliding are secreted through a gliding motility-associated protein secretion system, known as the type IX secretion system (T9SS). The fish pathogen Flavobacterium columnare produces spreading (rhizoid [Rz], soft [S]) and non-spreading (rough [R]) colony types, of which only the spreading Rz type is virulent. In this study, we explored the spreading behavior of these colony types by microscopic imaging and measured the expression of genes associated with gliding motility and T9SS (gldG, gldH, gldL, sprA, sprB, sprE, sprF, sprT, and porV) under high and low resource levels by using RT-qPCR (reverse transcription quantitative PCR). The spreading colony types responded to the low resource level with increased colony size. The non-spreading colony type, as well as the cells growing under high nutrient level expressed only moderate cell movements. Yet, a low nutrient level provoked more active gliding motility in individual cells and increased spreading by cooperative gliding. The gene expression survey demonstrated an increased expression level of sprA (a core component of T9SS) and sprF (needed for adhesin secretion) under low nutrient conditions. Surprisingly, the expression of gliding motility genes was not consistently associated with more active spreading behavior. Furthermore, no genetic differences were found between spreading and non-spreading colony types in the studied genes associated with gliding motility. Our study demonstrates that environmental nutrient level is an important regulator of both gliding motility and the expression of some of the associated genes. These results may help to understand the connections between nutrient concentration, gliding motility, and virulence of F. columnare.

The Type IX Secretion System Is Required for Virulence of the Fish Pathogen Flavobacterium columnare

Flavobacterium columnare, a member of the phylum Bacteroidetes, causes columnaris disease in wild and aquaculture-reared freshwater fish. The mechanisms responsible for columnaris disease are not known. Many members of the phylum Bacteroidetes use type IX secretion systems (T9SSs) to secrete enzymes, adhesins, and proteins involved in gliding motility. The F. columnare genome has all of the genes needed to encode a T9SS. gldN, which encodes a core component of the T9SS, was deleted in wild-type strains of F. columnare The F. columnare ΔgldN mutants were deficient in the secretion of several extracellular proteins and lacked gliding motility. The ΔgldN mutants exhibited reduced virulence in zebrafish, channel catfish, and rainbow trout, and complementation restored virulence. PorV is required for the secretion of a subset of proteins targeted to the T9SS. An F. columnare ΔporV mutant retained gliding motility but exhibited reduced virulence. Cell-free spent media from exponentially growing cultures of wild-type and complemented strains caused rapid mortality, but spent media from ΔgldN and ΔporV mutants did not, suggesting that soluble toxins are secreted by the T9SS.IMPORTANCE Columnaris disease, caused by F. columnare, is a major problem for freshwater aquaculture. Little is known regarding the virulence factors produced by F. columnare, and control measures are limited. Analysis of targeted gene deletion mutants revealed the importance of the type IX protein secretion system (T9SS) and of secreted toxins in F. columnare virulence. T9SSs are common in members of the phylum Bacteroidetes and likely contribute to the virulence of other animal and human pathogens.

Comparative genome analysis of fish pathogen Flavobacterium columnare reveals extensive sequence diversity within the species

Flavobacterium columnare is one of the deadliest fish pathogens causing devastating mortality in various freshwater fish species globally. To gain an insight into bacterial genomic contents and structures, comparative genome analyses were performed using the reference and newly sequenced genomes of F. columnare including genomovar I, II and I/II strains isolated from Thailand, Europe and the USA. Bacterial genomes varied in size from 3.09 to 3.39Mb (2714 to 3101 CDSs). The pan-genome analysis revealed open pan-genome nature of F. columnare strains, which possessed at least 4953 genes and tended to increase progressively with the addition of a new genome. Genomic islands (GIs) present in bacterial genomes were diverse, in which 65% (39 out of 60) of possible GIs were strain-specific. A CRISPR/cas investigation indicated at least two different CRISPR systems with varied spacer profiles. On the other hand, putative virulence genes, including those related to gliding motility, type IX secretion system (T9SS), outer membrane proteins (Omp), were equally distributed among F. columnare strains. The MLSA scheme categorized bacterial strains into nine different sequence types (ST 9-17). Phylogenetic analyses based on either 16S rRNA, MLSA and concatenated SNPs of core genome revealed the diversity of F. columnare strains. DNA homology analysis indicated that the estimated digital DNA-DNA hybridization (dDDH) between strains of genomovar I and II can be as low as 42.6%, while the three uniquely tilapia-originated strains from Thailand (1214, NK01 and 1215) were clearly dissimilar to other F. columnare strains as the dDDH values were only 27.7-30.4%. Collectively, this extensive diversity among bacterial strains suggested that species designation of F. columnare would potentially require re-emendation.