Diversity and ecology of protists revealed by metabarcoding

Nano- and microplastics drive the dynamic equilibrium of amoeba-associated bacteria and antibiotic resistance genes

As emerging pollutants, microplastics have become pervasive on a global scale, inflicting significant harm upon ecosystems. However, the impact of these microplastics on the symbiotic relationship between protists and bacteria remains poorly understood. In this study, we investigated the mechanisms through which nano- and microplastics of varying sizes and concentrations influence the amoeba-bacterial symbiotic system. The findings reveal that nano- and microplastics exert deleterious effects on the adaptability of the amoeba host, with the magnitude of these effects contingent upon particle size and concentration. Furthermore, nano- and microplastics disrupt the initial equilibrium in the symbiotic relationship between amoeba and bacteria, with nano-plastics demonstrating a reduced ability to colonize symbiotic bacteria within the amoeba host when compared to their microplastic counterparts. Moreover, nano- and microplastics enhance the relative abundance of antibiotic resistance genes and heavy metal resistance genes in the bacteria residing within the amoeba host, which undoubtedly increases the potential transmission risk of both human pathogens and resistance genes within the environment. In sum, the results presented herein provide a novel perspective and theoretical foundation for the study of interactions between microplastics and microbial symbiotic systems, along with the establishment of risk assessment systems for ecological environments and human health.

Microbial, holobiont, and Tree of Life eDNA/eRNA for enhanced ecological assessment

Microbial environmental DNA and RNA (collectively 'eNA') originate from a diverse and abundant array of microbes present in environmental samples. These eNA signals, largely representing whole organisms, serve as a powerful complement to signals derived from fragments or remnants of larger organisms. Integrating microbial data into the toolbox of ecosystem assessments and biotic indices therefore has the potential to transform how we use eNA data to understand biodiversity dynamics and ecosystem functions, and to inform the next generation of environmental monitoring. Incorporating holobiont and Tree of Life approaches into eNA analyses offers further holistic insight into the range of ecological interactions between microbes and other organisms, paving the way for advancing our understanding of, and ultimately manipulating ecosystem properties pertinent to environmental management, conservation, wildlife health, and food production.

Unveiling the microbial realm with VEBA 2.0: a modular bioinformatics suite for end-to-end genome-resolved prokaryotic, (micro)eukaryotic and viral multi-omics from either short- or long-read sequencing

The microbiome is a complex community of microorganisms, encompassing prokaryotic (bacterial and archaeal), eukaryotic, and viral entities. This microbial ensemble plays a pivotal role in influencing the health and productivity of diverse ecosystems while shaping the web of life. However, many software suites developed to study microbiomes analyze only the prokaryotic community and provide limited to no support for viruses and microeukaryotes. Previously, we introduced the Viral Eukaryotic Bacterial Archaeal (VEBA) open-source software suite to address this critical gap in microbiome research by extending genome-resolved analysis beyond prokaryotes to encompass the understudied realms of eukaryotes and viruses. Here we present VEBA 2.0 with key updates including a comprehensive clustered microeukaryotic protein database, rapid genome/protein-level clustering, bioprospecting, non-coding/organelle gene modeling, genome-resolved taxonomic/pathway profiling, long-read support, and containerization. We demonstrate VEBA's versatile application through the analysis of diverse case studies including marine water, Siberian permafrost, and white-tailed deer lung tissues with the latter showcasing how to identify integrated viruses. VEBA represents a crucial advancement in microbiome research, offering a powerful and accessible software suite that bridges the gap between genomics and biotechnological solutions.

The global genetic diversity of planktonic foraminifera reveals the structure of cryptic speciation in plankton

The nature and extent of diversity in the plankton has fascinated scientists for over a century. Initially, the discovery of many new species in the remarkably uniform and unstructured pelagic environment appeared to challenge the concept of ecological niches. Later, it became obvious that only a fraction of plankton diversity had been formally described, because plankton assemblages are dominated by understudied eukaryotic lineages with small size that lack clearly distinguishable morphological features. The high diversity of the plankton has been confirmed by comprehensive metabarcoding surveys, but interpretation of the underlying molecular taxonomies is hindered by insufficient integration of genetic diversity with morphological taxonomy and ecological observations. Here we use planktonic foraminifera as a study model and reveal the full extent of their genetic diversity and investigate geographical and ecological patterns in their distribution. To this end, we assembled a global data set of ~7600 ribosomal DNA sequences obtained from morphologically characterised individual foraminifera, established a robust molecular taxonomic framework for the observed diversity, and used it to query a global metabarcoding data set covering ~1700 samples with ~2.48 billion reads. This allowed us to extract and assign 1 million reads, enabling characterisation of the structure of the genetic diversity of the group across ~1100 oceanic stations worldwide. Our sampling revealed the existence of, at most, 94 distinct molecular operational taxonomic units (MOTUs) at a level of divergence indicative of biological species. The genetic diversity only doubles the number of formally described species identified by morphological features. Furthermore, we observed that the allocation of genetic diversity to morphospecies is uneven. Only 16 morphospecies disguise evolutionarily significant genetic diversity, and the proportion of morphospecies that show genetic diversity increases poleward. Finally, we observe that MOTUs have a narrower geographic distribution than morphospecies and that in some cases the MOTUs belonging to the same morphospecies (cryptic species) have different environmental preferences. Overall, our analysis reveals that even in the light of global genetic sampling, planktonic foraminifera diversity is modest and finite. However, the extent and structure of the cryptic diversity reveals that genetic diversification is decoupled from morphological diversification, hinting at different mechanisms acting at different levels of divergence.

Environmental DNA metabarcoding reveals diverse phytoplankton assemblages and potentially harmful algal distribution along the urban coasts of Türkiye

Marine phytoplankton are widely used to monitor the state of the water column due to their rapid changes in response to environmental conditions. In this study, we aimed to investigate the coastal phytoplankton assemblages, including bloom-forming species using high-throughput sequencing of 18S rRNA genes targeting the V4 region and their relationship with environmental variables along the Istanbul coasts of the Sea of Marmara. A total of 118 genera belonging to six phyla were detected. Among them, Dinoflagellata (36) and Bacillariophyta (26) were represented with the highest number of genera. According to the relative abundance of DNA reads, the most abundant taxa were Dinoflagellata_phylum (18.1%), Emiliania (8.4%), Biecheleria (8.4), and Noctiluca (8.1%). The ANOSIM test showed that there was a significant temporal difference in the assemblages, while the driving environmental factors were pH, water temperature, and salinity. According to the TRIX index, the trophic state of the coasts was highly mesotrophic and eutrophic. In addition, 45 bloom-forming and HAB taxa were detected and two species of Noctiluca and Emiliania, which frequently cause blooms in the area, were recorded in high abundance. Our results provide insight into the phytoplankton assemblages along the urbanized coastlines by analysing the V4 region of 18S rRNA. This data can support future studies that use both traditional methods and metabarcoding, employing various primers and targeting different genes and regions.

