Seasonal and interannual variability of the free-living and particle-associated bacteria of a coastal microbiome.

Marine microbial communities differ genetically, metabolically, and ecologically according to their lifestyle, and they may respond differently to environmental changes. In this study, we investigated the seasonal dynamics of bacterial assemblies in the free-living (FL) and particle-associated (PA) fractions across a span of 6 years in the Blanes Bay Microbial Observatory in the Northwestern Mediterranean. Both lifestyles showed marked seasonality. The trends in alpha diversity were similar, with lower values in spring-summer than in autumn-winter. Samples from both fractions were grouped seasonally and the percentage of community variability explained by the measured environmental variables was comparable (32% in FL and 31% in PA). Canonical analyses showed that biotic interactions were determinants of bacterioplankton dynamics and that their relevance varies depending on lifestyles. Time-decay curves confirmed a high degree of predictability in both fractions. Yet, 'seasonal' Amplicon Sequence Variants (ASVs) (as defined by Lomb Scargle time series analysis) in the PA communities represented 46% of the total relative abundance while these accounted for 30% in the FL fraction. These results demonstrate that bacteria inhabiting both fractions exhibit marked seasonality, highlighting the importance of accounting for both lifestyles to fully comprehend the dynamics of marine prokaryotic communities.

Impacts of a newly formed lava delta on the marine environment: Lava-induced upwelling and abrupt chlorophyll depletion during the Tajogaite eruption (La Palma, 2021).

The 2021 Tajogaite eruption in La Palma (Canary Islands, Spain) emitted vast volumes of lava during 85 days, which reached the ocean in several occasions at the western flank of the island. Most of these flows merged to create a primary lava delta, covering an area of 48 ha, with an additional 30 ha underwater. Here we characterize the effects of the lava-seawater interaction on the surrounding marine environment. The area was sampled during two multidisciplinary oceanographic cruises: the first one comprised the days before the lava reached the ocean and after the first contact; and the second took place a month later, when the lava delta was already formed but still receiving lava inputs. Physical-chemical anomalies were found in the whole water column at different depths up to 300 m in all measured parameters, such as turbidity (+9 NTU), dissolved oxygen concentration (-17.17 µmol kg-1), pHT25 (-0.1), and chlorophyll-a concentration (-0.33 mg m-3). Surface temperature increased up to +2.3 °C (28.5 °C) and surface salinity showed increases and decreases of -1.01 and +0.70, respectively, in a radius of 4 km around the lava delta. In the water column, the heated waters experimented a lava-induced upwelling, bringing deeper, nutrient-rich waters to shallower depths; however, this feature did not trigger any phytoplankton bloom. In fact, integrated chlorophyll-a showed an abrupt decrease of -41 % in just two days and -69 % a month later, compared to prior conditions. The chlorophyll-a depletion reached a distance larger than 2.5 km (not delimited).

Growth rates of marine prokaryotes are extremely diverse, even among closely related taxa.

Marine prokaryotes play crucial roles in ocean biogeochemical cycles, being their contribution strongly influenced by their growth rates. Hence, elucidating the variability and phylogenetic imprint of marine prokaryotes' growth rates are crucial for better determining the role of individual taxa in biogeochemical cycles. Here, we estimated prokaryotic growth rates at high phylogenetic resolution in manipulation experiments using water from the northwestern Mediterranean Sea. Experiments were run in the four seasons with different treatments that reduced growth limiting factors: predators, nutrient availability, viruses, and light. Single-amplicon sequence variants (ASVs)-based growth rates were calculated from changes in estimated absolute abundances using total prokaryotic abundance and the proportion of each individual ASV. The trends obtained for growth rates in the different experiments were consistent with other estimates based on total cell-counts, catalyzed reporter deposition fluorescence in situ hybridization subcommunity cell-counts or metagenomic-operational taxonomic units (OTUs). Our calculations unveil a broad range of growth rates (0.3-10 d-1) with significant variability even within closely related ASVs. Likewise, the impact of growth limiting factors changed over the year for individual ASVs. High numbers of responsive ASVs were shared between winter and spring seasons, as well as throughout the year in the treatments with reduced nutrient limitation and viral pressure. The most responsive ASVs were rare in the in situ communities, comprising a large pool of taxa with the potential to rapidly respond to environmental changes. Essentially, our results highlight the lack of phylogenetic coherence in the range of growth rates observed, and differential responses to the various limiting factors, even for closely related taxa.

