Homologous Recombination in Core Genomes Facilitates Marine Bacterial Adaptation

Critically evaluating the relative importance of phage in shaping microbial community composition

The ubiquity of bacteriophages (phages) and the major evolutionary and ecological impacts they can have on their microbial hosts has resulted in phages often cited as key drivers shaping microbial community composition (the relative abundances of species). However, the evidence for the importance of phages is mixed. Here, we critically review the theory and data exploring the role of phages in communities, identifying the conditions when phages are likely to be important drivers of community composition. At ecological scales, we conclude that phages are often followers rather than drivers of microbial population and community dynamics. While phages can affect strain diversity within species, there is yet to be strong evidence suggesting that fluctuations in species' strains affects community composition.

Ecological divergence of syntopic marine bacterial species is shaped by gene content and expression

Identifying mechanisms by which bacterial species evolve and maintain genomic diversity is particularly challenging for the uncultured lineages that dominate the surface ocean. A longitudinal analysis of bacterial genes, genomes, and transcripts during a coastal phytoplankton bloom revealed two co-occurring, highly related Rhodobacteraceae species from the deeply branching and uncultured NAC11-7 lineage. These have identical 16S rRNA gene amplicon sequences, yet their genome contents assembled from metagenomes and single cells indicate species-level divergence. Moreover, shifts in relative dominance of the species during dynamic bloom conditions over 7 weeks confirmed the syntopic species' divergent responses to the same microenvironment at the same time. Genes unique to each species and genes shared but divergent in per-cell inventories of mRNAs accounted for 5% of the species' pangenome content. These analyses uncover physiological and ecological features that differentiate the species, including capacities for organic carbon utilization, attributes of the cell surface, metal requirements, and vitamin biosynthesis. Such insights into the coexistence of highly related and ecologically similar bacterial species in their shared natural habitat are rare.

Effects of thermal desorption on ecotoxicological characteristics of heavy petroleum-contaminated soil

This study comprehensively evaluates the ecotoxicity of high-concentration heavy petroleum (HCHP)-contaminated soil before and after thermal desorption (TD) remediation at different temperatures and times. The results showed that the detoxification of contaminated soil was effectively achieved by extending the remediation duration at 400-600 °C. After treatment at 400 °C for 60 min, the toxicological indicators including bioluminescence EC₅₀ (acute toxicity), seed germination ratio (Gr) and plant biomass of Brassica juncea (subacute toxicity), and diversity of the microbial community (chronic toxicity) reached a maximum. The value of the SOS-Induction Factor (SOSIF), characterizing genotoxicity was below 1.5, indicating that it was non-toxic. Pearson's correlation analysis illustrated that the water-soluble fraction (WSF), ALK1-3 and ARO1-3 of petroleum hydrocarbons were the primary sources of ecotoxicity. Notably, although the total ratio of petroleum removed from the soil reached 87.26 ± 4.38 %-98.69 ± 1.61 % under high-temperature thermal desorption (HTTD, 500-600 °C), the ecotoxicity was not lower than that at 400 °C. The pyrolysis products of petroleum macromolecules and extreme changes in soil properties were the leading causes of soil ecotoxicity following HTTD. The inconsistency between the removal of petroleum pollutants and ecological health risks reveals the significance of soil ecotoxicological assessments for identifying TD remediation endpoints and process optimization.

Coastal Transient Niches Shape the Microdiversity Pattern of a Bacterioplankton Population with Reduced Genomes

Globally dominant marine bacterioplankton lineages are often limited in metabolic versatility, owing to their extensive genome reductions, and thus cannot take advantage of transient nutrient patches. It is therefore perplexing how the nutrient-poor bulk seawater sustains the pelagic streamlined lineages, each containing numerous populations. Here, we sequenced the genomes of 33 isolates of the recently discovered CHUG lineage (~2.6 Mbp), which have some of the smallest genomes in the globally abundant Roseobacter group (commonly over 4 Mbp). These genome-reduced bacteria were isolated from a transient habitat: seawater surrounding the brown alga, Sargassum hemiphyllum. Population genomic analyses showed that: (i) these isolates, despite sharing identical 16S rRNA genes, were differentiated into several genetically isolated populations through successive speciation events; (ii) only the first speciation event led to the genetic separation of both core and accessory genomes; and (iii) populations resulting from this event are differentiated at many loci involved in carbon utilization and oxygen respiration, corroborated by BiOLOG phenotype microarray assays and oxygen uptake kinetics experiments, respectively. These differentiated traits match well with the dynamic nature of the macroalgal seawater, in which the quantity and quality of carbon sources and the concentration of oxygen likely vary spatially and temporally, though other habitats, like fresh organic aggregates, cannot be ruled out. Our study implies that transient habitats in the overall nutrient-poor ocean can shape the microdiversity and population structure of genome-reduced bacterioplankton lineages.

