Pangenome Evolution in the Marine Bacterium Alteromonas

Dispersant addition, but not nutrients, stimulated blooms of multiple hydrocarbonoclastic genera in nutrient-replete coastal marine surface waters

The range of impacts of chemical dispersants on indigenous marine microbial communities and their activity remains poorly constrained. We tested the response of nearshore surface waters chronically exposed to oil leakage from a downed platform and supplied with nutrients by the Mississippi River to Corexit dispersant and nutrient additions. As assessed using 14C-labeled tracers, hexadecane mineralization potential was orders of magnitude higher in all unamended samples than in previously assessed bathypelagic communities. Nutrient additions stimulated microbial mortality but did not affect community composition and had no generalizable effect on hydrocarbon mineralization potential. By contrast, Corexit amendments caused a rapid shift in community composition and a drawdown of inorganic nitrogen and orthophosphate though no generalizable effect on hydrocarbon mineralization potential. The hydrocarbonoclastic community's response to dispersants is largely driven by the relative availability of organic substrates and nutrients, underscoring the role of environmental conditions and multiple interacting stressors on hydrocarbon degradation potential.

Flexible genomic island conservation across freshwater and marine Methylophilaceae

The evolutionary trajectory of Methylophilaceae includes habitat transitions from freshwater sediments to freshwater and marine pelagial that resulted in genome reduction (genome-streamlining) of the pelagic taxa. However, the extent of genetic similarities in the genomic structure and microdiversity of the two genome-streamlined pelagic lineages (freshwater "Ca. Methylopumilus" and the marine OM43 lineage) has so far never been compared. Here, we analyzed complete genomes of 91 "Ca. Methylopumilus" strains isolated from 14 lakes in Central Europe and 12 coastal marine OM43 strains. The two lineages showed a remarkable niche differentiation with clear species-specific differences in habitat preference and seasonal distribution. On the other hand, we observed a synteny preservation in their genomes by having similar locations and types of flexible genomic islands (fGIs). Three main fGIs were identified: a replacement fGI acting as phage defense, an additive fGI harboring metabolic and resistance-related functions, and a tycheposon containing nitrogen-, thiamine-, and heme-related functions. The fGIs differed in relative abundances in metagenomic datasets suggesting different levels of variability ranging from strain-specific to population-level adaptations. Moreover, variations in one gene seemed to be responsible for different growth at low substrate concentrations and a potential biogeographic separation within one species. Our study provides a first insight into genomic microdiversity of closely related taxa within the family Methylophilaceae and revealed remarkably similar dynamics involving mobile genetic elements and recombination between freshwater and marine family members.

Phosphate-related genomic islands as drivers of environmental adaptation in the streamlined marine alphaproteobacterial HIMB59

IMPORTANCE

These results shed light on the evolutionary strategies of microbes with streamlined genomes to adapt and survive in the oligotrophic conditions that dominate the surface waters of the global ocean. At the individual level, these microbes have been subjected to evolutionary constraints that have led to a more efficient use of nutrients, removing non-essential genes named as "streamlining theory." However, at the population level, they conserve a highly diverse gene pool in flexible genomic islands resulting in polyclonal populations on the same genomic background as an evolutionary response to environmental pressures. Localization of these islands at equivalent positions in the genome facilitates horizontal transfer between clonal lineages. This high level of environmental genomic heterogeneity could explain their cosmopolitan distribution. In the case of the order HIMB59 within the class Alphaproteobacteria, two factors exert evolutionary pressure and determine this intraspecific diversity: phages and the concentration of P in the environment.

Distantly related Alteromonas bacteriophages share tail fibers exhibiting properties of transient chaperone caps

The host recognition modules encoding the injection machinery and receptor binding proteins (RBPs) of bacteriophages are predisposed to mutation and recombination to maintain infectivity towards co-evolving bacterial hosts. In this study, we reveal how Alteromonas mediterranea schitovirus A5 shares its host recognition module, including tail fiber and cognate chaperone, with phages from distantly related families including Alteromonas myovirus V22. While the V22 chaperone is essential for producing active tail fibers, here we demonstrate production of functional A5 tail fibers regardless of chaperone co-expression. AlphaFold-generated models of tail fiber and chaperone pairs from phages A5, V22, and other Alteromonas phages reveal how amino acid insertions within both A5-like proteins results in a knob domain duplication in the tail fiber and a chaperone β-hairpin "tentacle" extension. These structural modifications are linked to differences in chaperone dependency between the A5 and V22 tail fibers. Structural similarity between the chaperones and intramolecular chaperone domains of other phage RBPs suggests an additional function of these chaperones as transient fiber "caps". Finally, our identification of homologous host recognition modules from morphologically distinct phages implies that horizontal gene transfer and recombination events between unrelated phages may be a more common process than previously thought among Caudoviricetes phages.

Exploring environmental intra-species diversity through non-redundant pangenome assemblies

At the genome level, microorganisms are highly adaptable both in terms of allele and gene composition. Such heritable traits emerge in response to different environmental niches and can have a profound influence on microbial community dynamics. As a consequence, any individual genome or population will contain merely a fraction of the total genetic diversity of any operationally defined "species", whose ecological potential can thus be only fully understood by studying all of their genomes and the genes therein. This concept, known as the pangenome, is valuable for studying microbial ecology and evolution, as it partitions genomes into core (present in all the genomes from a species, and responsible for housekeeping and species-level niche adaptation among others) and accessory regions (present only in some, and responsible for intra-species differentiation). Here we present SuperPang, an algorithm producing pangenome assemblies from a set of input genomes of varying quality, including metagenome-assembled genomes (MAGs). SuperPang runs in linear time and its results are complete, non-redundant, preserve gene ordering and contain both coding and non-coding regions. Our approach provides a modular view of the pangenome, identifying operons and genomic islands, and allowing to track their prevalence in different populations. We illustrate this by analysing intra-species diversity in Polynucleobacter, a bacterial genus ubiquitous in freshwater ecosystems, characterized by their streamlined genomes and their ecological versatility. We show how SuperPang facilitates the simultaneous analysis of allelic and gene content variation under different environmental pressures, allowing us to study the drivers of microbial diversification at unprecedented resolution.

Host range of strand-biased circularizing integrative elements: a new class of mobile DNA elements nesting in Gammaproteobacteria

BACKGROUND

The strand-biased circularizing integrative elements (SEs) are putatively non-mobilizable integrative elements for transmitting antimicrobial resistance genes. The transposition mode and the prevalence of SEs in prokaryotes remain vague.

RESULTS