Growth and mortality of aerobic anoxygenic phototrophs in the North Pacific Subtropical Gyre.

Aerobic anoxygenic phototrophic (AAP) bacteria harvest light energy using bacteriochlorophyll-containing reaction centers to supplement their mostly heterotrophic metabolism. While their abundance and growth have been intensively studied in coastal environments, much less is known about their activity in oligotrophic open ocean regions. Therefore, we combined in situ sampling in the North Pacific Subtropical Gyre, north of O'ahu island, Hawaii, with two manipulation experiments. Infra-red epifluorescence microscopy documented that AAP bacteria represented approximately 2% of total bacteria in the euphotic zone with the maximum abundance in the upper 50 m. They conducted active photosynthetic electron transport with maximum rates up to 50 electrons per reaction center per second. The in situ decline of bacteriochlorophyll concentration over the daylight period, an estimate of loss rates due to predation, indicated that the AAP bacteria in the upper 50 m of the water column turned over at rates of 0.75-0.90 d-1. This corresponded well with the specific growth rate determined in dilution experiments where AAP bacteria grew at a rate 1.05 ± 0.09 d-1. An amendment of inorganic nitrogen to obtain N:P = 32 resulted in a more than 10 times increase in AAP abundance over 6 days. The presented data document that AAP bacteria are an active part of the bacterioplankton community in the oligotrophic North Pacific Subtropical Gyre and that their growth was mostly controlled by nitrogen availability and grazing pressure.IMPORTANCEMarine bacteria represent a complex assembly of species with different physiology, metabolism, and substrate preferences. We focus on a specific functional group of marine bacteria called aerobic anoxygenic phototrophs. These photoheterotrophic organisms require organic carbon substrates for growth, but they can also supplement their metabolic needs with light energy captured by bacteriochlorophyll. These bacteria have been intensively studied in coastal regions, but rather less is known about their distribution, growth, and mortality in the oligotrophic open ocean. Therefore, we conducted a suite of measurements in the North Pacific Subtropical Gyre to determine the distribution of these organisms in the water column and their growth and mortality rates. A nutrient amendment experiment showed that aerobic anoxygenic phototrophs were limited by inorganic nitrogen. Despite this, they grew more rapidly than average heterotrophic bacteria, but their growth was balanced by intense grazing pressure.

Expanding success in the isolation of abundant marine bacteria after reduction in grazing and viral pressure and increase in nutrient availability.

Isolation of microorganisms is a useful approach to gathering knowledge about their genomic properties, physiology, and ecology, in addition to allowing the characterization of novel taxa. We performed an extensive isolation effort on samples from seawater manipulation experiments that were carried out during the four astronomical seasons in a coastal site of the northwest Mediterranean to evaluate the impact of grazing, viral mortality, resource competition reduction, and light presence/absence on bacterioplankton growth. Isolates were retrieved using two growth media, and their full 16S rRNA gene was sequenced to assess their identity and calculate their culturability across seasons and experimental conditions. A total of 1,643 isolates were obtained, mainly affiliated to the classes Gammaproteobacteria (44%), Alphaproteobacteria (26%), and Bacteroidia (17%). Isolates pertaining to class Gammaproteobacteria were the most abundant in all experiments, while Bacteroidia were preferentially enriched in the treatments with reduced grazing. Sixty-one isolates had a similarity below 97% to cultured taxa and are thus putatively novel. Comparison of isolate sequences with 16S rRNA gene amplicon sequences from the same samples showed that the percentage of reads corresponding to isolates was 21.4% within the whole data set, with dramatic increases in the summer virus-reduced (71%) and diluted (47%) treatments. In fact, we were able to isolate the top 10 abundant taxa in several experiments and from the whole data set. We also show that top-down and bottom-up controls differentially affect taxa in terms of culturability. Our results indicate that culturing marine bacteria using agar plates can be successful in certain ecological situations. IMPORTANCE Bottom-up and top-down controls greatly influence marine microbial community composition and dynamics, which in turn have effects on their culturability. We isolated a high amount of heterotrophic bacterial strains from experiments where seawater environmental conditions had been manipulated and found that decreasing grazing and viral pressure as well as rising nutrient availability are key factors increasing the success in culturing marine bacteria. Our data hint at factors influencing culturability and underpin bacterial cultures as a powerful way to discover new taxa.

