Probing the evolution, ecology and physiology of marine protists using transcriptomics

Distinct structure, assembly, and gene expression of microplankton in two Arctic estuaries with varied terrestrial inputs

The Laptev Sea is a major Marginal Sea in the Western Arctic Ocean. The Arctic amplification brought by global warming influences the hydrological properties of rivers passing through the permafrost zone, which would alter the biological community structure at continental margin. In this study, the structure, assembly, and gene expression of planktonic microbial communities in two estuaries (Protoka Ularovskaya River Estuary, PURE; Lena River Estuary, LRE) of Laptev Sea were examined to investigate the environmental effects of polar rivers. PURE and LRE exhibited distinct environmental characteristics: low temperature and high salinity for PURE, and high temperature and low salinity for LRE, influenced by runoff size. Salinity more closely influenced microbial communities in LRE, with freshwater species playing a significant role in community composition. The findings revealed differences between two estuaries in community composition and diversity. Prokaryotes and microeukaryotes had shown different assembly patterns in response to habitat changes caused by terrestrial freshwater input. Furthermore, compared with the PURE, the co-occurrence and inter-domain network of the LRE, which was more affected by terrestrial input, was more complex and stable. Functional gene prediction revealed a higher gene expression of methane metabolism in LRE than in PURE, particularly those related to methane oxidation, and this conclusion could help better explore the impact of global warming on the methane cycle in the Arctic Marginal Seas. This study explored the increased freshwater runoffs under the background of global warming dramatically affect Arctic microplankton communities from community structure, assembly and gene expression aspects.

Time-series metagenomics reveals changing protistan ecology of a temperate dimictic lake

BACKGROUND

Protists, single-celled eukaryotic organisms, are critical to food web ecology, contributing to primary productivity and connecting small bacteria and archaea to higher trophic levels. Lake Mendota is a large, eutrophic natural lake that is a Long-Term Ecological Research site and among the world's best-studied freshwater systems. Metagenomic samples have been collected and shotgun sequenced from Lake Mendota for the last 20 years. Here, we analyze this comprehensive time series to infer changes to the structure and function of the protistan community and to hypothesize about their interactions with bacteria.

RESULTS

Based on small subunit rRNA genes extracted from the metagenomes and metagenome-assembled genomes of microeukaryotes, we identify shifts in the eukaryotic phytoplankton community over time, which we predict to be a consequence of reduced zooplankton grazing pressures after the invasion of a invasive predator (the spiny water flea) to the lake. The metagenomic data also reveal the presence of the spiny water flea and the zebra mussel, a second invasive species to Lake Mendota, prior to their visual identification during routine monitoring. Furthermore, we use species co-occurrence and co-abundance analysis to connect the protistan community with bacterial taxa. Correlation analysis suggests that protists and bacteria may interact or respond similarly to environmental conditions. Cryptophytes declined in the second decade of the timeseries, while many alveolate groups (e.g., ciliates and dinoflagellates) and diatoms increased in abundance, changes that have implications for food web efficiency in Lake Mendota.

CONCLUSIONS

We demonstrate that metagenomic sequence-based community analysis can complement existing efforts to monitor protists in Lake Mendota based on microscopy-based count surveys. We observed patterns of seasonal abundance in microeukaryotes in Lake Mendota that corroborated expectations from other systems, including high abundance of cryptophytes in winter and diatoms in fall and spring, but with much higher resolution than previous surveys. Our study identified long-term changes in the abundance of eukaryotic microbes and provided context for the known establishment of an invasive species that catalyzes a trophic cascade involving protists. Our findings are important for decoding potential long-term consequences of human interventions, including invasive species introduction. Video Abstract.

Seasonal metabolic dynamics of microeukaryotic plankton: a year-long metatranscriptomic study in a temperate sea

Seasonal fluctuations profoundly affect marine microeukaryotic plankton composition and metabolism, but accurately tracking these changes has been a long-standing challenge. In this study, we present a year-long metatranscriptomic data set from the Southern Bight of the North Sea, shedding light on the seasonal dynamics in temperate plankton ecosystems. We observe distinct shifts in active plankton species and their metabolic processes in response to seasonal changes. We characterized the metabolic signatures of different seasonal phases in detail, thereby revealing the metabolic versatility of dinoflagellates, the heterotrophic dietary strategy of Phaeocystis during its late-stage blooms, and stark variations in summer and fall diatom abundance and metabolic activity across nearby sampling stations. Our data illuminate the varied contributions of microeukaryotic taxa to biomass production and nutrient cycling at different times of the year and allow delineation of their ecological niches.

IMPORTANCE

Ecosystem composition and metabolic functions of temperate marine microeukaryote plankton are strongly influenced by seasonal dynamics. Although monitoring of species composition of microeukaryotes has expanded recently, few methods also contain seasonally resolved information on ecosystem functioning. We generated a year-long spatially resolved metatranscriptomic data set to assess seasonal dynamics of microeukaryote species and their associated metabolic functions in the Southern Bight of the North Sea. Our study underscores the potential of metatranscriptomics as a powerful tool for advancing our understanding of marine ecosystem functionality and resilience in response to environmental changes, emphasizing its potential in continuous marine ecosystem monitoring to enhance our ecological understanding of the ocean's eukaryotic microbiome.

A brief history of metal recruitment in protozoan predation

Metals and metalloids are used as weapons for predatory feeding by unicellular eukaryotes on prokaryotes. This review emphasizes the role of metal(loid) bioavailability over the course of Earth's history, coupled with eukaryogenesis and the evolution of the mitochondrion to trace the emergence and use of the metal(loid) prey-killing phagosome as a feeding strategy. Members of the genera Acanthamoeba and Dictyostelium use metals such as zinc (Zn) and copper (Cu), and possibly metalloids, to kill their bacterial prey after phagocytosis. We provide a potential timeline on when these capacities first evolved and how they correlate with perceived changes in metal(loid) bioavailability through Earth's history. The origin of phagotrophic eukaryotes must have postdated the Great Oxidation Event (GOE) in agreement with redox-dependent modification of metal(loid) bioavailability for phagotrophic poisoning. However, this predatory mechanism is predicted to have evolved much later - closer to the origin of the multicellular metazoans and the evolutionary development of the immune systems.

Biogeographical and biodiversity patterns of planktonic microeukaryotes along the tropical western to eastern Pacific Ocean transect revealed by metabarcoding

Microeukaryotic plankton (0.2-200 µm), which are morphologically and genetically highly diverse, play a crucial role in ocean productivity and carbon consumption. The Pacific Ocean (PO), one of the world's largest oligotrophic regions, remains largely unexplored in terms of the biogeography and biodiversity of microeukaryotes based on large-scale sampling. We investigated the horizontal distribution of microeukaryotes along a 16,000 km transect from the west to the east of the PO. The alpha diversity indices showed a distinct decreasing trend from west to east, which was highly correlated with water temperature. The microeukaryotic community, which was clustered into the western, central, and eastern PO groups, displayed a significant distance-decay relationship. Syndiniales, a lineage of parasitic dinoflagellates, was ubiquitously distributed along the transect and dominated the community in terms of both sequence and zero-radius operational taxonomic unit (ZOTU) proportions. The prevailing dominance of Syndiniales-affiliated ZOTUs and their close associations with dinoflagellates, diatoms, and radiolarians, as revealed by SparCC correlation analysis, suggested that parasitism may be an important trophic strategy in the surface waters of the PO. Geographical distance and temperature were the most important environmental factors that significantly correlated with community structure. Overall, our study sheds more light on the distribution pattern of both alpha and beta diversities of microeukaryotic communities and highlighted the importance of parasitisms by Syndiniales across the tropical PO.IMPORTANCEUnderstanding the biogeographical and biodiversity patterns of microeukaryotic communities is essential to comprehending their roles in biogeochemical cycling. In this study, planktonic microeukaryotes were collected along a west-to-east Pacific Ocean transect (ca. 16,000 km). Our study revealed that the alpha diversity indices were highly correlated with water temperature, and the microeukaryotic communities displayed a distinct geographical distance-driven pattern. The predominance of the parasitic dinoflagellate lineage Syndiniales and their close relationship with other microeukaryotic groups suggest that parasitism may be a crucial survival strategy for microeukaryotes in the surface waters of the Pacific Ocean. Our findings expand our understanding of the biodiversity and biogeographical pattern of microeukaryotes and highlight the significance of parasitic Syndiniales in the surface ocean.