To corroborate the transposition mode and the prevalence of SEs, hypothetical transposition intermediates of an SE were searched for in genomic DNA fractions of an SE host. Then, the SE core genes were defined based on gene knockout experiments, and the synteny blocks of their distant homologs were searched for in the RefSeq complete genome sequence database using PSI-BLAST. A genomic DNA fractionation experiment revealed that SE copies are present in a double-stranded nicked circular form in vivo. Operonic structure of three conserved coding sequences (intA, tfp, intB) and srap located at the left end of SEs were identified as essential for attL × attR recombination. The synteny blocks of tfp and srap homologs were detected in 3.6% of the replicons of Gammaproteobacteria but not in other taxa, implying that SE movement is host-dependent. SEs have been discovered most frequently in the orders Vibrionales (19% of replicons), Pseudomonadales (18%), Alteromonadales (17%), and Aeromonadales (12%). Genomic comparisons revealed 35 new SE members with identifiable termini. SEs are present at 1 to 2 copies per replicon and have a median length of 15.7 kb. Three newly identified SE members carry antimicrobial resistance genes, like tmexCD-toprJ, mcr-9, and bla_GMA-1. Further experiments validated that three new SE members possess the strand-biased attL × attR recombination activity.

CONCLUSIONS

This study suggested that transposition intermediates of SEs are double-stranded circular DNA. The main hosts of SEs are a subset of free-living Gammaproteobacteria; this represents a rather narrow host range compared to those of mobile DNA element groups discovered to date. As the host range, genetic organization, and movements are unique among the mobile DNA elements, SEs provide a new model system for host-mobile DNA element coevolution studies.

Single-amplified genomes reveal most streamlined free-living marine bacteria

Evolutionary adaptations of prokaryotes to the environment sometimes result in genome reduction. Our knowledge of this phenomenon among free-living bacteria remains scarce. We address the dynamics and limits of genome reduction by examining one of the most abundant bacteria in the ocean, the SAR86 clade. Despite its abundance, comparative genomics has been limited by the absence of pure cultures and the poor representation in metagenome-assembled genomes. We co-assembled multiple previously available single-amplified genomes to obtain the first complete genomes from members of the four families. All families showed a convergent evolutionary trajectory with characteristic features of streamlined genomes, most pronounced in the TMED112 family. This family has a genome size of ca. 1 Mb and only 1 bp as median intergenic distance, exceeding values found in other abundant microbes such as SAR11, OM43 and Prochlorococcus. This genomic simplification led to a reduction in the biosynthesis of essential molecules, DNA repair-related genes, and the ability to sense and respond to environmental factors, which could suggest an evolutionary dependence on other co-occurring microbes for survival (Black Queen hypothesis). Therefore, these reconstructed genomes within the SAR86 clade provide new insights into the limits of genome reduction in free-living marine bacteria.

Phototroph-heterotroph interactions during growth and long-term starvation across Prochlorococcus and Alteromonas diversity

Due to their potential impact on ecosystems and biogeochemistry, microbial interactions, such as those between phytoplankton and bacteria, have been studied intensively using specific model organisms. Yet, to what extent interactions differ between closely related organisms, or how these interactions change over time, or culture conditions, remains unclear. Here, we characterize the interactions between five strains each of two globally abundant marine microorganisms, Prochlorococcus (phototroph) and Alteromonas (heterotroph), from the first encounter between individual strains and over more than a year of repeated cycles of exponential growth and long-term nitrogen starvation. Prochlorococcus-Alteromonas interactions had little effect on traditional growth parameters such as Prochlorococcus growth rate, maximal fluorescence, or lag phase, affecting primarily the dynamics of culture decline, which we interpret as representing cell mortality and lysis. The shape of the Prochlorococcus decline curve and the carrying capacity of the co-cultures were determined by the phototroph and not the heterotroph strains involved. Comparing various mathematical models of culture mortality suggests that Prochlorococcus death rate increases over time in mono-cultures but decreases in co-cultures, with cells potentially becoming more resistant to stress. Our results demonstrate intra-species differences in ecologically relevant co-culture outcomes. These include the recycling efficiency of N and whether the interactions are mutually synergistic or competitive. They also highlight the information-rich growth and death curves as a useful readout of the interaction phenotype.

Microbe Profile: Alteromonas macleodii - a widespread, fast-responding, 'interactive' marine bacterium

Alteromonas macleodii is a marine heterotrophic bacterium with widespread distribution - from temperate to tropical oceans, and from surface to deep waters. Strains of A. macleodii exhibit considerable genomic and metabolic variability, and can grow rapidly on diverse organic compounds. A. macleodii is a model organism for the study of population genomics, physiological adaptations and microbial interactions, with individual genomes encoding diverse phenotypic traits influenced by recombination and horizontal gene transfer.

mOTUpan: a robust Bayesian approach to leverage metagenome-assembled genomes for core-genome estimation

Recent advances in sequencing and bioinformatics have expanded the tree of life by providing genomes for uncultured environmentally relevant clades, either through metagenome-assembled genomes or through single-cell genomes. While this expanded diversity can provide novel insights into microbial population structure, most tools available for core-genome estimation are sensitive to genome completeness. Consequently, a major portion of the huge phylogenetic diversity uncovered by environmental genomic approaches remains excluded from such analyses. We present mOTUpan, a novel iterative Bayesian method for computing the core genome for sets of genomes of highly diverse completeness range. The likelihood for each gene cluster to belong to core or accessory genome is estimated by computing the probability of its presence/absence pattern in the target genome set. The core-genome prediction is computationally efficient and can be scaled up to thousands of genomes. It has shown comparable estimates to state-of-the-art tools Roary and PPanGGOLiN for high-quality genomes and is capable of using genomes at lower completeness thresholds. mOTUpan wraps a bootstrapping procedure to estimate the quality of a specific core-genome prediction, as the accuracy of each run will depend on the specific completeness distribution and the number of genomes in the dataset under scrutiny. mOTUpan is implemented in the mOTUlizer software package, and available at github.com/moritzbuck/mOTUlizer, under GPL 3.0 license.

Shallow-Water Hydrothermal Vents as Natural Accelerators of Bacterial Antibiotic Resistance in Marine Coastal Areas

Environmental contamination by heavy metals (HMs) poses several indirect risks to human health, including the co-spreading of genetic traits conferring resistance to both HMs and antibiotics among micro-organisms. Microbial antibiotic resistance (AR) acquisition is enhanced at sites anthropogenically polluted by HMs, as well as in remote systems naturally enriched in HMs, such as hydrothermal vents in the deep sea. However, to date, the possible role of hydrothermal vents at shallower water depths as hot spots of microbial AR gain and spreading has not been tested, despite the higher potential risks associated with the closer vicinity of such ecosystems to coasts and human activities. In this work, we collected waters and sediments at the Panarea shallow-water hydrothermal vents, testing the presence of culturable marine bacteria and their sensitivity to antibiotics and HMs. All of the bacterial isolates showed resistance to at least one antibiotic and one HM and, most notably, 80% of them displayed multi-AR on average to 12 (min 8, max 15) different antibiotics, as well as multi-HM tolerance. We show that our isolates displayed high similarity (≥99%) to common marine bacteria, affiliating with Actinobacteria, Gammaproteobacteria, Alphaproteobacteria and Firmicutes, and all displayed wide growth ranges for temperature and salinity during in vitro physiological tests. Notably, the analysis of the genomes available in public databases for their closest relatives highlighted the lack of genes for AR, posing new questions on the origin of multi-AR acquisition in this peculiar HM-rich environment. Overall, our results point out that shallow-water hydrothermal vents may contribute to enhance AR acquisition and spreading among common marine bacteria in coastal areas, highlighting this as a focus for future research.

