Evidence for environmental and ecological selection in a microbe with no geographic limits to gene flow

Intraspecific genetic diversity with unrestricted gene flow in the domoic acid-producing diatom Nitzschia navis-varingica (Bacillariophyceae) from the Western Pacific

The benthic pennate diatom Nitzschia navis-varingica, known for producing domoic acid (DA) and its isomers, is widely distributed in the Western Pacific (WP) region. To investigate the genetic differentiation and gene flow patterns among the populations in the WP, the genetic diversity of 354 strains of N. navis-varingica was analysed using two nuclear-encoded rDNA loci: the large subunit rDNA (LSU rDNA) and the internal transcribed spacer 2 (ITS2). Frustule morphology of each strain was examined by TEM. The LSU rDNA phylogeny revealed a monophyletic lineage encompassing all strains, with sequence divergences of <0.9 %. Phylogenetic and population genetic analyses of ITS2 identified eight distinct clades (designated as Groups A to H) with moderate to high genetic heterogeneity (0.5-19.7 %). The low genetic differentiations between the geographically separated populations (pairwise FST of <0.03) suggested high gene flow and lack of spatial genetic structuring. Molecular clock analysis of the ITS2 phylogeny traced the evolutionary history of N. navis-varingica to the Eocene Epoch, and the split between clades likely occurred from the mid-Miocene to Pleistocene Epochs (10.8-1.2 Ma). The population dispersal in the WP were likely influenced by historical events like the Quarternary glacial cycles during the period, contributing to its homogenous distributions in the region.

Genome-wide differentiation corresponds to climatic niches in two species of lichen-forming fungi

Lichens can withstand fluctuating environmental conditions such as hydration-desiccation cycles. Many species distribute across climate zones, suggesting population-level adaptations to conditions such as freezing and drought. Here, we aim to understand how climate affects population genomic patterns in lichenized fungi. We analysed population structure along elevational gradients in closely related Umbilicaria phaea (North American; two gradients) and Umbilicaria pustulata (European; three gradients). All gradients showed clear genomic breaks splitting populations into low-elevation (Mediterranean zone) and high-elevation (cold temperate zone). A total of 3301 SNPs in U. phaea and 138 SNPs in U. pustulata were driven to fixation between the two ends of the gradients. The difference between the species is likely due to differences in recombination rate: the sexually reproducing U. phaea has a higher recombination rate than the primarily asexually reproducing U. pustulata. Cline analysis revealed allele frequency transitions along all gradients at approximately 0°C, coinciding with the transition between the Mediterranean and cold temperate zones, suggesting freezing is a strong driver of population differentiation. Genomic scans further confirmed temperature-related selection targets. Both species showed similar differentiation patterns overall, but different selected alleles indicate convergent adaptation to freezing. Our results enrich our knowledge of fungal genomic functions related to temperature and climate, fungal population genomics, and species responses to environmental heterogeneity.

Population bottlenecks and sexual recombination shape diatom microevolution

Diatoms are single-celled organisms that contribute approximately 20% of the global primary production and play a crucial role in biogeochemical cycles and trophic chains. Despite their ecological importance, our knowledge of microevolution is limited. We developed a model using the SLiM evolutionary framework to address this knowledge gap. As a reference, we used the diatom Pseudo-nitzschia multistriata, which has been extensively studied in the Gulf of Naples. Our model recapitulates what we observe in natural populations, with microevolutionary processes that occur annually during a three-stage bloom phase. Interestingly, we found that non-bloom phases allow the population to maintain sex-generated diversity produced during blooms. This finding suggests that non-bloom phases are critical to counteract bloom-related pressures and mitigate genetic divergence at the species level. Moreover, our model showed that despite the consistent genetic differentiation during bloom phases, the population tends to return to pre-bloom states. While our model is limited to neutral dynamics, our study provides valuable insights into diatoms' microevolution, paving the way to explore the ecological implications of the life history dynamics of these organisms.

Connectivity mediates the spatial ecological impacts of a glyphosate-based herbicide in experimental metaecosystems

Metacommunity ecology has shown that connectivity is important for the persistence of a species locally and across connected ecosystems, however we do not know if ecological effects in freshwater ecosystems exposed to biocides leaking from agriculture depend on metaecosystem connectivity. We experimentally replicated metaecosystems in the laboratory using gradostats as a model system. We tested the effects of connectivity, in terms of node distance from the pollutant-source, flow rate, and a glyphosate-based herbicide, on phytoplankton productivity, diversity and stability. Gradostats were composed of interconnected equally spaced nodes where resources and phytoplankton move directionally along a gradient of increasing distance from the source of the polluting herbicide. We hypothesised that ecological effects would be stronger in the node situated closer to the point of herbicide input, but that flow would suppress phytoplankton populations in distant nodes. Overall, RoundUp impacted phytoplankton productivity and stability by reducing algal biomass and abundances. This occurred especially in the node closest to the diluted herbicide point-source and under high flow, where species abundances were heavily suppressed by the effects of the rapidly flowing herbicide. At low flow on the other hand, distant nodes where buffered from the effects of the slow-moving herbicide. No differences in beta and gamma diversity among replicate metaecosystems was found; however, a significant loss of alpha diversity in all metaecosystems occurred through time until the end of the experiment. Together, these results point to the importance of considering aquatic connectivity in management plans for monitoring and mitigating unintended ecological consequences of agrochemical runoff.

Multiannual patterns of genetic structure and mating type ratios highlight the complex bloom dynamics of a marine planktonic diatom

Understanding the genetic structure of populations and the processes responsible for its spatial and temporal dynamics is vital for assessing species' adaptability and survival in changing environments. We investigate the genetic fingerprinting of blooming populations of the marine diatom Pseudo-nitzschia multistriata in the Gulf of Naples (Mediterranean Sea) from 2008 to 2020. Strains were genotyped using microsatellite fingerprinting and natural samples were also analysed with Microsatellite Pool-seq Barcoding based on Illumina sequencing of microsatellite loci. Both approaches revealed a clonal expansion event in 2013 and a more stable genetic structure during 2017-2020 compared to previous years. The identification of a mating type (MT) determination gene allowed to assign MT to strains isolated over the years. MTs were generally at equilibrium with two notable exceptions, including the clonal bloom of 2013. The populations exhibited linkage equilibrium in most blooms, indicating that sexual reproduction leads to genetic homogenization. Our findings show that P. multistriata blooms exhibit a dynamic genetic and demographic composition over time, most probably determined by deeper-layer cell inocula. Occasional clonal expansions and MT imbalances can potentially affect the persistence and ecological success of planktonic diatoms.

Kalidium cuspidatum colonization changes the structure and function of salt crust microbial communities

The changes of soil moisture, salinity, and nutrients by halophyte colonization in high-salinity environment profoundly affect the assembly and structure of microbial communities. However, salt marshes in arid region have received little attention. This study was conducted in Lianhuachi Lake, a typical inland salt marsh wetland in China, to determine the physicochemical characteristics of salt crusts in [Kalidium cuspidatum (Ung.-Sternb.) Grub.] colonization areas and bulk soil, respectively, and to analyze the microbial community structure of salt crusts by high-throughput sequencing. Kalidium cuspidatum colonization significantly decreased total salinity, soil water content, and water-soluble ions of salt crusts and increased total carbon, total nitrogen, and total phosphorus content. At the same time, changes in physicochemical properties caused by Kalidium cuspidatum colonization affect the ecological processes of bacterial, fungal, and archaeal community assemblies in salt crusts. In addition, cross-kingdom network analysis showed that Kalidium cuspidatum colonization increased the complexity and stability of microbial networks in salt crust soils. Functional projections further showed that bacterial diversity had a potential driving effect on the nitrogen cycle function of salt crust. Our study further demonstrated the different ecological strategies of microorganisms for halophyte colonization in extreme environments and contributed to the understanding of restoration and management of salt marsh wetlands in arid region.

