Single-cell genomics reveals hundreds of coexisting subpopulations in wild Prochlorococcus

Large diversity in the O-chain biosynthetic cluster within populations of Pelagibacterales

Genomic diversity in prokaryotic species is largely due to the existence of extensive pangenomes, allowing different gene complements to be drawn depending on the strain. Here, we have studied the diversity of the O-chain polysaccharide biosynthesis cluster (OBC) in marine bacteria of the Pelagibacterales order as a proxy to measure such genetic diversity in a single population. The study of single-amplified genomes (SAGs) from the whole order found a pattern similar to that of other well-studied microbes, such as the Enterobacteriales or Alteromonas, where distinct OBCs represent strains containing different gene pools. We found that most of the OBC sharing happened among individuals of the same clonal frame (>99% average nucleotide identity). Moreover, given the parsimonious way this cluster changes, the diversity of the OBCs can be extrapolated to the size of the population's pangenome. This assumes that different OBCs correspond to lineages containing unique flexible gene pools, as seen in the aforementioned microbes. Through long-read metagenomics, we could detect 380 different OBCs at a single Mediterranean sampling site. Within a single population (single species and sample) of the endemic Ia.3/VII (gMED) genomospecies, we identified 158 OBCs, of which 130 were unique. These findings suggest that the gene pool within a single population might be substantial (several thousands). While this figure is large, it aligns with the complexity of the dissolved organic matter that these organisms can potentially degrade.IMPORTANCEDifferent strains of the same bacterial species contain very different gene pools. This has been long known by epidemiologists. However, it is unknown what gene pool is present in a single set of environmental conditions, i.e., the same time and place in free-living bacteria. Here, we have leveraged information from SAGs to analyze the diversity of the gene cluster coding for the O-chain polysaccharide, a typical component of the flexible gene pool classically used as a tool to differentiate strains in clinical microbiology. It evolves at a similar rate to the rest of the genome and does not seem to be affected by an arms race with phages. One single species of Pelagibacteriales (gMED) revealed an astounding diversity in one sample studied by long-read metagenomics. Our results point to a large gene pool (local pangenome) present in a single population, which is critical to interpreting the biological meaning of the pangenome, i.e., it provides intrapopulation diversity rather than characterizing strains with different distribution in time and/or space.

Intraspecific Diversity in Thermal Performance Determines Phytoplankton Ecological Niche

Temperature has a primary influence on phytoplankton physiology and ecology. We grew 12 strains of Gephyrocapsa huxleyi isolated from different-temperature regions for ~45 generations (2 months) and characterised acclimated thermal response curves across a temperature range. Even with similar temperature optima and overlapping cell size, strain growth rates varied between 0.45 and 1 day-1. Thermal niche widths varied from 16.7°C to 24.8°C, suggesting that strains use distinct thermal response mechanisms. We investigated the implications of this thermal intraspecific diversity using an ocean ecosystem simulation resolving phytoplankton thermal phenotypes. Model analogues of thermal 'generalists' and 'specialists' resulted in a distinctive global biogeography of thermal niche widths with a nonlinear latitudinal pattern. We leveraged model output to predict ranges of the 12 lab-reared strains and demonstrated how this approach could broadly refine geographic range predictions. Our combination of observations and modelled biogeography highlights the capacity of diverse groups to survive temperature shifts.

How Big Is Big? The Effective Population Size of Marine Bacteria

Genome-reduced bacteria constitute most of the cells in surface-ocean bacterioplankton communities. Their extremely large census population sizes (N c) have been unfoundedly translated to huge effective population sizes (N e)-the size of an ideal population carrying as much neutral genetic diversity as the actual population. As N e scales inversely with the strength of genetic drift, constraining the magnitude of N e is key to evaluating whether natural selection can overcome the power of genetic drift to drive evolutionary events. Determining the N e of extant species requires measuring the genomic mutation rate, a challenging step for most genome-reduced bacterioplankton lineages. Results for genome-reduced Prochlorococcus and CHUG are surprising-their N e values are an order of magnitude lower than those of less abundant lineages carrying large genomes, such as Ruegeria and Vibrio. As bacterioplankton genome reduction commonly occurred in the distant past, appreciating their population genetic mechanisms requires constraining their ancient N e values by other methods.

Isolation and Characterization of Multi-Trait Plant Growth-Promoting Endophytic Bacteria from Scots Pine Tissues

Scots pine (Pinus sylvestris) is a globally significant tree species with considerable economic importance in forestry. A major challenge in afforestation, particularly in stressful environments, is growing seedlings with high viability and stress resistance. Recent studies suggest that Pseudomonas strains can alleviate stress and promote growth in crops, though limited evidence exists for trees. This study aimed to assess the plant growth-promoting (PGP) properties of Pseudomonas strains isolated from Scots pine stems using in vitro assays, and to evaluate their potential as bioinoculants through a two-year long field trial. From over sixty bacterial isolates originating from Scots pine stem tissues, only four were selected as being similar to Pseudomonas bacteria. Through 16S rRNA gene sequencing, the isolates were identified as Pseudomonas putida P57, Pseudomonas lurida P88 and 10-1, and Stenotrophomonas maltophilia P77. All isolates inhibited fungal pathogens Botrytis cinerea and Fusarium sporotrichiella, and exhibited PGP activities including nitrogen fixation and production of IAA (1.24-17.74 mg/l), ammonia (4.06-12.71 μM/ml), and siderophores, with the highest value of 1.44 ± 0.19 for the P. lurida P88 strain. Additionally, the Pseudomonas strains demonstrated phosphate solubilization capacity. We revealed that bioinoculation with strains P57 and P88 enhanced field germination of seeds by 35-45% and increased aerial biomass of two-year-old seedlings by 80-140%. Both strains adhered to seed surface and colonized roots and stems at levels of 2.4-3.2 log CFU/g fresh tissue up to two years post-inoculation. These findings highlight the potential of these bacterial strains as effective bioinoculants for improving Scots pine seedling growth under natural conditions.

Novel isolates expand the physiological diversity of Prochlorococcus and illuminate its macroevolution

Prochlorococcus is a diverse picocyanobacterial genus and the most abundant phototroph on Earth. Its photosynthetic diversity divides it into high-light (HL)- or low-light (LL)-adapted groups representing broad phylogenetic grades-each composed of several monophyletic clades. Here, we physiologically characterize four new Prochlorococcus strains isolated from below the deep chlorophyll maximum in the North Pacific Ocean. We combine these physiological properties with genomic analyses to explore the evolution of photosynthetic antennae and discuss potential macroevolutionary implications. The isolates belong to deeply branching low-light-adapted clades that have no other cultivated representatives and display some unusual characteristics. For example, despite its otherwise low-light-adapted physiological characteristics, strain MIT1223 has low chl b2 content similar to high-light-adapted strains. Isolate genomes revealed that each strain contains a unique arsenal of pigment biosynthesis and binding alleles that have been horizontally acquired, contributing to the observed physiological diversity. Comparative genomic analysis of all picocyanobacteria reveals that Pcb, the major pigment carrying protein in Prochlorococcus, greatly increased in copy number and diversity per genome along a branch that coincides with the loss of facultative particle attachment. Collectively, these observations support a recently developed macroevolutionary model, in which niche-constructing radiations allowed ancestral lineages of picocyanobacteria to transition from a particle-attached to planktonic lifestyle and broadly colonize the euphotic zone.IMPORTANCEThe marine cyanobacterium, Prochlorococcus, is among the Earth's most abundant organisms, and much of its genetic and physiological diversity remains uncharacterized. Although field studies help reveal the scope of diversity, cultured isolates allow us to link genomic potential to physiological processes, illuminate eco-evolutionary feedbacks, and test theories arising from comparative genomics of wild cells. Here, we report the isolation and characterization of novel low-light (LL)-adapted Prochlorococcus strains that fill in multiple evolutionary gaps. These new strains are the first cultivated representatives of the LLVII and LLVIII paraphyletic grades of Prochlorococcus, which are broadly distributed in the lower regions of the ocean euphotic zone. Each of these grades is a unique, highly diverse section of the Prochlorococcus tree that separates distinct ecological groups: the LLVII grade branches between monophyletic clades that have facultatively particle-associated and constitutively planktonic lifestyles, whereas the LLVIII grade lies along the branch that leads to all high-light (HL)-adapted clades. Characterizing strains and genomes from these grades yields insights into the large-scale evolution of Prochlorococcus. The new LLVII and LLVIII strains are adapted to growth at very low irradiance levels and possess unique light-harvesting gene signatures and pigmentation. The LLVII strains represent the most basal Prochlorococcus group with a major expansion in photosynthetic antenna genes. Furthermore, a strain from the LLVIII grade challenges the paradigm that all LL-adapted Prochlorococcus exhibit high ratios of chl b:a2. These findings provide insights into the photophysiological evolution of Prochlorococcus and redefine what it means to be a low- vs high-light-adapted Prochlorococcus cell.