Succession of diversity, assembly mechanisms, and activities of the microeukaryotic community throughout Scrippsiella acuminata (Dinophyceae) bloom phases

Harmful algal bloom (HAB) is a rapidly expanding marine ecological hazard. Although numerous studies have been carried out about the ecological impact and the ecological mechanism of HAB outbreaks, few studies have comprehensively addressed the shifts of species composition, metabolic activity level, driving factors and community assembly mechanisms of microeukaryotic plankton in the course of the bloom event. To fill the gap of research, we conducted 18S ribosomal DNA and RNA sequencing during the initiation, development, sustenance and decline stages of a Scrippsiella acuminata (S. acuminata) bloom at the coastal sea of Fujian Province, China. We found that the bloom event caused a decrease in microeukaryotic plankton species diversity and increase in community homogeneity. Our results revealed that the RNA- and DNA-inferred communities were similar, but α-diversity was more dynamic in RNA- than in DNA-inferred communities. The main taxa with high projected metabolic activity (with RNA:DNA ratio as the proxy) during the bloom included dinoflagellates, Cercozoa, Chlorophyta, Protalveolata, and diatoms. The role of deterministic processes in microeukaryotic plankton community assembly increased during the bloom, but stochastic processes were always the dominant assembly mechanism throughout the bloom process. Our findings improve the understanding of temporal patterns, driving factors and assembly mechanisms underlying the microeukarytic plankton community in a dinoflagellate bloom.

Insight into planktonic protistan and fungal communities across the nutrient-depleted environment of the South Pacific Subtropical Gyre

Ocean microorganisms constitute ~70% of the marine biomass, contribute to ~50% of the Earth's primary production, and play a vital role in global biogeochemical cycles. The marine heterotrophic and mixotrophic protistan and fungal communities have often been overlooked mainly due to limitations in morphological species identification. Despite the accumulation of studies on biogeographic patterns observed in microbial communities, our understanding of the abundance and distribution patterns within the microbial community of the largest subtropical gyre, the South Pacific Gyre (SPG), remains incomplete. Here, we investigated the diversity and vertical composition of protistan and fungal communities in the water column of the ultra-oligotrophic SPG. Our results showed apparent differences in protistan community diversity in the photic and aphotic regions. The entire protistan community diversity was significantly affected by temperature, salinity, oxygen, and nutrient concentrations, while the parasitic community diversity was also affected by chlorophyll a concentration. The parasitic protists were assigned to the class Syndiniales accounting for over 98% of the total parasitic protists, exhibiting higher relative sequence abundance along the water depth and displaying consistent patterns among different sampling stations. In contrast to the protistan community, the fungal community along the SPG primarily clustered based on the sampling station and pelagic zones. In particular, our study reveals a significant presence of parasitic protists and functionally diverse fungi in SPG and their potential impact on carbon cycling in the gyre.IMPORTANCEOur findings carry important implications for understanding the distribution patterns of the previously unrecognized occurrence of parasitic protists and functionally diverse fungi in the nutrient-limited South Pacific Gyre. In particular, our study reveals a significant presence of parasitic Syndiniales, predominantly abundant in the upper 300 m of the aphotic zone in the gyre, and a distinct presence of fungal communities in the aphotic zone at the central part of the gyre. These findings strongly suggest that these communities play a substantial role in yet insufficiently described microbial food web. Moreover, our research enhances our understanding of their contribution to the dynamics of the food webs in oligotrophic gyres and is valuable for projecting the ecological consequences of future climate warming.

Chapter 6: The Breadth and Limits of Life on Earth

Scientific ideas about the potential existence of life elsewhere in the universe are predominantly informed by knowledge about life on Earth. Over the past ∼4 billion years, life on Earth has evolved into millions of unique species. Life now inhabits nearly every environmental niche on Earth that has been explored. Despite the wide variety of species and diverse biochemistry of modern life, many features, such as energy production mechanisms and nutrient requirements, are conserved across the Tree of Life. Such conserved features help define the operational parameters required by life and therefore help direct the exploration and evaluation of habitability in extraterrestrial environments. As new diversity in the Tree of Life continues to expand, so do the known limits of life on Earth and the range of environments considered habitable elsewhere. The metabolic processes used by organisms living on the edge of habitability provide insights into the types of environments that would be most suitable to hosting extraterrestrial life, crucial for planning and developing future astrobiology missions. This chapter will introduce readers to the breadth and limits of life on Earth and show how the study of life at the extremes can inform the broader field of astrobiology.

Chapter 1: The Astrobiology Primer 3.0

The Astrobiology Primer 3.0 (ABP3.0) is a concise introduction to the field of astrobiology for students and others who are new to the field of astrobiology. It provides an entry into the broader materials in this supplementary issue of Astrobiology and an overview of the investigations and driving hypotheses that make up this interdisciplinary field. The content of this chapter was adapted from the other 10 articles in this supplementary issue and thus represents the contribution of all the authors who worked on these introductory articles. The content of this chapter is not exhaustive and represents the topics that the authors found to be the most important and compelling in a dynamic and changing field.

Responses of microeukaryotic community structure to a Phaeocystis globosa bloom in a semi-enclosed subtropical bay

The occurrence of Phaeocystis globosa, a harmful algal bloom species in Chinese coastal waters, has significant impacts on marine organisms and poses a threat to the safety of coastal nuclear power plants. Although previous studies have established a close association between P. globosa blooms and the bacterial community, the relationship between the microeukaryotic community and P. globosa blooms remains poorly understood. In this study, the variations in the microeukaryotic community resulting from a P. globosa bloom were analyzed using 18S rRNA gene amplicon sequencing. The results indicated that the diversity of the microeukaryotic community during the bloom phase was significantly higher than that during the dissipation phase. The microeukaryotic community compositions varied significantly between the two phases of the P. globosa bloom. During the bloom phase, the dominant microeukaryotic was Viridiplantae, which was then replaced by Dinoflagellata during the dissipation phase. Co-occurrence network analysis showed that the relationship among the microeukaryotic community during the bloom phase was more complex than that during the dissipation phase, and the keystone taxa varied as the bloom progressed. Additionally, microeukaryotic community assembly was primarily driven by stochastic processes during the bloom phase based on the β-nearest taxon distance, whereas it was driven by both deterministic processes and stochastic processes during the dissipation phase. Overall, our findings provide novel insight into the mechanisms and interactions involved in microeukaryotic community dynamics in environments disturbed by P. globosa blooms.

Novel Plastid Genome Characteristics in Fugacium kawagutii and the Trend of Accelerated Evolution of Plastid Proteins in Dinoflagellates

Typical (peridinin-containing) dinoflagellates possess plastid genomes composed of small plasmids named "minicircles". Despite the ecological importance of dinoflagellate photosynthesis in corals and marine ecosystems, the structural characteristics, replication dynamics, and evolutionary forcing of dinoflagellate plastid genomes remain poorly understood. Here, we sequenced the plastid genome of the symbiodiniacean species Fugacium kawagutii and conducted comparative analyses. We identified psbT-coding minicircles, features previously not found in Symbiodiniaceae. The copy number of F. kawagutii minicircles showed a strong diel dynamics, changing between 3.89 and 34.3 copies/cell and peaking in mid-light period. We found that F. kawagutii minicircles are the shortest among all dinoflagellates examined to date. Besides, the core regions of the minicircles are highly conserved within genus in Symbiodiniaceae. Furthermore, the codon usage bias of the plastid genomes in Heterocapsaceae, Amphidiniaceae, and Prorocentraceae species are greatly influenced by selection pressure, and in Pyrocystaceae, Symbiodiniaceae, Peridiniaceae, and Ceratiaceae species are influenced by both natural selection pressure and mutation pressure, indicating a family-level distinction in codon usage evolution in dinoflagellates. Phylogenetic analysis using 12 plastid-encoded proteins and five nucleus-encoded plastid proteins revealed accelerated evolution trend of both plastid- and nucleus-encoded plastid proteins in peridinin- and fucoxanthin-dinoflagellate plastids compared to plastid proteins of nondinoflagellate algae. These findings shed new light on the structure and evolution of plastid genomes in dinoflagellates, which will facilitate further studies on the evolutionary forcing and function of the diverse dinoflagellate plastids. The accelerated evolution documented here suggests plastid-encoded sequences are potentially useful for resolving closely related dinoflagellates.

Microeukaryotic plankton evolutionary constraints in a subtropical river explained by environment and bacteria along differing taxonomic resolutions