Differentiated Evolutionary Strategies of Genetic Diversification in Atlantic and Pacific Thaumarchaeal Populations

Some marine microbes are seemingly “ubiquitous,” thriving across a wide range of environmental conditions. While the increased depth in metagenomic sequencing has led to a growing body of research on within-population heterogeneity in environmental microbial populations, there have been fewer systematic comparisons and characterizations of population-level genetic diversity over broader expanses of time and space. Here, we investigated the factors that govern the diversification of ubiquitous microbial taxa found within and between ocean basins. Specifically, we use mapped metagenomic paired reads to examine the genetic diversity of ammonia-oxidizing archaeal (“Candidatus Nitrosopelagicus brevis”) populations in the Pacific (Hawaii Ocean Time-series [HOT]) and Atlantic (Bermuda Atlantic Time Series [BATS]) Oceans sampled over 2 years. We observed higher nucleotide diversity in “Ca. N. brevis” at HOT, driven by a higher rate of homologous recombination. In contrast, “Ca. N. brevis” at BATS featured a more open pangenome with a larger set of genes that were specific to BATS, suggesting a history of dynamic gene gain and loss events. Furthermore, we identified highly differentiated genes that were regulatory in function, some of which exhibited evidence of recent selective sweeps. These findings indicate that different modes of genetic diversification likely incur specific adaptive advantages depending on the selective pressures that they are under. Within-population diversity generated by the environment-specific strategies of genetic diversification is likely key to the ecological success of “Ca. N. brevis.”

IMPORTANCE Ammonia-oxidizing archaea (AOA) are one of the most abundant chemolithoautotrophic microbes in the marine water column and are major contributors to global carbon and nitrogen cycling. Despite their ecological importance and geographical pervasiveness, there have been limited systematic comparisons and characterizations of their population-level genetic diversity over time and space. Here, we use metagenomic time series from two ocean observatories to address the fundamental questions of how abiotic and biotic factors shape the population-level genetic diversity and how natural microbial populations adapt across diverse habitats. We show that the marine AOA “Candidatus Nitrosopelagicus brevis” in different ocean basins exhibits distinct modes of genetic diversification in response to their selective regimes shaped by nutrient availability and patterns of environmental fluctuations. Our findings specific to “Ca. N. brevis” have broader implications, particularly in understanding the population-level responses to the changing climate and predicting its impact on biogeochemical cycles.

Cryptic niche differentiation of novel sediment ecotypes of Rugeria pomeroyi correlates with nitrate respiration

Marine intertidal sediments fluctuate in redox conditions and nutrient availability, and they are also known as an important sink of nitrogen mainly through denitrification, yet how denitrifying bacteria adapt to this dynamic habitat remains largely untapped. Here, we investigated novel intertidal benthic ecotypes of the model pelagic marine bacterium Ruegeria pomeroyi DSS-3 with a population genomic approach. While differing by only 1.3% at the 16S rRNA gene level, members of the intertidal benthic ecotypes are complete denitrifiers whereas the pelagic ecotype representative (DSS-3) is a partial denitrifier lacking a nitrate reductase. The intertidal benthic ecotypes are further differentiated by using non-homologous nitrate reductases and a different set of genes that allow alleviating oxidative stress and acquiring organic substrates. In the presence of nitrate, the two ecotypes showed contrasting growth patterns under initial oxygen concentrations at 1 vol% versus 7 vol% and supplemented with different carbon sources abundant in intertidal sediments. Collectively, this combination of evidence indicates that there are cryptic niches in coastal intertidal sediments that support divergent evolution of denitrifying bacteria. This knowledge will in turn help understand how these benthic environments operate to effectively remove nitrogen.