A Metagenomic and Amplicon Sequencing Combined Approach Reveals the Best Primers to Study Marine Aerobic Anoxygenic Phototrophs.

Studies based on protein-coding genes are essential to describe the diversity within bacterial functional groups. In the case of aerobic anoxygenic phototrophic (AAP) bacteria, the pufM gene has been established as the genetic marker for this particular functional group, although available primers are known to have amplification biases. We review here the existing primers for pufM gene amplification, design new ones, and evaluate their phylogenetic coverage. We then use samples from contrasting marine environments to evaluate their performance. By comparing the taxonomic composition of communities retrieved with metagenomics and with different amplicon approaches, we show that the commonly used PCR primers are biased towards the Gammaproteobacteria phylum and some Alphaproteobacteria clades. The metagenomic approach, as well as the use of other combinations of the existing and newly designed primers, show that these groups are in fact less abundant than previously observed, and that a great proportion of pufM sequences are affiliated to uncultured representatives, particularly in the open ocean. Altogether, the framework developed here becomes a better alternative for future studies based on the pufM gene and, additionally, serves as a reference for primer evaluation of other functional genes.

Phytoplankton dynamics in the Mar Menor, a Mediterranean coastal lagoon strongly impacted by eutrophication.

The Mar Menor hypersaline coastal lagoon has suffered serious degradation in the last three decades attributable to nutrient pollution. In 2015, the lagoon experienced an intensive bloom of cyanobacteria that triggered a drastic change of its ecosystem. Our analyses indicate that phytoplankton in 2016-2021 did not present a seasonal variability pattern; the community was mainly dominated by diatoms and punctually reached abundance peaks above 107 cell L-1 along with chlorophyll a concentrations exceeding 20 µg L-1. The predominant diatom genera during these blooms were different as well as the nutrient conditions under which they were produced. These high diatom abundances are unprecedented in the lagoon; in fact, our data indicate that the taxonomic composition, time variation patterns and cell abundance of phytoplankton in 2016-2021 differ notably in comparison to the data published before 2015. Consequently, our results support the finding that the trophic status of the lagoon has changed profoundly.

Seasonality of biogeochemically relevant microbial genes in a coastal ocean microbiome.

Microbes drive the biogeochemical cycles of marine ecosystems through their vast metabolic diversity. While we have a fairly good understanding of the spatial distribution of these metabolic processes in various ecosystems, less is known about their seasonal dynamics. We investigated the annual patterns of 21 biogeochemical relevant functions in an oligotrophic coastal ocean site by analysing the presence of key genes, analysing high-rank gene taxonomy and the dynamics of nucleotide variants. Most genes presented seasonality: photoheterotrophic processes were enriched during spring, phosphorous-related genes were dominant during summer, coinciding with potential phosphate limitation, and assimilatory nitrate reductases appeared mostly during summer and autumn, correlating negatively with nitrate availability. Additionally, we identified the main taxa driving each function at each season and described the role of underrecognized taxa such as Litoricolaceae in carbon fixation (rbcL), urea degradation (ureC), and CO oxidation (coxL). Finally, the seasonality of single variants of some families presented a decoupling between the taxonomic abundance patterns and the functional gene patterns, implying functional specialization of the different genera. Our study unveils the seasonality of key biogeochemical functions and the main taxonomic groups that harbour these relevant functions in a coastal ocean ecosystem.

Long-term patterns of an interconnected core marine microbiota.