Geographic population structure and distinct intra-population dynamics of globally abundant freshwater bacteria

Implications of geographic separation and temporal dynamics on the evolution of free-living bacterial species are widely unclear. However, the vast amount of metagenome sequencing data generated during the last decades from various habitats around the world provides an unprecedented opportunity for such investigations. Here, we exploited publicly available and new freshwater metagenomes in combination with the genomes of abundant freshwater bacteria to reveal geographic and temporal population structure. We focused on species that were detected across broad geographic ranges at high enough sequence coverage for meaningful population genomic analyses, associated with the predominant freshwater taxa acI, LD12, Polynucleobacter, and Candidatus Methylopumilus. Despite the broad geographic ranges, each species appeared as a sequence-discrete cluster, in contrast to abundant marine taxa, for which continuous diversity structures were reported on a global scale. Population differentiation increased significantly with spatial distance in all species, but notable dispersal barriers (e.g. oceanic) were not apparent. Yet, the different species showed contrasting rates of geographic divergence and strikingly different intra-population dynamics in time series within individual habitats. The change in an LD12 population over 7 years was minor (FST = 0.04) compared to differentiation between lakes, whereas a Polynucleobacter population displayed strong changes within merely 2 months (FST up to 0.54), similar in scale to differentiation between populations separated by thousands of kilometers. The slowly and steadily evolving LD12 population showed high strain diversity, whereas the dynamic Polynucleobacter population exhibited alternating clonal expansions of mostly two strains only. Based on the contrasting population structures, we propose distinct models of speciation.

Spatiotemporal dynamics of giant viruses within a deep freshwater lake reveal a distinct dark-water community

Giant viruses (GVs) significantly regulate the ecological dynamics of diverse ecosystems. Although metagenomics has expanded our understanding of their diversity and ecological roles played in marine environments, little is known about GVs of freshwater ecosystems. Most previous studies have employed short-read sequencing and therefore resulted in fragmented genomes, hampering accurate assessment of genetic diversity. We sought to bridge this knowledge gap and overcome previous technical limitations. We subjected spatiotemporal (2 depths × 12 months) samples from Lake Biwa to metagenome-assembled genome reconstruction enhanced by long-read metagenomics. This yielded 293 GV metagenome-assembled genomes. Of these, 285 included previously unknown species in five orders of nucleocytoviruses and the first representatives of freshwater mirusviruses, which exhibited marked divergence from marine-derived lineages. The good performance of our long-read metagenomic assembly was demonstrated by the detection of 42 (14.3%) genomes composed of single contigs with completeness values >90%. GVs were partitioned across water depths, with most species specific to either the sunlit epilimnion or the dark hypolimnion. Epilimnion-specific members tended to be transient and exhibit short and intense abundance peaks, in line with the fact that they regulate the surface algal blooms. During the spring bloom, mirusviruses and members of three nucleocytovirus families were among the most abundant viruses. In contrast, hypolimnion-specific ones, including a mirusvirus genome, were typically more persistent in the hypolimnion throughout the water-stratified period, suggesting that they infect hosts specific to the hypolimnion and play previously unexplored ecological roles in dark water microbial ecosystems.

A genus in the bacterial phylum Aquificota appears to be endemic to Aotearoa-New Zealand

Allopatric speciation has been difficult to examine among microorganisms, with prior reports of endemism restricted to sub-genus level taxa. Previous microbial community analysis via 16S rRNA gene sequencing of 925 geothermal springs from the Taupō Volcanic Zone (TVZ), Aotearoa-New Zealand, revealed widespread distribution and abundance of a single bacterial genus across 686 of these ecosystems (pH 1.2-9.6 and 17.4-99.8 °C). Here, we present evidence to suggest that this genus, Venenivibrio (phylum Aquificota), is endemic to Aotearoa-New Zealand. A specific environmental niche that increases habitat isolation was identified, with maximal read abundance of Venenivibrio occurring at pH 4-6, 50-70 °C, and low oxidation-reduction potentials. This was further highlighted by genomic and culture-based analyses of the only characterised species for the genus, Venenivibrio stagnispumantis CP.B2T, which confirmed a chemolithoautotrophic metabolism dependent on hydrogen oxidation. While similarity between Venenivibrio populations illustrated that dispersal is not limited across the TVZ, extensive amplicon, metagenomic, and phylogenomic analyses of global microbial communities from DNA sequence databases indicates Venenivibrio is geographically restricted to the Aotearoa-New Zealand archipelago. We conclude that geographic isolation, complemented by physicochemical constraints, has resulted in the establishment of an endemic bacterial genus.

Balanced biogeographic and local environmental effects determine the patterns of microbial diversity in biocrusts at multi-scales

Introduction

Biodiversity maintenance and its underlying mechanisms are central issues of ecology. However, predicting the composition turnovers of microbial communities at multiple spatial scales remains greatly challenging because they are obscured by the inconsistent impacts of climatic and local edaphic conditions on the assembly process.

Methods

Based on the Illumina MeSeq 16S/18S rRNA sequencing technology, we investigated soil bacterial and eukaryotic communities in biocrusts with different successional levels at a subcontinental scale of Northern China.

Results

Results showed that irrespective of spatial scale, bacterial α diversity increased but eukaryotic diversity decreased with the primary succession, whereas both β diversities decreased at the subcontinental scale compared with smaller scales, indicating that the biogeographic pattern of soil microorganisms was balanced by successional convergence and distance decay effect. We found that the convergence of bacterial and eukaryotic communities was attributed to the turnovers of generalist and specialist species, respectively. In this process, edaphic and climatic factors showed unique roles in the changes of diversity at local/subcontinental scales. Moreover, the taxonomic diversity tended to be more susceptible to climatic and edaphic conditions, while biotic factors (photosynthesis and pigments) were more important to phylogenetic diversity.

Conclusion

Taken together, our study provided comprehensive insights into understanding the pattern of microbial diversity at multiple spatial scales of drylands.

Genus-Wide Transcriptional Landscapes Reveal Correlated Gene Networks Underlying Microevolutionary Divergence in Diatoms

Marine microbes like diatoms make up the base of marine food webs and drive global nutrient cycles. Despite their key roles in ecology, biogeochemistry, and biotechnology, we have limited empirical data on how forces other than adaptation may drive diatom diversification, especially in the absence of environmental change. One key feature of diatom populations is frequent extreme reductions in population size, which can occur both in situ and ex situ as part of bloom-and-bust growth dynamics. This can drive divergence between closely related lineages, even in the absence of environmental differences. Here, we combine experimental evolution and transcriptome landscapes (t-scapes) to reveal repeated evolutionary divergence within several species of diatoms in a constant environment. We show that most of the transcriptional divergence can be captured on a reduced set of axes, and that repeatable evolution can occur along a single major axis of variation defined by core ortholog expression comprising common metabolic pathways. Previous work has associated specific transcriptional changes in gene networks with environmental factors. Here, we find that these same gene networks diverge in the absence of environmental change, suggesting these pathways may be central in generating phenotypic diversity as a result of both selective and random evolutionary forces. If this is the case, these genes and the functions they encode may represent universal axes of variation. Such axes that capture suites of interacting transcriptional changes during diversification improve our understanding of both global patterns in local adaptation and microdiversity, as well as evolutionary forces shaping algal cultivation.

How does speciation in marine plankton work?

Marine plankton species are ecologically important, yet, it remains unclear how they originate in the ocean, where few barriers are apparent to cause the most common type of speciation - divergence in isolation. Here I discuss the use of modern evolutionary genetic approaches to shed light on longstanding questions regarding their evolution.

