Metabarcoding reveals environmental factors influencing spatio‐temporal variation in pelagic micro‐eukaryotes

Different Microeukaryotic Trophic Groups Show Different Latitudinal Spatial Scale Dependences in Assembly Processes across the Continental Shelves of China

The relative role of stochasticity versus determinism is critically dependent on the spatial scale over which communities are studied. However, only a few studies have attempted to reveal how spatial scales influence the balance of different assembly processes. In this study, we investigated the latitudinal spatial scale dependences in assembly processes of microeukaryotic communities in surface water and sediment along the continental shelves of China. It was hypothesized that different microeukaryotic trophic groups (i.e., autotroph, heterotroph, mixotroph, and parasite) showed different latitudinal scale dependences in their assembly processes. Our results disclosed that the relative importance of different assembly processes depended on a latitudinal space scale for planktonic microeukaryotes. In surface water, as latitudinal difference increased, the relative contributions of homogenous selection and homogenizing dispersal decreased for the entire community, while those of heterogeneous selection and drift increased. The planktonic autotrophic and heterotrophic groups shifted from stochasticity-dominated processes to heterogeneous selection as latitudinal differences surpassed thresholds of 8° and 16°, respectively. For mixotrophic and parasitic groups, however, the assembly processes were always dominated by drift across different spatial scales. The balance of different assembly processes for the autotrophic group was mainly driven by temperature, whereas that of the heterotrophic group was driven by salinity and geographical distance. In sediment, neither the entire microeukaryotic community nor the four trophic groups showed remarkable spatial scale dependences in assembly processes; they were always overwhelmingly dominated by the drift. This work provides a deeper understanding of the distribution mechanisms of microeukaryotes along the continental shelves of China from the perspective of trophic groups.

New insights into microbial community coalescence in the land-sea continuum

The land-sea continuum constitutes a mixing zone where soil microbial communities encounter, via runoff, those inhabiting marine coastal sediment resulting in community coalescence. Here, we propose an experimental approach, mimicking the land-sea continuum, to study the microbial community coalescence events in different situations, by 16S and 18S rRNA genes metabarcoding. The microbial community structure of sediment diverged with the soil inputs. For prokaryotes, phylogenetic enrichment and amplicon sequence variants (ASVs) replacements characterized the community changes in sediment receiving soil inputs. For fungi, despite phylogenetic enrichment was not observed, the fungal ASVs richness was maintained by soil inputs. Comparison of microbial communities revealed ASVs specific to sediment receiving soil inputs, and also ASVs shared with soil and/or runoff. Among these specific ASVs, four bacterial and one fungal ASVs were identified as indicators of coalescence. Our study provides evidences that coalescence involves the mixing of microorganisms and of the environment.

Ecological Dynamics and Co-occurrences Among Prokaryotes and Microeukaryotes in a Diatom Bloom Process in Xiangshan Bay, China

Diatom blooms can significantly affect the succession of microbial communities, yet little is known about the assembly processes and interactions of microbial communities during autumn bloom events. In this study, we investigated the ecological effects of an autumn diatom bloom on prokaryotic communities (PCCs) and microeukaryotic communities (MECs), focusing on their assembly processes and interactions. The PCCs were largely dominated by Alphaproteobacteria, Gammaproteobacteria, Cyanobacteria, and Flavobacteria, while the MECs primarily included Diatomea, Dinoflagellata, and Chlorophyta. The succession of both PCCs and MECs was mainly driven by this diatom bloom and environmental factors, such as nitrate and silicate. Null modeling revealed that homogeneous selection had a more pronounced impact on the structure of PCCs compared with that of MECs. In particular, drift and dispersal limitation cannot be neglected in the assembly processes of MECs. Co-occurrence network analyses showed that Litorimicrobium, Cercozoa, Marine Group I (MGI), Cryptomonadales, Myrionecta, and Micromonas may affect the bloom process. In summary, these results elucidated the complex, robust interactions and obviously distinct assembly mechanisms of PCCs and MECs during a diatom bloom and extend our current comprehension of the ecological mechanisms and microbial interactions involved in an autumn diatom bloom process.