Microeukaryotic plankton communities are keystone components for keeping aquatic primary productivity. Currently, variations in microeukaryotic plankton diversity have often been explained by local ecological factors but not by evolutionary constraints. We used amplicon sequencing of 100 water samples across five years to investigate the ecological preferences of the microeukaryotic plankton community in a subtropical riverine ecosystem. We found that microeukaryotic plankton diversity was less associated with bacterial abundance (16S rRNA gene copy number) than bacterial diversity. Further, environmental effects exhibited a larger influence on microeukaryotic plankton community composition than bacterial community composition, especially at fine taxonomic levels. The evolutionary constraints of microeukaryotic plankton community increased with decreasing taxonomic resolution (from 97% to 91% similarity levels), but not significant change from 85% to 70% similarity levels. However, compared with the bacterial community, the evolutionary constraints were shown to be more affected by environmental variables. This study illustrated possible controlling environmental and bacterial drivers of microeukaryotic diversity and community assembly in a subtropical river, thereby indirectly reflecting on the quality status of the water environment by providing new clues on the microeukaryotic community assembly.

MarFERReT, an open-source, version-controlled reference library of marine microbial eukaryote functional genes

Metatranscriptomics generates large volumes of sequence data about transcribed genes in natural environments. Taxonomic annotation of these datasets depends on availability of curated reference sequences. For marine microbial eukaryotes, current reference libraries are limited by gaps in sequenced organism diversity and barriers to updating libraries with new sequence data, resulting in taxonomic annotation of about half of eukaryotic environmental transcripts. Here, we introduce Marine Functional EukaRyotic Reference Taxa (MarFERReT), a marine microbial eukaryotic sequence library designed for use with taxonomic annotation of eukaryotic metatranscriptomes. We gathered 902 publicly accessible marine eukaryote genomes and transcriptomes and assessed their sequence quality and cross-contamination issues, selecting 800 validated entries for inclusion in MarFERReT. Version 1.1 of MarFERReT contains reference sequences from 800 marine eukaryotic genomes and transcriptomes, covering 453 species- and strain-level taxa, totaling nearly 28 million protein sequences with associated NCBI and PR2 Taxonomy identifiers and Pfam functional annotations. The MarFERReT project repository hosts containerized build scripts, documentation on installation and use case examples, and information on new versions of MarFERReT.

Monitoring potentially pathogenic protists in sewage sludge using Metataxonomics

Intestinal parasites continue to pose a significant threat to human health worldwide, particularly among children. Contaminated water and soil serve as major transmission vehicles for these parasites and intestinal protists are among the most prevalent parasites in both developed and developing nations. Traditionally, parasites have been studied using human or animal fecal samples, while studying them in environmental samples has been challenging due to technical limitations. However, advancements in Next-Generation Sequencing (NGS) and bioinformatic approaches now enable the detection of parasite DNA in environmental samples. In this study, we applied a metataxonomic and phylogenetic strategy to detect and classify DNA of protists present in sewage sludge from two major cities in Colombia: Medellin and Cali. We successfully detected several human pathogenic parasites including Giardia intestinalis, Entamoeba histolytica, and Blastocystis sp., among other protists, in all sludge samples examined. We also investigated the entry and exit of parasite DNA from the San Fernando wastewater treatment plant (WWTP). We observed a higher number of parasite DNA sequences in the plant's influent wastewater, but we also detected the discharge of DNA from pathogenic parasites in both effluent waters and biosolids.

Pesticides in the Great Barrier Reef catchment area: Plausible risks to fish populations

Waterways that drain the Great Barrier Reef catchment area (GBRCA) transport pollutants to marine habitats, provide a critical corridor between freshwater and marine habitats for migratory fish species, and are of high socioecological value. Some of these waterways contain concentrations of pesticide active ingredients (PAIs) that exceed Australian ecotoxicity threshold values (ETVs) for ecosystem protection. In this article, we use a "pathway to harm" model with five key criteria to assess whether the available information supports the hypothesis that PAIs are or could have harmful effects on fish and arthropod populations. Strong evidence of the first three criteria and circumstantial weaker evidence of the fourth and fifth criteria are presented. Specifically, we demonstrate that exceedances of Australian and New Zealand ETVs for ecosystem protection are widespread in the GBRCA, that the PAI contaminated water occurs (spatially and temporally) in important habitats for fisheries, and that there are clear direct and indirect mechanisms by which PAIs could cause harmful effects. The evidence of individuals and populations of fish and arthropods being adversely affected species is more circumstantial but consistent with PAIs causing harmful effects in the freshwater ecosystems of Great Barrier Reef waterways. We advocate strengthening the links between PAI concentrations and fish health because of the cultural values placed on the freshwater ecosystems by relevant stakeholders and Traditional Owners, with the aim that stronger links between elevated PAI concentrations and changes in recreationally and culturally important fish species will inspire improvements in water quality. Integr Environ Assess Manag 2024;00:1-24. © 2023 Commonwealth of Australia and The Commonwealth Scientific and Industrial Research Organisation. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Eukaryotic genomes from a global metagenomic data set illuminate trophic modes and biogeography of ocean plankton

Metagenomics is a powerful method for interpreting the ecological roles and physiological capabilities of mixed microbial communities. Yet, many tools for processing metagenomic data are neither designed to consider eukaryotes nor are they built for an increasing amount of sequence data. EukHeist is an automated pipeline to retrieve eukaryotic and prokaryotic metagenome-assembled genomes (MAGs) from large-scale metagenomic sequence data sets. We developed the EukHeist workflow to specifically process large amounts of both metagenomic and/or metatranscriptomic sequence data in an automated and reproducible fashion. Here, we applied EukHeist to the large-size fraction data (0.8-2,000 µm) from Tara Oceans to recover both eukaryotic and prokaryotic MAGs, which we refer to as TOPAZ (Tara Oceans Particle-Associated MAGs). The TOPAZ MAGs consisted of >900 environmentally relevant eukaryotic MAGs and >4,000 bacterial and archaeal MAGs. The bacterial and archaeal TOPAZ MAGs expand upon the phylogenetic diversity of likely particle- and host-associated taxa. We use these MAGs to demonstrate an approach to infer the putative trophic mode of the recovered eukaryotic MAGs. We also identify ecological cohorts of co-occurring MAGs, which are driven by specific environmental factors and putative host-microbe associations. These data together add to a number of growing resources of environmentally relevant eukaryotic genomic information. Complementary and expanded databases of MAGs, such as those provided through scalable pipelines like EukHeist, stand to advance our understanding of eukaryotic diversity through increased coverage of genomic representatives across the tree of life.IMPORTANCESingle-celled eukaryotes play ecologically significant roles in the marine environment, yet fundamental questions about their biodiversity, ecological function, and interactions remain. Environmental sequencing enables researchers to document naturally occurring protistan communities, without culturing bias, yet metagenomic and metatranscriptomic sequencing approaches cannot separate individual species from communities. To more completely capture the genomic content of mixed protistan populations, we can create bins of sequences that represent the same organism (metagenome-assembled genomes [MAGs]). We developed the EukHeist pipeline, which automates the binning of population-level eukaryotic and prokaryotic genomes from metagenomic reads. We show exciting insight into what protistan communities are present and their trophic roles in the ocean. Scalable computational tools, like EukHeist, may accelerate the identification of meaningful genetic signatures from large data sets and complement researchers' efforts to leverage MAG databases for addressing ecological questions, resolving evolutionary relationships, and discovering potentially novel biodiversity.

Plastoquinone synthesis inhibition by tetrabromo biphenyldiol as a widespread algicidal mechanism of marine bacteria

Algae and bacteria have complex and intimate interactions in the ocean. Besides mutualism, bacteria have evolved a variety of molecular-based anti-algal strategies. However, limited by the unknown mechanism of synthesis and action of these molecules, these strategies and their global prevalence remain unknown. Here we identify a novel strategy through which a marine representative of the Gammaproteobacteria produced 3,3',5,5'-tetrabromo-2,2'-biphenyldiol (4-BP), that kills or inhibits diverse phytoplankton by inhibiting plastoquinone synthesis and its effect cascades to many other key metabolic processes of the algae. Through comparative genomic analysis between the 4-BP-producing bacterium and its algicidally inactive mutant, combined with gene function verification, we identified the gene cluster responsible for 4-BP synthesis, which contains genes encoding chorismate lyase, flavin-dependent halogenase and cytochrome P450. We demonstrated that in near in situ simulated algal blooming seawater, even low concentrations of 4-BP can cause changes in overall phytoplankton community structure with a decline in dinoflagellates and diatoms. Further analyses of the gene sequences from the Tara Oceans expeditions and 2750 whole genome sequences confirmed the ubiquitous presence of 4-BP synthetic genes in diverse bacterial members in the global ocean, suggesting that it is a bacterial tool potentially widely used in global oceans to mediate bacteria-algae antagonistic relationships.