Gene-network analysis predicts clinical response to immunotherapy in patients affected by NSCLC

OBJECTIVES

Predictive biomarkers of response to immune checkpoint inhibitors (ICIs) have been extensively studied in non-small cell lung cancer (NSCLC) with controversial results. Recently, gene-network analysis emerged as a new tool to address tumor biology and behavior, representing a potential tool to evaluate response to therapies.

METHODS

Clinical data and genetic profiles of 644 advanced NSCLCs were retrieved from cBioPortal and the Cancer Genome Atlas (TCGA); 243 ICI-treated NSCLCs were used to identify an immunotherapy response signatures via mutated gene network analysis and K-means unsupervised clustering. Signatures predictive values were tested in an external dataset of 242 cases and assessed versus a control group of 159 NSCLCs treated with standard chemotherapy.

RESULTS

At least two mutations in the coding sequence of genes belonging to the chromatin remodelling pathway (A signature), and/or at least two mutations of genes involved in cell-to-cell signalling pathways (B signature), showed positive prediction in ICI-treated advanced NSCLC. Signatures performed best when combined for patients undergoing first-line immunotherapy, and for those receiving combined ICIs.

CONCLUSIONS

Alterations in genes related to chromatin remodelling complexes and cell-to-cell crosstalk may force dysfunctional immune evasion, explaining susceptibility to immunotherapy. Therefore, exploring mutated gene networks could be valuable for determining essential biological interactions, contributing to treatment personalization.

Assessing the relative importance of stressors to the benthic index, M-AMBI: An example from U.S. estuaries

M-AMBI, a multivariate benthic index, has been used by European and American (U.S.) authorities to assess estuarine and coastal health and has been used in scientific studies throughout the world. It has been shown to be related to multiple pressures and stressors, but the relative importance of individual stressors within a multiple stressor context has not generally been assessed. In this study, we assembled data collected between 1999 and 2015 by the U.S. Environmental Protection Agency using consistent methods. These data included sediment and water quality measures and benthic invertebrate data which were used to calculate M-AMBI. We further assembled watersheds for all US estuaries with benthic data and calculated land use metrics. Random forest (RF) was used to identify those variables most strongly related to M-AMBI. Because RF is a compilation of multiple, nonlinear models, we then assessed which of these variables had a direct relationship with M-AMBI. The resulting variables were then assessed using RF to identify the subsets of variables that produced an effective and parsimonious model. This process was conducted at the national and ecoregional scale and the variables identified as being most important to predict M-AMBI were compared with literature reports of ecological patterns in a given area. At the national scale, better condition was correlated with clearer waters, lower amounts of agriculture in the watershed, and lower carbon and metal concentrations in estuarine sediments. Other stressors were identified as being important at the ecoregional scale, although sediment metal concentrations and watershed agriculture were identified as being important in most ecoregions. Our results suggest that this technique is useful to identify the most important variables impacting M-AMBI at broad spatial scales, even when the percentage of sites in Bad or Poor condition is low. This technique also provides an initial identification of important stressors that can be used to target more intensive local studies.

Strong population genomic structure of the toxic dinoflagellate Alexandrium minutum inferred from meta-transcriptome samples

Despite theoretical expectations, marine microeukaryote population are often highly structured and the mechanisms behind such patterns remain to be elucidated. These organisms display huge census population sizes, yet genotyping usually requires clonal strains originating from single cells, hindering proper population sampling. Estimating allelic frequency directly from population wide samples, without any isolation step, offers an interesting alternative. Here, we validate the use of meta-transcriptome environmental samples to determine the population genetic structure of the dinoflagellate Alexandrium minutum. Strain and meta-transcriptome based results both indicated a strong genetic structure for A. minutum in Western Europe, to the level expected between cryptic species. The presence of numerous private alleles, and even fixed polymorphism, would indicate ancient divergence and absence of gene flow between populations. Single nucleotide polymorphisms (SNPs) displaying strong allele frequency differences were distributed throughout the genome, which might indicate pervasive selection from standing genetic variation (soft selective sweeps). However, a few genomic regions displayed extremely low diversity that could result from the fixation of adaptive de novo mutations (hard selective sweeps) within the populations.

Whole-genome scanning reveals environmental selection mechanisms that shape diversity in populations of the epipelagic diatom Chaetoceros

Diatoms form a diverse and abundant group of photosynthetic protists that are essential players in marine ecosystems. However, the microevolutionary structure of their populations remains poorly understood, particularly in polar regions. Exploring how closely related diatoms adapt to different environments is essential given their short generation times, which may allow rapid adaptations, and their prevalence in marine regions dramatically impacted by climate change, such as the Arctic and Southern Oceans. Here, we address genetic diversity patterns in Chaetoceros, the most abundant diatom genus and one of the most diverse, using 11 metagenome-assembled genomes (MAGs) reconstructed from Tara Oceans metagenomes. Genome-resolved metagenomics on these MAGs confirmed a prevalent distribution of Chaetoceros in the Arctic Ocean with lower dispersal in the Pacific and Southern Oceans as well as in the Mediterranean Sea. Single-nucleotide variants identified within the different MAG populations allowed us to draw a landscape of Chaetoceros genetic diversity and revealed an elevated genetic structure in some Arctic Ocean populations. Gene flow patterns of closely related Chaetoceros populations seemed to correlate with distinct abiotic factors rather than with geographic distance. We found clear positive selection of genes involved in nutrient availability responses, in particular for iron (e.g., ISIP2a, flavodoxin), silicate, and phosphate (e.g., polyamine synthase), that were further supported by analysis of Chaetoceros transcriptomes. Altogether, these results highlight the importance of environmental selection in shaping diatom diversity patterns and provide new insights into their metapopulation genomics through the integration of metagenomic and environmental data.

Thermal trait variation may buffer Southern Ocean phytoplankton from anthropogenic warming

Despite the potential of standing genetic variation to rescue communities and shape future adaptation to climate change, high levels of uncertainty are associated with intraspecific trait variation in marine phytoplankton. Recent model intercomparisons have pointed to an urgent need to reduce uncertainty in the projected responses of marine ecosystems to climate change, including Southern Ocean (SO) surface waters, which are among the most rapidly warming habitats on Earth. Because SO phytoplankton growth responses to warming sea surface temperature (SST) are poorly constrained, we developed a high-throughput growth assay to simultaneously examine inter- and intra-specific thermal trait variation in a group of 43 taxonomically diverse and biogeochemically important SO phytoplankton called diatoms. We found significant differential growth performance among species across thermal traits, including optimum and maximum tolerated growth temperatures. Within species, coefficients of variation ranged from 3% to 48% among strains for those same key thermal traits. Using SO SST projections for 2100, we predicted biogeographic ranges that differed by up to 97% between the least and most tolerant strains for each species, illustrating the role that strain-specific differences in temperature response can play in shaping predictions of future phytoplankton biogeography. Our findings revealed the presence and scale of thermal trait variation in SO phytoplankton and suggest these communities may already harbour the thermal trait diversity required to withstand projected 21st-century SST change in the SO even under severe climate forcing scenarios.

Reduced representation sequencing accurately quantifies relative abundance and reveals population-level variation in Pseudo-nitzschia spp