Enhanced copper-resistance gene repertoire in Alteromonas macleodii strains isolated from copper-treated marine coatings

Copper is prevalent in coastal ecosystems due to its use as an algaecide and as an anti-fouling agent on ship hulls. Alteromonas spp. have previously been shown to be some of the early colonizers of copper-based anti-fouling paint but little is known about the mechanisms they use to overcome this initial copper challenge. The main models of copper resistance include the Escherichia coli chromosome-based Cue and Cus systems; the plasmid-based E. coli Pco system; and the plasmid-based Pseudomonas syringae Cop system. These were all elucidated from strains isolated from copper-rich environments of agricultural and/or enteric origin. In this work, copper resistance assays demonstrated the ability of Alteromonas macleodii strains CUKW and KCC02 to grow at levels lethal to other marine bacterial species. A custom database of Hidden Markov Models was designed based on proteins from the Cue, Cus, and Cop/Pco systems and used to identify potential copper resistance genes in CUKW and KCC02. Comparative genomic analyses with marine bacterial species and bacterial species isolated from copper-rich environments demonstrated that CUKW and KCC02 possess genetic elements of all systems, oftentimes with multiple copies, distributed throughout the chromosome and mega-plasmids. In particular, two copies of copA (the key player in cytoplasmic detoxification), each with its own apparent MerR-like transcriptional regulator, occur on a mega-plasmid, along with multiple copies of Pco homologs. Genes from both systems were induced upon exposure to elevated copper levels (100 μM- 3 mM). Genomic analysis identified one of the merR-copA clusters occurs on a genomic island (GI) within the plasmid, and comparative genomic analysis found that either of the merR-copA clusters, which also includes genes coding for a cupredoxin domain-containing protein and an isoprenylcysteine methyltransferase, occurs on a GI across diverse bacterial species. These genomic findings combined with the ability of CUKW and KCC02 to grow in copper-challenged conditions are couched within the context of the genome flexibility of the Alteromonas genus.

Enhanced Recovery of Microbial Genes and Genomes From a Marine Water Column Using Long-Read Metagenomics

Third-generation sequencing has penetrated little in metagenomics due to the high error rate and dependence for assembly on short-read designed bioinformatics. However, second-generation sequencing metagenomics (mostly Illumina) suffers from limitations, particularly in the assembly of microbes with high microdiversity and retrieval of the flexible (adaptive) fraction of prokaryotic genomes. Here, we have used a third-generation technique to study the metagenome of a well-known marine sample from the mixed epipelagic water column of the winter Mediterranean. We have compared PacBio Sequel II with the classical approach using Illumina Nextseq short reads followed by assembly to study the metagenome. Long reads allow for efficient direct retrieval of complete genes avoiding the bias of the assembly step. Besides, the application of long reads on metagenomic assembly allows for the reconstruction of much more complete metagenome-assembled genomes (MAGs), particularly from microbes with high microdiversity such as Pelagibacterales. The flexible genome of reconstructed MAGs was much more complete containing many adaptive genes (some with biotechnological potential). PacBio Sequel II CCS appears particularly suitable for cellular metagenomics due to its low error rate. For most applications of metagenomics, from community structure analysis to ecosystem functioning, long reads should be applied whenever possible. Specifically, for in silico screening of biotechnologically useful genes, or population genomics, long-read metagenomics appears presently as a very fruitful approach and can be analyzed from raw reads before a computationally demanding (and potentially artifactual) assembly step.

A Novel Broad Host Range Phage Infecting Alteromonas

Bacteriophages substantially contribute to bacterial mortality in the ocean and play critical roles in global biogeochemical processes. Alteromonas is a ubiquitous bacterial genus in global tropical and temperate waters, which can cross-protect marine cyanobacteria and thus has important ecological benefits. However, little is known about the biological and ecological features of Alteromonas phages (alterophages). Here, we describe a novel alterophage vB_AmeP-R8W (R8W), which belongs to the Autographiviridae family and infects the deep-clade Alteromonas mediterranea. R8W has an equidistant and icosahedral head (65 ± 1 nm in diameter) and a short tail (12 ± 2 nm in length). The genome size of R8W is 48,825 bp, with a G + C content of 40.55%. R8W possesses three putative auxiliary metabolic genes encoding proteins involved in nucleotide metabolism and DNA binding: thymidylate synthase, nucleoside triphosphate pyrophosphohydrolase, and PhoB. R8W has a rapid lytic cycle with a burst size of 88 plaque-forming units/cell. Notably, R8W has a wide host range, such that it can infect 35 Alteromonas strains; it exhibits a strong specificity for strains isolated from deep waters. R8W has two specific receptor binding proteins and a compatible holin-endolysin system, which contribute to its wide host range. The isolation of R8W will contribute to the understanding of alterophage evolution, as well as the phage-host interactions and ecological importance of alterophages.

Genomes of the "Candidatus Actinomarinales" Order: Highly Streamlined Marine Epipelagic Actinobacteria

Microbiology is in a new age in which sequence databases are primary sources of information about many microbes. However, in-depth analysis of environmental genomes thus retrieved is essential to substantiate the new knowledge.

“Candidatus Actinomarinales” was defined as a subclass of exclusively marine Actinobacteria with small cells and genomes. We have collected all the available genomes in databases to assess the diversity included in this group and analyzed it by comparative genomics. We have found the equivalent of five genera and 18 genomospecies. They have genome reduction parameters equal to those of freshwater actinobacterial “Candidatus Nanopelagicales” or marine alphaproteobacterial Pelagibacterales. Genome recruitment shows that they are found only in the photic zone and mainly in surface waters, with only one genus that is found preferentially at or below the deep chlorophyll maximum. “Ca. Actinomarinales” show a highly conserved core genome (80% of the gene families conserved for the whole order) with a saturation of genomic diversity of the flexible genome at the genomospecies level. We found only a flexible genomic island preserved throughout the order; it is related to the sugar decoration of the envelope and uses several tRNAs as hot spots to increase its genomic diversity. Populations had a discrete level of sequence diversity similar to other marine microbes but drastically different from the much higher levels found for Pelagibacterales. Genomic analysis suggests that they are all aerobic photoheterotrophs with one type 1 rhodopsin and a heliorhodopsin. Like other actinobacteria, they possess the F420 coenzyme biosynthesis pathway, and its lower reduction potential could provide access to an increased range of redox chemical transformations. Last, sequence analysis revealed the first “Ca. Actinomarinales” phages, including a prophage, with metaviromic islands related to sialic acid cleavage.