Yeast metagenomics: analytical challenges in the analysis of the eukaryotic microbiome

Even if their impact is often underestimated, yeasts and yeast-like fungi represent the most prevalent eukaryotic members of microbial communities on Earth. They play numerous roles in natural ecosystems and in association with their hosts. They are involved in the food industry and pharmaceutical production, but they can also cause diseases in other organisms, making the understanding of their biology mandatory. The ongoing loss of biodiversity due to overexploitation of environmental resources is a growing concern in many countries. Therefore, it becomes crucial to understand the ecology and evolutionary history of these organisms to systematically classify them. To achieve this, it is essential that our knowledge of the mycobiota reaches a level similar to that of the bacterial communities. To overcome the existing challenges in the study of fungal communities, the first step should be the establishment of standardized techniques for the correct identification of species, even from complex matrices, both in wet lab practices and in bioinformatic tools.

Transcriptional profile reveals the physiological responses to prey availability in the mixotrophic chrysophyte Poterioochromonas malhamensis

Mixotrophic flagellates, which have diverse nutritional modes and play important roles in connecting the microbial loop with the classical food chain, are ideal models to study the mechanisms of adaptation between different nutritional modes in protists. In their natural ecosystems, mixotrophic flagellates may encounter microalgal prey of different digestibility, which may affect the carbon flow. To date, a molecular biological view of the metabolic processes in the mixotrophic flagellate Poterioochromonas malhamensis during nutritional adaptation and feeding on microalgal prey of different digestibility is still lacking. Accordingly, this study focused on the gene expression differences in P. malhamensis under autotrophy, being fed by the digestible microalga Chlorella sorokiniana GT-1, and being fed by the indigestible microalga C. sorokiniana CMBB-146. Results showed that the growth rate of P. malhamensis under autotrophy was much lower than that when fed by digestible microalgae. Addition of C. sorokiniana CMBB-146 could only increase the growth rate of P. malhamensis in the first 3 days, but the cell concentration of P. malhamensis started to decrease gradually after 4 days. Compared to autotrophic P. malhamensis, total 6,583 and 3,510 genes were significantly and differentially expressed in P. malhamensis fed by digestible microalgae and indigestible microalgae, respectively. Compared to autotrophic cells, genes related to the ribosome, lysosome, glycolysis, gluconeogenesis, TCA cycle, β-oxidation, duplication, and β-1,3-glucan in P. malhamensis grazing on digestible prey were up-regulated, while genes related to light harvesting and key enzymes referring to chlorophyll were down-regulated. Genes related to apoptosis and necrosis in P. malhamensis were up-regulated after grazing on indigestible microalgae compared to the autotrophic group, which we suggest is associated with the up-regulation of genes related to lysosome enzymes. This study provides abundant information on the potential intracellular physiological responses of P. malhamensis during the process of nutritional adaptation.

Linking extreme seasonality and gene expression in Arctic marine protists

At high latitudes, strong seasonal differences in light availability affect marine organisms and regulate the timing of ecosystem processes. Marine protists are key players in Arctic aquatic ecosystems, yet little is known about their ecological roles over yearly cycles. This is especially true for the dark polar night period, which up until recently was assumed to be devoid of biological activity. A 12 million transcripts catalogue was built from 0.45 to 10 μm protist assemblages sampled over 13 months in a time series station in an Arctic fjord in Svalbard. Community gene expression was correlated with seasonality, with light as the main driving factor. Transcript diversity and evenness were higher during polar night compared to polar day. Light-dependent functions had higher relative expression during polar day, except phototransduction. 64% of the most expressed genes could not be functionally annotated, yet up to 78% were identified in Arctic samples from Tara Oceans, suggesting that Arctic marine assemblages are distinct from those from other oceans. Our study increases understanding of the links between extreme seasonality and biological processes in pico- and nanoplanktonic protists. Our results set the ground for future monitoring studies investigating the seasonal impact of climate change on the communities of microbial eukaryotes in the High Arctic.

Functional ecology of planktonic ciliates: Measuring mortality rates in response to starvation

Population dynamics of aquatic ciliates are controlled "bottom-up" via food supply and "top-down" by grazing and parasitism. While intrinsic growth rates of ciliates under saturating food conditions have been studied in some detail, mortality rates induced by starvation have received little attention thus far. To this end, we examined the response of three algivorous freshwater ciliate species to starvation using three different optical methods. Two of these methods, i.e. ciliate mortality rates (δ) estimated from (i) numerical response experiments and (ii) the rate of decline (ROD) in cell numbers, investigated the response of the ciliate population using conventional light microscopy. The third method, imaging cytometry using a FlowCAM instrument, monitored single cells during the starvation experiment. Like light microscopy, the FlowCAM approach estimated δ based on ROD in the experimental containers. However, imaging cytometry also measured the relative cellular chlorophyll a content in the ciliates' food vacuoles as a proxy for the nutritional status of the cells. The linear decline of the cellular chl. a yielded an independent estimate of δ that was similar to δ calculated from ROD. Additionally, the FlowCAM measurements revealed a high degree of phenotypic plasticity of the ciliates when exposed to starvation.

Gene expression dynamics of natural assemblages of heterotrophic flagellates during bacterivory

BACKGROUND

Marine heterotrophic flagellates (HF) are dominant bacterivores in the ocean, where they represent the trophic link between bacteria and higher trophic levels and participate in the recycling of inorganic nutrients for regenerated primary production. Studying their activity and function in the ecosystem is challenging since most of the HFs in the ocean are still uncultured. In the present work, we investigated gene expression of natural HF communities during bacterivory in four unamended seawater incubations.

RESULTS

The most abundant species growing in our incubations belonged to the taxonomic groups MAST-4, MAST-7, Chrysophyceae, and Telonemia. Gene expression dynamics were similar between incubations and could be divided into three states based on microbial counts, each state displaying distinct expression patterns. The analysis of samples where HF growth was highest revealed some highly expressed genes that could be related to bacterivory. Using available genomic and transcriptomic references, we identified 25 species growing in our incubations and used those to compare the expression levels of these specific genes. Video Abstract CONCLUSIONS: Our results indicate that several peptidases, together with some glycoside hydrolases and glycosyltransferases, are more expressed in phagotrophic than in phototrophic species, and thus could be used to infer the process of bacterivory in natural assemblages.

RNA outperforms DNA-based metabarcoding in assessing the diversity and response of microeukaryotes to environmental variables in the Arctic Ocean

The Arctic Ocean (AO) has a harsh environment characterized by low temperatures, extensive ice coverage, and periodic freezing and melting of sea ice, which has provided diverse habitats for microorganisms. Prior studies primarily focused on microeukaryote communities in the upper water or sea ice based on environmental DNA, leaving the composition of active microeukaryotes in the diverse AO environments largely unknown. This study provided a vertical assessment of microeukaryote communities in the AO from snow and ice to sea water at a depth of 1670 m using high-throughput sequencing of co-extracted DNA and RNA. RNA extracts depicted microeukaryote community structure and intergroup correlations more accurately and responded more sensitively to environmental conditions than those derived from DNA. Using RNA:DNA ratios as a proxy for relative activity of major taxonomic groups, the metabolic activities of major microeukaryote groups were determined along depth. Analysis of co-occurrence networks showed that parasitism between Syndiniales and dinoflagellates/ciliates in the deep ocean may be significant. This study increased our knowledge of the diversity of active microeukaryote communities and highlighted the importance of using RNA-based sequencing over DNA-based sequencing to examine the relationship between microeukaryote assemblages and the responses of microeukaryotes to environmental variables in the AO.

Reverse engineering environmental metatranscriptomes clarifies best practices for eukaryotic assembly

BACKGROUND

Diverse communities of microbial eukaryotes in the global ocean provide a variety of essential ecosystem services, from primary production and carbon flow through trophic transfer to cooperation via symbioses. Increasingly, these communities are being understood through the lens of omics tools, which enable high-throughput processing of diverse communities. Metatranscriptomics offers an understanding of near real-time gene expression in microbial eukaryotic communities, providing a window into community metabolic activity.

RESULTS