The shifts in microbial interactions and gene expression caused by temperature and nutrient loading influence Raphidiopsis raciborskii blooms

Climate change and the trophic status of water bodies are important factors in global occurrence of cyanobacterial blooms. The aim of this study was to explore the cyanobacteria‒bacterial interactions that occur during Raphidiopsis raciborskii (R. raciborskii) blooms by conducting microcosm simulation experiments at different temperatures (20 °C and 30 °C) and with different phosphorus concentrations (0.01 mg/L and 1 mg/L) using an ecological model of microbial behavior and by analyzing microbial self-regulatory strategies using weighted gene coexpression network analysis (WGCNA). Three-way ANOVA revealed significant effects of temperature and phosphorus on the growth of R. raciborskii (P < 0.001). The results of a metagenomics-based analysis of bacterioplankton revealed that the synergistic effects of both climate and trophic changes increased the ability of R. raciborskii to compete with other cyanobacteria for dominance in the cyanobacterial community. The antagonistic effects of climate and nutrient changes favored the occurrence of R. raciborskii blooms, especially in eutrophic waters at approximately 20 °C. The species diversity and richness indices differed between the eutrophication treatment group at 20 °C and the other treatment groups. The symbiotic bacterioplankton network revealed the complexity and stability of the symbiotic bacterioplankton network during blooms and identified the roles of key species in the network. The study also revealed a complex pattern of interactions between cyanobacteria and non-cyanobacteria dominated by altruism, as well as the effects of different behavioral patterns on R. raciborskii bloom occurrence. Furthermore, this study revealed self-regulatory strategies that are used by microbes in response to the dual pressures of temperature and nutrient loading. These results provide important insights into the adaptation of microbial communities in freshwater ecosystems to environmental change and provide useful theoretical support for aquatic environmental management and ecological restoration efforts.

Strain-specific metabarcoding reveals rapid evolution of copper tolerance in populations of the coastal diatom Skeletonema marinoi

Phytoplankton have short generation times, flexible reproduction strategies, large population sizes and high standing genetic diversity, traits that should facilitate rapid evolution under directional selection. We quantified local adaptation of copper tolerance in a population of the diatom Skeletonema marinoi from a mining-exposed inlet in the Baltic Sea and in a non-exposed population 100 km away. We hypothesized that mining pollution has driven evolution of elevated copper tolerance in the impacted population of S. marinoi. Assays of 58 strains originating from sediment resting stages revealed no difference in the average tolerance to copper between the two populations. However, variation within populations was greater at the mining site, with three strains displaying hyper-tolerant phenotypes. In an artificial evolution experiment, we used a novel intraspecific metabarcoding locus to track selection and quantify fitness of all 58 strains during co-cultivation in one control and one toxic copper treatment. As expected, the hyper-tolerant strains enabled rapid evolution of copper tolerance in the mining-exposed population through selection on available strain diversity. Within 42 days, in each experimental replicate a single strain dominated (30%-99% abundance) but different strains dominated the different treatments. The reference population developed tolerance beyond expectations primarily due to slowly developing plastic response in one strain, suggesting that different modes of copper tolerance are present in the two populations. Our findings provide novel empirical evidence that standing genetic diversity of phytoplankton resting stage allows populations to evolve rapidly (20-50 generations) and flexibly on timescales relevant for seasonal bloom progressions.

Modern microbiology: Embracing complexity through integration across scales

Microbes were the only form of life on Earth for most of its history, and they still account for the vast majority of life's diversity. They convert rocks to soil, produce much of the oxygen we breathe, remediate our sewage, and sustain agriculture. Microbes are vital to planetary health as they maintain biogeochemical cycles that produce and consume major greenhouse gases and support large food webs. Modern microbiologists analyze nucleic acids, proteins, and metabolites; leverage sophisticated genetic tools, software, and bioinformatic algorithms; and process and integrate complex and heterogeneous datasets so that microbial systems may be harnessed to address contemporary challenges in health, the environment, and basic science. Here, we consider an inevitably incomplete list of emergent themes in our discipline and highlight those that we recognize as the archetypes of its modern era that aim to address the most pressing problems of the 21st century.

Genomic representativeness and chimerism in large collections of SAGs and MAGs of marine prokaryoplankton

BACKGROUND

Single amplified genomes (SAGs) and metagenome-assembled genomes (MAGs) are the predominant sources of information about the coding potential of uncultured microbial lineages, but their strengths and limitations remain poorly understood. Here, we performed a direct comparison of two previously published collections of thousands of SAGs and MAGs obtained from the same, global environment.

RESULTS

We found that SAGs were less prone to chimerism and more accurately reflected the relative abundance and the pangenome content of microbial lineages inhabiting the epipelagic of the tropical and subtropical ocean, as compared to MAGs. SAGs were also better suited to link genome information with taxa discovered through 16S rRNA amplicon analyses. Meanwhile, MAGs had the advantage of more readily recovering genomes of rare lineages.

CONCLUSIONS

Our analyses revealed the relative strengths and weaknesses of the two most commonly used genome recovery approaches in environmental microbiology. These considerations, as well as the need for better tools for genome quality assessment, should be taken into account when designing studies and interpreting data that involve SAGs or MAGs. Video Abstract.

Decoding populations in the ocean microbiome

Understanding the characteristics and structure of populations is fundamental to comprehending ecosystem processes and evolutionary adaptations. While the study of animal and plant populations has spanned a few centuries, microbial populations have been under scientific scrutiny for a considerably shorter period. In the ocean, analyzing the genetic composition of microbial populations and their adaptations to multiple niches can yield important insights into ecosystem function and the microbiome's response to global change. However, microbial populations have remained elusive to the scientific community due to the challenges associated with isolating microorganisms in the laboratory. Today, advancements in large-scale metagenomics and metatranscriptomics facilitate the investigation of populations from many uncultured microbial species directly from their habitats. The knowledge acquired thus far reveals substantial genetic diversity among various microbial species, showcasing distinct patterns of population differentiation and adaptations, and highlighting the significant role of selection in structuring populations. In the coming years, population genomics is expected to significantly increase our understanding of the architecture and functioning of the ocean microbiome, providing insights into its vulnerability or resilience in the face of ongoing global change. Video Abstract.

The global speciation continuum of the cyanobacterium Microcoleus