BACKGROUND: Ocean microbes constitute ~ 70% of the marine biomass, are responsible for ~ 50% of the Earth's primary production and are crucial for global biogeochemical cycles. Marine microbiotas include core taxa that are usually key for ecosystem function. Despite their importance, core marine microbes are relatively unknown, which reflects the lack of consensus on how to identify them. So far, most core microbiotas have been defined based on species occurrence and abundance. Yet, species interactions are also important to identify core microbes, as communities include interacting species. Here, we investigate interconnected bacteria and small protists of the core pelagic microbiota populating a long-term marine-coastal observatory in the Mediterranean Sea over a decade. RESULTS: Core microbes were defined as those present in > 30% of the monthly samples over 10 years, with the strongest associations. The core microbiota included 259 Operational Taxonomic Units (OTUs) including 182 bacteria, 77 protists, and 1411 strong and mostly positive (~ 95%) associations. Core bacteria tended to be associated with other bacteria, while core protists tended to be associated with bacteria. The richness and abundance of core OTUs varied annually, decreasing in stratified warmers waters and increasing in colder mixed waters. Most core OTUs had a preference for one season, mostly winter, which featured subnetworks with the highest connectivity. Groups of highly associated taxa tended to include protists and bacteria with predominance in the same season, particularly winter. A group of 13 highly-connected hub-OTUs, with potentially important ecological roles dominated in winter and spring. Similarly, 18 connector OTUs with a low degree but high centrality were mostly associated with summer or autumn and may represent transitions between seasonal communities. CONCLUSIONS: We found a relatively small and dynamic interconnected core microbiota in a model temperate marine-coastal site, with potential interactions being more deterministic in winter than in other seasons. These core microbes would be essential for the functioning of this ecosystem over the year. Other non-core taxa may also carry out important functions but would be redundant and non-essential. Our work contributes to the understanding of the dynamics and potential interactions of core microbes possibly sustaining ocean ecosystem function.

Global diversity and distribution of aerobic anoxygenic phototrophs in the tropical and subtropical oceans.

The aerobic anoxygenic phototrophic (AAP) bacteria are common in most marine environments but their global diversity and biogeography remain poorly characterized. Here, we analyzed AAP communities across 113 globally-distributed surface ocean stations sampled during the Malaspina Expedition in the tropical and subtropical ocean. By means of amplicon sequencing of the pufM gene, a genetic marker for this functional group, we show that AAP communities along the surface ocean were mainly composed of members of the Halieaceae (Gammaproteobacteria), which were adapted to a large range of environmental conditions, and of different clades of the Alphaproteobacteria, which seemed to dominate under particular circumstances, such as in the oligotrophic gyres. AAP taxa were spatially structured within each of the studied oceans, with communities from adjacent stations sharing more taxonomic similarities. AAP communities were composed of a large pool of rare members and several habitat specialists. When compared to the surface ocean prokaryotic and picoeukaryotic communities, it appears that AAP communities display an idiosyncratic global biogeographical pattern, dominated by selection processes and less influenced by dispersal limitation. Our study contributes to the understanding of how AAP communities are distributed in the horizontal dimension and the mechanisms underlying their distribution across the global surface ocean.

Seasonal niche differentiation among closely related marine bacteria.

Bacteria display dynamic abundance fluctuations over time in marine environments, where they play key biogeochemical roles. Here, we characterized the seasonal dynamics of marine bacteria in a coastal oligotrophic time series station, tested how similar the temporal niche of closely related taxa is, and what are the environmental parameters modulating their seasonal abundance patterns. We further explored how conserved the niche is at higher taxonomic levels. The community presented recurrent patterns of seasonality for 297 out of 6825 amplicon sequence variants (ASVs), which constituted almost half of the total relative abundance (47%). For certain genera, niche similarity decreased as nucleotide divergence in the 16S rRNA gene increased, a pattern compatible with the selection of similar taxa through environmental filtering. Additionally, we observed evidence of seasonal differentiation within various genera as seen by the distinct seasonal patterns of closely related taxa. At broader taxonomic levels, coherent seasonal trends did not exist at the class level, while the order and family ranks depended on the patterns that existed at the genus level. This study identifies the coexistence of closely related taxa for some bacterial groups and seasonal differentiation for others in a coastal marine environment subjected to a strong seasonality.

Role of small-sized phytoplankton in triggering an ecosystem disruptive algal bloom in a Mediterranean hypersaline coastal lagoon.

Monthly samplings carried out in 2016-2019 and satellite color images from 2002 to 2019 have been combined to determine the onset and causative species of the ecosystem disruptive algal bloom (EDAB) that affects the Mar Menor coastal lagoon (Western Mediterranean Sea) since 2015. Substantial changes in satellite spectral reflectance attributable to increasing abundance of Synechococcus were registered in 2014. Furthermore, cell abundances of this species in 2016 were the largest ever obtained in the lagoon (6 106 cells mL-1), with values similar to those reported for other Mediterranean hypertrophic estuaries and coastal lagoons. These results suggest that the early changes leading to the EDAB started in 2014 and that Synechococcus played a relevant role in its development. Moreover, diatom and dinoflagellate abundances changed substantially in 2016-2019, ranging from 102 to more than 104 cells mL-1. Some of these changes were linked to flood, suggesting that EDAB has modified substantially the homeostatic capacity of the lagoon.