IMPORTANCE Microbiology is in a new age in which sequence databases are primary sources of information about many microbes. However, in-depth analysis of environmental genomes thus retrieved is essential to substantiate the new knowledge. Here, we study 182 genomes belonging to the only known exclusively marine pelagic group of the phylum Actinobacteria. The aquatic branch of this phylum is largely known from environmental sequencing studies (single-amplified genomes [SAGs] and metagenome-assembled genomes [MAGs]), and we have collected and analyzed the available information present in databases about the “Ca. Actinomarinales.” They are among the most streamlined microbes to live in the epipelagic zone of the ocean, and their study is critical to obtain a proper view of the diversity of Actinobacteria and their role in aquatic ecosystems.

The Evolutionary Success of the Marine Bacterium SAR11 Analyzed through a Metagenomic Perspective

As the most abundant bacteria in oceans, the Pelagibacterales order (here SAR11) plays an important role in the global carbon cycle, but the study of the evolutionary forces driving its evolution has lagged considerably due to the inherent difficulty of obtaining pure cultures. Multiple evolutionary models have been proposed to explain the diversification of distinct lineages within a population; however, the identification of many of these patterns in natural populations remains mostly enigmatic. We have used a metagenomic approach to explore microdiversity patterns in their natural habitats. Comparison with a collection of bacterial and archaeal groups from the same environments shows that SAR11 populations have a different evolutionary regime, where multiple genotypes coexist within the same population and remain stable over time. Widespread homologous recombination could be one of the main driving factors of this homogenization.

The SAR11 clade of Alphaproteobacteria is the most abundant group of planktonic cells in the near-surface epipelagic waters of the ocean, but the mechanisms underlying its exceptional success have not been fully elucidated. Here, we applied a metagenomic approach to explore microdiversity patterns by measuring the accumulation of synonymous and nonsynonymous mutations as well as homologous recombination in populations of SAR11 from different aquatic habitats (marine epipelagic, bathypelagic, and surface freshwater). The patterns of mutation accumulation and recombination were compared to those of other groups of representative marine microbes with multiple ecological strategies that share the same marine habitat, namely, Cyanobacteria (Prochlorococcus and Synechococcus), Archaea (“Candidatus Nitrosopelagicus” and Marine Group II Thalassoarchaea), and some heterotrophic marine bacteria (Alteromonas and Erythrobacter). SAR11 populations showed widespread recombination among distantly related members, preventing divergence leading to a genetically stable population. Moreover, their high intrapopulation sequence diversity with an enrichment in synonymous replacements supports the idea of a very ancient divergence and the coexistence of multiple different clones. However, other microbes analyzed seem to follow different evolutionary dynamics where processes of diversification driven by geographic and ecological instability produce a higher number of nonsynonymous replacements and lower intrapopulation sequence diversity. Together, these data shed light on some of the evolutionary and ecological processes that lead to the large genomic diversity in SAR11. Furthermore, this approach can be applied to other similar microbes that are difficult to culture in the laboratory, but abundant in nature, to investigate the underlying dynamics of their genomic evolution.

IMPORTANCE As the most abundant bacteria in oceans, the Pelagibacterales order (here SAR11) plays an important role in the global carbon cycle, but the study of the evolutionary forces driving its evolution has lagged considerably due to the inherent difficulty of obtaining pure cultures. Multiple evolutionary models have been proposed to explain the diversification of distinct lineages within a population; however, the identification of many of these patterns in natural populations remains mostly enigmatic. We have used a metagenomic approach to explore microdiversity patterns in their natural habitats. Comparison with a collection of bacterial and archaeal groups from the same environments shows that SAR11 populations have a different evolutionary regime, where multiple genotypes coexist within the same population and remain stable over time. Widespread homologous recombination could be one of the main driving factors of this homogenization.

Flexible genes establish widespread bacteriophage pan-genomes in cryoconite hole ecosystems

Bacteriophage genomes rapidly evolve via mutation and horizontal gene transfer to counter evolving bacterial host defenses; such arms race dynamics should lead to divergence between phages from similar, geographically isolated ecosystems. However, near-identical phage genomes can reoccur over large geographical distances and several years apart, conversely suggesting many are stably maintained. Here, we show that phages with near-identical core genomes in distant, discrete aquatic ecosystems maintain diversity by possession of numerous flexible gene modules, where homologous genes present in the pan-genome interchange to create new phage variants. By repeatedly reconstructing the core and flexible regions of phage genomes from different metagenomes, we show a pool of homologous gene variants co-exist for each module in each location, however, the dominant variant shuffles independently in each module. These results suggest that in a natural community, recombination is the largest contributor to phage diversity, allowing a variety of host recognition receptors and genes to counter bacterial defenses to co-exist for each phage.

Bacteriophages and their hosts are involved in a constant evolutionary arms race that should lead to divergence between phage genes over time. Here, the authors recruit metagenomic reads to virus reference genomes and genome fragments in samples from cryoconite holes and show that phages with near-identical core genomes maintain diversity by possession of numerous flexible gene modules, where homologous genes present in the pan-genome interchange to create new phage variants.

Chitosan Films Incorporated with Exopolysaccharides from Deep Seawater Alteromonas Sp

Two Alteromonas sp. strains isolated from deep seawater were grown to promote the production of exopolysaccharides (EPS, E611 and E805), which were incorporated into chitosan solutions to develop films. The combination of the major marine polysaccharides (chitosan and the isolated bacterial EPS) resulted in the formation of homogenous, transparent, colorless films, suggesting good compatibility between the two components of the film-forming formulation. With regards to optical properties, the films showed low values of gloss, in the range of 5-10 GU, indicating the formation of non-glossy and rough surfaces. In addition to the film surface, both showed hydrophobic character, with water contact angles higher than 100 º, regardless of EPS addition. Among the two EPS under analysis, chitosan films with E805 showed better mechanical performance, leading to resistant, flexible, easy to handle films.

Genomic Characteristics and Potential Metabolic Adaptations of Hadal Trench Roseobacter and Alteromonas Bacteria Based on Single-Cell Genomics Analyses