Speciation is a continuous process driven by genetic, geographic, and ecological barriers to gene flow. It is widely investigated in multicellular eukaryotes, yet we are only beginning to comprehend the relative importance of mechanisms driving the emergence of barriers to gene flow in microbial populations. Here, we explored the diversification of the nearly ubiquitous soil cyanobacterium Microcoleus. Our dataset consisted of 291 genomes, of which 202 strains and eight herbarium specimens were sequenced for this study. We found that Microcoleus represents a global speciation continuum of at least 12 lineages, which radiated during Eocene/Oligocene aridification and exhibit varying degrees of divergence and gene flow. The lineage divergence has been driven by selection, geographical distance, and the environment. Evidence of genetic divergence and selection was widespread across the genome, but we identified regions of exceptional differentiation containing candidate genes associated with stress response and biosynthesis of secondary metabolites.

Tools for microbial single-cell genomics for obtaining uncultured microbial genomes

The advent of next-generation sequencing technologies has facilitated the acquisition of large amounts of DNA sequence data at a relatively low cost, leading to numerous breakthroughs in decoding microbial genomes. Among the various genome sequencing activities, metagenomic analysis, which entails the direct analysis of uncultured microbial DNA, has had a profound impact on microbiome research and has emerged as an indispensable technology in this field. Despite its valuable contributions, metagenomic analysis is a "bulk analysis" technique that analyzes samples containing a wide diversity of microbes, such as bacteria, yielding information that is averaged across the entire microbial population. In order to gain a deeper understanding of the heterogeneous nature of the microbial world, there is a growing need for single-cell analysis, similar to its use in human cell biology. With this paradigm shift in mind, comprehensive single-cell genomics technology has become a much-anticipated innovation that is now poised to revolutionize microbiome research. It has the potential to enable the discovery of differences at the strain level and to facilitate a more comprehensive examination of microbial ecosystems. In this review, we summarize the current state-of-the-art in microbial single-cell genomics, highlighting the potential impact of this technology on our understanding of the microbial world. The successful implementation of this technology is expected to have a profound impact in the field, leading to new discoveries and insights into the diversity and evolution of microbes.

Single-cell genomics of a bloom-forming phytoplankton species reveals population genetic structure across continents

The study of microbial diversity over time and space is fundamental to the understanding of their ecology and evolution. The underlying processes driving these patterns are not fully resolved but can be studied using population genetic approaches. Here we investigated the population genetic structure of Gonyostomum semen, a bloom-forming phytoplankton species, across two continents. The species appears to be expanding in Europe, whereas similar trends are not observed in the USA. Our aim was to investigate if populations of Gonyostomum semen in Europe and in the USA are genetically differentiated, if there is population genetic structure within the continents, and what the potential drivers of differentiation are. To this end, we used a novel method based on single-amplified genomes combined with Restriction-site Associated DNA sequencing that allows de novo genotyping of natural single-cell isolates without the need for culturing. We amplified over 900 single-cell genomes from 25 lake populations across Europe and the USA and identified two distinct population clusters, one in Europe and another in the USA. Low genetic diversity in European populations supports the hypothesized recent expansion of Gonyostomum semen on this continent. Geographic population structure within each continent was associated with differences in environmental variables that may have led to ecological divergence of population clusters. Overall, our results show that single-amplified genomes combined with Restriction-site Associated DNA sequencing can be used to analyze microalgal population structure and differentiation based on single-cell isolates from natural, uncultured samples.

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.

Genomic basis of environmental adaptation in the widespread poly-extremophilic Exiguobacterium group

Delineating cohesive ecological units and determining the genetic basis for their environmental adaptation are among the most important objectives in microbiology. In the last decade, many studies have been devoted to characterizing the genetic diversity in microbial populations to address these issues. However, the impact of extreme environmental conditions, such as temperature and salinity, on microbial ecology and evolution remains unclear so far. In order to better understand the mechanisms of adaptation, we studied the (pan)genome of Exiguobacterium, a poly-extremophile bacterium able to grow in a wide range of environments, from permafrost to hot springs. To have the genome for all known Exiguobacterium type strains, we first sequenced those that were not yet available. Using a reverse-ecology approach, we showed how the integration of phylogenomic information, genomic features, gene and pathway enrichment data, regulatory element analyses, protein amino acid composition, and protein structure analyses of the entire Exiguobacterium pangenome allows to sharply delineate ecological units consisting of mesophilic, psychrophilic, halophilic-mesophilic, and halophilic-thermophilic ecotypes. This in-depth study clarified the genetic basis of the defined ecotypes and identified some key mechanisms driving the environmental adaptation to extreme environments. Our study points the way to organizing the vast microbial diversity into meaningful ecologically units, which, in turn, provides insight into how microbial communities adapt and respond to different environmental conditions in a changing world.

scMicrobe PTA: near complete genomes from single bacterial cells

Microbial genomes produced by standard single-cell amplification methods are largely incomplete. Here, we show that primary template-directed amplification (PTA), a novel single-cell amplification technique, generated nearly complete genomes from three bacterial isolate species. Furthermore, taxonomically diverse genomes recovered from aquatic and soil microbiomes using PTA had a median completeness of 81%, whereas genomes from standard multiple displacement amplification-based approaches were usually <30% complete. PTA-derived genomes also included more associated viruses and biosynthetic gene clusters.

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.

Dual thermal ecotypes coexist within a nearly genetically identical population of the unicellular marine cyanobacterium Synechococcus

The extent and ecological significance of intraspecific functional diversity within marine microbial populations is still poorly understood, and it remains unclear if such strain-level microdiversity will affect fitness and persistence in a rapidly changing ocean environment. In this study, we cultured 11 sympatric strains of the ubiquitous marine picocyanobacterium Synechococcus isolated from a Narragansett Bay (RI) phytoplankton community thermal selection experiment. Thermal performance curves revealed selection at cool and warm temperatures had subdivided the initial population into thermotypes with pronounced differences in maximum growth temperatures. Curiously, the genomes of all 11 isolates were almost identical (average nucleotide identities of >99.99%, with >99% of the genome aligning) and no differences in gene content or single nucleotide variants were associated with either cool or warm temperature phenotypes. Despite a very high level of genomic similarity, sequenced epigenomes for two strains showed differences in methylation on genes associated with photosynthesis. These corresponded to measured differences in photophysiology, suggesting a potential pathway for future mechanistic research into thermal microdiversity. Our study demonstrates that present-day marine microbial populations can harbor cryptic but environmentally relevant thermotypes which may increase their resilience to future rising temperatures.

The Bay of Bengal exposes abundant photosynthetic picoplankton and newfound diversity along salinity-driven gradients

The Bay of Bengal (BoB) is a 2,600,000 km2 expanse in the Indian Ocean upon which many humans rely. However, the primary producers underpinning food chains here remain poorly characterized. We examined phytoplankton abundance and diversity along strong BoB latitudinal and vertical salinity gradients-which have low temperature variation (27-29°C) between the surface and subsurface chlorophyll maximum (SCM). In surface waters, Prochlorococcus averaged 11.7 ± 4.4 × 104  cells ml-1 , predominantly HLII, whereas LLII and 'rare' ecotypes, HLVI and LLVII, dominated in the SCM. Synechococcus averaged 8.4 ± 2.3 × 104  cells ml-1 in the surface, declined rapidly with depth, and population structure of dominant Clade II differed between surface and SCM; Clade X was notable at both depths. Across all sites, Ostreococcus Clade OII dominated SCM eukaryotes whereas communities differentiated strongly moving from Arabian Sea-influenced high salinity (southerly; prasinophytes) to freshwater-influenced low salinity (northerly; stramenopiles, specifically, diatoms, pelagophytes, and dictyochophytes, plus the prasinophyte Micromonas) surface waters. Eukaryotic phytoplankton peaked in the south (1.9 × 104  cells ml-1 , surface) where a novel Ostreococcus was revealed, named here Ostreococcus bengalensis. We expose dominance of a single picoeukaryote and hitherto 'rare' picocyanobacteria at depth in this complex ecosystem where studies suggest picoplankton are replacing larger phytoplankton due to climate change.

