Abstract 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 communities of multiple species. Metatranscriptomics in particular offers an understanding of real-time gene expression in microbial eukaryotic communities, providing a window into community metabolic activity. However, these approaches are rarely validated. A systematic validation of metatranscriptome assembly and annotation methods is necessary to appropriately assess (1) the fidelity of our community composition measurements and (2) the confidence of novel taxonomic and functional content discovered with metatranscriptomics. Here, we present euk rhythmic , an automated and customizable multi-assembler workflow for eukaryotic metatranscriptome assembly, annotation, and analysis, and validate the ability of the pipeline to recapitulate simulated eukaryotic community-level expression data treated as a gold standard. Finally, we apply this pipeline to two previously published metatranscriptomic datasets, showing both recapitulation of previous results as well as enhanced taxonomic and functional gene discovery.

Abstract Meta-omics has become commonplace in the study of microbial eukaryotes. The explosion of available data has inspired large-scale analyses, including species or taxonomic group distribution mapping, gene catalog construction, and inference on the functional roles and activities of microbial eukaryotes in situ . However, genome and transcriptome databases are prone to misannotation biases, and meta-omic inventories may have no recoverable taxonomic annotation for more than half of assembled contigs or predicted proteins. Direct mapping solely to organisms of interest might introduce a problematic misattribution bias, while full databases can annotate any cataloged organism but may be imbalanced between taxa. Here, we explore the potential pitfalls of common approaches to taxonomic annotation of protistan meta-omic datasets. We argue that ongoing curation of genetic resources is critical in accurately annotating protists in situ in meta-omic datasets. Moreover, we propose that precise taxonomic annotation of meta-omic data is a clustering problem rather than a feasible alignment problem. We show that taxonomic membership of sequence clusters demonstrates more accurate estimated community composition than returning exact sequence labels, and overlap between clusters can address database shortcomings. Clustering approaches can be applied to diverse environments while continuing to exploit the wealth of annotation data collated in databases, and database selection and evaluation is a critical part of correctly annotating protistan taxonomy in environmental datasets. We re-analyze three environmental datasets at three levels of taxonomic hierarchy in order to illustrate the critical importance of both database completeness and curation in enabling accurate environmental interpretation.

Abstract The keystone marine nitrogen fixer Trichodesmium thrives in high dust environments, and while experimental observations suggest that Trichodesmium colonies can access the essential nutrient iron from dust particles, it is not known the extent to which this occurs in the field. Here we demonstrate that Trichodesmium colonies actively process mineral particles in nature with direct molecular impacts. Microscopy and synchrotron-based imaging demonstrated heterogeneous associations with particles consistent with iron oxide and iron silicate minerals. Metaproteomic analysis of individual colonies revealed enrichment of biogeochemically-relevant proteins including photosynthesis proteins and metalloproteins containing iron, nickel, copper and zinc when particles were present. The iron-storage protein ferritin was particularly enriched implying accumulation of particle-derived iron, and multiple iron acquisition pathways including Fe(II), Fe(III), and Fe-siderophore transporters were engaged, including evidence of superoxide-driven particle dissolution. While the particles clearly provided iron, there was also evidence that the concentrated metals had toxic effects. The molecular mechanisms allowing Trichodesmium to interact with particulate minerals are fundamental to its success and global impact on nitrogen biogeochemistry, and may contribute to the leaching of particulate trace metals with implications for global iron and carbon cycling.

Frequent blooms of phytoplankton occur in coastal upwelling zones creating hotspots of biological productivity in the ocean. As cold, nutrient-rich water is brought up to sunlit layers from depth, phytoplankton are also transported upwards to seed surface blooms that are often dominated by diatoms. The physiological response of phytoplankton to this process, commonly referred to as shift-up, is characterized by rapid growth rates and increases in nitrate assimilation. To examine the molecular underpinnings behind this phenomenon, metatranscriptomics was applied to a simulated upwelling experiment using natural phytoplankton communities from the California Upwelling Zone. An increase in diatom growth following five days of incubation was attributed to the genera Chaetoceros and Pseudo-nitzschia. Here we show that certain bloom-forming diatoms exhibit a distinct transcriptional response that coordinates shift-up where diatoms exhibited the greatest transcriptional change following upwelling; however, comparison of co-expressed genes exposed overrepresentation of distinct sets within each of the dominant phytoplankton groups. The analysis revealed that diatoms frontload genes involved in nitrogen assimilation likely in order to outcompete other groups for available nitrogen during upwelling events. We speculate that the evolutionary success of diatoms may be due, in part, to this proactive response to frequently encountered changes in their environment.

Microeukaryotes are key contributors to marine carbon cycling. Their physiology, ecology, and interactions with the chemical environment are poorly understood in offshore ecosystems, and especially in the deep ocean. Using the Autonomous Underwater Vehicle Clio, microbial communities along a 1050 km transect in the western North Atlantic Ocean were surveyed at 10–200 m vertical depth increments to capture metabolic signatures spanning oligotrophic, continental margin, and productive coastal ecosystems. Microeukaryotes were examined using a paired metatranscriptomic and metaproteomic approach. Here we show a diverse surface assemblage consisting of stramenopiles, dinoflagellates and ciliates represented in both the transcript and protein fractions, with foraminifera, radiolaria, picozoa, and discoba proteins enriched at >200 m, and fungal proteins emerging in waters >3000 m. In the broad microeukaryote community, nitrogen stress biomarkers were found at coastal sites, with phosphorus stress biomarkers offshore. This multi-omics dataset broadens our understanding of how microeukaryotic taxa and their functional processes are structured along environmental gradients of temperature, light, and nutrients. Using an autonomous underwater vehicle, this study presents an integrated biogeochemical and multiomic analysis of microbial eukaryotes from the North Atlantic Ocean. The work highlights diverse communities that shift through depth zones, with signatures of nutrient biomarkers changing across a coastal-offshore spatial gradient.