Certain species within the genus Pseudo-nitzschia are able to produce the neurotoxin domoic acid (DA), which can cause illness in humans, mass-mortality of marine animals, and closure of commercial and recreational shellfisheries during toxic events. Understanding and forecasting blooms of these harmful species is a primary management goal. However, accurately predicting the onset and severity of bloom events remains difficult, in part because the underlying drivers of bloom formation have not been fully resolved. Furthermore, Pseudo-nitzschia species often co-occur, and recent work suggests that the genetic composition of a Pseudo-nitzschia bloom may be a better predictor of toxicity than prevailing environmental conditions. We developed a novel next-generation sequencing assay using restriction site-associated DNA (2b-RAD) genotyping and applied it to mock Pseudo-nitzschia communities generated by mixing cultures of different species in known abundances. On average, 94% of the variance in observed species abundance was explained by the expected abundance. In addition, the false positive rate was low (0.45% on average) and unrelated to read depth, and false negatives were never observed. Application of this method to environmental DNA samples collected during natural Pseudo-nitzschia spp. bloom events in Southern California revealed that increases in DA were associated with increases in the relative abundance of P. australis. Although the absolute correlation across time-points was weak, an independent species fingerprinting assay (Automated Ribosomal Intergenic Spacer Analysis) supported this and identified other potentially toxic species. Finally, we assessed population-level genomic variation by mining SNPs from the environmental 2bRAD dataset. Consistent shifts in allele frequencies in P. pungens and P. subpacifica were detected between high and low DA years, suggesting that different intraspecific variants may be associated with prevailing environmental conditions or the presence of DA. Taken together, this method presents a potentially cost-effective and high-throughput approach for studies aiming to evaluate both population and species dynamics in mixed samples.

The microbiome of the dinoflagellate Prorocentrum cordatum in laboratory culture and its changes at higher temperatures

In the ocean, phytoplankton are dependent on communities of bacteria living in the phycosphere, a hot spot of metabolic and genetic exchange. Many types of interactions between phytoplankton and phycosphere bacteria have been shown, but it is unclear if the microbial communities associated with microalgae strains in culture collections are beneficial or harmful to the host strain. Here, we studied the microbial communities associated with four strains of the dinoflagellate Prorocentrum cordatum that had been isolated from distant geographical locations and maintained in culture collection for hundreds of generations. Community composition was determined by 16S rRNA gene amplicon sequencing. The dinoflagellate host strain was the strongest parameter separating communities, while growth phase, lifestyle (particle-attached versus free-living) and temperature had only a modulating effect. Although the strains had been isolated from distant locations in the Atlantic and Pacific Ocean, 14 ASVs were shared among all strains, the most abundant ones being Gilvibacter, Marivita, uncultivated Rhodobacteraceae, Marinobacter, Hyphomonadaceae, Cupriavidus, Variovorax, and Paucibacter. Adaptation to higher temperatures resulted in specific changes in each phycosphere microbiome, including increased abundance of rare community members. We then compared the growth of the four xenic cultures to that of the axenic P. cordatum CCMP1329. At 20°C, growth of the xenic cultures was similar or slower than that of CCMP1329. At 26°C, all four xenic cultures experienced a death phase, while the axenic culture stably remained in the stationary phase. At 30°C, only two of the xenic cultures were able to grow. A shift of dinoflagellate metabolism from autotrophy to mixotrophy and competition between dinoflagellate and bacteria for limiting nutrients, including essential vitamins, may contribute to these differences in growth patterns.

Phenological segregation suggests speciation by time in the planktonic diatom Pseudo-nitzschia allochrona sp. nov

The processes leading to the emergence of new species are poorly understood in marine plankton, where weak physical barriers and homogeneous environmental conditions limit spatial and ecological segregation. Here, we combine molecular and ecological information from a long-term time series and propose Pseudo-nitzschia allochrona, a new cryptic planktonic diatom, as a possible case of speciation by temporal segregation. The new species differs in several genetic markers (18S, 28S and ITS rDNA fragments and rbcL) from its closest relatives, which are morphologically very similar or identical, and is reproductively isolated from its sibling species P. arenysensis. Data from a long-term plankton time series show P. allochrona invariably occurring in summer-autumn in the Gulf of Naples, where its closely related species P. arenysensis, P. delicatissima, and P. dolorosa are instead found in winter-spring. Temperature and nutrients are the main factors associated with the occurrence of P. allochrona, which could have evolved in sympatry by switching its phenology and occupying a new ecological niche. This case of possible speciation by time shows the relevance of combining ecological time series with molecular information to shed light on the eco-evolutionary dynamics of marine microorganisms.

The rates of global bacterial and archaeal dispersal

The phylogenetic resolution at which microorganisms display geographic endemism, the rates at which they disperse at global scales, and the role of humans on global microbial dispersal are largely unknown. Answering these questions is necessary for interpreting microbial biogeography, ecology, and macroevolution and for predicting the spread of emerging pathogenic strains. To resolve these questions, I analyzed the geographic and evolutionary relationships between 36,795 bacterial and archaeal ("prokaryotic") genomes from ∼7000 locations around the world. I find clear signs of continental-scale endemism, including strong correlations between phylogenetic divergence and geographic distance. However, the phylogenetic scale at which endemism generally occurs is extremely small, and most "species" (defined by an average nucleotide identity ≥ 95%) and even closely related strains (average nucleotide identity ≥ 99.9%) are globally distributed. Human-associated lineages display faster dispersal rates than other terrestrial lineages; the average net distance between any two human-associated cell lineages diverging 50 years ago is roughly 580 km. These results suggest that many previously reported global-scale microbial biogeographical patterns are likely the result of recent or current environmental filtering rather than geographic endemism. For human-associated lineages, estimated transition rates between Europe and North America are particularly high, and much higher than for non-human associated terrestrial lineages, highlighting the role that human movement plays in global microbial dispersal. Dispersal was slowest for hot spring- and terrestrial subsurface-associated lineages, indicating that these environments may act as "isolated islands" of microbial evolution.

Host genotype structures the microbiome of a globally dispersed marine phytoplankton

Microscale interactions between marine phytoplankton and their bacterial microbiomes can influence ecosystem functioning and global biogeochemical cycling through complex exchanges of metabolites and sophisticated ecological processes. Previous investigation of the phytoplankton microbiome has not focused on the role of a host’s underlying genetic background. Through examination of a single phytoplankton species’ microbiome across the global ocean, we found that host genotype strongly influenced microbiome community composition, with associations that potentially persist across generations and ocean basins but assemble rapidly (within days). The long-term association of microbiomes with host genetic background could explain the evolution and maintenance of intricate phytoplankton–bacteria interactions.

Phytoplankton support complex bacterial microbiomes that rely on phytoplankton-derived extracellular compounds and perform functions necessary for algal growth. Recent work has revealed sophisticated interactions and exchanges of molecules between specific phytoplankton–bacteria pairs, but the role of host genotype in regulating those interactions is unknown. Here, we show how phytoplankton microbiomes are shaped by intraspecific genetic variation in the host using global environmental isolates of the model phytoplankton host Thalassiosira rotula and a laboratory common garden experiment. A set of 81 environmental T. rotula genotypes from three ocean basins and eight genetically distinct populations did not reveal a core microbiome. While no single bacterial phylotype was shared across all genotypes, we found strong genotypic influence of T. rotula, with microbiomes associating more strongly with host genetic population than with environmental factors. The microbiome association with host genetic population persisted across different ocean basins, suggesting that microbiomes may be associated with host populations for decades. To isolate the impact of host genotype on microbiomes, a common garden experiment using eight genotypes from three distinct host populations again found that host genotype influenced microbial community composition, suggesting that a process we describe as genotypic filtering, analogous to environmental filtering, shapes phytoplankton microbiomes. In both the environmental and laboratory studies, microbiome variation between genotypes suggests that other factors influenced microbiome composition but did not swamp the dominant signal of host genetic background. The long-term association of microbiomes with specific host genotypes reveals a possible mechanism explaining the evolution and maintenance of complex phytoplankton–bacteria chemical exchanges.

Marine phytoplankton functional types exhibit diverse responses to thermal change