Metabolite diversity among representatives of divergent Prochlorococcus ecotypes

IMPORTANCE

Approximately half of the annual carbon fixation on Earth occurs in the surface ocean through the photosynthetic activities of phytoplankton such as the ubiquitous picocyanobacterium Prochlorococcus. Ecologically distinct subpopulations (or ecotypes) of Prochlorococcus are central conduits of organic substrates into the ocean microbiome, thus playing important roles in surface ocean production. We measured the chemical profile of three cultured ecotype strains, observing striking differences among them that have implications for the likely chemical impact of Prochlorococcus subpopulations on their surroundings in the wild. Subpopulations differ in abundance along gradients of temperature, light, and nutrient concentrations, suggesting that these chemical differences could affect carbon cycling in different ocean strata and should be considered in models of Prochlorococcus physiology and marine carbon dynamics.

Microdiversity sustains the distribution of rhizosphere-associated bacterial species from the root surface to the bulk soil region in maize crop fields

Over the years, the microbial community of maize (Zea mays) rhizosphere has been extensively studied; however, the role of microdiversity sustain rhizosphere-associated microbial species distribution from root surface to bulk soil in mature maize is still unclear. Although operational taxonomic units (OTUs) have been used to classify species, amplicon sequence variants (ASVs) have been shown to be effective in representing microdiversity within OTUs at a finer genetic scale. Therefore, the aim of this study was to examine the role of microdiversity in influencing the distribution of rhizosphere-associated microbial species across environmental gradients from root surface to bulk soil at the OTU and ASV levels. Here, the microbial community structures of bulk, loosely bound, and tightly bound soil samples from maize rhizosphere were examined at OTU and ASV levels. The results showed that OTU and ASV methods exhibited similar microbial community structures in rhizosphere. Additionally, different ecotypes with varying distributions and habitat preferences were observed within the same bacterial OTU at the ASV level, indicating a rich bacterial microdiversity. In contrast, the fungal community exhibited low microdiversity, with no significant relationship between fungal microdiversity and persistence and variability. Moreover, the ecotypes observed within the bacterial OTUs were found to be positively or negatively associated with environmental factors, such as soil organic carbon (SOC), NO3 --N, NH4 +-N contents, and pH. Overall, the results showed that the rich microdiversity could sustain the distribution of rhizosphere-associated bacterial species across environmental gradients from root surface to bulk soil. Further genetic analyses of rhizosphere-associated bacterial species could have considerable implications for potential mediation of microdiversity for sustainable crop production.

Global scale phylogeography of functional traits and microdiversity in Prochlorococcus

Prochlorococcus is the most numerically abundant photosynthetic organism in the surface ocean. The Prochlorococcus high-light and warm-water adapted ecotype (HLII) is comprised of extensive microdiversity, but specific functional differences between microdiverse sub-clades remain elusive. Here we characterized both functional and phylogenetic diversity within the HLII ecotype using Bio-GO-SHIP metagenomes. We found widespread variation in gene frequency connected to local environmental conditions. Metagenome-assembled marker genes and genomes revealed a globally distributed novel HLII haplotype defined by adaptation to chronically low P conditions (HLII-P). Environmental correlation analysis revealed different factors were driving gene abundances verses phylogenetic differences. An analysis of cultured HLII genomes and metagenome-assembled genomes revealed a subclade within HLII, which corresponded to the novel HLII-P haplotype. This work represents the first global assessment of the HLII ecotype's phylogeography and corresponding functional differences. These findings together expand our understanding of how microdiversity structures functional differences and reveals the importance of nutrients as drivers of microdiversity in Prochlorococcus.

Uncultured prokaryotic genomes in the spotlight: An examination of publicly available data from metagenomics and single-cell genomics

Owing to the ineffectiveness of traditional culture techniques for the vast majority of microbial species, culture-independent analyses utilizing next-generation sequencing and bioinformatics have become essential for gaining insight into microbial ecology and function. This mini-review focuses on two essential methods for obtaining genetic information from uncultured prokaryotes, metagenomics and single-cell genomics. We analyzed the registration status of uncultured prokaryotic genome data from major public databases and assessed the advantages and limitations of both the methods. Metagenomics generates a significant quantity of sequence data and multiple prokaryotic genomes using straightforward experimental procedures. However, in ecosystems with high microbial diversity, such as soil, most genes are presented as brief, disconnected contigs, and lack association of highly conserved genes and mobile genetic elements with individual species genomes. Although technically more challenging, single-cell genomics offers valuable insights into complex ecosystems by providing strain-resolved genomes, addressing issues in metagenomics. Recent technological advancements, such as long-read sequencing, machine learning algorithms, and in silico protein structure prediction, in combination with vast genomic data, have the potential to overcome the current technical challenges and facilitate a deeper understanding of uncultured microbial ecosystems and microbial dark matter genes and proteins. In light of this, it is imperative that continued innovation in both methods and technologies take place to create high-quality reference genome databases that will support future microbial research and industrial applications.

High-resolution strain-level microbiome composition analysis from short reads

BACKGROUND

Bacterial strains under the same species can exhibit different biological properties, making strain-level composition analysis an important step in understanding the dynamics of microbial communities. Metagenomic sequencing has become the major means for probing the microbial composition in host-associated or environmental samples. Although there are a plethora of composition analysis tools, they are not optimized to address the challenges in strain-level analysis: highly similar strain genomes and the presence of multiple strains under one species in a sample. Thus, this work aims to provide a high-resolution and more accurate strain-level analysis tool for short reads.

RESULTS

In this work, we present a new strain-level composition analysis tool named StrainScan that employs a novel tree-based k-mers indexing structure to strike a balance between the strain identification accuracy and the computational complexity. We tested StrainScan extensively on a large number of simulated and real sequencing data and benchmarked StrainScan with popular strain-level analysis tools including Krakenuniq, StrainSeeker, Pathoscope2, Sigma, StrainGE, and StrainEst. The results show that StrainScan has higher accuracy and resolution than the state-of-the-art tools on strain-level composition analysis. It improves the F1 score by 20% in identifying multiple strains at the strain level.

CONCLUSIONS

By using a novel k-mer indexing structure, StrainScan is able to provide strain-level analysis with higher resolution than existing tools, enabling it to return more informative strain composition analysis in one sample or across multiple samples. StrainScan takes short reads and a set of reference strains as input and its source codes are freely available at https://github.com/liaoherui/StrainScan . Video Abstract.

Coexistence of Dominant Marine Phytoplankton Sustained by Nutrient Specialization

Prochlorococcus and Synechococcus are the two dominant picocyanobacteria in the low-nutrient surface waters of the subtropical ocean, but the basis for their coexistence has not been quantitatively demonstrated. Here, we combine in situ microcosm experiments and an ecological model to show that this coexistence can be sustained by specialization in the uptake of distinct nitrogen (N) substrates at low-level concentrations that prevail in subtropical environments. In field incubations, the response of both Prochlorococcus and Synechococcus to nanomolar N amendments demonstrates N limitation of growth in both populations. However, Prochlorococcus showed a higher affinity to ammonium, whereas Synechococcus was more adapted to nitrate uptake. A simple ecological model demonstrates that the differential nutrient preference inferred from field experiments with these genera may sustain their coexistence. It also predicts that as the supply of NO3- decreases, as expected under climate warming, the dominant genera should undergo a nonlinear shift from Synechococcus to Prochlorococcus, a pattern that is supported by subtropical field observations. Our study suggests that the evolution of differential nutrient affinities is an important mechanism for sustaining the coexistence of genera and that climate change is likely to shift the relative abundance of the dominant plankton genera in the largest biomes in the ocean. IMPORTANCE Our manuscript addresses the following fundamental question in microbial ecology: how do different plankton using the same essential nutrients coexist? Prochlorococcus and Synechococcus are the two dominant picocyanobacteria in the low-nutrient surface waters of the subtropical ocean, which support a significant amount of marine primary production. The geographical distributions of these two organisms are largely overlapping, but the basis for their coexistence in these biomes remains unclear. In this study, we combined in situ microcosm experiments and an ecosystem model to show that the coexistence of these two organisms can arise from specialization in the uptake of distinct nitrogen substrates; Prochlorococcus prefers ammonium, whereas Synechococcus prefers nitrate when these nutrients exist at low concentrations. Our framework can be used for simulating and predicting the coexistence in the future ocean and may provide hints toward understanding other similar types of coexistence.