Photoperiod-driven rhythms reveal multi-decadal stability of phytoplankton communities in a highly fluctuating coastal environment

Phytoplankton play a pivotal role in global biogeochemical and trophic processes and provide essential ecosystem services. However, there is still no broad consensus on how and to what extent their community composition responds to environmental variability. Here, high-frequency oceanographic and biological data collected over more than 25 years in a coastal Mediterranean site are used to shed light on the temporal patterns of phytoplankton species and assemblages in their environmental context. Because of the proximity to the coast and due to large-scale variations, environmental conditions showed variability on the short and long-term scales. Nonetheless, an impressive regularity characterised the annual occurrence of phytoplankton species and their assemblages, which translated into their remarkable stability over decades. Photoperiod was the dominant factor related to community turnover and replacement, which points at a possible endogenous regulation of biological processes associated with species-specific phenological patterns, in analogy with terrestrial plants. These results highlight the considerable stability and resistance of phytoplankton communities in response to different environmental pressures, which contrast the view of these organisms as passively undergoing changes that occur at different temporal scales in their habitat, and show how, under certain conditions, biological processes may prevail over environmental forcing.

Seasonal Diversity of Microeukaryotes in the Han River, Korea Through 18S rRNA Gene Metabarcoding

Freshwater ecosystems contain a large diversity of microeukaryotes that play important roles in maintaining their structure. Microeukaryote communities vary in composition and abundance on the basis of temporal and environmental variables and may serve as useful bioindicators of environmental changes. In the present study, 18S rRNA metabarcoding was employed to investigate the seasonal diversity of microeukaryote communities during four seasons in the Han River, Korea. In total, 882 unique operational taxonomic units (OTUs) were detected, including various diatoms, metazoans (e.g., arthropods and rotifers), chlorophytes, and fungi. Although alpha diversity revealed insignificant differences based on seasons, beta diversity exhibited a prominent variation in the community composition as per seasons. The analysis revealed that the diversity of microeukaryotes was primarily driven by seasonal changes in the prevailing conditions of environmental water temperature and dissolved oxygen. Moreover, potential indicator OTUs belonging to diatoms and chlorophytes were associated with seasonal and environmental factors. This analysis was a preliminary study that established a continuous monitoring system using metabarcoding. This approach could be an effective tool to manage the Han River along with other freshwater ecosystems in Korea.

Plankton biodiversity and species co-occurrence based on environmental DNA – a multiple marker study

Metabarcoding in combination with high-throughput sequencing (HTS) allows simultaneous detection of multiple taxa by targeting single or several taxonomically informative gene regions from environmental DNA samples. In this study, a multiple-marker HTS approach was applied to investigate the plankton diversity and seasonal succession in the Baltic Sea from winter to autumn. Four different markers targeting the 16S, 18S, and 28S ribosomal RNA genes were employed, including a marker for more efficient dinoflagellate detection. Typical seasonal changes were observed in phyto- and bacterioplankton communities. In phytoplankton, the appearance patterns of selected common, dominant, or harmful species followed the patterns also confirmed based on 20 years of phytoplankton monitoring data. In the case of zooplankton, both macro- and microzooplankton species were detected. However, no seasonal patterns were detected in their appearance. In total, 15 and 2 new zoo- and phytoplankton species were detected from the Baltic Sea. HTS approach was especially useful for detecting microzooplankton species as well as for investigating the co-occurrence and potential interactions of different taxa. The results of this study further exemplify the efficiency of metabarcoding for biodiversity monitoring and the advantage of employing multiple markers through the detection of species not identifiable based on a single marker survey and/or by traditional morphology-based methods.

Seven-year dynamics of testate amoeba communities driven more by stochastic than deterministic processes in two subtropical reservoirs