High above the rest: standing behaviors in the amoebae of Sappinia and Thecamoeba

Many terrestrial microbes have evolved cell behaviors that help them rise above their substrate, often to facilitate dispersal. One example of these behaviors is found in the amoebae of Sappinia pedata, which actively lift most of their cell mass above the substrate, known as standing. This standing behavior was first described in S. pedata in the 1890s from horse dung isolates but never molecularly characterized from dung. Our study expands this understanding, revealing the first molecularly confirmed S. pedata from herbivore dung in Mississippi, USA, and describing a new species, Sappinia dangeardi n. sp., with larger trophozoite cells. Additionally, we isolated another standing amoeba, Thecamoeba homeri n. sp., from soil, exhibiting a previously unreported "doughnut shape" transient behavior. In S. dangeardi n. sp., we discovered that standing is likely triggered by substrate drying, and that actin filaments actively localize in the "stalk" to support the standing cells, as observed through confocal microscopy. While the purpose of standing behaviors has not been investigated, we hypothesize it is energetically expensive and therefore a significant evolutionary strategy in these organisms. Overall, this study emphasizes behavioral adaptations to terrestrial environments within Amoebozoa, stressing the importance of diverse laboratory conditions that replicate natural habitats.

"Candidatus Intestinibacterium parameciiphilum"-member of the "Candidatus Paracaedibacteraceae" family (Alphaproteobacteria, Holosporales) inhabiting the ciliated protist Paramecium

Protists frequently host diverse bacterial symbionts, in particular those affiliated with the order Holosporales (Alphaproteobacteria). All characterised members of this bacterial lineage have been retrieved in obligate association with a wide range of eukaryotes, especially multiple protist lineages (e.g. amoebozoans, ciliates, cercozoans, euglenids, and nucleariids), as well as some metazoans (especially arthropods and related ecdysozoans). While the genus Paramecium and other ciliates have been deeply investigated for the presence of symbionts, known members of the family "Candidatus Paracaedibacteraceae" (Holosporales) are currently underrepresented in such hosts. Herein, we report the description of "Candidatus Intestinibacterium parameciiphilum" within the family "Candidatus Paracaedibacteraceae", inhabiting the cytoplasm of Paramecium biaurelia. This novel bacterium is almost twice as big as its relative "Candidatus Intestinibacterium nucleariae" from the opisthokont Nuclearia and does not present a surrounding halo. Based on phylogenetic analyses of 16S rRNA gene sequences, we identified six further potential species-level lineages within the genus. Based on the provenance of the respective samples, we investigated the environmental distribution of the representatives of "Candidatus Intestinibacterium" species. Obtained results are consistent with an obligate endosymbiotic lifestyle, with protists, in particular freshwater ones, as hosts. Thus, available data suggest that association with freshwater protists could be the ancestral condition for the members of the "Candidatus Intestinibacterium" genus.

Arctic cyanobacterial mat community diversity decreases with latitude across the Canadian Arctic

Cyanobacterial mats are commonly reported as hotspots of microbial diversity across polar environments. These thick, multilayered microbial communities provide a refuge from extreme environmental conditions, with many species able to grow and coexist despite the low allochthonous nutrient inputs. The visibly dominant phototrophic biomass is dependent on internal nutrient recycling by heterotrophic organisms within the mats; however, the specific contribution of heterotrophic protists remains little explored. In this study, mat community diversity was examined along a latitudinal gradient (55-83°N), spanning subarctic taiga, tundra, polar desert, and the High Arctic ice shelves. The prokaryotic and eukaryotic communities were targeted, respectively, by V4 16S ribosomal RNA (rRNA) and V9 18S rRNA gene amplicon high-throughput sequencing. Prokaryotic and eukaryotic richness decreased, in tandem with decreasing temperatures and shorter seasons of light availability, from the subarctic to the High Arctic. Taxonomy-based annotation of the protist community revealed diverse phototrophic, mixotrophic, and heterotrophic genera in all mat communities, with fewer parasitic taxa in High Arctic communities. Co-occurrence network analysis identified greater heterogeneity in eukaryotic than prokaryotic community structure among cyanobacterial mats across the Canadian Arctic. Our findings highlight the sensitivity of microbial eukaryotes to environmental gradients across northern high latitudes.

eDNA metabarcoding reveals biodiversity and depth stratification patterns of dinoflagellate assemblages within the epipelagic zone of the western Coral Sea

BACKGROUND

Dinoflagellates play critical roles in the functioning of marine ecosystems but also may pose a hazard to human and ecosystem health by causing harmful algal blooms (HABs). The Coral Sea is a biodiversity hotspot, but its dinoflagellate assemblages in pelagic waters have not been studied by modern sequencing methods. We used metabarcoding of the 18 S rRNA V4 amplicon to assess the diversity and structure of dinoflagellate assemblages throughout the water column to a depth of 150 m at three stations in the Western Coral Sea. Additionally, at one station we compared metabarcoding with morphological methods to optimise identification and detection of dinoflagellates.

RESULTS

Stratification of dinoflagellate assemblages was evident in depth-specific relative abundances of taxonomic groups; the greatest difference was between the 5-30 m assemblages and the 130-150 m assemblages. The relative abundance of Dinophyceae (photosynthetic and heterotrophic) decreased with increasing depth, whereas that of Syndiniales (parasitic) increased with increasing depth. The composition of major taxonomic groups was similar among stations. Taxonomic richness and diversity of amplicon sequence variants (ASVs) were similar among depths and stations; however, the abundance of dominant taxa was highest within 0-30 m, and the abundance of rare taxa was highest within 130-150 m, indicating adaptations to specific depth strata. The number of unclassified ASVs at the family and species levels was very high, particularly for Syndinian representatives.

CONCLUSIONS