Linking niche size and phylogenetic signals to predict future soil microbial relative abundances

Bacteria provide ecosystem services (e.g., biogeochemical cycling) that regulate climate, purify water, and produce food and other commodities, yet their distribution and likely responses to change or intervention are difficult to predict. Using bacterial 16S rRNA gene surveys of 1,381 soil samples from the Biomes of Australian Soil Environment (BASE) dataset, we were able to model relative abundances of soil bacterial taxonomic groups and describe bacterial niche space and optima. Hold out sample validated hypothetical causal networks (structural equation models; SEM) were able to predict the relative abundances of bacterial taxa from environmental data and elucidate soil bacterial niche space. By using explanatory SEM properties as indicators of microbial traits, we successfully predicted soil bacterial response, and in turn potential ecosystem service response, to near-term expected changes in the Australian climate. The methods developed enable prediction of continental-scale changes in bacterial relative abundances, and demonstrate their utility in predicting changes in bacterial function and thereby ecosystem services. These capabilities will be strengthened in the future with growing genome-level data.

Metapangenomic investigation provides insight into niche differentiation of methanogenic populations from the subsurface serpentinizing environment, Samail Ophiolite, Oman

Serpentinization reactions produce highly reduced waters that have hyperalkaline pH and that can have high concentrations of H₂ and CH₄. Putatively autotrophic methanogenic archaea have been identified in the subsurface waters of the Samail Ophiolite, Sultanate of Oman, though the strategies to overcome hyperalkaline pH and dissolved inorganic carbon limitation remain to be fully understood. Here, we recovered metagenome assembled genomes (MAGs) and applied a metapangenomic approach to three different Methanobacterium populations to assess habitat-specific functional gene distribution. A Type I population was identified in the fluids with neutral pH, while a Type II and "Mixed" population were identified in the most hyperalkaline fluids (pH 11.63). The core genome of all Methanobacterium populations highlighted potential DNA scavenging techniques to overcome phosphate or nitrogen limitation induced by environmental conditions. With particular emphasis on the Mixed and Type II population found in the most hyperalkaline fluids, the accessory genomes unique to each population reflected adaptation mechanisms suggesting lifestyles that minimize niche overlap. In addition to previously reported metabolic capability to utilize formate as an electron donor and generate intracellular CO₂, the Type II population possessed genes relevant to defense against antimicrobials and assimilating potential osmoprotectants to provide cellular stability. The accessory genome of the Mixed population was enriched in genes for multiple glycosyltransferases suggesting reduced energetic costs by adhering to mineral surfaces or to other microorganisms, and fostering a non-motile lifestyle. These results highlight the niche differentiation of distinct Methanobacterium populations to circumvent the challenges of serpentinization impacted fluids through coexistence strategies, supporting our ability to understand controls on methanogenic lifestyles and adaptations within the serpentinizing subsurface fluids of the Samail Ophiolite.

Intra-population genomic diversity of the bloom-forming cyanobacterium, Aphanizomenon gracile, at low spatial scale

Cyanobacteria are oxygenic photosynthetic bacteria that perform a substantial part of the global primary production. Some species are responsible for catastrophic environmental events, called blooms, which have become increasingly common in lakes and freshwater bodies as a consequence of global changes. Genotypic diversity is considered essential for marine cyanobacterial population, allowing it to cope with spatio-temporal environmental variations and to adapt to specific micro-niches in the ecosystem. This aspect is underestimated in the study of bloom development, however, and given little notice in studies of the ecology of harmful cyanobacteria. Here we compared the genomes of four strains of Aphanizomenon gracile, a species of filamentous toxinogenic cyanobacteria (Nostocales) found worldwide in fresh and brackish water. Millimeter-sized fascicles were isolated from a single water sample and have been maintained in culture since 2010. A comparative study revealed extensive heterogeneity in gene contents, despite similar genome size and high similarity indices. These variations were mainly associated with mobile genetic elements and biosynthetic gene clusters. For some of the latter, metabolomic analysis confirmed the production of related secondary metabolites, such as cyanotoxins and carotenoids, which are thought to play a fundamental role in the cyanobacterial fitness. Altogether, these results demonstrated that an A. gracile bloom could be a highly diverse population at low spatial scale and raised questions about potential exchanges of essential metabolites between individuals.

A compendium of bacterial and archaeal single-cell amplified genomes from oxygen deficient marine waters

Oxygen-deficient marine waters referred to as oxygen minimum zones (OMZs) or anoxic marine zones (AMZs) are common oceanographic features. They host both cosmopolitan and endemic microorganisms adapted to low oxygen conditions. Microbial metabolic interactions within OMZs and AMZs drive coupled biogeochemical cycles resulting in nitrogen loss and climate active trace gas production and consumption. Global warming is causing oxygen-deficient waters to expand and intensify. Therefore, studies focused on microbial communities inhabiting oxygen-deficient regions are necessary to both monitor and model the impacts of climate change on marine ecosystem functions and services. Here we present a compendium of 5,129 single-cell amplified genomes (SAGs) from marine environments encompassing representative OMZ and AMZ geochemical profiles. Of these, 3,570 SAGs have been sequenced to different levels of completion, providing a strain-resolved perspective on the genomic content and potential metabolic interactions within OMZ and AMZ microbiomes. Hierarchical clustering confirmed that samples from similar oxygen concentrations and geographic regions also had analogous taxonomic compositions, providing a coherent framework for comparative community analysis.

Quantifying Cell Heterogeneity and Subpopulations Using Single Cell Metabolomics

Mass spectrometry (MS) has become an indispensable tool for metabolomics studies. However, due to the lack of applicable experimental platforms, suitable algorithm, software, and quantitative analyses of cell heterogeneity and subpopulations, investigating global metabolomics profiling at the single cell level remains challenging. We combined the Single-probe single cell MS (SCMS) experimental technique with a bioinformatics software package, SinCHet-MS (Single Cell Heterogeneity for Mass Spectrometry), to characterize changes of tumor heterogeneity, quantify cell subpopulations, and prioritize the metabolite biomarkers of each subpopulation. As proof of principle studies, two melanoma cancer cell lines, the primary (WM115; with a lower drug resistance) and the metastatic (WM266-4; with a higher drug resistance), were used as models. Our results indicate that after the treatment of the anticancer drug vemurafenib, a new subpopulation emerged in WM115 cells, while the proportion of the existing subpopulations was changed in the WM266-4 cells. In addition, metabolites for each subpopulation can be prioritized. Combining the SCMS experimental technique with a bioinformatics tool, our label-free approach can be applied to quantitatively study cell heterogeneity, prioritize markers for further investigation, and improve the understanding of cell metabolism in human diseases and response to therapy.

Evaluation of the rbcL marker for metabarcoding of marine diatoms and inference of population structure of selected genera