Testate amoebae are widely distributed in natural ecosystems and play an important role in the material cycle and energy flow. However, community assembly of testate amoebae is not well understood, especially with regard to the relative importance of the stochastic and deterministic processes over time. In this study, we used Illumina high-throughput sequencing to explore the community assembly of testate amoebae from surface waters in two reservoirs of subtropical China over a seven-year period. Majority of testate amoebae belonged to the rare taxa because their relative abundances were typically lower than 0.01% of the total eukaryotic plankton community. The testate amoeba community dynamics exhibited a stronger interannual than seasonal variation in both reservoirs. Further, species richness, rather than species turnover, accounted for the majority of community variation. Environmental variables explained less than 20% of the variation in community composition of testate amoebae, and the community assembly appeared to be strongly driven by stochastic processes. Based on the Sloan neutral community model, it was found that neutral processes explained more than 65% of community variation. More importantly, the Stegen null model analysis showed that the stochastic processes (e.g., ecological drift) explained a significantly higher percentage of community assembly than deterministic processes over seven years, although deterministic processes were more influential in certain years. Our results provide new perspectives for understanding the ecological patterns, processes and mechanisms of testate amoeba communities in freshwater ecosystems at temporal scale, and have important implications for monitoring plankton diversity and protecting drinking-water resources.

Marine biomonitoring with eDNA: Can metabarcoding of water samples cut it as a tool for surveying benthic communities?

In the marine realm, biomonitoring using environmental DNA (eDNA) of benthic communities requires destructive direct sampling or the setting-up of settlement structures. Comparatively much less effort is required to sample the water column, which can be accessed remotely. In this study we assess the feasibility of obtaining information from the eukaryotic benthic communities by sampling the adjacent water layer. We studied two different rocky-substrate benthic communities with a technique based on quadrat sampling. We also took replicate water samples at four distances (0, 0.5, 1.5, and 20 m) from the benthic habitat. Using broad range primers to amplify a ca. 313 bp fragment of the cytochrome oxidase subunit I gene, we obtained a total of 3,543 molecular operational taxonomic units (MOTUs). The structure obtained in the two environments was markedly different, with Metazoa, Archaeplastida and Stramenopiles being the most diverse groups in benthic samples, and Hacrobia, Metazoa and Alveolata in the water. Only 265 MOTUs (7.5%) were shared between benthos and water samples and, of these, 180 (5.1%) were identified as benthic taxa that left their DNA in the water. Most of them were found immediately adjacent to the benthos, and their number decreased as we moved apart from the benthic habitat. It was concluded that water eDNA, even in the close vicinity of the benthos, was a poor proxy for the analysis of benthic structure, and that direct sampling methods are required for monitoring these complex communities via metabarcoding.

Temporal Variability and Ecological Interactions of Parasitic Marine Syndiniales in Coastal Protist Communities

