Diel metabolic patterns revealed by in situ transcriptome and proteome in a vertically migratory copepod

Zooplankton undergo a diel vertical migration (DVM) which exposes them to gradients of light, temperature, oxygen, and food availability on a predictable daily schedule. Disentangling the co-varying and potentially synergistic interactions on metabolic rates has proven difficult, despite the importance of this migration for the delivery of metabolic waste products to the distinctly different daytime (deep) and nighttime (surface) habitats. This study examines the transcriptomic and proteomic profiles of the circumglobal migratory copepod, Pleuromamma xiphias, over the diel cycle. The transcriptome showed that 96% of differentially expressed genes were upregulated during the middle of the day - the period often considered to be of lowest zooplankton activity. The changes in protein abundance were more spread out over time, peaking (42% of comparisons) in the early evening. Between 9:00 and 15:00, both the transcriptome and proteome datasets showed increased expression related to chitin synthesis and degradation. Additionally, at 09:00 and 22:00, there were increases in myosin and vitellogenin proteins, potentially linked to the stress of migration and/or reproductive investment. Based on protein abundances detected, there is an inferred switch in broad metabolic processes, shifting from electron transport system in the day to glycolysis and glycogen mobilization in the afternoon/evening. These observations provide evidence of the diel impact of DVM on transcriptomic and proteomic pathways that likely influence metabolic processes and subsequent excretion products, and clarify how this behaviour results in the direct rapid transport of waste metabolites from the surface to the deep ocean.

Proteomic response of early juvenile Pacific oysters (Crassostrea gigas) to temperature

Pacific oysters (Crassostrea gigas) are a valuable aquaculture product that provides important ecosystem benefits. Among other threats, climate-driven changes in ocean temperature can impact oyster metabolism, survivorship, and immune function. We investigated how elevated temperature impacts larval oysters during settlement (19-33 days post-fertilization), using shotgun proteomics with data-independent acquisition to identify proteins present in the oysters after 2 weeks of exposure to 23 °C or 29 °C. Oysters maintained at elevated temperatures were larger and had a higher settlement rate, with 86% surviving to the end of the experiment; these oysters also had higher abundance trends of proteins related to metabolism and growth. Oysters held at 23 °C were smaller, had a decreased settlement rate, displayed 100% mortality, and had elevated abundance trends of proteins related to immune response. This novel use of proteomics was able to capture characteristic shifts in protein abundance that hint at important differences in the phenotypic response of Pacific oysters to temperature regimes. Additionally, this work has produced a robust proteomic product that will be the basis for future research on bivalve developmental processes.

Coupled microbiome analyses highlights relative functional roles of bacteria in a bivalve hatchery

BACKGROUND

Microbial communities are ubiquitous throughout ecosystems and are commensal with hosts across taxonomic boundaries. Environmental and species-specific microbiomes are instrumental in maintaining ecosystem and host health, respectively. The introduction of pathogenic microbes that shift microbiome community structure can lead to illness and death. Understanding the dynamics of microbiomes across a diversity of environments and hosts will help us to better understand which taxa forecast survival and which forecast mortality events.

RESULTS

We characterized the bacterial community microbiome in the water of a commercial shellfish hatchery in Washington state, USA, where the hatchery has been plagued by recurring and unexplained larval mortality events. By applying the complementary methods of metagenomics and metaproteomics we were able to more fully characterize the bacterial taxa in the hatchery at high (pH 8.2) and low (pH 7.1) pH that were metabolically active versus present but not contributing metabolically. There were shifts in the taxonomy and functional profile of the microbiome between pH and over time. Based on detected metagenomic reads and metaproteomic peptide spectral matches, some taxa were more metabolically active than expected based on presence alone (Deltaproteobacteria, Alphaproteobacteria) and some were less metabolically active than expected (e.g., Betaproteobacteria, Cytophagia). There was little correlation between potential and realized metabolic function based on Gene Ontology analysis of detected genes and peptides.

CONCLUSION

The complementary methods of metagenomics and metaproteomics contribute to a more full characterization of bacterial taxa that are potentially active versus truly metabolically active and thus impact water quality and inter-trophic relationships.

Larval Geoduck (Panopea generosa) Proteomic Response to Ciliates