Seasonal impact of grazing, viral mortality, resource availability and light on the group-specific growth rates of coastal Mediterranean bacterioplankton.

Estimation of prokaryotic growth rates is critical to understand the ecological role and contribution of different microbes to marine biogeochemical cycles. However, there is a general lack of knowledge on what factors control the growth rates of different prokaryotic groups and how these vary between sites and along seasons at a given site. We carried out several manipulation experiments during the four astronomical seasons in the coastal NW Mediterranean in order to evaluate the impact of grazing, viral mortality, resource competition and light on the growth and loss rates of prokaryotes. Gross and net growth rates of different bacterioplankton groups targeted by group-specific CARD-FISH probes and infrared microscopy (for aerobic anoxygenic phototrophs, AAP), were calculated from changes in cell abundances. Maximal group-specific growth rates were achieved when both predation pressure and nutrient limitation were experimentally minimized, while only a minimal effect of viral pressure on growth rates was observed; nevertheless, the response to predation removal was more remarkable in winter, when the bacterial community was not subjected to nutrient limitation. Although all groups showed increases in their growth rates when resource competition as well as grazers and viral pressure were reduced, Alteromonadaceae consistently presented the highest rates in all seasons. The response to light availability was generally weaker than that to the other factors, but it was variable between seasons. In summer and spring, the growth rates of AAP were stimulated by light whereas the growth of the SAR11 clade (likely containing proteorhodopsin) was enhanced by light in all seasons. Overall, our results set thresholds on bacterioplankton group-specific growth and mortality rates and contribute to estimate the seasonally changing contribution of various bacterioplankton groups to the function of microbial communities. Our results also indicate that the least abundant groups display the highest growth rates, contributing to the recycling of organic matter to a much greater extent than what their abundances alone would predict.

Assessment of microbial plankton diversity as an ecological indicator in the NW Mediterranean coast.

High-throughput sequencing of microbial assemblages has been proposed as an alternative methodology to the traditional ones used in marine monitoring and environmental assessment. Here, we evaluated pico- and nanoplankton diversity as ecological indicators in NW Mediterranean coastal waters by comparing their diversity in samples subjected to varying degrees of continental pressures. Using metabarcoding of the 16S and 18S rRNA genes, we explored whether alphadiversity indices, abundance of Operational Taxonomic Units and taxonomic groups (and their ratios) provide information on the ecological quality of coastal waters. Our results revealed that only eukaryotic diversity metrics and a limited number of prokaryotic and eukaryotic taxa displayed potential in assessing continental influences in our surveyed area, resulting thus in a restrained potential of microbial plankton diversity as an ecological indicator. Therefore, incorporating microbial plankton diversity in environmental assessment could not always result in a significant improvement of current marine monitoring strategies.

Assessing Viral Abundance and Community Composition in Four Contrasting Regions of the Southern Ocean.

We explored how changes of viral abundance and community composition among four contrasting regions in the Southern Ocean relied on physicochemical and microbiological traits. During January-February 2015, we visited areas north and south of the South Orkney Islands (NSO and SSO) characterized by low temperature and salinity and high inorganic nutrient concentration, north of South Georgia Island (NSG) and west of Anvers Island (WA), which have relatively higher temperatures and lower inorganic nutrient concentrations. Surface viral abundance (VA) was highest in NSG (21.50 ± 10.70 × 106 viruses mL-1) and lowest in SSO (2.96 ± 1.48 × 106 viruses mL-1). VA was positively correlated with temperature, prokaryote abundance and prokaryotic heterotrophic production, chlorophyll a, diatoms, haptophytes, fluorescent organic matter, and isoprene concentration, and was negatively correlated with inorganic nutrients (NO3-, SiO42-, PO43-), and dimethyl sulfide (DMS) concentrations. Viral communities determined by randomly amplified polymorphic DNA-polymerase chain reaction (RAPD-PCR) were grouped according to the sampling location, being more similar within them than among regions. The first two axes of a canonical correspondence analysis, including physicochemical (temperature, salinity, inorganic nutrients-NO3-, SiO42-, and dimethyl sulfoniopropionate -DMSP- and isoprene concentrations) and microbiological (chlorophyll a, haptophytes and diatom, and prokaryote abundance and prokaryotic heterotrophic production) factors accounted for 62.9% of the variance. The first axis, temperature-related, accounted for 33.8%; the second one, salinity-related, accounted for 29.1%. Thus, different environmental situations likely select different hosts for viruses, leading to distinct viral communities.