Here we present a workflow for eukaryotic metatranscriptome assembly, and validate the ability of the pipeline to recapitulate real and manufactured eukaryotic community-level expression data. We also include an open-source tool for simulating environmental metatranscriptomes for testing and validation purposes. We reanalyze previously published metatranscriptomic datasets using our metatranscriptome analysis approach.

CONCLUSION

We determined that a multi-assembler approach improves eukaryotic metatranscriptome assembly based on recapitulated taxonomic and functional annotations from an in-silico mock community. The systematic validation of metatranscriptome assembly and annotation methods provided here is a necessary step to assess the fidelity of our community composition measurements and functional content assignments from eukaryotic metatranscriptomes.

Full-Length Transcriptome Analysis of the Ichthyotoxic Harmful Alga Heterosigma akashiwo (Raphidophyceae) Using Single-Molecule Real-Time Sequencing

The raphidophyte Heterosigma akashiwo is a harmful algal species. The bloom of this organism has been associated with the massive mortality of fish in many coastal waters. To investigate the molecular mechanism of H. akashiwo blooms, having a reliable reference transcriptome of this species is essential. Therefore, in this study, a full-length transcriptome of H. akashiwo was obtained by single-molecule real-time sequencing. In total, 45.44 Gb subread bases were generated, and 16,668 unigenes were obtained after the sequencing data processing. A total of 8666 (52.00%) unigenes were successfully annotated using seven public databases. Among them, mostly phosphorus and nitrogen metabolism genes were detected. Moreover, there were 300 putative transcription factors, 4392 putative long non-coding RNAs, and 7851 simple sequence repeats predicted. This study provides a valuable reference transcriptome for understanding how H. akashiwo blooms at a molecular level.

Community assembly patterns and processes of microbiome responses to habitats and Mytilopsis sallei invasion in the tidal zones of the Pearl River Estuary

The sustainability of estuarine ecosystem functions depends on the stabilization of microbial ecological processes. However, due to the unique and variable habitat characteristics of estuarine areas, in-depth studies on ecological processes such as the spatial distribution and assembly patterns of microbial community structure are lacking. As methods to elucidate this structure, we used 16S rDNA, 18S rDNA and ITS sequencing technologies to study the composition, diversity, spatial pattern and aggregation mechanism of the bacterial, protist and fungal communities in the tidal zones of the Pearl River Estuary (PRETZ). The abundance of bacterial communities was much higher than that of protists and fungi, and the spatial pattern was obvious in PRETZ. The application of neutral and null models revealed the assembly process of three microbial communities dominated by stochastic processes. Among the stochastic processes, undominated processes (64.03 %, 62.45 %, and 59.29 %) were the most critical processes in the assembly of bacterial, fungal and protist communities. Meanwhile, environmental variables, geographic locations, and biological factors were associated with the composition and assembly of bacterial, protist, and fungal communities. Among the environmental variables, dissolved oxygen and salinity were the main predictors that jointly affected the differences in the community structure of the three microorganisms, and geographic location was the second predictor affecting the community structure of the three microorganisms and had a more pronounced effect on the diversity and network structure of the bacterial and fungal communities. However, biological factors exerted a weaker effect on the microbial community structure than spatial factors and only affected bacteria and protists; the invasive species Mytilopsis sallei only affected the process of protist community assembly. In addition, environmental variables affected the relative importance of stochastic processes. In summary, the formation of microbial communities in the PRETZ was affected by random processes, environmental variables, geographic location, and invasive species.

Characterizing potential pathogens from intracellular bacterial community of protists in wastewater treatment plants

Protists are a trophically diverse and biogeochemically significant component of water environments and are widely reported as hosts of bacteria. However, the potential role of protists in wastewater treatment plants (WWTPs) as reservoirs for human pathogens does not appear to have received adequate attention. Here, a combination of fluorescence-activated cell sorting and Illumina sequencing was applied to characterize the dynamics of the internalized bacterial community of the enriched protists from the influents and effluents of five WWTPs. The results showed that Proteobacteria (mainly Betaproteobacteria) dominate the intracellular bacterial communities of protists in both influents and effluents of WWTPs, accounting for 72.6% of the total intracellular bacterial communities. The most frequently detected genus was Sulfuricurvum in the influent samples, Chryseobacterium and Pseudomonas were most prevalent in the effluent samples. Compared with the influents, a more diverse and abundant intracellular bacterial community was observed in the effluents. Moreover, the potential intracellular bacterial pathogens were 26 times higher in effluents than in influents, with Pseudomonas fluorescens and Pseudomonas putida significantly enriched in effluents. This work provides insights into the dynamics of bacterial communities and potential pathogens harbored by protists in the influents and effluents from WWTPs, contributing to the improved evaluation of biosafety in WWTPs.

Eukaryotic plankton community assembly and influencing factors between continental shelf and slope sites in the northern South China Sea

Eukaryotic plankton are pivotal members of marine ecosystems playing crucial roles in marine food webs and biogeochemical cycles. However, understanding the patterns and drivers of their community assembly remains a grand challenge. A study was conducted in the northern South China Sea (SCS) to address this issue. Here, 49 samples were collected and size-fractionated from discrete depths at continental shelf and continental slope in the northern SCS over a diel cycle. From high throughput sequencing of the 18S rDNA gene V4 region, 2463 operational taxonomic units (OTUs) were retrieved. Alveolata and Opisthokonta overwhelmingly dominated the assemblages in the abundance (44.76%, 31.08%) and species richness (59%, 12%). Biodiversity was higher in the slope than the shelf and increased with depth. Temperature and salinity appeared to be the most important deterministic drivers of taxon composition. Community structure was influenced by multiple factors in the importance order of: environmental factors (temperature + salinity) > spatial factor > water depth > sampling time. Furthermore, the neutral model explained more variations in the smaller-sized (0.22-3 μm) community (24%) than larger-sized (3-200 μm) community (16%) but generally explained less variations than did deterministic processes. Additionally, our data indicated that the larger plankton might be more environmentally filtered and less plastic whereas the smaller plankton had stronger dispersal ability. This study sheds light on the differential contributions of the deterministic process and stochastic process and complexities of assembly mechanisms in shaping the community assembly of micro-nano and pico-eukaryotic biospheres in a subtropical ocean.

Long-Read Sequencing Improves Recovery of Picoeukaryotic Genomes and Zooplankton Marker Genes from Marine Metagenomes

Long-read sequencing offers the potential to improve metagenome assemblies and provide more robust assessments of microbial community composition and function than short-read sequencing. We applied Pacific Biosciences (PacBio) CCS (circular consensus sequencing) HiFi shotgun sequencing to 14 marine water column samples and compared the results with those for short-read metagenomes from the corresponding environmental DNA samples. We found that long-read metagenomes varied widely in quality and biological information. The community compositions of the corresponding long- and short-read metagenomes were frequently dissimilar, suggesting higher stochasticity and/or bias associated with PacBio sequencing. Long reads provided few improvements to the assembly qualities, gene annotations, and prokaryotic metagenome-assembled genome (MAG) binning results. However, only long reads produced high-quality eukaryotic MAGs and contigs containing complete zooplankton marker gene sequences. These results suggest that high-quality long-read metagenomes can improve marine community composition analyses and provide important insight into eukaryotic phyto- and zooplankton genetics, but the benefits may be outweighed by the inconsistent data quality. IMPORTANCE Ocean microbes provide critical ecosystem services, but most remain uncultivated. Their communities can be studied through shotgun metagenomic sequencing and bioinformatic analyses, including binning draft microbial genomes. However, most sequencing to date has been done using short-read technology, which rarely yields genome sequences of key microbes like SAR11. Long-read sequencing can improve metagenome assemblies but is hampered by technological shortcomings and high costs. In this study, we compared long- and short-read sequencing of marine metagenomes. We found a wide range of long-read metagenome qualities and minimal improvements to microbiome analyses. However, long reads generated draft genomes of eukaryotic algal species and provided full-length marker gene sequences of zooplankton species, including krill and copepods. These results suggest that long-read sequencing can provide greater genetic insight into the wide diversity of eukaryotic phyto- and zooplankton that interact as part of and with the marine microbiome.

Genomic differentiation of three pico-phytoplankton species in the Mediterranean Sea