Syndiniales are a ubiquitous group of protist parasites that infect and kill a wide range of hosts, including harmful bloom-forming dinoflagellates. Despite the importance of parasitism as an agent of plankton mortality, parasite-host dynamics remain poorly understood, especially over time, hindering the inclusion of parasitism in food web and ecosystem models. For a full year in the Skidaway River Estuary (Georgia), we employed weekly 18S rRNA sampling and co-occurrence network analysis to characterize temporal parasite-host infection dynamics of Syndiniales. Over the year, Syndiniales exhibited strong temporal variability, with higher relative abundance from June to October (7 to 28%) than other months in the year (0.01% to 6%). Nonmetric dimensional scaling of Syndiniales composition revealed tight clustering in June to October that coincided with elevated temperatures (23 to 31°C), though in general, abiotic factors poorly explained composition (canonical correspondence analysis [CCA] and partial least-squares [PLS]) and were less important in the network than biotic relationships. Syndiniales amplicon sequence variants (ASVs) were well represented in the co-occurrence network (20% of edges) and had significant positive associations (Spearman r > 0.7), inferred to be putative parasite-host relationships, with known dinoflagellate hosts (e.g., Akashiwo and Gymnodinium) and other protist groups (e.g., ciliates, radiolarians, and diatoms). Positive associations rarely involved a single Syndiniales and dinoflagellate species, implying flexible parasite-host infection dynamics. These findings provide insight into the temporal dynamics of Syndiniales over a full year and reinforce the importance of single-celled parasites in driving plankton population dynamics. Further empirical work is needed to confirm network interactions and to incorporate parasitism within the context of ecosystem models.IMPORTANCE Protist parasites in the marine alveolate group, Syndiniales, have been observed within infected plankton host cells for decades, and recently, global-scale efforts (Tara Ocean exploration) have confirmed their importance within microbial communities. Yet, protist parasites remain enigmatic, particularly with respect to their temporal dynamics and parasite-host interactions. We employed weekly 18S amplicon surveys over a full year in a coastal estuary, revealing strong temporal shifts in Syndiniales parasites, with highest relative abundance during warmer summer to fall months. Though influenced by temperature, Syndiniales population dynamics were also driven by a high frequency of biological interactions with other protist groups, as determined through co-occurrence network analysis. Parasitic interactions implied by the network highlighted a range of confirmed (dinoflagellates) and putative (diatoms) interactions and suggests parasites may be less selective in their preferred hosts. Understanding parasite-host dynamics over space and time will improve our ability to include parasitism as a loss term in microbial food web models.

Variability and Community Composition of Marine Unicellular Eukaryote Assemblages in a Eutrophic Mediterranean Urban Coastal Area with Marked Plankton Blooms and Red Tides

The Thessaloniki Bay is a eutrophic coastal area which has been characterized in recent years by frequent and intense phytoplankton blooms and red tides. The aim of the study was to investigate the underexplored diversity of marine unicellular eukaryotes in four different sampling sites in Thessaloniki Bay during a year of plankton blooms, red tides, and mucilage aggregates. High-Throughput Sequencing (HTS) was applied in extracted DNA from weekly water samples targeting the 18S rRNA gene. In almost all samples, phytoplankton blooms and/or red tides and mucilage aggregates were observed. The metabarcoding analysis has detected the known unicellular eukaryotic groups frequently observed in the Bay, dominated by Bacillariophyta and Dinoflagellata, and revealed taxonomic groups previously undetected in the study area (MALVs, MAST, and Cercozoa). The dominant OTUs were closely related to species known to participate in red tides, harmful blooms, and mucilage aggregates. Other OTUs, present also during the blooms in low abundance (number of reads), were closely related to known harmful species, suggesting the occurrence of rare taxa with potential negative impacts on human health not detectable with classical microscopy. Overall, the unicellular eukaryote assemblages showed temporal patterns rather than small-scale spatial separation responding to the variability of physical and chemical factors.

Multi-marker eDNA metabarcoding survey to assess the environmental impact of three offshore gas platforms in the North Adriatic Sea (Italy)

The environmental DNA (eDNA) metabarcoding represents a new promising tool for biomonitoring and environmental impact assessment. One of the main advantages of eDNA metabarcoding, compared to the traditional morphotaxonomy-based methods, is to provide a more holistic biodiversity information that includes inconspicuous morphologically non-identifiable taxa. Here, we use eDNA metabarcoding to survey marine biodiversity in the vicinity of the three offshore gas platforms in North Adriatic Sea (Italy). We isolated eDNA from 576 water and sediment samples collected at 32 sampling sites situated along four axes at increasing distances from the gas platforms. We obtained about 46 million eDNA sequences for 5 markers from nuclear 18S V1V2, 18S V4, 18S 37F and mitochondrial 16S and COI genes that cover a wide diversity of benthic and planktonic eukaryotes. Our results showed some impact of platform activities on benthic and pelagic communities at very close distance (<50 m), while communities for intermediate (125 m, 250 m, 500 m) and reference (1000 m, 2000 m) sites did not show any particular biodiversity changes that could be related to platforms activities. The most significant community change along the distance gradient was obtained with the 18S V1V2 marker targeting benthic eukaryotes, even though other markers showed similar trends, but to a lesser extent. These results were congruent with the AMBI index inferred from the eDNA sequences assigned to benthic macrofauna. We finally explored the relation between various physicochemical parameters, including hydrocarbons, on benthic community in the case of one of the platforms. Our results showed that these communities were not significantly impacted by most of hydrocarbons, but rather by macro-elements and sediment texture.