Marine phytoplankton generate half of global primary production, making them essential to ecosystem functioning and biogeochemical cycling. Though phytoplankton are phylogenetically diverse, studies rarely designate unique thermal traits to different taxa, resulting in coarse representations of phytoplankton thermal responses. Here we assessed phytoplankton functional responses to temperature using empirically derived thermal growth rates from four principal contributors to marine productivity: diatoms, dinoflagellates, cyanobacteria, and coccolithophores. Using modeled sea surface temperatures for 1950-1970 and 2080-2100, we explored potential alterations to each group's growth rates and geographical distribution under a future climate change scenario. Contrary to the commonly applied Eppley formulation, our data suggest phytoplankton functional types may be characterized by different temperature coefficients (Q₁₀), growth maxima thermal dependencies, and thermal ranges which would drive dissimilar responses to each degree of temperature change. These differences, when applied in response to global simulations of future temperature, result in taxon-specific projections of growth and geographic distribution, with low-latitude coccolithophores facing considerable decreases and cyanobacteria substantial increases in growth rates. These results suggest that the singular effect of changing temperature may alter phytoplankton global community structure, owing to the significant variability in thermal response between phytoplankton functional types.

The biogeographic differentiation of algal microbiomes in the upper ocean from pole to pole

Eukaryotic phytoplankton are responsible for at least 20% of annual global carbon fixation. Their diversity and activity are shaped by interactions with prokaryotes as part of complex microbiomes. Although differences in their local species diversity have been estimated, we still have a limited understanding of environmental conditions responsible for compositional differences between local species communities on a large scale from pole to pole. Here, we show, based on pole-to-pole phytoplankton metatranscriptomes and microbial rDNA sequencing, that environmental differences between polar and non-polar upper oceans most strongly impact the large-scale spatial pattern of biodiversity and gene activity in algal microbiomes. The geographic differentiation of co-occurring microbes in algal microbiomes can be well explained by the latitudinal temperature gradient and associated break points in their beta diversity, with an average breakpoint at 14 °C ± 4.3, separating cold and warm upper oceans. As global warming impacts upper ocean temperatures, we project that break points of beta diversity move markedly pole-wards. Hence, abrupt regime shifts in algal microbiomes could be caused by anthropogenic climate change.

Genomic adaptations enabling Acidithiobacillus distribution across wide-ranging hot spring temperatures and pHs

BACKGROUND

Terrestrial hot spring settings span a broad spectrum of physicochemistries. Physicochemical parameters, such as pH and temperature, are key factors influencing differences in microbial composition across diverse geothermal areas. Nonetheless, analysis of hot spring pools from the Taupo Volcanic Zone (TVZ), New Zealand, revealed that some members of the bacterial genus, Acidithiobacillus, are prevalent across wide ranges of hot spring pHs and temperatures. To determine the genomic attributes of Acidithiobacillus that inhabit such diverse conditions, we assembled the genomes of 19 uncultivated hot spring Acidithiobacillus strains from six geothermal areas and compared these to 37 publicly available Acidithiobacillus genomes from various habitats.

RESULTS

Analysis of 16S rRNA gene amplicons from 138 samples revealed that Acidithiobacillus comprised on average 11.4 ± 16.8% of hot spring prokaryotic communities, with three Acidithiobacillus amplicon sequence variants (ASVs) (TVZ_G1, TVZ_G2, TVZ_G3) accounting for > 90% of Acidithiobacillus in terms of relative abundance, and occurring in 126 out of 138 samples across wide ranges of temperature (17.5-92.9 °C) and pH (1.0-7.5). We recovered 19 environmental genomes belonging to each of these three ASVs, as well as a fourth related group (TVZ_G4). Based on genome average nucleotide identities, the four groups (TVZ_G1-TVZ_G4) constitute distinct species (ANI < 96.5%) of which three are novel Acidithiobacillus species (TVZ_G2-TVZ_G4) and one belongs to Acidithiobacillus caldus (TVZ_G1). All four TVZ Acidithiobacillus groups were found in hot springs with temperatures above the previously known limit for the genus (up to 40 °C higher), likely due to significantly higher proline and GC contents than other Acidithiobacillus species, which are known to increase thermostability. Results also indicate hot spring-associated Acidithiobacillus have undergone genome streamlining, likely due to thermal adaptation. Moreover, our data suggest that Acidithiobacillus prevalence across varied hot spring pHs is supported by distinct strategies, whereby TVZ_G2-TVZ_G4 regulate pH homeostasis mostly through Na⁺/H⁺ antiporters and proton-efflux ATPases, whereas TVZ_G1 mainly relies on amino acid decarboxylases.

CONCLUSIONS

This study provides insights into the distribution of Acidithiobacillus species across diverse hot spring physichochemistries and determines genomic features and adaptations that potentially enable Acidithiobacillus species to colonize a broad range of temperatures and pHs in geothermal environments. Video Abstract.

Long-Term m5C Methylome Dynamics Parallel Phenotypic Adaptation in the Cyanobacterium Trichodesmium

A major challenge in modern biology is understanding how the effects of short-term biological responses influence long-term evolutionary adaptation, defined as a genetically determined increase in fitness to novel environments. This is particularly important in globally important microbes experiencing rapid global change, due to their influence on food webs, biogeochemical cycles, and climate. Epigenetic modifications like methylation have been demonstrated to influence short-term plastic responses, which ultimately impact long-term adaptive responses to environmental change. However, there remains a paucity of empirical research examining long-term methylation dynamics during environmental adaptation in nonmodel, ecologically important microbes. Here, we show the first empirical evidence in a marine prokaryote for long-term m5C methylome modifications correlated with phenotypic adaptation to CO2, using a 7-year evolution experiment (1,000+ generations) with the biogeochemically important marine cyanobacterium Trichodesmium. We identify m5C methylated sites that rapidly changed in response to high (750 µatm) CO2 exposure and were maintained for at least 4.5 years of CO2 selection. After 7 years of CO2 selection, however, m5C methylation levels that initially responded to high-CO2 returned to ancestral, ambient CO2 levels. Concurrently, high-CO2 adapted growth and N2 fixation rates remained significantly higher than those of ambient CO2 adapted cell lines irrespective of CO2 concentration, a trend consistent with genetic assimilation theory. These data demonstrate the maintenance of CO2-responsive m5C methylation for 4.5 years alongside phenotypic adaptation before returning to ancestral methylation levels. These observations in a globally distributed marine prokaryote provide critical evolutionary insights into biogeochemically important traits under global change.

Continental-Scale Gene Flow Prevents Allopatric Divergence of Pelagic Freshwater Bacteria

Allopatric divergence is one of the principal mechanisms for speciation of macro-organisms. Microbes by comparison are assumed to disperse more freely and to be less limited by dispersal barriers. However, thermophilic prokaryotes restricted to geothermal springs have shown clear signals of geographic isolation, but robust studies on this topic for microbes with less strict habitat requirements are scarce. Furthermore, it has only recently been recognized that homologous recombination among conspecific individuals provides species coherence in a wide range of prokaryotes. Recombination barriers thus may define prokaryotic species boundaries, yet, the extent to which geographic distance between populations gives rise to such barriers is an open question. Here, we investigated gene flow and population structure in a widespread species of pelagic freshwater bacteria, Polynucleobacter paneuropaeus. Through comparative genomics of 113 conspecific strains isolated from freshwater lakes and ponds located across a North–South range of more than 3,000 km, we were able to reconstruct past gene flow events. The species turned out to be highly recombinogenic as indicated by significant signs of gene transfer and extensive genome mosaicism. Although genomic differences increased with spatial distance on a regional scale (<170 km), such correlations were mostly absent on larger scales up to 3,400 km. We conclude that allopatric divergence in European P. paneuropaeus is minor, and that effective gene flow across the sampled geographic range in combination with a high recombination efficacy maintains species coherence.