The innate immune response is active in invertebrate larvae from early development. Induction of immune response pathways may occur as part of the natural progression of larval development, but an up-regulation of pathways can also occur in response to a pathogen. Here, we took advantage of a protozoan ciliate infestation of a larval geoduck clam culture in a commercial hatchery to investigate the molecular underpinnings of the innate immune response of the larvae to the pathogen. Larval proteomes were analyzed on days 4-10 post-fertilization; ciliates were present on days 8 and 10 post-fertilization. Through comparisons with larval cultures that did not encounter ciliates, proteins implicated in the response to ciliate presence were identified using mass spectrometry-based proteomics. Ciliate response proteins included many associated with ribosomal synthesis and protein translation, suggesting the importance of protein synthesis during the larval immune response. There was also an increased abundance of proteins typically associated with the stress and immune responses during ciliate exposure, such as heat shock proteins, glutathione metabolism, and the reactive oxygen species response. These findings provide a basic understanding of the bivalve molecular response to a mortality-inducing ciliate and improved characterization of the ontogenetic development of the innate immune response.

Growth phase proteomics of the heterotrophic marine bacterium Ruegeria pomeroyi

The heterotrophic marine bacterium, Ruegeria pomeroyi, was experimentally cultured under environmentally realistic carbon conditions and with a tracer-level addition of 13C-labeled leucine to track bacterial protein biosynthesis through growth phases. A combination of methods allowed observation of real-time bacterial protein production to understand metabolic priorities through the different growth phases. Over 2000 proteins were identified in each experimental culture from exponential and stationary growth phases. Within two hours of the 13C-labeled leucine addition, R. pomeroyi significantly assimilated the newly encountered substrate into new proteins. This dataset provides a fundamental baseline for understanding growth phase differences in molecular physiology of a cosmopolitan marine bacterium.

Pacific geoduck (Panopea generosa) resilience to natural pH variation

Pacific geoduck aquaculture is a growing industry, however, little is known about how geoduck respond to varying environmental conditions, or how the industry will fare under projected climate conditions. To understand how geoduck production may be impacted by low pH associated with ocean acidification, multi-faceted environmental heterogeneity needs to be included to understand species and community responses. In this study, eelgrass habitats and environmental heterogeneity across four estuarine bays were leveraged to examine low pH effects on geoduck under different natural regimes, using targeted proteomics to assess physiology. Juvenile geoduck were deployed in eelgrass and adjacent unvegetated habitats for 30 days while pH, temperature, dissolved oxygen, and salinity were monitored. Across the four bays, pH was lower in unvegetated habitats compared to eelgrass habitats. However this did not impact geoduck growth, survival, or proteomic abundance patterns in gill tissue. Temperature and dissolved oxygen differences across all locations corresponded to differences in growth and targeted protein abundance patterns. Specifically, three protein abundance levels (trifunctional-enzyme β-subunit, puromycin-sensitive aminopeptidase, and heat shock protein 90-α) and shell growth positively correlated with dissolved oxygen variability and inversely correlated with mean temperature. These results demonstrate that geoduck may be resilient to low pH in a natural setting, but other abiotic factors (i.e. temperature, dissolved oxygen variability) may have a greater influence on geoduck physiology. In addition this study contributes to the understanding of how eelgrass patches influences water chemistry.

Variations in Copepod Proteome and Respiration Rate in Association with Diel Vertical Migration and Circadian Cycle

The diel vertical migration of zooplankton is a process during which individuals spend the night in surface waters and retreat to depth during the daytime, with substantial implications for carbon transport and the ecology of midwater ecosystems. The physiological consequences of this daily pattern have, however, been poorly studied beyond investigations of speed and the energetic cost of swimming. Many other processes are likely influenced, such as fuel use, energetic trade-offs, underlying diel (circadian) rhythms, and antioxidant responses. Using a new reference transcriptome, proteomic analyses were applied to compare the physiological state of a migratory copepod, Pleuromamma xiphias, immediately after arriving to the surface at night and six hours later. Oxygen consumption was monitored semi-continuously to explore underlying cyclical patterns in metabolic rate under dark-dark conditions. The proteomic analysis suggests a distinct shift in physiology that reflects migratory exertion and changes in metabolism. These proteomic analyses are supported by the respiration experiments, which show an underlying cycle in metabolic rate, with a peak at dawn. This project generates molecular tools (transcriptome and proteome) that will allow for more detailed understanding of the underlying physiological processes that influence and are influenced by diel vertical migration. Further, these studies suggest that P. xiphias is a tractable model for continuing investigations of circadian and diel vertical migration influences on plankton physiology. Previous studies did not account for this cyclic pattern of respiration and may therefore have unrepresented respiratory carbon fluxes from copepods by about 24%.