For more than a decade, high-throughput sequencing has transformed the study of marine planktonic communities and has highlighted the extent of protist diversity in these ecosystems. Nevertheless, little is known relative to their genomic diversity at the species-scale as well as their major speciation mechanisms. An increasing number of data obtained from global scale sampling campaigns is becoming publicly available, and we postulate that metagenomic data could contribute to deciphering the processes shaping protist genomic differentiation in the marine realm. As a proof of concept, we developed a findable, accessible, interoperable and reusable (FAIR) pipeline and focused on the Mediterranean Sea to study three a priori abundant protist species: Bathycoccus prasinos, Pelagomonas calceolata and Phaeocystis cordata. We compared the genomic differentiation of each species in light of geographic, environmental and oceanographic distances. We highlighted that isolation-by-environment shapes the genomic differentiation of B. prasinos, whereas P. cordata is impacted by geographic distance (i.e. isolation-by-distance). At present time, the use of metagenomics to accurately estimate the genomic differentiation of protists remains challenging since coverages are lower compared to traditional population surveys. However, our approach sheds light on ecological and evolutionary processes occurring within natural marine populations and paves the way for future protist population metagenomic studies.

Upwelling and nutrient dynamics in the Arabian Gulf and sea of Oman

This study demonstrates the vertical and horizontal distribution of nutrients and the seasonal response of nutrients to upwelling in the Arabian Gulf and the Sea of Oman. Thus, monthly data on nitrate, phosphate, and silicate are obtained from the World Ocean Atlas 2018 (WOA), as well as estimates of coastal and curl driven upwelling in both regions. The results of the study indicate that the Sea of Oman's surface and deep waters contained higher concentrations of nutrients than the Arabian Gulf by 80%. In addition, both regions have exhibited a general increase in the vertical distribution of nutrients as the depth increases. Among the aforementioned nutrients, nitrate is found to be a more limiting nutrient for phytoplankton growth than phosphate as the nitrate-to-phosphate ratios (N:P) in surface waters are lower (≤ 4.6:1) than the Redfield ratio (16:1). As for the upwelling, curl-driven upwelling accounts for more than half of the total upwelling in both regions, and both play an important role in nutrient transport. Thus, nutrients are upwelled from the subsurface to the mixed layer at a rate of 50% in the Oman Sea from 140 m to 20 m during the summer and to 40 m during the winter. Similarly, the Arabian Gulf shows 50% transport for nitrates, but 32% for phosphates, from 20 m to 5-10 m. However, due to the abundance of diatoms at the surface of the Arabian Gulf, the surface silicate content is 30% higher than that of the deeper waters.

Microbial Dark Matter: from Discovery to Applications

With the rapid increase of the microbiome samples and sequencing data, more and more knowledge about microbial communities has been gained. However, there is still much more to learn about microbial communities, including billions of novel species and genes, as well as countless spatiotemporal dynamic patterns within the microbial communities, which together form the microbial dark matter. In this work, we summarized the dark matter in microbiome research and reviewed current data mining methods, especially artificial intelligence (AI) methods, for different types of knowledge discovery from microbial dark matter. We also provided case studies on using AI methods for microbiome data mining and knowledge discovery. In summary, we view microbial dark matter not as a problem to be solved but as an opportunity for AI methods to explore, with the goal of advancing our understanding of microbial communities, as well as developing better solutions to global concerns about human health and the environment.

Protist Diversity and Metabolic Strategy in Freshwater Lakes Are Shaped by Trophic State and Watershed Land Use on a Continental Scale

Protists play key roles in aquatic food webs as primary producers, predators, nutrient recyclers, and symbionts. However, a comprehensive view of protist diversity in freshwaters has been challenged by the immense environmental heterogeneity among lakes worldwide. We assessed protist diversity in the surface waters of 366 freshwater lakes across a north temperate to subarctic range covering nearly 8.4 million km² of Canada. Sampled lakes represented broad gradients in size, trophic state, and watershed land use. Hypereutrophic lakes contained the least diverse and most distinct protist communities relative to nutrient-poor lakes. Greater taxonomic variation among eutrophic lakes was mainly a product of heterotroph and mixotroph diversity, whereas phototroph assemblages were more similar under high-nutrient conditions. Overall, local physicochemical factors, particularly ion and nutrient concentrations, elicited the strongest responses in community structure, far outweighing the effects of geographic gradients. Despite their contrasting distribution patterns, obligate phototroph and heterotroph turnover was predicted by an overlapping set of environmental factors, while the metabolic plasticity of mixotrophs may have made them less predictable. Notably, protist diversity was associated with variation in watershed soil pH and agricultural crop coverage, pointing to human impact on the land-water interface that has not been previously identified in studies on smaller scales. Our study exposes the importance of both within-lake and external watershed characteristics in explaining protist diversity and biogeography, critical information for further developing an understanding of how freshwater lakes and their watersheds are impacted by anthropogenic stressors.

IMPORTANCE Freshwater lakes are experiencing rapid changes under accelerated anthropogenic stress and a warming climate. Microorganisms underpin aquatic food webs, yet little is known about how freshwater microbial communities are responding to human impact. Here, we assessed the diversity of protists and their myriad ecological roles in lakes varying in size across watersheds experiencing a range of land use pressures by leveraging data from a continental-scale survey of Canadian lakes. We found evidence of human impact on protist assemblages through an association with lake trophic state and extending to agricultural activity and soil characteristics in the surrounding watershed. Furthermore, trophic state appeared to explain the distributions of phototrophic and heterotrophic protists in contrasting ways. Our findings highlight the vulnerability of lake ecosystems to increased land use and the importance of assessing terrestrial interfaces to elucidate freshwater ecosystem dynamics.

Development of transcriptomics-based growth rate indices in two model eukaryotes and relevance to metatranscriptomic datasets

Growth rate estimation is important to understand the flow of energy and nutrient elements in an ecosystem, but it has remained challenging, especially on microscopic organisms. In this study, we propose four growth rate indices that use mRNA abundance ratios between nuclear and mitochondrial genes: (1) total nuclear and mitochondrial mRNA ratio (Nuc:Mito‐TmRNA); (2) nuclear and mitochondrial ribosomal protein mRNA ratio (Nuc:Mito‐RPmRNA); (3) gene ontology (GO) terms and total mitochondrial mRNA ratios; and (4) nuclear and mitochondrial specific gene mRNA ratio. We examine these proposed ratios using RNA‐Seq datasets of Daphnia magna, and Saccharomyces cerevisiae retrieved from the NCBI Short Read Archive. The results showed that both Nuc:Mito‐TmRNA and Nuc:Mito‐RPmRNA ratio indices showed significant correlations with the growth rate for both species. A large number of GO terms mRNA ratios showed significant correlations with the growth rate of S. cerevisiae. Lastly, we identified mRNA ratios of several specific nuclear and mitochondrial gene pairs that showed significant correlations. We foresee future implications for the proposed mRNA ratios used in metatranscriptome analyses to estimate the growth rate of communities and species.

Evolution of Phytoplankton as Estimated from Genetic Diversity

Phytoplankton are photosynthetic, single-celled organisms producing almost half of all oxygen on Earth and play a central role as prey for higher organisms, making them irreplaceable in the marine food web. As Global Change proceeds, imposing rapidly intensifying selection pressures, phytoplankton are forced to undergo evolution, local extinction, or redistribution, with potentially cascading effects throughout the marine ecosystem. Recent results from the field of population genetics display high levels of standing genetic diversity in natural phytoplankton populations, providing ample ‘evolutionary options’ and implying high adaptive potential to changing conditions. This potential for adaptive evolution is realized in several studies of experimental evolution, even though most of these studies investigate the evolution of only single strains. This, however, shows that phytoplankton not only evolve from standing genetic diversity, but also rely on de novo mutations. Recent global sampling campaigns show that the immense intraspecific diversity of phytoplankton in the marine ecosystem has been significantly underestimated, meaning we are only studying a minor portion of the relevant variability in the context of Global Change and evolution. An increased understanding of genomic diversity is primarily hampered by the low number of ecologically representative reference genomes of eukaryotic phytoplankton and the functional annotation of these. However, emerging technologies relying on metagenome and transcriptome data may offer a more realistic understanding of phytoplankton diversity.

TreeTuner: A pipeline for minimizing redundancy and complexity in large phylogenetic datasets

Various bioinformatics protocols have been developed for trimming the number of operational taxonomic units (OTUs) in phylogenetic datasets, but they typically require significant manual intervention. Here we present TreeTuner, a semiautomated pipeline that allows both coarse and fine-scale tuning of large protein sequence phylogenetic datasets via the minimization of OTU redundancy. TreeTuner facilitates preliminary investigation of such datasets as well as more rigorous downstream analysis of specific subsets of OTUs.