Population differentiation of Rhodobacteraceae along with coral compartments

Coral mucus, tissue, and skeleton harbor compositionally different microbiota, but how these coral compartments shape the microbial evolution remains unexplored. Here, we sampled bacteria inhabiting a prevalent coral species Platygyra acuta and sequenced genomes of 234 isolates comprising two populations in Rhodobacteraceae, an alphaproteobacterial lineage representing a significant but variable proportion (5-50%) of the coral microbiota. The Ruegeria population (20 genomes) contains three clades represented by eight, six, and six isolates predominantly sampled from the skeleton (outgroup), mucus (clade-M), and skeleton (clade-S), respectively. The clade-M possesses functions involved in the utilization of coral osmolytes abundant in the mucus (e.g., methylamines, DMSP, taurine, and L-proline), whereas the clade-S uniquely harbors traits that may promote adaptation to the low-energy and diurnally anoxic skeleton (e.g., sulfur oxidation and swimming motility). These between-clade genetic differences were largely supported by physiological assays. Expanded analyses by including genomes of 24 related isolates (including seven new genomes) from other marine environments suggest that clade-M and clade-S may have diversified in non-coral habitats, but they also consolidated a key role of distinct coral compartments in diversifying many of the above-mentioned traits. The unassigned Rhodobacteraceae population (214 genomes) varies only at a few dozen nucleotide sites across the whole genomes, but the number of between-compartment migration events predicted by the Slatkin-Maddison test supported that dispersal limitation between coral compartments is another key mechanism diversifying microbial populations. Collectively, our results suggest that different coral compartments represent ecologically distinct and microgeographically separate habitats that drive the evolution of the coral microbiota.

Comparative genomics reveals broad genetic diversity, extensive recombination and nascent ecological adaptation in Micrococcus luteus

Background

Micrococcus luteus is a group of actinobacteria that is widely used in biotechnology and is being thought as an emerging nosocomial pathogen. With one of the smallest genomes of free-living actinobacteria, it is found in a wide range of environments, but intraspecies genetic diversity and adaptation strategies to various environments remain unclear. Here, comparative genomics, phylogenomics, and genome-wide association studies were used to investigate the genomic diversity, evolutionary history, and the potential ecological differentiation of the species.

Results

High-quality genomes of 66 M. luteus strains were downloaded from the NCBI GenBank database and core and pan-genome analysis revealed a considerable intraspecies heterogeneity. Phylogenomic analysis, gene content comparison, and average nucleotide identity calculation consistently indicated that the species has diverged into three well-differentiated clades. Population structure analysis further suggested the existence of an unknown ancestor or the fourth, yet unsampled, clade. Reconstruction of gene gain/loss events along the evolutionary history revealed both early events that contributed to the inter-clade divergence and recent events leading to the intra-clade diversity. We also found convincing evidence that recombination has played a key role in the evolutionary process of the species, with upto two-thirds of the core genes having been affected by recombination. Furthermore, distribution of mammal-associated strains (including pathogens) on the phylogenetic tree suggested that the last common ancestor had a free-living lifestyle, and a few recently diverged lineages have developed a mammal-associated lifestyle separately. Consistently, genome-wide association analysis revealed that mammal-associated strains from different lineages shared genes functionally relevant to the host-associated lifestyle, indicating a recent ecological adaption to the new host-associated habitats.

Conclusions

These results revealed high intraspecies genomic diversity of M. luteus and highlighted that gene gain/loss events and extensive recombination events played key roles in the genome evolution. Our study also indicated that, as a free-living species, some lineages have recently developed or are developing a mammal-associated lifestyle. This study provides insights into the mechanisms that drive the genome evolution and adaption to various environments of a bacterial species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07432-5.

Microbial community responses to soil parameters and their effects on petroleum degradation during bio-electrokinetic remediation