Taxonomic differences shape the responses of freshwater aerobic anoxygenic phototrophic bacterial communities to light and predation.

Aerobic anoxygenic phototrophic (AAP) bacteria are a phylogenetically diverse and ubiquitous group of prokaryotes that use organic matter but can harvest light using bacteriochlorophyll a. Although the factors regulating AAP ecology have long been investigated through field surveys, the few available experimental studies have considered AAPs as a group, thus disregarding the potential differential responses between taxonomically distinct AAP assemblages. Here, we used sequencing of the pufM gene to describe the diversity of AAPs in 10 environmentally distinct temperate lakes, and to investigate the taxonomic responses of AAP communities in these lakes when subjected to similar experimental manipulations of light and predator removal. The studied communities were clearly dominated by Limnohabitans AAP but presented a clear taxonomic segregation between lakes presumably driven by local conditions, which was maintained after experimental manipulations. Predation reduction (but not light exposure) caused significant compositional shifts across most assemblages, but the magnitude of these changes could not be clearly related to changes in bulk AAP abundances or taxonomic richness of AAP assemblages during experiments. Only a few operational taxonomic units, which differed taxonomically between lakes, were found to respond positively during experimental treatments. Our results highlight that different freshwater AAP communities respond differently to similar control mechanisms, highlighting that in-depth knowledge on AAP diversity is essential to understand the ecology and potential role of these photoheterotrophs.

Impact of grazing, resource availability and light on prokaryotic growth and diversity in the oligotrophic surface global ocean.

The impact of grazing, resource competition and light on prokaryotic growth and taxonomic composition in subtropical and tropical surface waters were studied through 10 microcosm experiments conducted between 30°N and 30°S in the Atlantic, Pacific and Indian oceans. Under natural sunlight conditions, significant changes in taxonomic composition were only observed after the reduction of grazing by sample filtration in combination with a decrease in resource competition by sample dilution. Sunlight exposure significantly reduced prokaryote growth (11 ± 6%) and community richness (14 ± 4%) compared to continuous darkness but did not significantly change community composition. The largest growth inhibition after sunlight exposure occurred at locations showing deep mixed layers. The reduction of grazing had an expected and significant positive effect on growth, but caused a significant decrease in community richness (16 ± 6%), suggesting that the coexistence of many different OTUs is partly promoted by the presence of predators. Dilution of the grazer-free prokaryotic community significantly enhanced growth at the level of community, but consistently and sharply reduced the abundance of Prochlorococcus and SAR11 populations. The decline of these oligotrophic bacterial taxa following an increase in resource availability is consistent with their high specialization for exploiting the limited resources available in the oligotrophic warm ocean.

Delineation of ecologically distinct units of marine Bacteroidetes in the Northwestern Mediterranean Sea.

Bacteroidetes is one of the dominant phyla of ocean bacterioplankton, yet its diversity and population structure is poorly understood. To advance in the delineation of ecologically meaningful units within this group, we constructed near full-length 16S rRNA gene clone libraries from contrasting marine environments in the NW Mediterranean. Based on phylogeny and the associated ecological variables (depth and season), 24 different Bacteroidetes clades were delineated. By considering their relative abundance (from iTag amplicon sequencing studies), we described the distribution patterns of each of these clades, delimiting them as Ecologically Significant Taxonomic Units (ESTUs). Spatially, there was almost no overlap among ESTUs at different depths. In deep waters there was predominance of Owenweeksia, Leeuwenhoekiella, Muricauda-related genera, and some depth-associated ESTUs within the NS5 and NS2b marine clades. Seasonally, multi-annual dynamics of recurring ESTUs were present with dominance of some ESTUs within the NS4, NS5 and NS2b marine clades along most of the year, but with variable relative frequencies between months. A drastic change towards the predominance of Formosa-related ESTUs and one ESTU from the NS5 marine clade was typically present after the spring bloom. Even though there are no isolates available for these ESTUs to determine their physiology, correlation analyses identified the environmental preference of some of them. Overall, our results suggest that there is a high degree of niche specialisation within these closely related clades. This work constitutes a step forward in disentangling the ecology of marine Bacteroidetes, which are essential players in organic matter processing in the oceans.