For complete details on the use and execution of this protocol, please refer to Maruyama et al. (2013) and Sibbald et al. (2019).

•

Minimizes sequence redundancy in large phylogenetic datasets

•

Trims thousands of operational taxonomic units (OTUs) from a preliminary tree

•

Maintains the desired minimal taxonomic diversity and retains specific OTUs for analysis

•

Available coarse- and fine-tuning options depending on the phylogenetic question at hand

Population and molecular responses to warming in Netzelia tuberspinifera - An endemic and sensitive protist from East Asia

With the effects of global warming becoming ever more obvious, biodiversity conservation is facing severe challenges. Currently, a deeper understanding the mechanisms of the effects of warming on sensitive species has become an important topic in aquatic biodiversity and ecological management. Our study first overcame the "challenge" for a sensitive indicator species (Netzelia tuberspinifera, an endemic testate amoeba species in East Asia) of culturing under laboratory conditions, and then explored its molecular response mechanisms to warming using transcriptomic analysis. Our data indicate that temperature mainly drove the geographical and seasonal variation of N. tuberspinifera populations. Transcriptomic results indicate that when the temperature is <25 °C, rising temperature triggers the biosynthesis of ribosomes; while the temperature is >25 °C, it triggers molecular processes related with cell division, test formation and general biomass increase. However, once the temperature exceeds 40 °C, N. tuberspinifera is unable to survive. Following from these results, the distribution of N. tuberspinifera might expand towards higher altitude or latitude regions under global warming. For the first time, our study showed direct evidence for sensitive protozoa species that presents a very narrow adaptation mechanism to local climate. Our work provides fundamental data for regional biodiversity conservation and scientific reference in subtropical and tropical waterbodies.

Ecology of planktonic ciliates in a changing world: Concepts, methods, and challenges

Plankton ecologists ultimately focus on forecasting, both applied and environmental outcomes. We review how appreciating planktonic ciliates has become central to these predictions. We explore the 350‐year‐old canon on planktonic ciliates and examine its steady progression, which has been punctuated by conceptual insights and technological breakthroughs. By reflecting on this process, we offer suggestions as to where future leaps are needed, with an emphasis on predicting outcomes of global warming. We conclude that in terms of climate change research: (i) climatic hotspots (e.g. polar oceans) require attention; (ii) simply adding ciliate measurements to zooplankton/phytoplankton‐based sampling programs is inappropriate; (iii) elucidating the rare biosphere's functional ecology requires culture‐independent genetic methods; (iv) evaluating genetic adaptation (microevolution) and population composition shifts is required; (v) contrasting marine and freshwaters needs attention; (vi) mixotrophy needs attention; (vii) laboratory and field studies must couple automated measurements and molecular assessment of functional gene expression; (viii) ciliate trophic diversity requires appreciation; and (ix) marrying gene expression and function, coupled with climate change scenarios is needed. In short, continued academic efforts and financial support are essential to achieve the above; these will lead to understanding how ciliates will respond to climate change, providing tools for forecasting.

Resolving the microalgal gene landscape at the strain level: A novel hybrid transcriptome of Emiliania huxleyi CCMP3266

Microalgae are key ecological players with a complex evolutionary history. Genomic diversity, in addition to limited availability of high-quality genomes, challenge studies that aim to elucidate molecular mechanisms underlying microalgal ecophysiology. Here, we present a novel and comprehensive transcriptomic hybrid approach to generate a reference for genetic analyses, and resolve the microalgal gene landscape at the strain level. The approach is demonstrated for a strain of the coccolithophore microalga Emiliania huxleyi, which is a species complex with considerable genome variability. The investigated strain is commonly studied as a model for algal-bacterial interactions, and was therefore sequenced in the presence of bacteria to elicit the expression of interaction-relevant genes. We applied complementary PacBio Iso-Seq full-length cDNA, and poly(A)-independent Illumina total RNA sequencing, which resulted in a de novo assembled, near complete hybrid transcriptome. In particular, hybrid sequencing improved the reconstruction of long transcripts and increased the recovery of full-length transcript isoforms. To use the resulting hybrid transcriptome as a reference for genetic analyses, we demonstrate a method that collapses the transcriptome into a genome-like dataset, termed "synthetic genome" (sGenome). We used the sGenome as a reference to visually confirm the robustness of the CCMP3266 gene assembly, to conduct differential gene expression analysis, and to characterize novel E. huxleyi genes. The newly-identified genes contribute to our understanding of E. huxleyi genome diversification, and are predicted to play a role in microbial interactions. Our transcriptomic toolkit can be implemented in various microalgae to facilitate mechanistic studies on microalgal diversity and ecology. Importance Microalgae are key players in the ecology and biogeochemistry of our oceans. Efforts to implement genomic and transcriptomic tools in laboratory studies involving microalgae suffer from the lack of published genomes. In the case of coccolithophore microalgae, the problem has long been recognized; the model species Emiliania huxleyi is a species complex with genomes composed of a core, and a large variable portion. To study the role of the variable portion in niche adaptation, and specifically in microbial interactions, strain-specific genetic information is required. Here we present a novel transcriptomic hybrid approach, and generated strain-specific genome-like information. We demonstrate our approach on an E. huxleyi strain that is co-cultivated with bacteria. By constructing a "synthetic genome", we generated comprehensive gene annotations that enabled accurate analyses of gene expression patterns. Importantly, we unveiled novel genes in the variable portion of E. huxleyi that play putative roles in microbial interactions.

A cell-cell atlas approach for understanding symbiotic interactions between microbes

Natural environments are composed of a huge diversity of microorganisms interacting with each other to form complex functional networks. Our understanding of the operative nature of host-symbiont associations is limited because propagating such associations in a laboratory is challenging. The advent of single-cell technologies applied to, for example, animal cells and apicomplexan parasites has revolutionized our understanding of development and disease. Such cell atlas approaches generate maps of cell-specific processes and variations within cellular populations. These methods can now be combined with cellular-imaging so that interaction stage versus transcriptome state can be quantized for microbe-microbe interactions. We predict that the combination of these methods applied to the study of symbioses will transform our understanding of many ecological interactions, including those sampled directly from natural environments.

Diversity and ecology of protists revealed by metabarcoding

Protists are the dominant eukaryotes in the biosphere where they play key functional roles. While protists have been studied for over a century, it is the high-throughput sequencing of molecular markers from environmental samples - the approach of metabarcoding - that has revealed just how diverse, and abundant, these small organisms are. Metabarcoding is now routine to survey environmental diversity, so data have rapidly accumulated from a multitude of environments and at different sampling scales. This mass of data has provided unprecedented opportunities to study the taxonomic and functional diversity of protists, and how this diversity is organised in space and time. Here, we use metabarcoding as a common thread to discuss the state of knowledge in protist diversity research, from technical considerations of the approach to important insights gained on diversity patterns and the processes that might have structured this diversity. In addition to these insights, we conclude that metabarcoding is on the verge of an exciting added dimension thanks to the maturation of high-throughput long-read sequencing, so that a robust eco-evolutionary framework of protist diversity is within reach.

Changes in protist communities in drainages across the Pearl River Delta under anthropogenic influence

Drainages in the Pearl River Delta urban agglomeration (PRDUA) host vital aquatic ecosystems and face enormous pressures from human activities in one of the largest urban agglomerations in the world. Despite being crucial components of aquatic ecosystems, the interactions and assembly processes of the protistan community are rarely explored in areas with serious anthropogenic disturbance. To elucidate the mechanisms of these processes, we used environmental DNA sequencing of 18S rDNA to investigate the influence of environmental factors and species interactions on the protistan community and its assembly in drainages of the PRDUA during summer. The protistan community showed a high level of diversity and a marked spatial pattern in this region. Community assembly was driven primarily by stochastic processes based on the Sloan neutral community model, explaining 74.28%, 75.82%, 73.67%, 74.40% and 51.24% of community variations in the BJ (Beijiang), XJ (Xijiang), PRD (Pearl River Delta), PRE (Pearl River Estuary) areas and in total, respectively. Meanwhile, environmental variables including temperature, pH, dissolved oxygen, transparency, nutrients and land use were strongly correlated with the composition and assembly of the protistan community, explaining 40.40% of variation in the protistan community. Furthermore, the bacterial community was simultaneously analysed by the 16S rDNA sequencing. Co-occurrence network analysis revealed that species interactions within bacteria (81.41% positive) or protists (82.80% positive), and those between bacteria and protists (50% positive and 50% negative) impacted the protistan community assembly. In summary, stochastic processes dominated, whereas species interactions and environmental factors also played important roles in shaping the protistan communities in drainages across the PRDUA. This study provides insights into the ecological patterns, assembly processes and species interactions underlying protistan dynamics in urban aquatic ecosystems experiencing serious anthropogenic disturbance.