Diatoms are one of the most important phytoplankton groups in the world's oceans. There are responsible for up to 40% of the photosynthetic activity in the Ocean, and they play an important role in the silicon and carbon cycles by decoupling carbon from atmospheric interactions through sinking and export. These processes are strongly influenced by the taxonomic composition of diatom assemblages. Traditionally, these have been assessed using microscopy, which in some cases is not reliable or reproducible. Next-generation sequencing enabled us to study diversity in a high-throughput manner and uncover new distribution patterns and diversity. However, phylogenetic markers used for this purpose, such as various 18S rDNA regions, are often insufficient because they cannot distinguish between some taxa. In this work, we demonstrate the performance of the chloroplast-encoded rbcL marker for metabarcoding marine diatoms compared to microscopy and 18S-V9 metabarcoding using a series of monthly samples from the Gulf of Trieste (GoT), northern Adriatic Sea. We demonstrate that rbcL is able to detect more taxa compared to 18S-V9 metabarcoding or microscopy, while the overall structure of the diatom assemblage was comparable to the other two methods with some variations, that were taxon dependent. In total, 6 new genera and 22 new diatom species for the study region were identified. We were able to spot misidentification of genera obtained with microscopy such as Pseudo-nitzschia galaxiae, which was mistaken for Cylindrotheca closterium, as well as genera that were completely overlooked, such as Minidiscus and several genera from the Cymatosiraceae family. Furthermore, on the example of two well-studied genera in the region, namely Chaetoceros and particularly Pseudo-nitzschia, we show how the rbcL method can be used to infer even deeper phylogenetic and ecologically significant differences at the species population level. Despite a very thorough community analysis obtained by rbcL the incompleteness of reference databases was still evident, and we shed light on possible improvements. Our work has further implications for studies dealing with taxa distribution and population structure, as well as carbon and silica flux models and networks.

Adaptive genetic traits in pelagic freshwater microbes

Pelagic microbes have adopted distinct strategies to inhabit the pelagial of lakes and oceans and can be broadly categorized in two groups: free-living, specialized oligotrophs and patch-associated generalists or copiotrophs. In this review, we aim to identify genomic traits that enable pelagic freshwater microbes to thrive in their habitat. To do so, we discuss the main genetic differences of pelagic marine and freshwater microbes that are both dominated by specialized oligotrophs and the difference to freshwater sediment microbes, where copiotrophs are more prevalent. We phylogenomically analysed a collection of >7700 metagenome-assembled genomes, classified habitat preferences on different taxonomic levels, and compared the metabolic traits of pelagic freshwater, marine, and freshwater sediment microbes. Metabolic differences are mainly associated with transport functions, environmental information processing, components of the electron transport chain, osmoregulation and the isoelectric point of proteins. Several lineages with known habitat transitions (Nitrososphaeria, SAR11, Methylophilaceae, Synechococcales, Flavobacteriaceae, Planctomycetota) and the underlying mechanisms in this process are discussed in this review. Additionally, the distribution, ecology and genomic make-up of the most abundant freshwater prokaryotes are described in details in separate chapters for Actinobacteriota, Bacteroidota, Burkholderiales, Verrucomicrobiota, Chloroflexota, and 'Ca. Patescibacteria'.

Differential global distribution of marine picocyanobacteria gene clusters reveals distinct niche-related adaptive strategies

The ever-increasing number of available microbial genomes and metagenomes provides new opportunities to investigate the links between niche partitioning and genome evolution in the ocean, especially for the abundant and ubiquitous marine picocyanobacteria Prochlorococcus and Synechococcus. Here, by combining metagenome analyses of the Tara Oceans dataset with comparative genomics, including phyletic patterns and genomic context of individual genes from 256 reference genomes, we show that picocyanobacterial communities thriving in different niches possess distinct gene repertoires. We also identify clusters of adjacent genes that display specific distribution patterns in the field (eCAGs) and are thus potentially involved in the same metabolic pathway and may have a key role in niche adaptation. Several eCAGs are likely involved in the uptake or incorporation of complex organic forms of nutrients, such as guanidine, cyanate, cyanide, pyrimidine, or phosphonates, which might be either directly used by cells, for example for the biosynthesis of proteins or DNA, or degraded to inorganic nitrogen and/or phosphorus forms. We also highlight the enrichment of eCAGs involved in polysaccharide capsule biosynthesis in Synechococcus populations thriving in both nitrogen- and phosphorus-depleted areas vs. low-iron (Fe) regions, suggesting that the complexes they encode may be too energy-consuming for picocyanobacteria thriving in the latter areas. In contrast, Prochlorococcus populations thriving in Fe-depleted areas specifically possess an alternative respiratory terminal oxidase, potentially involved in the reduction of Fe(III) to Fe(II). Altogether, this study provides insights into how phytoplankton communities populate oceanic ecosystems, which is relevant to understanding their capacity to respond to ongoing climate change.

Draft genomes of three closely related low light-adapted Prochlorococcus

OBJECTIVES

The marine cyanobacterium Prochlorococcus is a critical part of warm ocean ecosystems and a model for studying microbial evolution and ecology. To expand the representation of this organism's vast wild diversity in sequence collections, we performed a set of isolation efforts targeting low light-adapted Prochlorococcus. Three genomes resulting from this larger body of work are described here.

DATA DESCRIPTION

We present draft-quality Prochlorococcus genomes from enrichment cultures P1344, P1361, and P1363, sampled in the North Pacific. The genomes were built from Illumina paired reads assembled de novo. Supporting datasets of raw reads, assessments, and sequences from co-enriched heterotrophic marine bacteria are also provided. These three genomes represent members of the low light-adapted LLIV Prochlorococcus clade that are closely related, with 99.9% average nucleotide identity between pairs, yet vary in gene content. Expanding the powerful toolkit of Prochlorococcus genomes, these sequences provide an opportunity to study fine-scale variation and microevolutionary processes.

Ten years of marine evolutionary biology-Challenges and achievements of a multidisciplinary research initiative

The Centre for Marine Evolutionary Biology (CeMEB) at the University of Gothenburg, Sweden, was established in 2008 through a 10-year research grant of 8.7 m€ to a team of senior researchers. Today, CeMEB members have contributed >500 scientific publications, 30 PhD theses and have organised 75 meetings and courses, including 18 three-day meetings and four conferences. What are the footprints of CeMEB, and how will the centre continue to play a national and international role as an important node of marine evolutionary research? In this perspective article, we first look back over the 10 years of CeMEB activities and briefly survey some of the many achievements of CeMEB. We furthermore compare the initial goals, as formulated in the grant application, with what has been achieved, and discuss challenges and milestones along the way. Finally, we bring forward some general lessons that can be learnt from a research funding of this type, and we also look ahead, discussing how CeMEB's achievements and lessons can be used as a springboard to the future of marine evolutionary biology.

Global scale phylogeography of functional traits and microdiversity in Prochlorococcus

Prochlorococcus is the most numerically abundant photosynthetic organism in the surface ocean. The Prochlorococcus high-light and warm-water adapted ecotype (HLII) is comprised of extensive microdiversity, but specific functional differences between microdiverse sub-clades remain elusive. Here we characterized both functional and phylogenetic diversity within the HLII ecotype using Bio-GO-SHIP metagenomes. We found widespread variation in gene frequency connected to local environmental conditions. Metagenomically assembled marker genes and genomes revealed a globally distributed novel HLII haplotype defined by adaptation to chronically low P conditions (HLII-P). Environmental correlation analysis revealed different factors were driving gene abundances verses phylogenetic differences. An analysis of cultured HLII genomes and metagenomically assembled genomes revealed a subclade within HLII, which corresponded to the novel HLII-P haplotype. This work represents the first global assessment of the HLII ecotype’s phylogeography and corresponding functional differences. These findings together expand our understanding of how microdiversity structures functional differences and reveals the importance of nutrients as drivers of microdiversity in Prochlorococcus .

Novel integrative elements and genomic plasticity in ocean ecosystems