Heterotrophic bacteria such as those from the Roseobacter group and genus Alteromonas dominate the hadal zones of oceans; however, we know little about the genomic characteristics and potential metabolic adaptations of hadal trench-dwelling bacteria. Here, we report multiple single amplified genomes (SAGs) belonging to Roseobacter and Alteromonas, recovered from the hadal zone of the Mariana Trench. While phylogenetic analyses show that these hadal SAGs cluster with their surface relatives, an analysis of genomic recruitment indicates that they have higher relative abundances in the hadal zone of the Mariana Trench. Comparative genomic analyses between the hadal SAGs and reference genomes of closely related shallow-water relatives indicate that genes involved in the mobilome (prophages and transposons) are overrepresented among the unique genes of the hadal Roseobacter and Alteromonas SAGs; the functional proteins encoded by this category of genes also shows higher amino acid sequence variation than those encoded by other gene sets within the Roseobacter SAGs. We also found that genes involved in cell wall/membrane/envelope biogenesis, transcriptional regulation, and metal transport may be important for the adaptation of hadal Roseobacter and Alteromonas lineages. These results imply that the modification of cell surface-related proteins and transporters is the major direction of genomic evolution in Roseobacter and Alteromonas bacteria adapting to the hadal environment, and that prophages and transposons may be the key factors driving this process.

Differential Responses of a Coastal Prokaryotic Community to Phytoplanktonic Organic Matter Derived from Cellular Components and Exudates

The phytoplanktonic production and prokaryotic consumption of organic matter significantly contribute to marine carbon cycling. Organic matter released from phytoplankton via three processes (exudation of living cells, cell disruption through grazing, and viral lysis) shows distinct chemical properties. We herein investigated the effects of phytoplanktonic whole-cell fractions (WF) (representing cell disruption by grazing) and extracellular fractions (EF) (representing exudates) prepared from Heterosigma akashiwo, a bloom-forming Raphidophyceae, on prokaryotic communities using culture-based experiments. We analyzed prokaryotic community changes for two weeks. The shift in cell abundance by both treatments showed similar dynamics, reaching the first peak (~4.1×10⁶‍ ‍cells‍ ‍mL^–1) on day 3 and second peak (~1.1×10⁶‍ ‍cells‍ ‍mL^–1) on day 13. We classified the sequences obtained into operational taxonomic units (OTUs). A Bray-Curtis dissimilarity analysis revealed that the OTU-level community structure changed distinctively with the two treatments. Ten and 13 OTUs were specifically abundant in the WF and EF treatments, respectively. These OTUs were assigned as heterotrophic bacteria mainly belonging to the Alteromonadales (Gammaproteobacteria) and Bacteroidetes clades and showed successive dynamics following the addition of organic matter. We also analyzed the dynamics of these OTUs in the ocean using publicly available metagenomic data from a natural coastal bloom in Monterey Bay, USA. At least two WF treatment OTUs showed co-occurrence with H. akashiwo, indicating that the blooms of H. akashiwo also affect these OTUs in the ocean. The present results strongly suggest that the thriving and dead cells of uninfected phytoplankton differentially influence the marine prokaryotic community.

Alteromonas Myovirus V22 Represents a New Genus of Marine Bacteriophages Requiring a Tail Fiber Chaperone for Host Recognition

Marine phages play a variety of critical roles in regulating the microbial composition of our oceans. Despite constituting the majority of genetic diversity within these environments, there are relatively few isolates with complete genome sequences or in-depth analyses of their host interaction mechanisms, such as characterization of their receptor binding proteins (RBPs). Here, we present the 92,760-bp genome of the Alteromonas-targeting phage V22. Genomic and morphological analyses identify V22 as a myovirus; however, due to a lack of sequence similarity to any other known myoviruses, we propose that V22 be classified as the type phage of a new Myoalterovirus genus within the Myoviridae family. V22 shows gene homology and synteny with two different subfamilies of phages infecting enterobacteria, specifically within the structural region of its genome. To improve our understanding of the V22 adsorption process, we identified putative RBPs (gp23, gp24, and gp26) and tested their ability to decorate the V22 propagation strain, Alteromonas mediterranea PT11, as recombinant green fluorescent protein (GFP)-tagged constructs. Only GFP-gp26 was capable of bacterial recognition and identified as the V22 RBP. Interestingly, production of functional GFP-gp26 required coexpression with the downstream protein gp27. GFP-gp26 could be expressed alone but was incapable of host recognition. By combining size-exclusion chromatography with fluorescence microscopy, we reveal how gp27 is not a component of the final RBP complex but instead is identified as a new type of phage-encoded intermolecular chaperone that is essential for maturation of the gp26 RBP.IMPORTANCE Host recognition by phage-encoded receptor binding proteins (RBPs) constitutes the first step in all phage infections and the most critical determinant of host specificity. By characterizing new types of RBPs and identifying their essential chaperones, we hope to expand the repertoire of known phage-host recognition machineries. Due to their genetic plasticity, studying RBPs and their associated chaperones can shed new light onto viral evolution affecting phage-host interactions, which is essential for fields such as phage therapy or biotechnology. In addition, since marine phages constitute one of the most important reservoirs of noncharacterized genetic diversity on the planet, their genomic and functional characterization may be of paramount importance for the discovery of novel genes with potential applications.

Transcriptomic Study of Substrate-Specific Transport Mechanisms for Iron and Carbon in the Marine Copiotroph Alteromonas macleodii

As the major facilitators of the turnover of organic matter in the marine environment, the ability of heterotrophic bacteria to acquire specific compounds within the diverse range of dissolved organic matter will affect the regeneration of essential nutrients such as iron and carbon. TonB-dependent transporters are a prevalent cellular tool in Gram-negative bacteria that allow a relatively high-molecular-weight fraction of organic matter to be directly accessed. However, these transporters are not well characterized in marine bacteria, limiting our understanding of the flow of specific substrates through the marine microbial loop. Here, we characterize the TonB-dependent transporters responsible for iron and carbon acquisition in a representative marine copiotroph and examine their distribution across the genus Alteromonas. We provide evidence that substrate-specific bioavailability is niche specific, particularly for iron complexes, indicating that transport capacity may serve as a significant control on microbial community dynamics and the resultant cycling of organic matter.

Iron is an essential micronutrient for all microbial growth in the marine environment, and in heterotrophic bacteria, iron is tightly linked to carbon metabolism due to its central role as a cofactor in enzymes of the respiratory chain. Here, we present the iron- and carbon-regulated transcriptomes of a representative marine copiotroph, Alteromonas macleodii ATCC 27126, and characterize its cellular transport mechanisms. ATCC 27126 has distinct metabolic responses to iron and carbon limitation and, accordingly, uses distinct sets of TonB-dependent transporters for the acquisition of iron and carbon. These distinct sets of TonB-dependent transporters were of a similar number, indicating that the diversity of carbon and iron substrates available to ATCC 27126 is of a similar scale. For the first time in a marine bacterium, we have also identified six characteristic inner membrane permeases for the transport of siderophores via an ATPase-independent mechanism. An examination of the distribution of specific TonB-dependent transporters in 31 genomes across the genus Alteromonas points to niche specialization in transport capacity, particularly for iron. We conclude that the substrate-specific bioavailability of both iron and carbon in the marine environment will likely be a key control on the processing of organic matter through the microbial loop.