Dinoflagellate assemblages in open water of the Coral Sea are highly diverse and taxonomically stratified by depth; patterns of relative abundance along the depth gradient reflect environmental factors and ecological processes. Metabarcoding detects more species richness than does traditional microscopical methods of sample analysis, yet the methods are complementary, with morphological analysis revealing additional richness. The large number of unclassified dinoflagellate-ASVs indicates a need for improved taxonomic reference databases and suggests presence of dinoflagellate-crypto and-morphospecies.

Protist communities of microbial mats from the extreme environments of five saline Andean lagoons at high altitudes in the Atacama Desert

Introduction

Heterotrophic protists colonizing microbial mats have received little attention over the last few years, despite their importance in microbial food webs. A significant challenge originates from the fact that many protists remain uncultivable and their functions remain poorly understood.

Methods

Metabarcoding studies of protists in microbial mats across high-altitude lagoons of different salinities (4.3-34 practical salinity units) were carried out to provide insights into their vertical stratification at the millimeter scale. DNA and cDNA were analyzed for selected stations.

Results

Sequence variants classified as the amoeboid rhizarian Rhogostoma and the ciliate Euplotes were found to be common members of the heterotrophic protist communities. They were accompanied by diatoms and kinetoplastids. Correlation analyses point to the salinity of the water column as a main driver influencing the structure of the protist communities at the five studied microbial mats. The active part of the protist communities was detected to be higher at lower salinities (<20 practical salinity units).

Discussion

We found a restricted overlap of the protist community between the different microbial mats indicating the uniqueness of these different aquatic habitats. On the other hand, the dominating genotypes present in metabarcoding were similar and could be isolated and sequenced in comparative studies (Rhogostoma, Euplotes, Neobodo). Our results provide a snapshot of the unculturable protist diversity thriving the benthic zone of five athalossohaline lagoons across the Andean plateau.

Insights into freshwater ciliate diversity through high throughput DNA metabarcoding

The freshwater bodies of India are highly biodiverse but still understudied, especially concerning ciliates. Ciliates constitute a significant portion of eukaryotic diversity and play crucial roles in microbial loops, nutrient recycling, and ecosystem maintenance. The present study aimed to elucidate ciliate diversity in three freshwater sites in the Delhi region of India: Okhla Bird Sanctuary (OBS), Sanjay Lake (SL), and Raj Ghat pond (RJ). This study represents the first investigation into the taxonomic diversity and richness of freshwater ciliates in India using a high-throughput DNA metabarcoding approach. For the analysis, total environmental DNA was extracted from the three freshwater samples, followed by sequencing of the 18S V4 barcode region and subsequent phylogenetic analyses. Operational taxonomic units (OTU) analyses revealed maximum species diversity in OBS (106), followed by SL (104) and RJ (99) sites. Ciliates from the classes Oligohymenophorea, Prostomatea, and Spirotrichea were dominant in the three sites. The study discusses the ability of the metabarcoding approach to uncover unknown and rare species. The study highlights the need for refined reference databases and cautious interpretation of the high-throughput sequencing-generated data while emphasizing the complementary nature of molecular and morphological approaches in studying ciliate diversity.

The high molecular diversity in Noctiluca scintillans is dominated by intra-genomic variations revealed by single cell high-throughput sequencing of 18S rDNA V4

The application of high-throughput sequencing (HTS) technologies has revolutionized research on phytoplankton biodiversity by generating an unprecedented amount of molecular data in marine ecosystem surveys. However, high-level of molecular diversity uncovered in HTS-based metabarcoding analyses may lead to overinterpretation of phytoplankton diversity due to excessive intra-genomic variations (IGVs). The aims in this study are to explore the nature of phytoplankton molecular diversity and to test the hypothesis. We carried out single-cell metabarcoding analysis of 18S rDNA V4 sequences obtained in single Noctiluca scintillans cells isolated from various sites in coastal waters of China. Results showed that each single N. scintillans cell harbored a high level of IGVs with about 100 amplicon sequence variants (ASVs). The large numbers of non-dominant ASVs identified in N. scintillans cells, which might correspond to the larger numbers of ASVs annotated as N. scintillans and showed similar temporal dynamics in metabarcoding analyses, could inflate the inter-species diversity or intra-species genetic diversity. In addition, there were large numbers of additional ASVs that were not annotated as N. scintillans. These non-N. scintillans ASVs might represent diverse preys for N. scintillans, consistent with previous reports that N. scintillans may act as chance predator of a broad-spectrum preys. This single-cell study has unambiguously demonstrated that the existence of high levels of IGVs in N. scintillans and most likely many other phytoplankton species, demonstrating that the majority of the molecular diversity revealed in metabarcoding analysis, which were generally interpreted as the sum of inter-species diversity and intra-species diversity, actually included high levels of IGVs and should be interpreted with caution.

VESPA: an optimized protocol for accurate metabarcoding-based characterization of vertebrate eukaryotic endosymbiont and parasite assemblages

Protocols for characterizing taxonomic assemblages by deep sequencing of short DNA barcode regions (metabarcoding) have revolutionized our understanding of microbial communities and are standardized for bacteria, archaea, and fungi. Unfortunately, comparable methods for host-associated eukaryotes have lagged due to technical challenges. Despite 54 published studies, issues remain with primer complementarity, off-target amplification, and lack of external validation. Here, we present VESPA (Vertebrate Eukaryotic endoSymbiont and Parasite Analysis) primers and optimized metabarcoding protocol for host-associated eukaryotic community analysis. Using in silico prediction, panel PCR, engineered mock community standards, and clinical samples, we demonstrate VESPA to be more effective at resolving host-associated eukaryotic assemblages than previously published methods and to minimize off-target amplification. When applied to human and non-human primate samples, VESPA enables reconstruction of host-associated eukaryotic endosymbiont communities more accurately and at finer taxonomic resolution than microscopy. VESPA has the potential to advance basic and translational science on vertebrate eukaryotic endosymbiont communities, similar to achievements made for bacterial, archaeal, and fungal microbiomes.

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.

Sea-ice melt determines seasonal phytoplankton dynamics and delimits the habitat of temperate Atlantic taxa as the Arctic Ocean atlantifies