Oceanic dispersal barriers in a holoplanktonic gastropod

Pteropods, a group of holoplanktonic gastropods, are regarded as bioindicators of the effects of ocean acidification on open ocean ecosystems, because their thin aragonitic shells are susceptible to dissolution. While there have been recent efforts to address their capacity for physiological acclimation, it is also important to gain predictive understanding of their ability to adapt to future ocean conditions. However, little is known about the levels of genetic variation and large-scale population structuring of pteropods, key characteristics enabling local adaptation. We examined the spatial distribution of genetic diversity in the mitochondrial cytochrome c oxidase I (COI) and nuclear 28S gene fragments, as well as shell shape variation, across a latitudinal transect in the Atlantic Ocean (35°N-36°S) for the pteropod Limacina bulimoides. We observed high levels of genetic variability (COI π = 0.034, 28S π = 0.0021) and strong spatial structuring (COI ΦST  = 0.230, 28S ΦST  = 0.255) across this transect. Based on the congruence of mitochondrial and nuclear differentiation, as well as differences in shell shape, we identified a primary dispersal barrier in the southern Atlantic subtropical gyre (15-18°S). This barrier is maintained despite the presence of expatriates, a gyral current system, and in the absence of any distinct oceanographic gradients in this region, suggesting that reproductive isolation between these populations must be strong. A secondary dispersal barrier supported only by 28S pairwise ΦST comparisons was identified in the equatorial upwelling region (between 15°N and 4°S), which is concordant with barriers observed in other zooplankton species. Both oceanic dispersal barriers were congruent with regions of low abundance reported for a similar basin-scale transect that was sampled 2 years later. Our finding supports the hypothesis that low abundance indicates areas of suboptimal habitat that result in barriers to gene flow in widely distributed zooplankton species. Such species may in fact consist of several populations or (sub)species that are adapted to local environmental conditions, limiting their potential for adaptive responses to ocean changes. Future analyses of genome-wide diversity in pteropods could provide further insight into the strength, formation and maintenance of oceanic dispersal barriers.

The Microbial Conveyor Belt: Connecting the Globe through Dispersion and Dormancy

Despite the recent increase in knowledge concerning microorganisms, the processes determining their global distribution and functioning have not been disentangled. Microbial dormant stages are adapted to endure specific adverse conditions related to their dispersion path, suggesting that dispersion is not entirely a stochastic process. Long-term dormancy enhances microbial dispersion, promoting the ubiquity of microorganisms. The evidence leads us to propose that there is a global, recurrent, and spatially cyclical dispersion of microorganisms that we have called the Microbial Conveyor Belt. These dispersion cycles directly influence the distribution of microorganisms, the global cycling of inorganic and organic matter, and thus the Earth system's functioning.

Phenotypic trait variability as an indication of adaptive capacity in a cosmopolitan marine diatom

Determining the adaptive capacity of marine phytoplankton is important in predicting changes in phytoplankton responses to ocean warming. Phytoplankton may consist of high levels of standing phenotypic and genetic variability, the basis of rapid evolution; however, few studies have quantified trait variability within and amongst closely related diatom species. Using 35 clonal cultures of the ubiquitous marine diatom Leptocylindrus isolated from six locations, spanning 2000 km of the south-eastern Australian coastline, we found evidence of significant intraspecific morphological and metabolic trait variability, which for 8 of 9 traits (growth rate, biovolume, C:N, silica deposition, silica incorporation rate, chl-a, and photosynthetic efficiency under dark adapted, growth irradiance, and high-light adaptation) were greater within a species than between species. Moreover, only two traits revealed a latitudinal trend with strains isolated from lower latitudes showing significantly higher silicification rates and protein:lipid content compared to their higher latitude counterparts. These data mirror recent studies on diatom intraspecific genetic diversity, which has found comparable levels of genetic diversity at a single site to those thousands of kilometres apart, and provide evidence of a functional role of diatom diversity that will allow for rapid adaptation via ecological selection on standing variation in response to changing conditions.

Geospatial immune variability illuminates differential evolution of lung adenocarcinoma

Remarkable progress in molecular analyses has improved our understanding of the evolution of cancer cells toward immune escape1-5. However, the spatial configurations of immune and stromal cells, which may shed light on the evolution of immune escape across tumor geographical locations, remain unaddressed. We integrated multiregion exome and RNA-sequencing (RNA-seq) data with spatial histology mapped by deep learning in 100 patients with non-small cell lung cancer from the TRACERx cohort6. Cancer subclones derived from immune cold regions were more closely related in mutation space, diversifying more recently than subclones from immune hot regions. In TRACERx and in an independent multisample cohort of 970 patients with lung adenocarcinoma, tumors with more than one immune cold region had a higher risk of relapse, independently of tumor size, stage and number of samples per patient. In lung adenocarcinoma, but not lung squamous cell carcinoma, geometrical irregularity and complexity of the cancer-stromal cell interface significantly increased in tumor regions without disruption of antigen presentation. Decreased lymphocyte accumulation in adjacent stroma was observed in tumors with low clonal neoantigen burden. Collectively, immune geospatial variability elucidates tumor ecological constraints that may shape the emergence of immune-evading subclones and aggressive clinical phenotypes.

Disentangling the mechanisms shaping the surface ocean microbiota

BACKGROUND

The ocean microbiota modulates global biogeochemical cycles and changes in its configuration may have large-scale consequences. Yet, the underlying ecological mechanisms structuring it are unclear. Here, we investigate how fundamental ecological mechanisms (selection, dispersal and ecological drift) shape the smallest members of the tropical and subtropical surface-ocean microbiota: prokaryotes and minute eukaryotes (picoeukaryotes). Furthermore, we investigate the agents exerting abiotic selection on this assemblage as well as the spatial patterns emerging from the action of ecological mechanisms. To explore this, we analysed the composition of surface-ocean prokaryotic and picoeukaryotic communities using DNA-sequence data (16S- and 18S-rRNA genes) collected during the circumglobal expeditions Malaspina-2010 and TARA-Oceans.

RESULTS

We found that the two main components of the tropical and subtropical surface-ocean microbiota, prokaryotes and picoeukaryotes, appear to be structured by different ecological mechanisms. Picoeukaryotic communities were predominantly structured by dispersal-limitation, while prokaryotic counterparts appeared to be shaped by the combined action of dispersal-limitation, selection and drift. Temperature-driven selection appeared as a major factor, out of a few selected factors, influencing species co-occurrence networks in prokaryotes but not in picoeukaryotes, indicating that association patterns may contribute to understand ocean microbiota structure and response to selection. Other measured abiotic variables seemed to have limited selective effects on community structure in the tropical and subtropical ocean. Picoeukaryotes displayed a higher spatial differentiation between communities and a higher distance decay when compared to prokaryotes, consistent with a scenario of higher dispersal limitation in the former after considering environmental heterogeneity. Lastly, random dynamics or drift seemed to have a more important role in structuring prokaryotic communities than picoeukaryotic counterparts.

CONCLUSIONS

The differential action of ecological mechanisms seems to cause contrasting biogeography, in the tropical and subtropical ocean, among the smallest surface plankton, prokaryotes and picoeukaryotes. This suggests that the idiosyncrasy of the main constituents of the ocean microbiota should be considered in order to understand its current and future configuration, which is especially relevant in a context of global change, where the reaction of surface ocean plankton to temperature increase is still unclear. Video Abstract.

Tracing the Origin of Planktonic Protists in an Ancient Lake