Rhythmicity of coastal marine picoeukaryotes, bacteria and archaea despite irregular environmental perturbations

Seasonality in marine microorganisms has been classically observed in phytoplankton blooms, and more recently studied at the community level in prokaryotes, but rarely investigated at the scale of individual microbial taxa. Here we test if specific marine eukaryotic phytoplankton, bacterial and archaeal taxa display yearly rhythms at a coastal site impacted by irregular environmental perturbations. Our seven-year study in the Bay of Banyuls (North Western Mediterranean Sea) shows that despite some fluctuating environmental conditions, many microbial taxa displayed significant yearly rhythms. The robust rhythmicity was found in both autotrophs (picoeukaryotes and cyanobacteria) and heterotrophic prokaryotes. Sporadic meteorological events and irregular nutrient supplies did, however, trigger the appearance of less common non-rhythmic taxa. Among the environmental parameters that were measured, the main drivers of rhythmicity were temperature and day length. Seasonal autotrophs may thus be setting the pace for rhythmic heterotrophs. Similar environmental niches may be driving seasonality as well. The observed strong association between Micromonas and SAR11, which both need thiamine precursors for growth, could be a first indication that shared nutritional niches may explain some rhythmic patterns of co-occurrence.

River Flow Impacts Bacterial and Archaeal Community Structure in Surface Sediments in the Northern Gulf of Mexico

Meiobenthic community structure in the northern Gulf of Mexico has been shown to be driven by geographical differences due to inshore-offshore gradients and location relative to river discharge. Samples collected along three transects spanning Mobile Bay, Alabama, showed significant differences in meiobenthic communities east of the bay compared to those sampled from the west. In contrast, analysis of bacterial and archaeal communities from the same sediment samples shows that the inshore-offshore gradient has minimal impact on their community structure. Significant differences in community structure were observed for Bacteria and Archaea between the east and west samples, but there was no difference in richness or diversity. Grouped by sediment type, higher richness was observed in silty samples compared to sandy samples. Significant differences were also observed among sediment types for community structure with bacteria communities in silty samples having more anaerobic sulfate reducers compared to aerobic heterotrophs, which had higher abundances in sandy sediments. This is likely due to increased organic matter in the silty sediments from the overlying river leading to low oxygen habitats. Most archaeal sequences represented poorly characterized high-level taxa, limiting interpretation of their distributions. Overlap between groups based on transect and sediment characteristics made determining which factor is more important in structuring bacterial and archaeal communities difficult. However, both factors are driven by discharge from the Mobile River. Although inshore-offshore gradients do not affect Bacteria or Archaea to the same extent as the meiobenthic communities, all three groups are strongly affected by sediment characteristics.

The future of biotic indices in the ecogenomic era: Integrating (e)DNA metabarcoding in biological assessment of aquatic ecosystems

The bioassessment of aquatic ecosystems is currently based on various biotic indices that use the occurrence and/or abundance of selected taxonomic groups to define ecological status. These conventional indices have some limitations, often related to difficulties in morphological identification of bioindicator taxa. Recent development of DNA barcoding and metabarcoding could potentially alleviate some of these limitations, by using DNA sequences instead of morphology to identify organisms and to characterize a given ecosystem. In this paper, we review the structure of conventional biotic indices, and we present the results of pilot metabarcoding studies using environmental DNA to infer biotic indices. We discuss the main advantages and pitfalls of metabarcoding approaches to assess parameters such as richness, abundance, taxonomic composition and species ecological values, to be used for calculation of biotic indices. We present some future developments to fully exploit the potential of metabarcoding data and improve the accuracy and precision of their analysis. We also propose some recommendations for the future integration of DNA metabarcoding to routine biomonitoring programs.