IMPORTANCE As the major facilitators of the turnover of organic matter in the marine environment, the ability of heterotrophic bacteria to acquire specific compounds within the diverse range of dissolved organic matter will affect the regeneration of essential nutrients such as iron and carbon. TonB-dependent transporters are a prevalent cellular tool in Gram-negative bacteria that allow a relatively high-molecular-weight fraction of organic matter to be directly accessed. However, these transporters are not well characterized in marine bacteria, limiting our understanding of the flow of specific substrates through the marine microbial loop. Here, we characterize the TonB-dependent transporters responsible for iron and carbon acquisition in a representative marine copiotroph and examine their distribution across the genus Alteromonas. We provide evidence that substrate-specific bioavailability is niche specific, particularly for iron complexes, indicating that transport capacity may serve as a significant control on microbial community dynamics and the resultant cycling of organic matter.

Metagenome Mining Reveals Hidden Genomic Diversity of Pelagimyophages in Aquatic Environments

SAR11 clade members are among the most abundant bacteria on Earth. Their study is complicated by their great diversity and difficulties in being grown and manipulated in the laboratory. On the other hand, and due to their extraordinary abundance, metagenomic data sets provide enormous richness of information about these microbes. Given the major role played by phages in the lifestyle and evolution of prokaryotic cells, the contribution of several new bacteriophage genomes preying on this clade opens windows into the infection strategies and life cycle of its viruses. Such strategies could provide models of attack of large-genome phages preying on streamlined aquatic microbes.

The SAR11 clade is one of the most abundant bacterioplankton groups in surface waters of most of the oceans and lakes. However, only 15 SAR11 phages have been isolated thus far, and only one of them belongs to the Myoviridae family (pelagimyophages). Here, we have analyzed 26 sequences of myophages that putatively infect the SAR11 clade. They have been retrieved by mining ca. 45 Gbp aquatic assembled cellular metagenomes and viromes. Most of the myophages were obtained from the cellular fraction (0.2 μm), indicating a bias against this type of virus in viromes. We have found the first myophages that putatively infect Candidatus Fonsibacter (freshwater SAR11) and another group putatively infecting bathypelagic SAR11 phylogroup Ic. The genomes have similar sizes and maintain overall synteny in spite of low average nucleotide identity values, revealing high similarity to marine cyanomyophages. Pelagimyophages recruited metagenomic reads widely from several locations but always much more from cellular metagenomes than from viromes, opposite to what happens with pelagipodophages. Comparing the genomes resulted in the identification of a hypervariable island that is related to host recognition. Interestingly, some genes in these islands could be related to host cell wall synthesis and coinfection avoidance. A cluster of curli-related proteins was widespread among the genomes, although its function is unclear.

IMPORTANCE SAR11 clade members are among the most abundant bacteria on Earth. Their study is complicated by their great diversity and difficulties in being grown and manipulated in the laboratory. On the other hand, and due to their extraordinary abundance, metagenomic data sets provide enormous richness of information about these microbes. Given the major role played by phages in the lifestyle and evolution of prokaryotic cells, the contribution of several new bacteriophage genomes preying on this clade opens windows into the infection strategies and life cycle of its viruses. Such strategies could provide models of attack of large-genome phages preying on streamlined aquatic microbes.

Genomic, metabolic and phenotypic variability shapes ecological differentiation and intraspecies interactions of Alteromonas macleodii

Ecological differentiation between strains of bacterial species is shaped by genomic and metabolic variability. However, connecting genotypes to ecological niches remains a major challenge. Here, we linked bacterial geno- and phenotypes by contextualizing pangenomic, exometabolomic and physiological evidence in twelve strains of the marine bacterium Alteromonas macleodii, illuminating adaptive strategies of carbon metabolism, microbial interactions, cellular communication and iron acquisition. In A. macleodii strain MIT1002, secretion of amino acids and the unique capacity for phenol degradation may promote associations with Prochlorococcus cyanobacteria. Strain 83-1 and three novel Pacific isolates, featuring clonal genomes despite originating from distant locations, have profound abilities for algal polysaccharide utilization but without detrimental implications for Ecklonia macroalgae. Degradation of toluene and xylene, mediated via a plasmid syntenic to terrestrial Pseudomonas, was unique to strain EZ55. Benzoate degradation by strain EC673 related to a chromosomal gene cluster shared with the plasmid of A. mediterranea EC615, underlining that mobile genetic elements drive adaptations. Furthermore, we revealed strain-specific production of siderophores and homoserine lactones, with implications for nutrient acquisition and cellular communication. Phenotypic variability corresponded to different competitiveness in co-culture and geographic distribution, indicating linkages between intraspecific diversity, microbial interactions and biogeography. The finding of "ecological microdiversity" helps understanding the widespread occurrence of A. macleodii and contributes to the interpretation of bacterial niche specialization, population ecology and biogeochemical roles.

Ecological and functional capabilities of an uncultured Kordia sp

Cultivable bacteria represent only a fraction of the diversity in microbial communities. However, the official procedures for classification and characterization of a novel prokaryotic species still rely on isolates. Nevertheless, due to single cell genomics, it is possible to retrieve genomes from environmental samples by sequencing them individually, and to assign specific genes to a specific taxon, regardless of their ability to grow in culture. In this study, a complete description was performed for uncultured Kordia sp. TARA_039_SRF, a proposed novel species within the genus Kordia, using culture-independent techniques. The type material was a high-quality draft genome (94.97% complete, 4.65% gene redundancy) co-assembled using ten nearly identical single amplified genomes (SAGs) from surface seawater in the North Indian Ocean during the Tara Oceans Expedition. The assembly process was optimized to obtain the best possible assembly metrics and a less fragmented genome. The closest relative of the species was Kordia periserrulae, which shared 97.56% similarity of the 16S rRNA gene, 75% orthologs and 89.13% average nucleotide identity. The functional potential of the proposed novel species included proteorhodopsin, the ability to incorporate nitrate, cytochrome oxidases with high affinity for oxygen, and CAZymes that were unique features within the genus. Its abundance at different depths and size fractions was also evaluated together with its functional annotation, revealing that its putative ecological niche could be particles of phytoplanktonic origin. It could putatively attach to these particles and consume them while sinking to the deeper and oxygen depleted layers of the North Indian Ocean.

High genomic variability in the plant pathogenic bacterium Pectobacterium parmentieri deciphered from de novo assembled complete genomes

BACKGROUND