MS analysis of a dilution series of bacteria:phytoplankton to improve detection of low abundance bacterial peptides

Assigning links between microbial activity and biogeochemical cycles in the ocean is a primary objective for ecologists and oceanographers. Bacteria represent a small ecosystem component by mass, but act as the nexus for both nutrient transformation and organic matter recycling. There are limited methods to explore the full suite of active bacterial proteins largely responsible for degradation. Mass spectrometry (MS)-based proteomics now has the potential to document bacterial physiology within these complex systems. Global proteome profiling using MS, known as data dependent acquisition (DDA), is limited by the stochastic nature of ion selection, decreasing the detection of low abundance peptides. The suitability of MS-based proteomics methods in revealing bacterial signatures outnumbered by phytoplankton proteins was explored using a dilution series of pure bacteria (Ruegeria pomeroyi) and diatoms (Thalassiosira pseudonana). Two common acquisition strategies were utilized: DDA and selected reaction monitoring (SRM). SRM improved detection of bacterial peptides at low bacterial cellular abundance that were undetectable with DDA from a wide range of physiological processes (e.g. amino acid synthesis, lipid metabolism, and transport). We demonstrate the benefits and drawbacks of two different proteomic approaches for investigating species-specific physiological processes across relative abundances of bacteria that vary by orders of magnitude.

Differential impacts of individual and combined exposures of deoxynivalenol and zearalenone on the HepaRG human hepatic cell proteome

Numerous surveys have highlighted the natural co-occurrence of deoxynivalenol (DON) and zearalenone (ZEA) mycotoxins in food and feed. Nevertheless, data regarding cellular mechanisms involved in response to their individual and simultaneous exposures are lacking. In this study, in order to analyze how low mycotoxin doses could impact cellular physiology and homeostasis, proteomic profiles of proliferating human hepatic cells (HepaRG) exposed for 1h and 24h to low DON and ZEA cytotoxicity levels (0.2 and 20μM respectively), alone or in combination, were analyzed by LC-MS/MS. Proteome analyses of mycotoxin-treated cells identified 4000 proteins with about 1.4% and 3.7% of these proteins exhibiting a significantly modified abundance compared to controls after 1h or 24h, respectively. Analysis of the Gene Ontology biological process annotations showed that cell cycle, proliferation and/or development as well as on DNA metabolic processes were affected for most treatments. Overall, different proteins, and thus biological processes, were impacted depending on the considered mycotoxin and exposure duration. Finally, despite the important proteome changes observed following 24h exposure to both mycotoxins, only the uptake of ZEA by the cells was suggested by the mycotoxin quantification in cell supernatants.

BIOLOGICAL SIGNIFICANCE

This study investigated the proteomic changes that occurred after DON and ZEA (individually and in combination) short exposures at low cytotoxicity levels in proliferating HepaRG cells using LC-MS/MS. The obtained results showed that the cellular response is time- and mycotoxin or mixture-dependent. In particular, after 1h exposure, the DON+ZEA combination led to more proteomic changes than DON or ZEA alone, whereas the opposite was observed after 24h. In addition, the significant cellular response to stress induced by ZEA after 24h exposure seemed to be reduced when combined with DON. Thus, these results supported a possible mitigation by the hepatocytes when exposed to the mycotoxin mixture for a long duration. These findings represent an essential step to further explore adaptive cell response to mycotoxin exposure using with more complex incubation kinetics and combining different "omics" tools. Moreover, as mycotoxin quantification in cell supernatants showed different behaviors for DON and ZEA, this also raises the question about how mycotoxins actually trigger the cell response.

Human-mediated evolution in a threatened species? Juvenile life-history changes in Snake River salmon