The Arctic Ocean is one of the regions where anthropogenic environmental change is progressing most rapidly and drastically. The impact of rising temperatures and decreasing sea ice on Arctic marine microbial communities is yet not well understood. Microbes form the basis of food webs in the Arctic Ocean, providing energy for larger organisms. Previous studies have shown that Atlantic taxa associated with low light are robust to more polar conditions. We compared to which extent sea ice melt influences light-associated phytoplankton dynamics and biodiversity over two years at two mooring locations in the Fram Strait. One mooring is deployed in pure Atlantic water, and the second in the intermittently ice-covered Marginal Ice Zone. Time-series analysis of amplicon sequence variants abundance over a 2-year period, allowed us to identify communities of co-occurring taxa that exhibit similar patterns throughout the annual cycle. We then examined how alterations in environmental conditions affect the prevalence of species. During high abundance periods of diatoms, polar phytoplankton populations dominated, while temperate taxa were weakly represented. Furthermore, we found that polar pelagic and ice-associated taxa, such as Fragilariopsis cylindrus and Melosira arctica, were more common in Atlantic conditions, while temperate taxa, such as Odontella aurita and Proboscia alata, were less abundant under polar conditions. This suggests that sea ice melt may act as a barrier to the northward expansion of temperate phytoplankton, preventing their dominance in regions still strongly influenced by polar conditions. Our findings highlight the complex interactions between sea ice melt, phytoplankton dynamics, and biodiversity in the Arctic.

Evaluation of sequential filtration and centrifugation to capture environmental DNA and survey microbial eukaryotic communities in aquatic environments

Sequential membrane filtration of water samples is commonly used to monitor the diversity of aquatic microbial eukaryotes. This capture method is efficient to focus on specific taxonomic groups within a size fraction, but it is time-consuming. Centrifugation, often used to collect microorganisms from pure culture, could be seen as an alternative to capture microbial eukaryotic communities from environmental samples. Here, we compared the two capture methods to assess diversity and ecological patterns of eukaryotic communities in the Thau lagoon, France. Water samples were taken twice a month over a full year and sequential filtration targeting the picoplankton (0.2-3 μm) and larger organisms (>3 μm) was used in parallel to centrifugation. The microbial eukaryotic community in the samples was described using an environmental DNA approach targeting the V4 region of the 18S rRNA gene. The most abundant divisions in the filtration fractions and the centrifugation pellet were Dinoflagellata, Metazoa, Ochrophyta, Cryptophyta. Chlorophyta were dominant in the centrifugation pellet and the picoplankton fraction but not in the larger fraction. Diversity indices and structuring patterns of the community in the two size fractions and the centrifugation pellet were comparable. Twenty amplicon sequence variants were significantly differentially abundant between the two size fractions and the centrifugation pellet, and their temporal patterns of abundance in the two fractions combined were similar to those obtained with centrifugation. Overall, centrifugation led to similar ecological conclusions as the two filtrated fractions combined, thus making it an attractive time-efficient alternative to sequential filtration.

Increased contribution of parasites in microbial eukaryotic communities of different Aegean Sea coastal systems

Background-Aim

Protistan communities have a major contribution to biochemical processes and food webs in coastal ecosystems. However, related studies are scarce and usually limited in specific groups and/or sites. The present study examined the spatial structure of the entire protistan community in seven different gulfs and three different depths in a regional Mediterranean Sea, aiming to define taxa that are important for differences detected in the marine microbial network across the different gulfs studied as well as their trophic interactions.

Methods

Protistan community structure analysis was based on the diversity of the V2-V3 hypervariable region of the 18S rRNA gene. Operational taxonomic units (OTUs) were identified using a 97% sequence identity threshold and were characterized based on their taxonomy, trophic role, abundance and niche specialization level. The differentially abundant, between gulfs, OTUs were considered for all depths and interactions amongst them were calculated, with statistic and network analysis.

Results

It was shown that Dinophyceae, Bacillariophyta and Syndiniales were the most abundant groups, prevalent in all sites and depths. Gulfs separation was more striking at surface corroborating with changes in environmental factors, while it was less pronounced in higher depths. The study of differentially abundant, between gulfs, OTUs revealed that the strongest biotic interactions in all depths occurred between parasite species (mainly Syndiniales) and other trophic groups. Most of these species were generalists but not abundant highlighting the importance of rare species in protistan community assemblage.

Conclusion

Overall this study revealed the emergence of parasites as important contributors in protistan network regulation regardless of depth.

Parasites of firebugs in Austria with focus on the "micro"-diversity of the cosmopolitan trypanosomatid Leptomonas pyrrhocoris

In this work, we investigated parasites of the firebug Pyrrhocoris apterus in Austria and demonstrated that in addition to the extensively studied Leptomonas pyrrhocoris, it can also be infected by Blastocrithidia sp. and by a mermithid, which for the first time has been characterized using molecular methods. This diversity can be explained by the gregarious lifestyle, as well as the coprophagous and cannibalistic behavior of the insect hosts that makes them susceptible to various parasites. In addition, we showed no tight association of the L. pyrrhocoris haplotypes and geographical locations (at least, considering the relatively small scale of locations in Austria) implying that the natural populations of L. pyrrhocoris are mixed due to the mobility of their firebug hosts.

Mesoscale oceanographic meanders influence protist community function and structure in the southern Indian Ocean

The interface between the nutrient-rich Southern Ocean and oligotrophic Indian Ocean creates unique environmental conditions that can strongly influence biological processes. We investigated protist communities across a mesoscale meander of the Subtropical Front within the Southern Indian Ocean. 18S V9 rDNA metabarcoding suggests a diverse protist community in which the dinoflagellates and parasitic Syndiniales were abundant. Diversity was highest in frontal waters of the mesoscale meander, with differences in community structure inside and outside the meander. While the overall community was dominated by mixotrophic taxa, the frontal boundary of the meander had increased abundances of heterotrophic taxa, with potential implications for net atmospheric CO2 drawdown. Pulse amplitude modulated (PAM) fluorimetry revealed significant differences in the photophysiology of phytoplankton communities inside and outside the meander. By using single-cell PAM microscopy, we identified physiological differences between dinoflagellate and coccolithophore taxa, which may have contributed to changes in photophysiology observed at community level. Overall, our results demonstrate that frontal areas have a strong impact on the composition of protist communities in the Southern Ocean with important implications for understanding biological processes in this region.

Microbiome predators in changing soils