TRAPID 2.0: a web application for taxonomic and functional analysis of de novo transcriptomes

Advances in high-throughput sequencing have resulted in a massive increase of RNA-Seq transcriptome data. However, the promise of rapid gene expression profiling in a specific tissue, condition, unicellular organism or microbial community comes with new computational challenges. Owing to the limited availability of well-resolved reference genomes, de novo assembled (meta)transcriptomes have emerged as popular tools for investigating the gene repertoire of previously uncharacterized organisms. Yet, despite their potential, these datasets often contain fragmented or contaminant sequences, and their analysis remains difficult. To alleviate some of these challenges, we developed TRAPID 2.0, a web application for the fast and efficient processing of assembled transcriptome data. The initial processing phase performs a global characterization of the input data, providing each transcript with several layers of annotation, comprising structural, functional, and taxonomic information. The exploratory phase enables downstream analyses from the web application. Available analyses include the assessment of gene space completeness, the functional analysis and comparison of transcript subsets, and the study of transcripts in an evolutionary context. A comparison with similar tools highlights TRAPID’s unique features. Finally, analyses performed within TRAPID 2.0 are complemented by interactive data visualizations, facilitating the extraction of new biological insights, as demonstrated with diatom community metatranscriptomes.

Diel transcriptional oscillations of light-sensitive regulatory elements in open-ocean eukaryotic plankton communities

Most organisms coordinate key biological events to coincide with the day/night cycle. These diel oscillations are entrained through the activity of light-sensitive photoreceptors that allow organisms to respond rapidly to changes in light exposure. In the ocean, the plankton community must additionally contend with dramatic changes in the quantity and quality of light over depth. Here, we show that the predominantly blue-light field in the open-ocean environment may have driven expansion of blue light-sensitive regulatory elements in open-ocean eukaryotic plankton derived from secondary and tertiary endosymbiosis. The diel transcription of genes encoding light-sensitive elements indicate that photosynthetic and heterotrophic marine protists respond to and anticipate fluctuating light conditions in the dynamic marine environment.

The 24-h cycle of light and darkness governs daily rhythms of complex behaviors across all domains of life. Intracellular photoreceptors sense specific wavelengths of light that can reset the internal circadian clock and/or elicit distinct phenotypic responses. In the surface ocean, microbial communities additionally modulate nonrhythmic changes in light quality and quantity as they are mixed to different depths. Here, we show that eukaryotic plankton in the North Pacific Subtropical Gyre transcribe genes encoding light-sensitive proteins that may serve as light-activated transcription factors, elicit light-driven electrical/chemical cascades, or initiate secondary messenger-signaling cascades. Overall, the protistan community relies on blue light-sensitive photoreceptors of the cryptochrome/photolyase family, and proteins containing the Light-Oxygen-Voltage (LOV) domain. The greatest diversification occurred within Haptophyta and photosynthetic stramenopiles where the LOV domain was combined with different DNA-binding domains and secondary signal-transduction motifs. Flagellated protists utilize green-light sensory rhodopsins and blue-light helmchromes, potentially underlying phototactic/photophobic and other behaviors toward specific wavelengths of light. Photoreceptors such as phytochromes appear to play minor roles in the North Pacific Subtropical Gyre. Transcript abundance of environmental light-sensitive protein-encoding genes that display diel patterns are found to primarily peak at dawn. The exceptions are the LOV-domain transcription factors with peaks in transcript abundances at different times and putative phototaxis photoreceptors transcribed throughout the day. Together, these data illustrate the diversity of light-sensitive proteins that may allow disparate groups of protists to respond to light and potentially synchronize patterns of growth, division, and mortality within the dynamic ocean environment.

De novo transcriptome assembly of the green alga Ankistrodesmus falcatus

Ankistrodesmus falcatus is a globally distributed freshwater chlorophyte that is a candidate for biofuel production, is used to study the effects of toxins on aquatic communities, and is used as food in zooplankton research. Each of these research fields is transitioning to genomic tools. We created a reference transcriptome for of A. falcatus using NextGen sequencing and de novo assembly methods including Trinity, Velvet-Oases, and EvidentialGene. The assembled transcriptome has a total of 17,997 contigs, an N50 value of 2,462, and a GC content of 64.8%. BUSCO analysis recovered 83.3% of total chlorophyte BUSCOs and 82.5% of the eukaryotic BUSCOs. A portion (7.9%) of these supposedly single-copy genes were found to have transcriptionally active, distinct duplicates. We annotated the assembly using the dammit annotation pipeline, resulting in putative functional annotation for 68.89% of the assembly. Using available rbcL sequences from 16 strains (10 species) of Ankistrodesmus, we constructed a neighbor-joining phylogeny to illustrate genetic distances of our A. falcatus strain to other members of the genus. This assembly will be valuable for researchers seeking to identify Ankistrodesmus sequences in metatranscriptomic and metagenomic field studies and in experiments where separating expression responses of zooplankton and their algal food sources through bioinformatics is important.

High and specific diversity of protists in the deep-sea basins dominated by diplonemids, kinetoplastids, ciliates and foraminiferans

Heterotrophic protists (unicellular eukaryotes) form a major link from bacteria and algae to higher trophic levels in the sunlit ocean. Their role on the deep seafloor, however, is only fragmentarily understood, despite their potential key function for global carbon cycling. Using the approach of combined DNA metabarcoding and cultivation-based surveys of 11 deep-sea regions, we show that protist communities, mostly overlooked in current deep-sea foodweb models, are highly specific, locally diverse and have little overlap to pelagic communities. Besides traditionally considered foraminiferans, tiny protists including diplonemids, kinetoplastids and ciliates were genetically highly diverse considerably exceeding the diversity of metazoans. Deep-sea protists, including many parasitic species, represent thus one of the most diverse biodiversity compartments of the Earth system, forming an essential link to metazoans.

Soil microbial legacies differ following drying-rewetting and freezing-thawing cycles

Climate change alters frequencies and intensities of soil drying-rewetting and freezing-thawing cycles. These fluctuations affect soil water availability, a crucial driver of soil microbial activity. While these fluctuations are leaving imprints on soil microbiome structures, the question remains if the legacy of one type of weather fluctuation (e.g., drying-rewetting) affects the community response to the other (e.g., freezing-thawing). As both phenomenons give similar water availability fluctuations, we hypothesized that freezing-thawing and drying-rewetting cycles have similar effects on the soil microbiome. We tested this hypothesis by establishing targeted microcosm experiments. We created a legacy by exposing soil samples to a freezing-thawing or drying-rewetting cycle (phase 1), followed by an additional drying-rewetting or freezing-thawing cycle (phase 2). We measured soil respiration and analyzed soil microbiome structures. Across experiments, larger CO₂ pulses and changes in microbiome structures were observed after rewetting than thawing. Drying-rewetting legacy affected the microbiome and CO₂ emissions upon the following freezing-thawing cycle. Conversely, freezing-thawing legacy did not affect the microbial response to the drying-rewetting cycle. Our results suggest that drying-rewetting cycles have stronger effects on soil microbial communities and CO₂ production than freezing-thawing cycles and that this pattern is mediated by sustained changes in soil microbiome structures.

New paradigm in diatom omics and genetic manipulation

Diatoms are one of the most heterogeneous eukaryotic plankton known for regulating earth's biogeochemical cycles and maintaining the marine ecosystems ever since the late Eocene epoch. The advent of multidisciplinary omics approach has both epitomized and revolutionized the nature of their chimeric genetic toolkit, ecophysiology, and metabolic adaptability as well as their interaction with other communities. In addition, advanced functional annotation of transcriptomic and proteomic data using cutting edge bioinformatics tools together with high-resolution genome-scale mathematical modeling has effectively proven as the catapult in solving genetic bottlenecks in microbial as well as diatom exploration. In this review, a corroborative summation of the robust work done in manipulating, engineering, and sequencing of the diatom genomes besides underpinning the holistic application of omics in transcription and translation has been discussed in order to shrewd their multifarious novel potential in the field of biotechnology and provide an insight into their dynamic evolutionary relevance.