Pectobacterium parmentieri is a newly established species within the plant pathogenic family Pectobacteriaceae. Bacteria belonging to this species are causative agents of diseases in economically important crops (e.g. potato) in a wide range of different environmental conditions, encountered in Europe, North America, Africa, and New Zealand. Severe disease symptoms result from the activity of P. parmentieri virulence factors, such as plant cell wall degrading enzymes. Interestingly, we observe significant phenotypic differences among P. parmentieri isolates regarding virulence factors production and the abilities to macerate plants. To establish the possible genomic basis of these differences, we sequenced 12 genomes of P. parmentieri strains (10 isolated in Poland, 2 in Belgium) with the combined use of Illumina and PacBio approaches. De novo genome assembly was performed with the use of SPAdes software, while annotation was conducted by NCBI Prokaryotic Genome Annotation Pipeline.

RESULTS

The pan-genome study was performed on 15 genomes (12 de novo assembled and three reference strains: P. parmentieri CFBP 8475T, P. parmentieri SCC3193, P. parmentieri WPP163). The pan-genome includes 3706 core genes, a high number of accessory (1468) genes, and numerous unique (1847) genes. We identified the presence of well-known genes encoding virulence factors in the core genome fraction, but some of them were located in the dispensable genome. A significant fraction of horizontally transferred genes, virulence-related gene duplications, as well as different CRISPR arrays were found, which can explain the observed phenotypic differences. Finally, we found also, for the first time, the presence of a plasmid in one of the tested P. parmentieri strains isolated in Poland.

CONCLUSIONS

We can hypothesize that a large number of the genes in the dispensable genome and significant genomic variation among P. parmentieri strains could be the basis of the potential wide host range and widespread diffusion of P. parmentieri. The obtained data on the structure and gene content of P. parmentieri strains enabled us to speculate on the importance of high genomic plasticity for P. parmentieri adaptation to different environments.

Genome Variation in the Model Halophilic Bacterium Salinibacter ruber

The halophilic bacterium Salinibacter ruber is an abundant and ecologically important member of halophilic communities worldwide. Given its broad distribution and high intraspecific genetic diversity, S. ruber is considered one of the main models for ecological and evolutionary studies of bacterial adaptation to hypersaline environments. However, current insights on the genomic diversity of this species is limited to the comparison of the genomes of two co-isolated strains. Here, we present a comparative genomic analysis of eight S. ruber strains isolated at two different time points in each of two different Mediterranean solar salterns. Our results show an open pangenome with contrasting evolutionary patterns in the core and accessory genomes. We found that the core genome is shaped by extensive homologous recombination (HR), which results in limited sequence variation within population clusters. In contrast, the accessory genome is modulated by horizontal gene transfer (HGT), with genomic islands and plasmids acting as gateways to the rest of the genome. In addition, both types of genetic exchange are modulated by restriction and modification (RM) or CRISPR-Cas systems. Finally, genes differentially impacted by such processes reveal functional processes potentially relevant for environmental interactions and adaptation to extremophilic conditions. Altogether, our results support scenarios that conciliate “Neutral” and “Constant Diversity” models of bacterial evolution.

Homologous Recombination in Core Genomes Facilitates Marine Bacterial Adaptation

Acquisition of ecologically relevant genes is common among ocean bacteria, but whether it has a major impact on genome evolution in marine environments remains unknown. Here, we analyzed the core genomes of 16 phylogenetically diverse and ecologically relevant bacterioplankton lineages, each consisting of up to five genomes varying at the strain level. Statistical approaches identified from each lineage up to ∼50 loci showing anomalously high divergence at synonymous sites, which is best explained by recombination with distantly related organisms. The enriched gene categories in these outlier loci match well with the characteristics previously identified as the key phenotypes of these lineages. Examples are antibiotic synthesis and detoxification in Phaeobacter inhibens, exopolysaccharide production in Alteromonas macleodii, hydrocarbon degradation in Marinobacter hydrocarbonoclasticus, and cold adaptation in Pseudoalteromonas haloplanktis Intriguingly, the outlier loci feature polysaccharide catabolism in Cellulophaga baltica but not in Cellulophaga lytica, consistent with their primary habitat preferences in macroalgae and beach sands, respectively. Likewise, analysis of Prochlorococcus showed that photosynthesis-related genes listed in the outlier loci are found only in the high-light-adapted ecotype and not in the low-light adapted ecotype. These observations strongly suggest that recombination with distant relatives is a key mechanism driving the ecological diversification among marine bacterial lineages.IMPORTANCE Acquisition of new metabolic genes has been known as an important mechanism driving bacterial evolution and adaptation in the ocean, but acquisition of novel alleles of existing genes and its potential ecological role have not been examined. Guided by population genetic theories, our genomic analysis showed that divergent allele acquisition is prevalent in phylogenetically diverse marine bacterial lineages and that the affected loci often encode metabolic functions that underlie the known ecological roles of the lineages under study.

Description of Alteromonas abrolhosensis sp. nov., isolated from sea water of Abrolhos Bank, Brazil

Two Gram-negative, motile, aerobic bacteria isolated from waters of the Abrolhos Bank were classified through a whole genome-based taxonomy. Strains PEL67E^T and PEL68C shared 99% 16S rRNA and dnaK sequence identity with Alteromonas marina SW-47^T and Alteromonas macleodii ATCC 27126^T. In silico DNA-DNA Hybridization, i.e. genome-to-genome distance (GGD), average amino acid identity (AAI) and average nucleotide identity (ANI) showed that PEL67E^T and PEL68C had identity values between 33-36, 86-88 and 83-84%, and 85-86 and 83%, respectively, towards their close neighbors A. macleodii ATCC 27126^T and A. marina SW-47^T. The DNA G + C contents of PEL67E^T and PEL68C were 44.5%. The phenotypic features that differentiate PEL67E^T and PEL68C strains from their close neighbors were assimilation of galactose and activity of phosphatase, and lack of mannitol, maltose, acetate, xylose and glycerol assimilation and lack of lipase, α and β-glucosidase activity. The new species Alteromonas abrolhosensis is proposed. The type strain is PEL67E^T (= CBAS 610^T = CAIM 1925^T).

Benefit from decline: the primary transcriptome of Alteromonas macleodii str. Te101 during Trichodesmium demise