Horizontal gene transfer accelerates microbial evolution. The marine picocyanobacterium Prochlorococcus exhibits high genomic plasticity, yet the underlying mechanisms are elusive. Here, we report a novel family of DNA transposons-"tycheposons"-some of which are viral satellites while others carry cargo, such as nutrient-acquisition genes, which shape the genetic variability in this globally abundant genus. Tycheposons share distinctive mobile-lifecycle-linked hallmark genes, including a deep-branching site-specific tyrosine recombinase. Their excision and integration at tRNA genes appear to drive the remodeling of genomic islands-key reservoirs for flexible genes in bacteria. In a selection experiment, tycheposons harboring a nitrate assimilation cassette were dynamically gained and lost, thereby promoting chromosomal rearrangements and host adaptation. Vesicles and phage particles harvested from seawater are enriched in tycheposons, providing a means for their dispersal in the wild. Similar elements are found in microbes co-occurring with Prochlorococcus, suggesting a common mechanism for microbial diversification in the vast oligotrophic oceans.

Associations between picocyanobacterial ecotypes and cyanophage host genes across ocean basins and depth

Background

Cyanophages, viruses that infect cyanobacteria, are globally abundant in the ocean's euphotic zone and are a potentially important cause of mortality for marine picocyanobacteria. Viral host genes are thought to increase viral fitness by either increasing numbers of genes for synthesizing nucleotides for virus replication, or by mitigating direct stresses imposed by the environment. The encoding of host genes in viral genomes through horizontal gene transfer is a form of evolution that links viruses, hosts, and the environment. We previously examined depth profiles of the proportion of cyanophage containing various host genes in the Eastern Tropical North Pacific Oxygen Deficient Zone (ODZ) and at the subtropical North Atlantic (BATS). However, cyanophage host genes have not been previously examined in environmental depth profiles across the oceans.

Methodology

We examined geographical and depth distributions of picocyanobacterial ecotypes, cyanophage, and their viral-host genes across ocean basins including the North Atlantic, Mediterranean Sea, North Pacific, South Pacific, and Eastern Tropical North and South Pacific ODZs using phylogenetic metagenomic read placement. We determined the proportion of myo and podo-cyanophage containing a range of host genes by comparing to cyanophage single copy core gene terminase (terL). With this large dataset (22 stations), network analysis identified statistical links between 12 of the 14 cyanophage host genes examined here with their picocyanobacteria host ecotypes.

Results

Picyanobacterial ecotypes, and the composition and proportion of cyanophage host genes, shifted dramatically and predictably with depth. For most of the cyanophage host genes examined here, we found that the composition of host ecotypes predicted the proportion of viral host genes harbored by the cyanophage community. Terminase is too conserved to illuminate the myo-cyanophage community structure. Cyanophage cobS was present in almost all myo-cyanophage and did not vary in proportion with depth. We used the composition of cobS phylotypes to track changes in myo-cyanophage composition.

Conclusions

Picocyanobacteria ecotypes shift with changes in light, temperature, and oxygen and many common cyanophage host genes shift concomitantly. However, cyanophage phosphate transporter gene pstS appeared to instead vary with ocean basin and was most abundant in low phosphate regions. Abundances of cyanophage host genes related to nutrient acquisition may diverge from host ecotype constraints as the same host can live in varying nutrient concentrations. Myo-cyanophage community in the anoxic ODZ had reduced diversity. By comparison to the oxic ocean, we can see which cyanophage host genes are especially abundant (nirA, nirC, and purS) or not abundant (myo psbA) in ODZs, highlighting both the stability of conditions in the ODZ and the importance of nitrite as an N source to ODZ endemic LLV Prochlorococcus.

Fine-scale diversity of microbial communities due to satellite niches in boom and bust environments

Recent observations have revealed that closely related strains of the same microbial species can stably coexist in natural and laboratory settings subject to boom and bust dynamics and serial dilutions, respectively. However, the possible mechanisms enabling the coexistence of only a handful of strains, but not more, have thus far remained unknown. Here, using a consumer-resource model of microbial ecosystems, we propose that by differentiating along Monod parameters characterizing microbial growth rates in high and low nutrient conditions, strains can coexist in patterns similar to those observed. In our model, boom and bust environments create satellite niches due to resource concentrations varying in time. These satellite niches can be occupied by closely related strains, thereby enabling their coexistence. We demonstrate that this result is valid even in complex environments consisting of multiple resources and species. In these complex communities, each species partitions resources differently and creates separate sets of satellite niches for their own strains. While there is no theoretical limit to the number of coexisting strains, in our simulations, we always find between 1 and 3 strains coexisting, consistent with known experiments and observations.

Widespread yet persistent low abundance of TIM5-like cyanophages in the oceans

Ocean ecosystems are inhabited by a diverse set of viruses that impact microbial mortality and evolution. However, the distribution and abundances of specific viral lineages, particularly those from the large bank of rare viruses, remains largely unknown. Here, we assessed the diversity and abundance of the TIM5-like cyanophages. The sequencing of three new TIM5-like cyanophage genomes and environmental amplicons of a signature gene from the Red Sea revealed highly conserved gene content and sequence similarity. We adapted the polony method, a solid-phase polymerase chain reaction assay, to quantify TIM5-like cyanophages during three 2000 km expeditions in the Pacific Ocean and four annual cycles in the Red Sea. TIM5-like cyanophages were widespread, detected at all latitudes and seasons surveyed throughout the photic zone. Yet they were generally rare, ranging between <100 and 4000 viruses·ml^-1 . Occasional peaks in abundance of 10- to 100-fold were observed, reaching 71,000 viruses·ml^-1 . These peaks were ephemeral and seasonally variable in the Red Sea. Infection levels, quantified during one such peak, were very low. These characteristics of low diversity and abundance, as well as variable outbreaks, distinguishes the TIM5-like lineage from other major cyanophage lineages and illuminates that rare virus lineages can be persistent and widespread in the oceans.

Shifts in evolutionary balance of phenotypes under environmental changes

Environments shape communities by driving individual interactions and the evolutionary outcome of competition. In static, homogeneous environments a robust, evolutionary stable, outcome is sometimes reachable. However, inherently stochastic, this evolutionary process need not stabilize, resulting in a dynamic ecological state, often observed in microbial communities. We use evolutionary games to study the evolution of phenotypic competition in dynamic environments. Under the assumption that phenotypic expression depends on the environmental shifts, existing periodic relationships may break or result in formation of new periodicity in phenotypic interactions. The exact outcome depends on the environmental shift itself, indicating the importance of understanding how environments influence affected systems. Under periodic environmental fluctuations, a stable state preserving dominant phenotypes may exist. However, rapid environmental shifts can lead to critical shifts in the phenotypic evolutionary balance. This might lead to environmentally favoured phenotypes dominating making the system vulnerable. We suggest that understanding of the robustness of the system's current state is necessary to anticipate when it will shift to a new equilibrium via understanding what level of perturbations the system can take before its equilibrium changes. Our results provide insights in how microbial communities can be steered to states where they are dominated by desired phenotypes.

A Unique Set of Auxiliary Metabolic Genes Found in an Isolated Cyanophage Sheds New Light on Marine Phage-Host Interactions