Microbiome predators shape the soil microbiome and thereby soil functions. However, this knowledge has been obtained from small-scale observations in fundamental rather than applied settings and has focused on a few species under ambient conditions. Therefore, there are several unaddressed questions on soil microbiome predators: (1) What is the role of microbiome predators in soil functioning? (2) How does global change affect microbiome predators and their functions? (3) How can microbiome predators be applied in agriculture? We show that there is sufficient evidence for the vital role of microbiome predators in soils and stress that global changes impact their functions, something that urgently needs to be addressed to better understand soil functioning as a whole. We are convinced that there is a potential for the application of microbiome predators in agricultural settings, as they may help to sustainably increase plant growth. Therefore, we plea for more applied research on microbiome predators.

Spatiotemporal variation of phytoplankton communities and water quality among seaweed, shellfish and cage fish culture systems

Various marine aquaculture systems have different impacts on the environment, but few assessments were focused on the environmental impact by different systems in the same region. To study the effects of various aquaculture systems on phytoplankton community structure and water properties, 5 surveys were carried out in seaweed (Gracilaria lemaneiformis, GL), shellfish (Mytilus coruscus, MC) and cage fish (Larimichthys crocea, LC) mariculture areas in Dongji island, Zhejiang, China from June to September 2020. Significant differences were observed in some environmental parameters and phytoplankton communities among three aquaculture systems. The dissolved oxygen concentrations and Secchi depth in the surface waters in GL area were relatively higher than in the blank and other areas. As for nutrients concentration, LC and MC areas had higher concentrations than blank area, while GL area was the lowest. Though Chlorophyll-a concentration displayed fluctuations, relatively lower concentrations were found in GL area. Shannon diversity index was found to be relatively constant and higher in GL area. The Non-metric multidimensional scaling results revealed that phytoplankton composition had a distinct pattern among sampling times. The correlations and Redundancy analysis showed that total nitrogen, salinity and transparency were primary environmental factors associated with phytoplankton composition. Our study confirmed that different marine aquaculture systems can cause environmental fluctuations. Among the three systems, seaweed cultivation can bring multiple positive effects by improving surrounding water quality and increasing the phytoplankton composition. G. lemaneiformis culture in summer has great positive impact on seawater environment as it can maintain the ecological balance and reduce the risk of harmful algal blooms (HABs), and therefore, it is strongly recommended more G. lemaneiformis cultivation in this region in summer.

High Microeukaryotic Diversity in the Cold-Seep Sediment

Microeukaryotic diversity, community structure, and their regulating mechanisms remain largely unclear in chemosynthetic ecosystems. Here, using high-throughput sequencing data of 18S rRNA genes, we explored microeukaryotic communities from the Haima cold seep in the northern South China Sea. We compared three distinct habitats: active, less active, and non-seep regions, with vertical layers (0-25 cm) from sediment cores. The results showed that seep regions harbored more abundant and diverse parasitic microeukaryotes (e.g., Apicomplexa and Syndiniales) as indicator species, compared to nearby non-seep region. Microeukaryotic community heterogeneity was larger between habitats than within habitat, and greatly increased when considering molecular phylogeny, suggesting the local diversification in cold-seep sediments. Microeukaryotic α-diversity at cold seeps was positively increased by metazoan richness and dispersal rate of microeukaryotes, while its β-diversity was promoted by heterogeneous selection mainly from metazoan communities (as potential hosts). Their combined effects led to the significant higher γ-diversity (i.e., total diversity in a region) at cold seeps than non-seep regions, suggesting cold-seep sediment as a hotspot for microeukaryotic diversity. Our study highlights the importance of microeukaryotic parasitism in cold-seep sediment and has implications for the roles of cold seep in maintaining and promoting marine biodiversity.

Response of planktonic microbial assemblages to disturbance in an urban sub-tropical estuary

Microbes are sensitive indicators of estuarine processes because they respond rapidly to dynamic disturbance events. As most of the world's population lives in urban areas and climate change-related disturbance events are becoming more frequent, estuaries bounded by cities are experiencing increasing stressors, at the same time that their ecosystem services are required more than ever. Here, using a multidisciplinary approach, we determined the response of planktonic microbial assemblages in response to seasonality and a rainfall disturbance in an urban estuary bounded by Australia's largest city, Sydney. We used molecular barcoding (16S, 18S V4 rRNA) and microscopy-based identification to compare microbial assemblages at locations with differing characteristics and urbanisation histories. Across 142 samples, we identified 8,496 unique free-living bacterial zOTUs, 8,175 unique particle associated bacterial zOTUs, and 1,920 unique microbial eukaryotic zOTUs. Using microscopy, we identified only the top <10% abundant, larger eukaryotic taxa (>10 µm), however quantification was possible. The site with the greater history of anthropogenic impact showed a more even community of associated bacteria and eukaryotes, and a significant increase in dissolved inorganic nitrogen following rainfall, when compared to the more buffered site. This coincided with a reduced proportional abundance of Actinomarina and Synechococcus spp., a change in SAR 11 clades, and an increase in the eukaryotic microbial groups Dinophyceae, Mediophyceae and Bathyoccocaceae, including a temporary dominance of the harmful algal bloom dinoflagellate Prorocentrum cordatum (syn. P. minimum). Finally, a validated hydrodynamic model of the estuary supported these results, showing that the more highly urbanised and upstream location consistently experienced a higher magnitude of salinity reduction in response to rainfall events during the study period. The best abiotic variables to explain community dissimilarities between locations were TDP, PN, modelled temperature and salinity (r = 0.73) for the free living bacteria, TP for the associated bacteria (r = 0.43), and modelled temperature (r = 0.28) for the microbial eukaryotic communities. Overall, these results show that a minor disturbance such as a brief rainfall event can significantly shift the microbial assemblage of an anthropogenically impacted area within an urban estuary to a greater degree than a seasonal change, but may result in a lesser response to the same disturbance at a buffered, more oceanic influenced location. Fine scale research into the factors driving the response of microbial communities in urban estuaries to climate related disturbances will be necessary to understand and implement changes to maintain future estuarine ecosystem services.

Extreme environments offer an unprecedented opportunity to understand microbial eukaryotic ecology, evolution, and genome biology

Research in extreme environments has substantially expanded our understanding of the ecology and evolution of life on Earth, but a major group of organisms has been largely overlooked: microbial eukaryotes (i.e., protists). In this Perspective, we summarize data from over 80 studies of protists in extreme environments and identify focal lineages that are of significant interest for further study, including clades within Echinamoebida, Heterolobosea, Radiolaria, Haptophyta, Oomycota, and Cryptophyta. We argue that extreme environments are prime sampling targets to fill gaps in the eukaryotic tree of life and to increase our understanding of the ecology, metabolism, genome architecture, and evolution of eukaryotic life.