A comparative analysis of metabarcoding and morphology-based identification of benthic communities across different regional seas

In a world of declining biodiversity, monitoring is becoming crucial. Molecular methods, such as metabarcoding, have the potential to rapidly expand our knowledge of biodiversity, supporting assessment, management, and conservation. In the marine environment, where hard substrata are more difficult to access than soft bottoms for quantitative ecological studies, Artificial Substrate Units (ASUs) allow for standardized sampling. We deployed ASUs within five regional seas (Baltic Sea, Northeast Atlantic Ocean, Mediterranean Sea, Black Sea, and Red Sea) for 12-26 months to measure the diversity and community composition of macroinvertebrates. We identified invertebrates using a traditional approach based on morphological characters, and by metabarcoding of the mitochondrial cytochrome oxidase I (COI) gene. We compared community composition and diversity metrics obtained using the two methods. Diversity was significantly correlated between data types. Metabarcoding of ASUs allowed for robust comparisons of community composition and diversity, but not all groups were successfully sequenced. All locations were significantly different in taxonomic composition as measured with both kinds of data. We recovered previously known regional biogeographical patterns in both datasets (e.g., low species diversity in the Black and Baltic Seas, affinity between the Bay of Biscay and the Mediterranean). We conclude that the two approaches provide complementary information and that metabarcoding shows great promise for marine monitoring. However, until its pitfalls are addressed, the use of metabarcoding in monitoring of rocky benthic assemblages should be used in addition to classical approaches rather than instead of them.

Under the canopy: Community-wide effects of invasive algae in Marine Protected Areas revealed by metabarcoding

We analysed with multigene (18S and COI) metabarcoding the effects of the proliferation of invasive seaweeds on rocky littoral communities in two Spanish Marine Protected Areas. The invasive algae studied were Caulerpa cylindracea, Lophocladia lallemandii and Asparagopsis armata. They are canopy-forming, landscape-dominant seaweeds, and we were interested in their effects on the underlying communities of meiobenthos and macrobenthos, separated in two size fractions through sieving. A new semiquantitative treatment of metabarcoding data is introduced. The results for both markers showed that the presence of the invasive seaweed had a significant effect on the understory communities for Lophocladia lallemandii and Asparagopsis armata but not for Caulerpa cylindracea. Likewise, changes in MOTU richness and diversity with invasion status varied in magnitude and direction depending on the alga considered. Our results showed that metabarcoding allows monitoring of the less conspicuous, but not least important, effects of the presence of dominant invasive seaweeds.

How are species interactions structured in species-rich communities? A new method for analysing time-series data

Estimation of intra- and interspecific interactions from time-series on species-rich communities is challenging due to the high number of potentially interacting species pairs. The previously proposed sparse interactions model overcomes this challenge by assuming that most species pairs do not interact. We propose an alternative model that does not assume that any of the interactions are necessarily zero, but summarizes the influences of individual species by a small number of community-level drivers. The community-level drivers are defined as linear combinations of species abundances, and they may thus represent e.g. the total abundance of all species or the relative proportions of different functional groups. We show with simulated and real data how our approach can be used to compare different hypotheses on community structure. In an empirical example using aquatic microorganisms, the community-level drivers model clearly outperformed the sparse interactions model in predicting independent validation data.

Spatial and temporal variation of intertidal nematodes in the northern Gulf of Mexico after the Deepwater Horizon oil spill