Evaluations of human impacts on Earth's ecosystems often ignore evolutionary changes in response to altered selective regimes. Freshwater habitats for Snake River fall Chinook salmon (SRFCS), a threatened species in the US, have been dramatically changed by hydropower development and other watershed modifications. Associated biological changes include a shift in juvenile life history: Historically essentially 100% of juveniles migrated to sea as subyearlings, but a substantial fraction have migrated as yearlings in recent years. In contemplating future management actions for this species should major Snake River dams ever be removed (as many have proposed), it will be important to understand whether evolution is at least partially responsible for this life-history change. We hypothesized that if this trait is genetically based, parents who migrated to sea as subyearlings should produce faster-growing offspring that would be more likely to reach a size threshold to migrate to sea in their first year. We tested this with phenotypic data for over 2,600 juvenile SRFCS that were genetically matched to parents of hatchery and natural origin. Three lines of evidence supported our hypothesis: (i) the animal model estimated substantial heritability for juvenile growth rate for three consecutive cohorts; (ii) linear modeling showed an association between juvenile life history of parents and offspring growth rate; and (iii) faster-growing juveniles migrated at greater speeds, as expected if they were more likely to be heading to sea. Surprisingly, we also found that parents reared a full year in a hatchery produced the fastest growing offspring of all-apparently an example of cross-generational plasticity associated with artificial propagation. We suggest that SRFCS is an example of a potentially large class of species that can be considered to be "anthro-evolutionary"-signifying those whose evolutionary trajectories have been profoundly shaped by altered selective regimes in human-dominated landscapes.

An Alignment-Free "Metapeptide" Strategy for Metaproteomic Characterization of Microbiome Samples Using Shotgun Metagenomic Sequencing

In principle, tandem mass spectrometry can be used to detect and quantify the peptides present in a microbiome sample, enabling functional and taxonomic insight into microbiome metabolic activity. However, the phylogenetic diversity constituting a particular microbiome is often unknown, and many of the organisms present may not have assembled genomes. In ocean microbiome samples, with particularly diverse and uncultured bacterial communities, it is difficult to construct protein databases that contain the bulk of the peptides in the sample without losing detection sensitivity due to the overwhelming number of candidate peptides for each tandem mass spectrum. We describe a method for deriving "metapeptides" (short amino acid sequences that may be represented in multiple organisms) from shotgun metagenomic sequencing of microbiome samples. In two ocean microbiome samples, we constructed site-specific metapeptide databases to detect more than one and a half times as many peptides as by searching against predicted genes from an assembled metagenome and roughly three times as many peptides as by searching against the NCBI environmental proteome database. The increased peptide yield has the potential to enrich the taxonomic and functional characterization of sample metaproteomes.

Proteomics of Colwellia psychrerythraea at subzero temperatures - a life with limited movement, flexible membranes and vital DNA repair

The mechanisms that allow psychrophilic bacteria to remain metabolically active at subzero temperatures result from form and function of their proteins. We present first proteomic evidence of physiological changes of the marine psychrophile Colwellia psychrerythraea 34H (Cp34H) after exposure to subzero temperatures (-1, and -10°C in ice) through 8 weeks. Protein abundance was compared between different treatments to understand the effects of temperature and time, independently and jointly, within cells transitioning to, and being maintained in ice. Parallel [3H]-leucine and [3H]-thymidine incubations indicated active protein and DNA synthesis to -10°C. Mass spectrometry-based proteomics identified 1763 proteins across four experimental treatments. Proteins involved in osmolyte regulation and polymer secretion were found constitutively present across all treatments, suggesting that they are required for metabolic success below 0°C. Differentially abundant protein groups indicated a reallocation of resources from DNA binding to DNA repair and from motility to chemo-taxis and sensing. Changes to iron and nitrogen metabolism, cellular membrane structures, and protein synthesis and folding were also revealed. By elucidating vital strategies during life in ice, this study provides novel insight into the extensive molecular adaptations that occur in cold-adapted marine organisms to sustain cellular function in their habitat.

Shotgun proteomics reveals physiological response to ocean acidification in Crassostrea gigas

BACKGROUND