Ancient lakes are among the most interesting models for evolution studies because their biodiversity is the result of a complex combination of migration and speciation. Here, we investigate the origin of single celled planktonic eukaryotes from the oldest lake in the world—Lake Baikal (Russia). By using 18S rDNA metabarcoding, we recovered 1414 Operational Taxonomic Units (OTUs) belonging to protists populating surface waters (1–50 m) and representing pico/nano-sized cells. The recovered communities resembled other lacustrine freshwater assemblages found elsewhere, especially the taxonomically unclassified protists. However, our results suggest that a fraction of Baikal protists could belong to glacial relicts and have close relationships with marine/brackish species. Moreover, our results suggest that rapid radiation may have occurred among some protist taxa, partially mirroring what was already shown for multicellular organisms in Lake Baikal. We found 16% of the OTUs belonging to potential species flocks in Stramenopiles, Alveolata, Opisthokonta, Archaeplastida, Rhizaria, and Hacrobia. Putative flocks predominated in Chrysophytes, which are highly diverse in Lake Baikal. Also, the 18S rDNA of a number of species (7% of the total) differed >10% from other known sequences. These taxa as well as those belonging to the flocks may be endemic to Lake Baikal. Overall, our study points to novel diversity of planktonic protists in Lake Baikal, some of which may have emerged in situ after evolutionary diversification.

Microbial evolutionary strategies in a dynamic ocean

Marine microbes form the base of ocean food webs and drive ocean biogeochemical cycling. Yet little is known about the ability of microbial populations to adapt as they are advected through changing conditions. Here, we investigated the interplay between physical and biological timescales using a model of adaptation and an eddy-resolving ocean circulation climate model. Two criteria were identified that relate the timing and nature of adaptation to the ratio of physical to biological timescales. Genetic adaptation was impeded in highly variable regimes by nongenetic modifications but was promoted in more stable environments. An evolutionary trade-off emerged where greater short-term nongenetic transgenerational effects (low-γ strategy) enabled rapid responses to environmental fluctuations but delayed genetic adaptation, while fewer short-term transgenerational effects (high-γ strategy) allowed faster genetic adaptation but inhibited short-term responses. Our results demonstrate that the selective pressures for organisms within a single water mass vary based on differences in generation timescales resulting in different evolutionary strategies being favored. Organisms that experience more variable environments should favor a low-γ strategy. Furthermore, faster cell division rates should be a key factor in genetic adaptation in a changing ocean. Understanding and quantifying the relationship between evolutionary and physical timescales is critical for robust predictions of future microbial dynamics.

A genomics approach reveals the global genetic polymorphism, structure, and functional diversity of ten accessions of the marine model diatom Phaeodactylum tricornutum

Diatoms emerged in the Mesozoic period and presently constitute one of the main primary producers in the world's ocean and are of a major economic importance. In the current study, using whole genome sequencing of ten accessions of the model diatom Phaeodactylum tricornutum, sampled at broad geospatial and temporal scales, we draw a comprehensive landscape of the genomic diversity within the species. We describe strong genetic subdivisions of the accessions into four genetic clades (A-D) with constituent populations of each clade possessing a conserved genetic and functional makeup, likely a consequence of the limited dispersal of P. tricornutum in the open ocean. We further suggest dominance of asexual reproduction across all the populations, as implied by high linkage disequilibrium. Finally, we show limited yet compelling signatures of genetic and functional convergence inducing changes in the selection pressure on many genes and metabolic pathways. We propose these findings to have significant implications for understanding the genetic structure of diatom populations in nature and provide a framework to assess the genomic underpinnings of their ecological success and impact on aquatic ecosystems where they play a major role. Our work provides valuable resources for functional genomics and for exploiting the biotechnological potential of this model diatom species.

Single-amino acid variants reveal evolutionary processes that shape the biogeography of a global SAR11 subclade

Members of the SAR11 order Pelagibacterales dominate the surface oceans. Their extensive diversity challenges emerging operational boundaries defined for microbial 'species' and complicates efforts of population genetics to study their evolution. Here, we employed single-amino acid variants (SAAVs) to investigate ecological and evolutionary forces that maintain the genomic heterogeneity within ubiquitous SAR11 populations we accessed through metagenomic read recruitment using a single isolate genome. Integrating amino acid and protein biochemistry with metagenomics revealed that systematic purifying selection against deleterious variants governs non-synonymous variation among very closely related populations of SAR11. SAAVs partitioned metagenomes into two main groups matching large-scale oceanic current temperatures, and six finer proteotypes that connect distant oceanic regions. These findings suggest that environmentally-mediated selection plays a critical role in the journey of cosmopolitan surface ocean microbial populations, and the idea 'everything is everywhere but the environment selects' has credence even at the finest resolutions.

Coherence of Microcystis species revealed through population genomics

Microcystis is a genus of freshwater cyanobacteria, which causes harmful blooms in ecosystems worldwide. Some Microcystis strains produce harmful toxins such as microcystin, impacting drinking water quality. Microcystis colony morphology, rather than genetic similarity, is often used to classify Microcystis into morphospecies. Yet colony morphology is a plastic trait, which can change depending on environmental and laboratory culture conditions, and is thus an inadequate criterion for species delineation. Furthermore, Microcystis populations are thought to disperse globally and constitute a homogeneous gene pool. However, this assertion is based on relatively incomplete characterization of Microcystis genomic diversity. To better understand these issues, we performed a population genomic analysis of 33 newly sequenced genomes mainly from Canada and Brazil. We identified 17 Microcystis clusters of genomic similarity, five of which correspond to monophyletic clades containing at least three newly sequenced genomes. Four out of these five clades match to named morphospecies. Notably, M. aeruginosa is paraphyletic, distributed across 12 genomic clusters, suggesting it is not a coherent species. A few clades of closely related isolates are specific to a unique geographic location, suggesting biogeographic structure over relatively short evolutionary time scales. Higher homologous recombination rates within than between clades further suggest that monophyletic groups might adhere to a Biological Species-like concept, in which barriers to gene flow maintain species distinctness. However, certain genes-including some involved in microcystin and micropeptin biosynthesis-are recombined between monophyletic groups in the same geographic location, suggesting local adaptation.

High dispersal levels and lake warming are emergent drivers of cyanobacterial community assembly in peri-Alpine lakes

Disentangling the relative importance of deterministic and stochastic processes in shaping natural communities is central to ecology. Studies about community assembly over broad temporal and spatial scales in aquatic microorganisms are scarce. Here, we used 16S rDNA sequence data from lake sediments to test for community assembly patterns in cyanobacterial phylogenies across ten European peri-Alpine lakes and over a century of eutrophication and climate warming. We studied phylogenetic similarity in cyanobacterial assemblages over spatial and temporal distance, and over environmental gradients, comparing detected patterns with theoretical expectations from deterministic and stochastic processes. We found limited evidence for deviation of lake communities from a random assembly model and no significant effects of geographic distance on phylogenetic similarity, suggesting no dispersal limitation and high levels of stochastic assembly. We detected a weak influence of phosphorus, but no significant effect of nitrogen levels on deviation of community phylogenies from random. We found however a significant decay of phylogenetic similarity for non-random communities over a gradient of air temperature and water column stability. We show how phylogenetic data from sedimentary archives can improve our understanding of microbial community assembly processes, and support previous evidence that climate warming has been the strongest environmental driver of cyanobacterial community assembly over the past century.

Global pattern of phytoplankton diversity driven by temperature and environmental variability

Despite their importance to ocean productivity, global patterns of marine phytoplankton diversity remain poorly characterized. Although temperature is considered a key driver of general marine biodiversity, its specific role in phytoplankton diversity has remained unclear. We determined monthly phytoplankton species richness by using niche modeling and >540,000 global phytoplankton observations to predict biogeographic patterns of 536 phytoplankton species. Consistent with metabolic theory, phytoplankton richness in the tropics is about three times that in higher latitudes, with temperature being the most important driver. However, below 19°C, richness is lower than expected, with ~8°- 14°C waters (~35° to 60° latitude) showing the greatest divergence from theoretical predictions. Regions of reduced richness are characterized by maximal species turnover and environmental variability, suggesting that the latter reduces species richness directly, or through enhancing competitive exclusion. The nonmonotonic relationship between phytoplankton richness and temperature suggests unanticipated complexity in responses of marine biodiversity to ocean warming.