Groundwater environmental DNA metabarcoding reveals hidden diversity and reflects land-use and geology

Despite being the most important source of liquid freshwater on the planet, groundwater is severely threatened by climate change, agriculture, or industrial mining. It is thus extensively monitored for pollutants and declines in quantity. The organisms living in groundwater, however, are rarely the target of surveillance programmes and little is known about the fauna inhabiting underground habitats. The difficulties accessing groundwater, the lack of expertise, and the apparent scarcity of these organisms challenge sampling and prohibit adequate knowledge on groundwater fauna. Environmental DNA (eDNA) metabarcoding provides an approach to overcome these limitations but is largely unexplored. Here, we sampled water in 20 communal spring catchment boxes used for drinking water provisioning in Switzerland, with a high level of replication at both filtration and amplification steps. We sequenced a portion of the COI mitochondrial gene, which resulted in 4917 ASVs, yet only 3% of the reads could be assigned to a species, genus, or family with more than 90% identity. Careful evaluation of the unassigned reads corroborated that these sequences were true COI sequences belonging mostly to diverse eukaryotic groups, not present in the reference databases. Principal component analyses showed a strong correlation of the community composition with the surface land-use (agriculture vs. forest) and geology (fissured rock vs. unconsolidated sediment). While incomplete reference databases limit the assignment of taxa in groundwater eDNA metabarcoding, we showed that taxonomy-free approaches can reveal large hidden diversity and couple it with major land-use drivers, revealing their imprint on chemical and biological properties of groundwater.

An Environmental DNA Primer for Microbial and Restoration Ecology

Environmental DNA (eDNA) sequencing-DNA collected from the environment from living cells or shed DNA-was first developed for working with microbes and has greatly benefitted microbial ecologists for decades since. These tools have only become increasingly powerful with the advent of metabarcoding and metagenomics. Most new studies that examine diverse assemblages of bacteria, archaea, protists, fungi, and viruses lean heavily into eDNA using these newer technologies, as the necessary sequencing technology and bioinformatic tools have become increasingly affordable and user friendly. However, eDNA methods are rapidly evolving, and sometimes it can feel overwhelming to simply keep up with the basics. In this review, we provide a starting point for microbial ecologists who are new to DNA-based methods by detailing the eDNA methods that are most pertinent, including study design, sample collection and storage, selecting the right sequencing technology, lab protocols, equipment, and a few bioinformatic tools. Furthermore, we focus on how eDNA work can benefit restoration and what modifications are needed when working in this subfield.

Beyond the limits of the unassigned protist microbiome: inferring large-scale spatio-temporal patterns of Syndiniales marine parasites

Marine protists are major components of the oceanic microbiome that remain largely unrepresented in culture collections and genomic reference databases. The exploration of this uncharted protist diversity in oceanic communities relies essentially on studying genetic markers from the environment as taxonomic barcodes. Here we report that across 6 large scale spatio-temporal planktonic surveys, half of the genetic barcodes remain taxonomically unassigned at the genus level, preventing a fine ecological understanding for numerous protist lineages. Among them, parasitic Syndiniales (Dinoflagellata) appear as the least described protist group. We have developed a computational workflow, integrating diverse 18S rDNA gene metabarcoding datasets, in order to infer large-scale ecological patterns at 100% similarity of the genetic marker, overcoming the limitation of taxonomic assignment. From a spatial perspective, we identified 2171 unassigned clusters, i.e., Syndiniales sequences with 100% similarity, exclusively shared between the Tropical/Subtropical Ocean and the Mediterranean Sea among all Syndiniales orders and 25 ubiquitous clusters shared within all the studied marine regions. From a temporal perspective, over 3 time-series, we highlighted 39 unassigned clusters that follow rhythmic patterns of recurrence and are the best indicators of parasite community's variation. These clusters withhold potential as ecosystem change indicators, mirroring their associated host community responses. Our results underline the importance of Syndiniales in structuring planktonic communities through space and time, raising questions regarding host-parasite association specificity and the trophic mode of persistent Syndiniales, while providing an innovative framework for prioritizing unassigned protist taxa for further description.

Diversity of protist genera in periphyton of tufa-depositing karstic river

Purpose

In aquatic ecosystems, protists play a crucial role and cover numerous ecological functions. The karstic Krka River (Croatia) is a unique hotspot for high diversity of aquatic organisms, especially protists. The main objective of the present study was to obtain a detailed overview of the protist community structure in the periphyton of the Krka River and to determine the differences in protist diversity along the river.

Methods

Protist diversity was detected by amplicon sequencing of the hypervariable region V9 of the 18S rRNA gene, using the universal eukaryotic primer pair.

Results

The three main groups of protists were as follows: Ciliophora, Cercozoa, and Bacillariophyta. In terms of abundance of protist OTUs, the shade plot revealed an evident difference from the upstream to downstream river section, which increased between locations from Krka spring to Skradinski buk. Diversity was explored using measures of alpha and beta diversity. Alpha diversity showed an increasing trend in the downstream direction of the river. The location effect, or clustering/grouping of samples by location, was confirmed by the PERMANOVA permutation test of beta diversity.

Conclusion

The combination of alpha and beta diversity can help provide deeper insight into the study of diversity patterns, but also point out to decline in species diversity and allow for effective ways to protect aquatic karst habitats in future management.

Metabarcoding of the Three Domains of Life in Aquatic Saline Ecosystems

High-throughput amplicon sequencing, known as metabarcoding, is a powerful technique to decipher exhaustive microbial diversity considering specific gene markers. While most of the studies investigating ecosystem functioning through microbial diversity targeted only one domain of life, either bacteria, or archaea or microeukaryotes, the remaining challenge in microbial ecology is to uncover the integrated view of microbial diversity occurring in ecosystems. Indeed, interactions occurring between the different microbial counterparts are now recognized having a great impact on stability and resilience of ecosystems. Here, we summarize protocols describing sampling, molecular, and simultaneous metabarcoding of bacteria, archaea, and microeukaryotes, as well as a bioinformatic pipeline allowing the study of exhaustive microbial diversity in natural aquatic saline samples.

metaPR2 : A database of eukaryotic 18S rRNA metabarcodes with an emphasis on protists