Cyanophages, viruses that infect cyanobacteria, are abundant and widely distributed in aquatic ecosystems, playing important roles in regulating the abundance, activity, diversity, and evolution of cyanobacteria. A T4-like cyanophage, S-SCSM1, infecting Synechococcus and Prochlorococcus strains of different ecotypes, was isolated from the South China Sea in this study. For the first time, a mannose-6-phosphate isomerase (MPI) gene was identified in the cultured cyanophage. At least 11 phylogenetic clusters of cyanophage MPIs were retrieved and identified from the marine metagenomic data sets, indicating that cyanophage MPIs in the marine environment are extremely diverse. The existence of 24 genes encoding 2-oxoglutarate (2OG)-Fe(II) oxygenase superfamily proteins in the S-SCSM1 genome emphasizes their potential importance and diverse functions in reprogramming host metabolism during phage infection. Novel cell wall synthesis and modification genes found in the S-SCSM1 genome indicate that diverse phenotypic modifications imposed by phages on cyanobacterial hosts remain to be discovered. Two noncoding RNAs of cis-regulatory elements in the S-SCSM1 genome were predicted to be associated with host exopolysaccharide metabolism and photosynthesis. The isolation and genomic characterization of cyanophage S-SCSM1 provide more information on the genetic diversity of cyanophages and phage-host interactions in the marine environment. IMPORTANCE Cyanophages play important ecological roles in aquatic ecosystems. Genomic and proteomic characterizations of the T4-like cyanophage S-SCSM1 indicate that novel and diverse viral genes and phage-host interactions in the marine environment remain unexplored. The first identified mannose-6-phosphate isomerase (MPI) gene from a cultured cyanophage was found in the S-SCSM1 genome, although MPIs were previously found in viral metagenomes at high frequencies similar to those of the cyanophage photosynthetic gene psbA. The presence of 24 genes encoding 2-oxoglutarate (2OG)-Fe(II) oxygenase superfamily proteins, novel cell wall synthesis and modification genes, a nonbleaching protein A gene, and 2 noncoding RNAs of cis-regulatory elements in the S-SCSM1 genome as well as the presence of a virion-associated regulatory protein indicate the diverse functions that cyanophages have in reprogramming the metabolism and modifying the phenotypes of hosts during infection.

Pan-genome evolution and its association with divergence of metabolic functions in Bifidobacterium genus

Previous studies were mainly focused on genomic evolution and diversity of type species of Bifidobacterium genus due to their health-promoting effect on host. However, those studies were mainly based on species-level taxonomic resolution, adaptation, and characterization of carbohydrate metabolic features of the bifidobacterial species. Here, a comprehensive analysis of the type strain genome unveils the association of pan-genome evolution with the divergence of metabolic function of the Bifidobacterium genus. This study has also demonstrated that horizontal gene transfer, as well as genome expansion and reduction events, leads to the divergence of metabolic functions in Bifidobacterium genus. Furthermore, the genome-based search of probiotic traits among all the available bifidobacterial type strains gives hints on type species, that could confer health benefits to nutrient-deficient individuals. Altogether, the present study provides insight into the developments of genomic evolution, functional divergence, and potential probiotic type species of the Bifidobacterium genus.

High genomic differentiation and limited gene flow indicate recent cryptic speciation within the genus Laspinema (cyanobacteria)

The sympatric occurrence of closely related lineages displaying conserved morphological and ecological traits is often characteristic of free-living microbes. Gene flow, recombination, selection, and mutations govern the genetic variability between these cryptic lineages and drive their differentiation. However, sequencing conservative molecular markers (e.g., 16S rRNA) coupled with insufficient population-level sampling hindered the study of intra-species genetic diversity and speciation in cyanobacteria. We used phylogenomics and a population genomic approach to investigate the extent of local genomic diversity and the mechanisms underlying sympatric speciation of Laspinema thermale. We found two cryptic lineages of Laspinema. The lineages were highly genetically diverse, with recombination occurring more frequently within than between them. That suggests the existence of a barrier to gene flow, which further maintains divergence. Genomic regions of high population differentiation harbored genes associated with possible adaptations to high/low light conditions and stress stimuli, although with a weak diversifying selection. Overall, the diversification of Laspinema species might have been affected by both genomic and ecological processes.

Single-cell genomics reveals population structures from in vitro evolutionary studies of Salmonella

Single-cell DNA sequencing has the potential to reveal detailed hierarchical structures in evolving populations of cells. Single cell approaches are increasingly used to study clonal evolution in human ageing and cancer but have not yet been deployed to study evolving clonal microbial populations. Here, we present an approach for single bacterial genomic analysis for in vitro evolution experiments using FACS isolation of individual bacteria followed by whole-genome amplification and sequencing. We apply this to the experimental evolution of a hypermutator strain of Salmonella in response to antibiotic stress (ciprofloxacin). By analysing sequence polymorphisms in individual cells from populations we identified the presence and prevalence of sub-populations which have acquired polymorphisms in genes previously demonstrated to be associated with ciprofloxacin susceptibility. We were also able to identify that the population exposed to antibiotic stress was able to develop resistance whilst maintaining diversity. This population structure could not be resolved from bulk sequence data, and our results show how high-throughput single-cell sequencing can enhance experimental studies of bacterial evolution.

Cyanophages from a less virulent clade dominate over their sister clade in global oceans

Environmental virus communities are highly diverse. However, the infection physiology underlying the evolution of diverse phage lineages and their ecological consequences are largely unknown. T7-like cyanophages are abundant in nature and infect the marine unicellular cyanobacteria, Synechococcus and Prochlorococcus, important primary producers in the oceans. Viruses belonging to this genus are divided into two distinct phylogenetic clades: clade A and clade B. These viruses have narrow host-ranges with clade A phages primarily infecting Synechococcus genotypes, while clade B phages are more diverse and can infect either Synechococcus or Prochlorococcus genotypes. Here we investigated infection properties (life history traits) and environmental abundances of these two clades of T7-like cyanophages. We show that clade A cyanophages have more rapid infection dynamics, larger burst sizes and greater virulence than clade B cyanophages. However, clade B cyanophages were at least 10-fold more abundant in all seasons, and infected more cyanobacteria, than clade A cyanophages in the Red Sea. Models predicted that steady-state cyanophage abundances, infection frequency, and virus-induced mortality, peak at intermediate virulence values. Our findings indicate that differences in infection properties are reflected in virus phylogeny at the clade level. They further indicate that infection properties, together with differences in subclade diversity and host repertoire, have important ecological consequences with the less aggressive, more diverse virus clade having greater ecological impacts.

The microbiome of a bacterivorous marine choanoflagellate contains a resource-demanding obligate bacterial associate

Microbial predators such as choanoflagellates are key players in ocean food webs. Choanoflagellates, which are the closest unicellular relatives of animals, consume bacteria and also exhibit marked biological transitions triggered by bacterial compounds, yet their native microbiomes remain uncharacterized. Here we report the discovery of a ubiquitous, uncultured bacterial lineage we name Candidatus Comchoanobacterales ord. nov., related to the human pathogen Coxiella and physically associated with the uncultured marine choanoflagellate Bicosta minor. We analyse complete ‘Comchoano’ genomes acquired after sorting single Bicosta cells, finding signatures of obligate host-dependence, including reduction of pathways encoding glycolysis, membrane components, amino acids and B-vitamins. Comchoano encode the necessary apparatus to import energy and other compounds from the host, proteins for host-cell associations and a type IV secretion system closest to Coxiella’s that is expressed in Pacific Ocean metatranscriptomes. Interactions between choanoflagellates and their microbiota could reshape the direction of energy and resource flow attributed to microbial predators, adding complexity and nuance to marine food webs.

Choanoflagellates are the closest living unicellular relatives of animals and are important bacterivorous predators in the ocean. Here the authors show that the microbiome of this predator includes an obligate, host resource-dependent bacterial associate.

Diversity strengthens competing teams

How does the composition of a collection of individuals affect its outcome in competition with other collections of individuals? Assuming that individuals can be different, we develop a model to interpolate between individual-level interactions and collective-level consequences. Rooted in theoretical mathematics, the model is not constrained to any specific context. Potential applications include research, education, sports, politics, ecology, agriculture, algorithms and finance. Our first main contribution is a game theoretic model that interpolates between the internal composition of an ensemble of individuals and the repercussions for the ensemble as a whole in competition with others. The second main contribution is the rigorous identification of all equilibrium points and strategies. These equilibria suggest a mechanistic underpinning for biological and physical systems to tend towards increasing diversity due to the strength it imparts to the system in competition with others.