This study evaluated the interactions among total petroleum hydrocarbons (TPH), soil parameters, and microbial communities during the bio-electrokinetic (BIO-EK) remediation process. The study was conducted on a petroleum-contaminated saline-alkali soil inoculated with petroleum-degrading bacteria with a high saline-alkali resistance. The results showed that the degradation of TPH was better explained by second-order kinetics, and the efficacy and sustainability of the BIO-EK were closely related to soil micro-environmental factors and microbial community structures. During a 98-d remediation process, the removal rate of TPH was highest in the first 35 d, and then decreased gradually in the later period, which was concurrent with changes in the soil physicochemical properties (conductivity, inorganic ions, pH, moisture, and temperature) and subsequent shifts in the microbial community structures. According to the redundancy analysis (RDA), TPH, soil temperature, and electric conductivity, as well as SO₄^2-, Cl^-, and K⁺ played a better role in explaining the changes in the microbial community at 0-21 d. However, pH and NO₃^- better explained the changes in the microbial community at 63-98 d. In particular, the dominant genera, Marinobacter and Bacillus, showed a positive correlation with TPH, conductivity, and SO₄^2-, Cl^-, and K⁺, but a negative relationship with pH and NO₃. Rhodococcus was positively correlated with soil temperature. The efficacy and sustainability of the BIO-EK remediation process is likely to be improved by controlling these properties.

A hydrogen-oxidizing bacterium enriched from the open ocean resembling a symbiont

A new autotrophic hydrogen-oxidizing Chromatiaceae bacterium, namely bacterium CTD079, was enriched from a water column sample at 1500 m water depth in the southern Pacific Ocean. Based on the phylogeny of 16S rRNA genes, it was closely related to a scaly snail endosymbiont (99.2% DNA sequence identity) whose host so far is only known to colonize hydrothermal vents along the Indian ridge. The average nucleotide identity between the genomes of CTD079 and the snail endosymbiont was 91%. The observed differences likely reflect adaptations to their specific habitats. For example, CTD079 encodes additional enzymes like the formate dehydrogenase increasing the organism's spectrum of energy generation pathways. Other additional physiological features of CTD079 included the increase of viral defence strategies, secretion systems and specific transporters for essential elements. These important genome characteristics suggest an adaptation to life in the open ocean.

Multiple Megaplasmids Confer Extremely High Levels of Metal Tolerance in Alteromonas Strains

Alteromonas is a widely distributed genus of marine Gammaproteobacteria, with representatives shown to be key players in diverse processes, including biogeochemical cycling and biofouling of marine substrata. While Alteromonas spp. are early colonizers of copper-based antifouling paints on marine vessels, their mechanism of tolerance is poorly understood. PacBio whole-genome sequencing of Alteromonas macleodii strains CUKW and KCC02, isolated from Cu/Ni alloy test coupons submerged in oligotrophic coastal waters, indicated the presence of multiple megaplasmids (ca. 200 kb) in both. A pulsed-field gel electrophoresis method was developed and used to confirm the presence of multiple megaplasmids in these two strains; it was then used to screen additional Alteromonas strains for which little to no sequencing data exist. Plasmids were not detected in any of the other strains. Bioinformatic analysis of the CUKW and KCC02 plasmids identified numerous genes associated with metal resistance. Copper resistance orthologs from both the Escherichia coli Cue and Cus and Pseudomonas syringae Cop systems were present, at times as multiple copies. Metal growth assays in the presence of copper, cobalt, manganese, and zinc performed with 10 Alteromonas strains demonstrated the ability of CUKW and KCC02 to grow at metal concentrations inhibitory to all the other strains tested. This study reports multiple megaplasmids in Alteromonas strains. Bioinformatic analysis of the CUKW and KCC02 plasmids indicate that they harbor elements of the Tra system conjugation apparatus, although their type of mobility remains to be experimentally verified.IMPORTANCE Copper is commonly used as an antifouling agent on ship hulls. Alteromonas spp. are early colonizers of copper-based antifouling paint, but their mechanism of tolerance is poorly understood. Sequencing of A. macleodii strains isolated from copper test materials for marine ships indicated the presence of multiple megaplasmids. Plasmids serve as key vectors in horizontal gene transfer and confer traits such as metal resistance, detoxification, ecological interaction, and antibiotic resistance. Bioinformatic analysis identified many metal resistance genes and genes associated with mobility. Understanding the molecular mechanisms and capacity for gene transfer within marine biofilms provides a platform for the development of novel antifouling solutions targeting genes involved in copper tolerance and biofilm formation.

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.