Nematodes are an abundant and diverse interstitial component of sedimentary habitats that have been reported to serve as important bioindicators. Though the 2010 Deepwater Horizon (DWH) disaster occurred 60 km offshore in the Gulf of Mexico (GOM) at a depth of 1525 m, oil rose to the surface and washed ashore, subjecting large segments of coastline in the northern GOM to contamination. Previous metabarcoding work shows intertidal nematode communities were negatively affected by the oil spill. Here we examine the subsequent recovery of nematode community structure at five sites along the Alabama coast over a two-year period. The latter part of the study (July 2011-July 2012) also included an examination of nematode vertical distribution in intertidal sediments. Results showed nematode composition within this region was more influenced by sample locality than time and depth. The five sampling sites were characterized by distinct nematode assemblages that varied by sampling dates. Nematode diversity decreased four months after the oil spill but increased after one year, returning to previous levels at all sites except Bayfront Park (BP). There was no significant difference among nematode assemblages in reference to vertical distribution. Although the composition of nematode assemblages changed, the feeding guilds they represented were not significantly different even though some variation was noted. Data from morphological observations integrated with metabarcoding data indicated similar spatial variation in nematode distribution patterns, indicating the potential of using these faster approaches to examine overall disturbance impact trends within communities. Heterogeneity of microhabitats in the intertidal zone indicates that future sampling and fine-scale studies of nematodes are needed to examine such anthropogenic effects.

Genetic assessment of meiobenthic community composition and spatial distribution in coastal sediments along northern Gulf of Mexico

Meiobenthic (meiofauna and micro-eukaryotes) organisms are important contributors to ecosystem functioning in aquatic environments through their roles in nutrient transport, sediment stability, and food web interactions. Despite their ecological importance, information pertaining to variation of these communities at various spatial and temporal scales is not widely known. Many studies in the Gulf of Mexico (GOM) have focused either on deep sea or continental shelf areas, while little attention has been paid to bays and coastal regions. Herein, we take a holistic approach by using high-throughput sequencing approaches to examine spatial variation in meiobenthic communities within Alabama bays and the coastal northern GOM region. Sediment samples were collected along three transects (Mississippi Sound: MS, FOCAL: FT, and Orange Beach: OB) from September 2010 to April 2012 and community composition was determined by metabarcoding the V9 hypervariable region of the nuclear18S rRNA gene. Results showed that Stramenopiles (diatoms), annelids, arthropods (copepods), and nematodes were the dominate groups within samples, while there was presence of other phyla throughout the dataset. Location played a larger role than time sampled in community composition. However, samples were collected over a short temporal scale. Samples clustered in reference to transect, with the most eastern transect (OB) having a distinct community composition in comparison to the other two transects (MS and FT). Communities also differed in reference to region (Bay versus Shelf). Bulk density and percent inorganic carbon were the only measured environmental factors that were correlated with community composition.

Spatial and temporal patterns in the Pelagibacteraceae across an estuarine gradient

Marine bacterial communities show strong spatial and seasonal patterns, often characterized by changes at high taxonomic levels. The Pelagibacteraceae are common members of bacterial communities, with well-documented biogeography at the subclade level. To identify patterns within the subclades, the abundance and diversity of Pelagibacteraceae were analyzed over a two-year period at four stations across an estuarine gradient. Pelagibacteraceae was the most abundant bacterial family, averaging 27% of the community, but varying from 1% to 57% in any one sample. Highest abundances were detected in autumn and winter. Pelagibacteraceae richness was lowest at the most inshore site, and highest in autumn and winter at all sites. Shannon diversity decreased in winter, when a few OTUs dominated the community. Dissolved oxygen, dissolved silicate and prokaryote abundance explained most of the variability in the Pelagibacteraceae communities, with salinity differentiating low salinity communities. The 10 most abundant OTUs included OTUs that varied across sites, with little seasonality as well as those with small site effects, but strong seasonal patterns indicating differences in the niches of individual OTUs. While salinity was important in structuring low salinity communities, higher salinity communities appear to be responding to additional environmental parameters including oxygen, nutrients and other organisms.