Abstract Protists (microeukaryotes) are key contributors to marine carbon cycling, influencing the transfer of energy to higher trophic levels and the vertical movement of carbon to the ocean interior. Their physiology, ecology, and interactions with the chemical environment are still poorly understood in offshore ecosystems, and especially in the deep ocean. Using the Autonomous Underwater Vehicle (AUV) Clio , the microbial community along a 1,050 km transect in the western North Atlantic Ocean was surveyed at 10-200 m vertical depth increments to capture metabolic microeukaryote signatures spanning a gradient of oligotrophic, continental margin, and productive coastal ecosystems. Plankton biomass was collected along the surface of this transect and across depth features, and taxonomy and metabolic function were examined using a paired metatranscriptomic and metaproteomic approach. A shift in the microeukaryote community composition was observed from the euphotic zone through the mesopelagic and into the bathypelagic ocean. A diverse surface assemblage consisting of haptophytes, stramenopiles, dinoflagellates and ciliates was represented in both the transcript and protein fractions, with foraminifera, radiolaria, picozoa, and discoba proteins enriched at >200 m depth, and fungal proteins emerging in waters >3,000 m depth. In the broad microeukaryote community, nitrogen stress biomarkers were found in productive coastal sites, with phosphorus stress biomarkers in offshore waters where Saharan dust input is thought to supply iron and nitrogen. This multi-omics dataset broadens our understanding of how microeukaryotic taxa and their functional processes are structured along environmental gradients of temperature, light, macronutrients, and trace metals.

ABSTRACT Animals integrate developmental and nutritional signals before committing crucial resources to growth and reproduction; however, the pathways that perceive and respond to these inputs remain poorly understood. Here, we demonstrate that DRL-1 and FLR-4, which share similarity with mammalian mitogen-activated protein kinases, maintain lipid homeostasis in the C. elegans intestine. DRL-1 and FLR-4 function in a protein complex at the plasma membrane to promote development, as mutations in drl-1 or flr-4 confer slow growth, small body size, and impaired lipid homeostasis. To identify factors that oppose DRL-1/FLR-4, we performed a forward genetic screen for suppressors of the drl-1 mutant phenotypes and identified mutations in flr-2 and fshr-1 , which encode the orthologues of follicle stimulating hormone and its putative G protein-coupled receptor, respectively. In the absence of DRL-1/FLR-4, neuronal FLR-2 acts through intestinal FSHR-1 and Protein Kinase A signaling to restrict growth. Furthermore, we show that opposing signaling through DRL-1 and FLR-2 coordinates TIR-1 phase transition, which modulates downstream p38/PMK-1 activity, lipid homeostasis, and development. Finally, we identify a surprising noncanonical role for the developmental transcription factor PHA-4/FOXA in the intestine where it restricts growth in response to impaired DRL-1 signaling. Our work uncovers a complex multi-tissue signaling network that converges on p38 signaling to maintain homeostasis during development.

ABSTRACT Ion channels are necessary for proper water and nutrient absorption in the intestine, which supports cellular metabolism and organismal growth. While a role for Na + co-transporters and pumps in intestinal nutrient absorption is well defined, how individual K + uniporters function to maintain ion homeostasis is poorly understood. Using Caenorhabditis elegans , we show that a gain-of-function mutation in twk-26 , which encodes a two-pore domain K + ion channel orthologous to human KCNK3, facilitates nutrient absorption and suppresses the metabolic and developmental defects displayed by impaired intestinal MAP Kinase (MAPK) signaling. Mutations in drl-1 and flr-4, which encode two components of this MAPK pathway, cause severe growth defects, reduced lipid storage, and a dramatic increase in autophagic lysosomes, which mirror dietary restriction phenotypes. Additionally, these MAPK mutants display structural defects of the intestine and an impaired defecation motor program. We find that activation of TWK-26 reverses the dietary restriction-like state of the MAPK mutants by restoring intestinal nutrient absorption without correcting the intestinal bloating or defecation defects. This study provides unique insight into the mechanisms by which intestinal K + ion channels support intestinal metabolic homeostasis.

Abstract The continental shelf of the Western Antarctic Peninsula (WAP) is a highly variable system characterized by strong cross-shelf gradients, rapid regional change and large blooms of phytoplankton, notably diatoms. Rapid environmental changes coincide with shifts in plankton community composition and productivity, food web dynamics and biogeochemistry. Despite progress in identifying important environmental factors influencing plankton community composition in the WAP, the molecular basis for their survival in this oceanic region, as well as variations in species abundance, metabolism and distribution remain largely unresolved. Across a gradient of physicochemical parameters, we analyzed the metabolic profiles of phytoplankton as assessed through metatranscriptomic sequencing. Distinct phytoplankton communities and metabolisms closely mirrored the strong gradients in oceanographic parameters that existed from coastal to offshore regions. Diatoms were abundant in coastal, southern regions, where colder and fresher waters were conducive to a bloom of the centric diatom, Actinocyclus . Members of this genus invested heavily in growth and energy production; carbohydrate, amino acid and nucleotide biosynthesis pathways; and coping with oxidative stress, resulting in uniquely expressed metabolic profiles compared to other diatoms. We observed strong molecular evidence for iron limitation in shelf and slope regions of the WAP, where diatoms in these regions employed iron-starved induced proteins, a geranylgeranyl reductase, aquaporins, and urease, among other strategies, while limiting the use of iron-containing proteins. The metatranscriptomic survey performed here revealed functional differences in diatom communities and provides further insight into the environmental factors influencing the growth of diatoms and their predicted response to changes in ocean conditions.