Shared genomic outliers across two divergent population clusters of a highly threatened seagrass

The seagrass, Zostera capensis, occurs across a broad stretch of coastline and wide environmental gradients in estuaries and sheltered bays in southern and eastern Africa. Throughout its distribution, habitats are highly threatened and poorly protected, increasing the urgency of assessing the genomic variability of this keystone species. A pooled genomic approach was employed to obtain SNP data and examine neutral genomic variation and to identify potential outlier loci to assess differentiation across 12 populations across the ∼9,600 km distribution of Z. capensis. Results indicate high clonality and low genomic diversity within meadows, which combined with poor protection throughout its range, increases the vulnerability of this seagrass to further declines or local extinction. Shared variation at outlier loci potentially indicates local adaptation to temperature and precipitation gradients, with Isolation-by-Environment significantly contributing towards shaping spatial variation in Z. capensis. Our results indicate the presence of two population clusters, broadly corresponding to populations on the west and east coasts, with the two lineages shaped only by frequency differences of outlier loci. Notably, ensemble modelling of suitable seagrass habitat provides evidence that the clusters are linked to historical climate refugia around the Last Glacial Maxi-mum. Our work suggests a complex evolutionary history of Z. capensis in southern and eastern Africa that will require more effective protection in order to safeguard this important ecosystem engineer into the future.

Community structural differences shape microbial responses to high molecular weight organic matter

The extent to which differences in microbial community structure result in variations in organic matter (OM) degradation is not well understood. Here, we tested the hypothesis that distinct marine microbial communities from North Atlantic surface and bottom waters would exhibit varying compositional succession and functional shifts in response to the same pool of complex high molecular weight (HMW-OM). We also hypothesized that microbial communities would produce a broader spectrum of enzymes upon exposure to HMW-OM, indicating a greater potential to degrade these compounds than reflected by initial enzymatic activities. Our results show that community succession in amended mesocosms was congruent with cell growth, increased bacterial production and most notably, with substantial shifts in enzymatic activities. In all amended mesocosms, closely related taxa that were initially rare became dominant at time frames during which a broader spectrum of active enzymes were detected compared to initial timepoints, indicating a similar response among different communities. However, succession on the whole-community level, and the rates, spectra and progression of enzymatic activities, reveal robust differences among distinct communities from discrete water masses. These results underscore the crucial role of rare bacterial taxa in ocean carbon cycling and the importance of bacterial community structure for HMW-OM degradation.

The role of intraspecific variation in the ecological and evolutionary success of diatoms in changing environments

Intraspecific variation in diatoms has been shown to play a key role in species' responses to several important environmental factors such as light, salinity, temperature and nutrients. Furthermore, modelling efforts indicate that this variation within species extends bloom periods, and likely provides sufficient variability in competitive interactions between species under hydrographically variable conditions. The intraspecific variation most likely corresponds to optimal fitness in temporary microhabitats and may help to explain the paradox of the plankton. Here, we examine the implications of intraspecific variation for the ecology and success of diatoms in general and emphasize the potential implications for our understanding of carbon metabolism in these important organisms. Additionally, data from palaeoecological studies have the potential for evaluating genetic variation through past climate changes, going thousands of years back in time. We suggest pathways for future research including the adoption of multiple strains of individual species into studies of diatom carbon metabolism, to refine our understanding of the variation within and between species, and the inclusion of experimental evolution as a tool for understanding potential evolutionary responses of diatom carbon metabolism to climate change.This article is part of the themed issue 'The peculiar carbon metabolism in diatoms'.

Morphological diversity and phylogeny of the diatom genus Entomoneis (Bacillariophyta) in marine plankton: six new species from the Adriatic Sea

The diatom genus Entomoneis is known from the benthos and plankton of marine, brackish, and freshwaters. Entomoneis includes diatoms with a bilobate keel elevated above the valve surface, a sigmoid canal raphe, and numerous girdle bands. Owing mostly to the scarcity of molecular data for a diverse set of species, the phylogeny of Entomoneis has not been investigated in depth. The few previous studies that included Entomoneis were focused on broader questions and the available data were from a small number of either unidentified Entomoneis or well-known species (e.g., E. paludosa). Since the first description of new species combining both molecular and morphological characters (E. tenera), we have continued to cultivate and investigate Entomoneis in the plankton of the Adriatic Sea. Combined multigene phylogeny (SSU rDNA sequences, rbcL, and psbC genes) and morphological observations (LM, SEM and TEM) revealed six new Entomoneis species supported by phylogenetic and morphological data: E. pusilla, E. gracilis, E. vilicicii, E. infula, E. adriatica, and E. umbratica. The most important morphological features for species delineation were cell shape, the degree and mode of torsion, valve apices, the appearance and structure of the transition between keel and valve body, the ultrastructure and the shape of the girdle bands, and the arrangement and density of perforations along the valve and valvocopulae. Our results highlight the underappreciated diversity of Entomoneis and call for a more in-depth morphological and molecular investigation of this genus especially in planktonic habitats.

The 'Neglected' Soil Virome - Potential Role and Impact

Bacteriophages are among the most abundant and diverse biological units in the biosphere. They have contributed to our understanding of the central dogma of biology and have been instrumental in the evolutionary success of bacterial pathogens. In contrast to our current understanding of marine viral communities, the soil virome and its function in terrestrial ecosystems has remained relatively understudied. Here, we examine, in a comparative fashion, the knowledge gathered from studies performed in soil versus marine settings. We address the information with respect to the abundance, diversity, ecological significance, and effects of, in particular, bacteriophages on their host's evolutionary trajectories. We also identify the main challenges that soil virology faces and the studies that are required to accompany the current developments in marine settings.

Cryptic diversity within the harmful dinoflagellate Akashiwo sanguinea in coastal Chinese waters is related to differentiated ecological niches

Blooms of the harmful dinoflagellate Akashiwo sanguinea are responsible for the mass mortality of fish and invertebrates in coastal waters. This cosmopolitan species includes several genetically differentiated clades. Four clonal cultures were established by isolating single cells from Xiamen Harbour (the East China Sea) for morphological and genetic analyses. The cultures displayed identical morphology but were genetically different, thus revealing presence of cryptic diversity in the study area. New details of the apical structure complex of Akashiwo sanguinea were also found. To investigate whether the observed cryptic diversity was related to environmental differentiation, 634 cells were obtained from seasonal water samples collected from 2008 to 2012. These cells were sequenced by single-cell PCR. For comparison with Chinese material, additional large subunit ribosomal DNA sequences were obtained for three established strains from Malaysian and French waters. To examine potential ecological differentiation of the distinct genotypes, growth responses of the studied strains were tested under laboratory conditions at temperatures of 12°C to 33°C. These experiments showed four distinct ribotypes of A. sanguinea globally, with the ribotypes A and B co-occuring in Xiamen Harbour. Ribotype A of A. sanguinea was present year-round in Xiamen Harbour, but it only bloomed in the winter and spring, thus corresponding to the winter type. In contrast, A. sanguinea ribotype B bloomed only in the summer, corresponding to the summer type. This differentiation supports the temperature optimum conditions that were established for these two ribotypes in the laboratory. Ribotype A grew better at lower temperatures compared to ribotype B which preferred higher temperatures. These findings support the idea that various ribotypes of A. sanguinea correspond to distinct ecotypes and allopatric speciation occurred in different climatic regions followed by dispersal.