Long-term seasonal and interannual variability of marine aerobic anoxygenic photoheterotrophic bacteria.

We studied the long-term temporal dynamics of the aerobic anoxygenic phototrophic (AAP) bacteria, a relevant functional group in the coastal marine microbial food web, using high-throughput sequencing of the pufM gene coupled with multivariate, time series and co-occurrence analyses at the Blanes Bay Microbial Observatory (NW Mediterranean). Additionally, using metagenomics, we tested whether the used primers captured accurately the seasonality of the most relevant AAP groups. Phylogroup K (Gammaproteobacteria) was the greatest contributor to community structure over all seasons, with phylogroups E and G (Alphaproteobacteria) being prevalent in spring. Diversity indices showed a clear seasonal trend, with maximum values in winter, which was inverse to that of AAP abundance. Multivariate analyses revealed sample clustering by season, with a relevant proportion of the variance explained by day length, temperature, salinity, phototrophic nanoflagellate abundance, chlorophyll a, and silicate concentration. Time series analysis showed robust rhythmic patterns of co-occurrence, but distinct seasonal behaviors within the same phylogroup, and even within different amplicon sequence variants (ASVs) conforming the same operational taxonomic unit (OTU). Altogether, our results picture the AAP assemblage as highly seasonal and recurrent but containing ecotypes showing distinctive temporal niche partitioning, rather than being a cohesive functional group.

Primer Design for an Accurate View of Picocyanobacterial Community Structure by Using High-Throughput Sequencing.

High-throughput sequencing (HTS) of the 16S rRNA gene has been used successfully to describe the structure and dynamics of microbial communities. Picocyanobacteria are important members of bacterioplankton communities, and, so far, they have predominantly been targeted using universal bacterial primers, providing a limited resolution of the picocyanobacterial community structure and dynamics. To increase such resolution, the study of a particular target group is best approached with the use of specific primers. Here, we aimed to design and evaluate specific primers for aquatic picocyanobacterial genera to be used with high-throughput sequencing. Since the various regions of the 16S rRNA gene have different degrees of conservation in different bacterial groups, we therefore first determined which hypervariable region of the 16S rRNA gene provides the highest taxonomic and phylogenetic resolution for the genera Synechococcus, Prochlorococcus, and Cyanobium An in silico analysis showed that the V5, V6, and V7 hypervariable regions appear to be the most informative for this group. We then designed primers flanking these hypervariable regions and tested them in natural marine and freshwater communities. We successfully detected that most (97%) of the obtained reads could be assigned to picocyanobacterial genera. We defined operational taxonomic units as exact sequence variants (zero-radius operational taxonomic units [zOTUs]), which allowed us to detect higher genetic diversity and infer ecologically relevant information about picocyanobacterial community composition and dynamics in different aquatic systems. Our results open the door to future studies investigating picocyanobacterial diversity in aquatic systems.IMPORTANCE The molecular diversity of the aquatic picocyanobacterial community cannot be accurately described using only the available universal 16S rRNA gene primers that target the whole bacterial and archaeal community. We show that the hypervariable regions V5, V6, and V7 of the 16S rRNA gene are better suited to study the diversity, community structure, and dynamics of picocyanobacterial communities at a fine scale using Illumina MiSeq sequencing. Due to its variability, it allows reconstructing phylogenies featuring topologies comparable to those generated when using the complete 16S rRNA gene sequence. Further, we successfully designed a new set of primers flanking the V5 to V7 region whose specificity for picocyanobacterial genera was tested in silico and validated in several freshwater and marine aquatic communities. This work represents a step forward for understanding the diversity and ecology of aquatic picocyanobacteria and sets the path for future studies on picocyanobacterial diversity.