Ocean acidification as a result of increased anthropogenic CO2 emissions is occurring in marine and estuarine environments worldwide. The coastal ocean experiences additional daily and seasonal fluctuations in pH that can be lower than projected end-of-century open ocean pH reductions. In order to assess the impact of ocean acidification on marine invertebrates, Pacific oysters (Crassostrea gigas) were exposed to one of four different p CO2 levels for four weeks: 400 μatm (pH 8.0), 800 μatm (pH 7.7), 1000 μatm (pH 7.6), or 2800 μatm (pH 7.3).

RESULTS

At the end of the four week exposure period, oysters in all four p CO2 environments deposited new shell, but growth rate was not different among the treatments. However, micromechanical properties of the new shell were compromised by elevated p CO2. Elevated p CO2 affected neither whole body fatty acid composition, nor glycogen content, nor mortality rate associated with acute heat shock. Shotgun proteomics revealed that several physiological pathways were significantly affected by ocean acidification, including antioxidant response, carbohydrate metabolism, and transcription and translation. Additionally, the proteomic response to a second stress differed with p CO2, with numerous processes significantly affected by mechanical stimulation at high versus low p CO2 (all proteomics data are available in the ProteomeXchange under the identifier PXD000835).

CONCLUSIONS

Oyster physiology is significantly altered by exposure to elevated p CO2, indicating changes in energy resource use. This is especially apparent in the assessment of the effects of p CO2 on the proteomic response to a second stress. The altered stress response illustrates that ocean acidification may impact how oysters respond to other changes in their environment. These data contribute to an integrative view of the effects of ocean acidification on oysters as well as physiological trade-offs during environmental stress.

Shotgun proteomics as a viable approach for biological discovery in the Pacific oyster

Shotgun proteomics offers an efficient means to characterize proteins in a complex mixture, particularly when sufficient genomic resources are available. In order to assess the practical application of shotgun proteomics in the Pacific oyster (Crassostrea gigas), liquid chromatography coupled with tandem mass spectrometry was used to characterize the gill proteome. Using information from the recently published Pacific oyster genome, 1043 proteins were identified. Biological samples (n = 4) and corresponding technical replicates (three) were similar in both specific proteins identified and expression, as determined by normalized spectral abundance factor. A majority of the proteins identified (703) were present in all biological samples. Functional analysis of the protein repertoire illustrates that these proteins represent a wide range of biological processes, supporting the dynamic function of the gill. These insights are important for understanding environmental influences on the oyster, because the gill tissue acts as the interface between the oyster and its environment. In silico analysis indicated that this sequencing effort identified a large proportion of the complete gill proteome. Together, these data demonstrate that shotgun sequencing is a viable approach for biological discovery and will play an important role in future studies of oyster physiology.

Characterization of genes involved in ceramide metabolism in the Pacific oyster (Crassostrea gigas)

Background

The lipid signaling molecule, ceramide, is a key component of the vertebrate stress response, however, there is limited information concerning its role in invertebrate species. In order to identify genes involved in ceramide metabolism in bivalve molluscs, Pacific oyster genomic resources were examined for genes associated with ceramide metabolism and signaling.

Results

Several genes were identified including full-length sequences characterized for serine palmitoyltransferase-1, 3-ketodihydrosphingosine reductase, acid ceramidase, and ceramide glucosyltransferase. Genes involved in ceramide synthesis and metabolism are conserved across taxa in both form and function. Expression analysis as assessed by quantitative PCR indicated all genes were expressed at high levels in gill tissue. The role of the ceramide pathway genes in the invertebrate stress response was also explored by measuring expression levels in adult oysters exposed to Vibrio vulnificus. Two genes demonstrated increased expression during the bacterial challenge: a gene involved in hydrolytic breakdown of ceramide (acid ceramidase) and a gene involved in de novo generation of ceramide (3-ketodihydrosphingosine reductase), suggesting a possible role of ceramide in the invertebrate stress and immune responses.

Conclusions

In silico and laboratory results support that Pacific oysters have the basic components of the ceramide metabolism pathway. These results also indicate that ceramide may have analogous functions in vertebrates and invertebrates. The gene expression pattern of acid ceramidase and 3-kethodihydrosphingosine reductase in response to bacterial exposure especially supports that ceramide and sphingolipid metabolism may be involved in the oyster’s stress and/or immune responses.