In recent years, metabarcoding has become the method of choice for investigating the composition and assembly of microbial eukaryotic communities. The number of environmental data sets published has increased very rapidly. Although unprocessed sequence files are often publicly available, processed data, in particular clustered sequences, are rarely available in a usable format. Clustered sequences are reported as operational taxonomic units (OTUs) with different similarity levels or more recently as amplicon sequence variants (ASVs). This hampers comparative studies between different environments and data sets, for example examining the biogeographical patterns of specific groups/species, as well analysing the genetic microdiversity within these groups. Here, we present a newly-assembled database of processed 18S rRNA metabarcodes that are annotated with the PR² reference sequence database. This database, called metaPR² , contains 41 data sets corresponding to more than 4000 samples and 90,000 ASVs. The database, which is accessible through both a web-based interface (https://shiny.metapr2.org) and an R package, should prove very useful to all researchers working on protist diversity in a variety of systems.

Estuarine microbial networks and relationships vary between environmentally distinct communities

Microbial interactions have profound impacts on biodiversity, biogeochemistry, and ecosystem functioning, and yet, they remain poorly understood in the ocean and with respect to changing environmental conditions. We applied hierarchical clustering of an annual 16S and 18S amplicon dataset in the Skidaway River Estuary, which revealed two similar clusters for prokaryotes (Bacteria and Archaea) and protists: Cluster 1 (March-May and November-February) and Cluster 2 (June-October). We constructed co-occurrence networks from each cluster to explore how microbial networks and relationships vary between environmentally distinct periods in the estuary. Cluster 1 communities were exposed to significantly lower temperature, sunlight, NO₃, and SiO₄; only NH₄ was higher at this time. Several network properties (e.g., edge number, degree, and centrality) were elevated for networks constructed with Cluster 1 vs. 2 samples. There was also evidence that microbial nodes in Cluster 1 were more connected (e.g., higher edge density and lower path length) compared to Cluster 2, though opposite trends were observed when networks considered Prokaryote-Protist edges only. The number of Prokaryote-Prokaryote and Prokaryote-Protist edges increased by >100% in the Cluster 1 network, mainly involving Flavobacteriales, Rhodobacterales, Peridiniales, and Cryptomonadales associated with each other and other microbial groups (e.g., SAR11, Bacillariophyta, and Strombidiida). Several Protist-Protist associations, including Bacillariophyta correlated with Syndiniales (Dino-Groups I and II) and an Unassigned Dinophyceae group, were more prevalent in Cluster 2. Based on the type and sign of associations that increased in Cluster 1, our findings indicate that mutualistic, competitive, or predatory relationships may have been more representative among microbes when conditions were less favorable in the estuary; however, such relationships require further exploration and validation in the field and lab. Coastal networks may also be driven by shifts in the abundance of certain taxonomic or functional groups. Sustained monitoring of microbial communities over environmental gradients, both spatial and temporal, is critical to predict microbial dynamics and biogeochemistry in future marine ecosystems.

Defensive symbiosis against giant viruses in amoebae

Protists are important regulators of microbial communities and key components in food webs with impact on nutrient cycling and ecosystem functioning. In turn, their activity is shaped by diverse intracellular parasites, including bacterial symbionts and viruses. Yet, bacteria-virus interactions within protists are poorly understood. Here, we studied the role of bacterial symbionts of free-living amoebae in the establishment of infections with nucleocytoplasmic large DNA viruses (Nucleocytoviricota). To investigate these interactions in a system that would also be relevant in nature, we first isolated and characterized a giant virus (Viennavirus, family Marseilleviridae) and a sympatric potential Acanthamoeba host infected with bacterial symbionts. Subsequently, coinfection experiments were carried out, using the fresh environmental isolates as well as additional amoeba laboratory strains. Employing fluorescence in situ hybridization and qPCR, we show that the bacterial symbiont, identified as Parachlamydia acanthamoebae, represses the replication of the sympatric Viennavirus in both recent environmental isolates as well as Acanthamoeba laboratory strains. In the presence of the symbiont, virions are still taken up, but viral factory maturation is inhibited, leading to survival of the amoeba host. The symbiont also suppressed the replication of the more complex Acanthamoeba polyphaga mimivirus and Tupanvirus deep ocean (Mimiviridae). Our work provides an example of an intracellular bacterial symbiont protecting a protist host against virus infections. The impact of virus-symbiont interactions on microbial population dynamics and eventually ecosystem processes requires further attention.

DNA sequence and taxonomic gap analyses to quantify the coverage of aquatic cyanobacteria and eukaryotic microalgae in reference databases: Results of a survey in the Alpine region

The taxonomic identification of organisms based on the amplification of specific genetic markers (metabarcoding) implicitly requires adequate discriminatory information and taxonomic coverage of environmental DNA sequences in taxonomic databases. These requirements were quantitatively examined by comparing the determination of cyanobacteria and microalgae obtained by metabarcoding and light microscopy. We used planktic and biofilm samples collected in 37 lakes and 22 rivers across the Alpine region. We focused on two of the most used and best represented genetic markers in the reference databases, namely the 16S rRNA and 18S rRNA genes. A sequence gap analysis using blastn showed that, in the identity range of 99-100%, approximately 30% (plankton) and 60% (biofilm) of the sequences did not find any close counterpart in the reference databases (NCBI GenBank). Similarly, a taxonomic gap analysis showed that approximately 50% of the cyanobacterial and eukaryotic microalgal species identified by light microscopy were not represented in the reference databases. In both cases, the magnitude of the gaps differed between the major taxonomic groups. Even considering the species determined under the microscope and represented in the reference databases, 22% and 26% were still not included in the results obtained by the blastn at percentage levels of identity ≥95% and ≥97%, respectively. The main causes were the absence of matching sequences due to amplification and/or sequencing failure and potential misidentification in the microscopy step. Our results quantitatively demonstrated that in metabarcoding the main obstacles in the classification of 16S rRNA and 18S rRNA sequences and interpretation of high-throughput sequencing biomonitoring data were due to the existence of important gaps in the taxonomic completeness of the reference databases and the short length of reads. The study focused on the Alpine region, but the extent of the gaps could be much greater in other less investigated geographic areas.

Composition and spatial-temporal dynamics of phytoplankton community shaped by environmental selection and interactions in the Jiaozhou Bay

The Jiaozhou Bay as a model marine ecosystem in China has been intensively investigated over the last 90 years. However, detailed phytoplankton community composition, spatial-temporal dynamics, and its assembly mechanism were still unclear. To address these, we systematically examined the composition and spatial-temporal dynamics of phytoplankton in the Jiaozhou Bay through high-throughput sequencing of 18S rDNA V4. Analysis of 468 samples from 12 sampling sites over one full year revealed much higher phytoplankton diversity than previous reports, and strong seasonal succession patterns. Some phytoplankton also showed spatial variations, although the phytoplankton community didn't show significant distance-decay pattern. Environmental factors (especially temperature), species-species interrelationships and unique resting stages were uncovered to be the main contributors instead of stochastic process in shaping the phytoplankton community assembly. The overwhelming positive correlations between phytoplankton and other protists suggested that coevolution might be critical in this marine ecosystem. Complementary distributions of different amplicon sequence variants (ASVs) of same genera, such as Skeletonema marinoi (ASV_2) and Skeletonema tropicum (ASV_263) of the genus Skeletonema, suggested that phytoplankton have evolved differentially to exploit a wide range of ecological niches. This study laid a solid foundation for asertaining phytoplankton composition and spatial-temporal dynamics in temperate seawaters and mechanisms underlying phytoplankton community assembly, allowing in-depth studies of marine ecology.

Effects of operational taxonomic unit inference methods on soil microeukaryote community analysis using long-read metabarcoding

Long amplicon metabarcoding has opened the door for phylogenetic analysis of the largely unknown communities of microeukaryotes in soil. Here, we amplified and sequenced the ITS and LSU regions of the rDNA operon (around 1500 bp) from grassland soils using PacBio SMRT sequencing. We tested how three different methods for generation of operational taxonomic units (OTUs) effected estimated richness and identified taxa, and how well large-scale ecological patterns associated with shifting environmental conditions were recovered in data from the three methods. The field site at Kungsängen Nature Reserve has drawn frequent visitors since Linnaeus's time, and its species rich vegetation includes the largest population of Fritillaria meleagris in Sweden. To test the effect of different OTU generation methods, we sampled soils across an abrupt moisture transition that divides the meadow community into a Carex acuta dominated plant community with low species richness in the wetter part, which is visually distinct from the mesic-dry part that has a species rich grass-dominated plant community including a high frequency of F. meleagris. We used the moisture and plant community transition as a framework to investigate how detected belowground microeukaryotic community composition was influenced by OTU generation methods. Soil communities in both moisture regimes were dominated by protists, a large fraction of which were taxonomically assigned to Ciliophora (Alveolata) while 30%-40% of all reads were assigned to kingdom Fungi. Ecological patterns were consistently recovered irrespective of OTU generation method used. However, different methods strongly affect richness estimates and the taxonomic and phylogenetic resolution of the characterized community with implications for how well members of the microeukaryotic communities can be recognized in the data.

Evolving perspective on the origin and diversification of cellular life and the virosphere

The tree of life (TOL) is a powerful framework to depict the evolutionary history of cellular organisms through time, from our microbial origins to the diversification of multicellular eukaryotes that shape the visible biosphere today. During the past decades, our perception of the TOL has fundamentally changed, in part, due to profound methodological advances, which allowed a more objective approach to studying organismal and viral diversity and led to the discovery of major new branches in the TOL as well as viral lineages. Phylogenetic and comparative genomics analyses of these data have, among others, revolutionized our understanding of the deep roots and diversity of microbial life, the origin of the eukaryotic cell, eukaryotic diversity, as well as the origin, and diversification of viruses. In this review, we provide an overview of some of the recent discoveries on the evolutionary history of cellular organisms and their viruses and discuss a variety of complementary techniques that we consider crucial for making further progress in our understanding of the TOL and its interconnection with the virosphere.

Size-Fractionated Microbiome Structure in Subarctic Rivers and a Coastal Plume Across DOC and Salinity Gradients

Little is known about the microbial diversity of rivers that flow across the changing subarctic landscape. Using amplicon sequencing (rRNA and rRNA genes) combined with HPLC pigment analysis and physicochemical measurements, we investigated the diversity of two size fractions of planktonic Bacteria, Archaea and microbial eukaryotes along environmental gradients in the Great Whale River (GWR), Canada. This large subarctic river drains an extensive watershed that includes areas of thawing permafrost, and discharges into southeastern Hudson Bay as an extensive plume that gradually mixes with the coastal marine waters. The microbial communities differed by size-fraction (separated with a 3-μm filter), and clustered into three distinct environmental groups: (1) the GWR sites throughout a 150-km sampling transect; (2) the GWR plume in Hudson Bay; and (3) small rivers that flow through degraded permafrost landscapes. There was a downstream increase in taxonomic richness along the GWR, suggesting that sub-catchment inputs influence microbial community structure in the absence of sharp environmental gradients. Microbial community structure shifted across the salinity gradient within the plume, with changes in taxonomic composition and diversity. Rivers flowing through degraded permafrost had distinct physicochemical and microbiome characteristics, with allochthonous dissolved organic carbon explaining part of the variation in community structure. Finally, our analyses of the core microbiome indicated that while a substantial part of all communities consisted of generalists, most taxa had a more limited environmental range and may therefore be sensitive to ongoing change.

Biodiversity and Spatial-Temporal Dynamics of Margalefidinium Species in Jiaozhou Bay, China

Many Margalefidinium species are cosmopolitan harmful algal bloom (HAB) species that have caused huge economic and ecological damage. Despite extensive research on Margalefidinium species, the biodiversity and spatial-temporal dynamics of these species remain obscure. Jiaozhou Bay is an ideal area for HAB research, being one of the earliest marine survey areas in China. In this study, we carried out the first metabarcoding study on the temporal and spatial dynamics of Margalefidinium species using the 18S rDNA V4 region as the molecular marker and samples collected monthly at 12 sampling sites in Jiaozhou Bay in 2019. Two harmful Margalefidinium species (M. polykrikoides and M. fulvescens) were identified with potentially high genetic diversity (although we cannot rule out the possibility of intra-genome sequence variations). Both M. polykrikoides and M. fulvescens demonstrated strong temporal preference with a sharp peak of abundance in early autumn (September), but without showing strong location preference in Jiaozhou Bay. Our results revealed that temperature might be the main driver for their temporal dynamics. Knowledge of biodiversity and spatial-temporal dynamics of the Margalefidinium species may shed light on the understanding of mechanisms underlying strongly biased occurrences of Margalefidinium blooms recorded globally.