Interactions between co-existing microorganisms deeply affect the physiology of the involved organisms and, ultimately, the function of the ecosystem as a whole. Copiotrophic Alteromonas are marine gammaproteobacteria that thrive during the late stages of phytoplankton blooms in the marine environment and in laboratory co-cultures with cyanobacteria such as Trichodesmium. The response of this heterotroph to the sometimes rapid and transient changes in nutrient supply when the phototroph crashes is not well understood. Here, we isolated and sequenced the strain Alteromonas macleodii str. Te101 from a laboratory culture of Trichodesmium erythraeum IMS101, yielding a chromosome of 4.63 Mb and a single plasmid of 237 kb. Increasing salinities to ≥43 ppt inhibited the growth of Trichodesmium but stimulated growth of the associated Alteromonas. We characterized the transcriptomic responses of both microorganisms and identified the complement of active transcriptional start sites in Alteromonas at single-nucleotide resolution. In replicate cultures, a similar set of genes became activated in Alteromonas when growth rates of Trichodesmium declined and mortality was high. The parallel activation of fliA, rpoS and of flagellar assembly and growth-related genes indicated that Alteromonas might have increased cell motility, growth, and multiple biosynthetic activities. Genes with the highest expression in the data set were three small RNAs (Aln1a-c) that were identified as analogs of the small RNAs CsrB-C in E. coli or RsmX-Z in pathogenic bacteria. Together with the carbon storage protein A (CsrA) homolog Te101_05290, these RNAs likely control the expression of numerous genes in responding to changes in the environment.

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.

Genomic divergence and cohesion in a species of pelagic freshwater bacteria

Background

In many prokaryotic genera a clustered phylogeny is observed, akin to the occurrence of species in sexually reproducing organisms. For some taxa, homologous recombination has been invoked as the underlying mechanism providing genomic cohesion among conspecific individuals. Whether this mechanism is applicable to prokaryotes in freshwaters with low habitat connectivity – i.e. elevated geographic barriers to gene flow – is unclear. To investigate further we studied genomic trends within the globally abundant PnecC cluster (genus Polynucleobacter, Betaproteobacteria) and analyzed homologous recombination within the affiliated species P. asymbioticus.

Results

Comparisons among 20 PnecC genomes revealed a clearly discontinuous distribution of nucleotide sequence similarities. Among the nine conspecific individuals (P. asymbioticus) all average nucleotide identity (ANI) values were greater than 97%, whereas all other comparisons exhibited ANI values lower than 85%. The reconstruction of recombination and mutation events for the P. asymbioticus core genomes yielded an r/m ratio of 7.4, which is clearly above estimated thresholds for recombination to act as a cohesive force. Hotspots of recombination were found to be located in the flanking regions of genomic islands. Even between geographically separated habitats a high flux of recombination was evident. While a biogeographic population structure was suggested from MLST data targeting rather conserved loci, such a structure was barely visible when whole genome data was considered. However, both MLST and whole genome data showed evidence of differentiation between two lineages of P. asymbioticus. The ratios of non-synonymous to synonymous substitution rates as well as growth rates in transplantation experiments suggested that this divergence was not selectively neutral.

Conclusions

The high extent of homologous recombination among P. asymbioticus bacteria can act as a cohesive force that effectively counteracts genetic divergence. At least on a regional scale, homologous recombination can act across geographically separated ecosystems and therefore plays an important role in the evolution and consistency of bacterial freshwater species. A species model akin to the biological species concept may be applicable for P. asymbioticus. Nonetheless, two genetically distinct lineages have emerged and further research may clarify if their divergence has been initiated by reinforced geographical barriers or has been evolving in sympatry.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4199-z) contains supplementary material, which is available to authorized users.

Microdiversification of a Pelagic Polynucleobacter Species Is Mainly Driven by Acquisition of Genomic Islands from a Partially Interspecific Gene Pool

Microdiversification of a planktonic freshwater bacterium was studied by comparing 37 Polynucleobacter asymbioticus strains obtained from three geographically separated sites in the Austrian Alps. Genome comparison of nine strains revealed a core genome of 1.8 Mb, representing 81% of the average genome size. Seventy-five percent of the remaining flexible genome is clustered in genomic islands (GIs). Twenty-four genomic positions could be identified where GIs are potentially located. These positions are occupied strain specifically from a set of 28 GI variants, classified according to similarities in their gene content. One variant, present in 62% of the isolates, encodes a pathway for the degradation of aromatic compounds, and another, found in 78% of the strains, contains an operon for nitrate assimilation. Both variants were shown in ecophysiological tests to be functional, thus providing the potential for microniche partitioning. In addition, detected interspecific horizontal exchange of GIs indicates a large gene pool accessible to Polynucleobacter species. In contrast to core genes, GIs are spread more successfully across spatially separated freshwater habitats. The mobility and functional diversity of GIs allow for rapid evolution, which may be a key aspect for the ubiquitous occurrence of Polynucleobacter bacteria.

IMPORTANCE

Assessing the ecological relevance of bacterial diversity is a key challenge for current microbial ecology. The polyphasic approach which was applied in this study, including targeted isolation of strains, genome analysis, and ecophysiological tests, is crucial for the linkage of genetic and ecological knowledge. Particularly great importance is attached to the high number of closely related strains which were investigated, represented by genome-wide average nucleotide identities (ANI) larger than 97%. The extent of functional diversification found on this narrow phylogenetic scale is compelling. Moreover, the transfer of metabolically relevant genomic islands between more distant members of the Polynucleobacter community provides important insights toward a better understanding of the evolution of these globally abundant freshwater bacteria.

Networking in microbes: conjugative elements and plasmids in the genus Alteromonas

Background

To develop evolutionary models for the free living bacterium Alteromonas the genome sequences of isolates of the genus have been extensively analyzed. However, the main genetic exchange drivers in these microbes, conjugative elements (CEs), have not been considered in detail thus far. In this work, CEs have been searched in several complete Alteromonas genomes and their sequence studied to understand their role in the evolution of this genus. Six genomes are reported here for the first time.

Results

We have found nine different plasmids of sizes ranging from 85 to 600 Kb, most of them were found in a single strain. Networks of gene similarity could be established among six of the plasmids that were also connected with another cluster of plasmids found in Shewanella strains. The cargo genes found in these plasmids included cassettes found before in chromosome flexible genomic islands of Alteromonas strains. We describe also the plasmids pAMCP48-600 and pAMCP49-600, the largest found in Alteromonas thus far (ca. 600 Kb) and containing all the hallmarks to be classified as chromids. We found in them some housekeeping genes and a cluster that code for an exocellular polysaccharide. They could represent the transport vectors for the previously described replacement flexible genomic islands. Integrative and conjugative elements (ICEs) were more common than plasmids and showed similar patterns of variation with cargo genes coding for components of additive flexible genomic islands. A nearly identical ICE was found in A. mediterranea MED64 and Vibrio cholera AHV1003 isolated from a human pathogen, indicating the potential exchange of these genes across phylogenetic distances exceeding the family threshold.

Conclusion

We have seen evidence of how CEs can be vectors to transfer gene cassettes acquired in the chromosomal flexible genomic islands, both of the additive and replacement kind. These CEs showed evidence of how genetic material is exchanged among members of the same species but also (albeit less frequently) across genus and family barriers. These gradients of exchange frequency are probably one of the main drivers of species origin and maintenance in prokaryotes and also provide these taxa with large genetic diversity.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3461-0) contains supplementary material, which is available to authorized users.