Proteogenomics of Gammarus fossarum to document the reproductive system of amphipods

Electron-activated dissociation (EAD) for the complementary annotation of metabolites and lipids through data-dependent acquisition analysis and feature-based molecular networking, applied to the sentinel amphipod Gammarus fossarum

The past decades have marked the rise of metabolomics and lipidomics as the -omics sciences which reflect the most phenotypes in living systems. Mass spectrometry-based approaches are acknowledged for both quantification and identification of molecular signatures, the latter relying primarily on fragmentation spectra interpretation. However, the high structural diversity of biological small molecules poses a considerable challenge in compound annotation. Feature-based molecular networking (FBMN) combined with database searches currently sets the gold standard for annotation of large datasets. Nevertheless, FBMN is usually based on collision-induced dissociation (CID) data, which may lead to unsatisfying information. The use of alternative fragmentation methods, such as electron-activated dissociation (EAD), is undergoing a re-evaluation for the annotation of small molecules, as it gives access to additional fragmentation routes. In this study, we apply the performances of data-dependent acquisition mass spectrometry (DDA-MS) under CID and EAD fragmentation along with FBMN construction, to perform extensive compound annotation in the crude extracts of the freshwater sentinel organism Gammarus fossarum. We discuss the analytical aspects of the use of the two fragmentation modes, perform a general comparison of the information delivered, and compare the CID and EAD fragmentation pathways for specific classes of compounds, including previously unstudied species. In addition, we discuss the potential use of FBMN constructed with EAD fragmentation spectra to improve lipid annotation, compared to the classic CID-based networks. Our approach has enabled higher confidence annotations and finer structure characterization of 823 features, including both metabolites and lipids detected in G. fossarum extracts.

GroEL-Proteotyping of Bacterial Communities Using Tandem Mass Spectrometry

Profiling bacterial populations in mixed communities is a common task in microbiology. Sequencing of 16S small subunit ribosomal-RNA (16S rRNA) gene amplicons is a widely accepted and functional approach but relies on amplification primers and cannot quantify isotope incorporation. Tandem mass spectrometry proteotyping is an effective alternative for taxonomically profiling microorganisms. We suggest that targeted proteotyping approaches can complement traditional population analyses. Therefore, we describe an approach to assess bacterial community compositions at the family level using the taxonomic marker protein GroEL, which is ubiquitously found in bacteria, except a few obligate intracellular species. We refer to our method as GroEL-proteotyping. GroEL-proteotyping is based on high-resolution tandem mass spectrometry of GroEL peptides and identification of GroEL-derived taxa via a Galaxy workflow and a subsequent Python-based analysis script. Its advantage is that it can be performed with a curated and extendable sample-independent database and that GroEL can be pre-separated by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) to reduce sample complexity, improving GroEL identification while simultaneously decreasing the instrument time. GroEL-proteotyping was validated by employing it on a comprehensive raw dataset obtained through a metaproteome approach from synthetic microbial communities as well as real human gut samples. Our data show that GroEL-proteotyping enables fast and straightforward profiling of highly abundant taxa in bacterial communities at reasonable taxonomic resolution.

Development of a multi-omics extraction method for ecotoxicology: investigation of the reproductive cycle of Gammarus fossarum

Omics study exemplified by proteomics, lipidomics or metabolomics, provides the opportunity to get insight of the molecular modifications occurring in living organisms in response to contaminants or in different physiological conditions. However, individual omics discloses only a single layer of information leading to a partial image of the biological complexity. Multiplication of samples preparation and processing can generate analytical variations resulting from several extractions and instrumental runs. To get all the -omics information at the proteins, metabolites and lipids level coming from a unique sample, a specific sample preparation must be optimized. In this study, we streamlined a biphasic extraction procedure based on a MTBE/Methanol mixture to provide the simultaneous extraction of polar (proteins, metabolites) and apolar compounds (lipids) for multi-omics analyses from a unique biological sample by a liquid chromatography (LC)/mass spectrometry (MS)/MS-based targeted approach. We applied the methodology for the study of female amphipod Gammarus fossarum during the reproductive cycle. Multivariate data analyses including Partial Least Squares Discriminant Analysis and multiple factor analysis were applied for the integration of the multi-omics data sets and highlighted molecular signatures, specific to the different stages.

Interest of a multispecies approach in active biomonitoring: Application in the Meuse watershed

A biomonitoring approach based on a single model species cannot be representative of the contaminations impacts on the ecosystem overall. As part of the Interreg DIADeM program ("Development of an integrated approach for the diagnosis of the water quality of the River Meuse"), a study was conducted to establish the proof of concept that the use of a multispecies active biomonitoring approach improves diagnostic of aquatic systems. The complementarity of the biomarker responses was tested in four model species belonging to various ecological compartments: the bryophyte Fontinalis antipyretica, the bivalve Dreissena polymorpha, the amphipod Gammarus fossarum and the fish Gasterosteus aculeatus. The species have been caged upstream and downstream from five wastewater treatment plants (WWTPs) in the Meuse watershed. After the exposure, a battery of biomarkers was measured and results were compiled in an Integrated Biomarker Response (IBR) for each species. A multispecies IBR value was then proposed to assess the quality of the receiving environment upstream the WWTPs. The effluent toxicity was variable according to the caged species and the WWTP. However, the calculated IBR were high for all species and upstream sites, suggesting that the water quality was already downgraded upstream the WWTP. This contamination of the receiving environment was confirmed by the multispecies IBR which has allowed to rank the rivers from the less to the most contaminated. This study has demonstrated the interest of the IBR in the assessment of biological impacts of a point-source contamination (WWTP effluent) but also of the receiving environment, thanks to the use of independent references. Moreover, this study has highlighted the complementarity between the different species and has emphasized the interest of this multispecies approach to consider the variability of the species exposition pathway and sensibility as well as the mechanism of contaminants toxicity in the final diagnosis.

Subcellular Distribution of Dietary Methyl-Mercury in Gammarus fossarum and Its Impact on the Amphipod Proteome

The transfer of methyl-Hg (MeHg) from food is central for its effects in aquatic animals, but we still lack knowledge concerning its impact on invertebrate primary consumers. In aquatic environments, cell walls of plants are particularly recalcitrant to degradation and as such remain available as a food source for long periods. Here, the impact at the proteomic level of dietary MeHg in Gammarus fossarum was established and linked to subcellular distribution of Hg. Individuals of G. fossarum were fed with MeHg in cell wall or intracellular compartments of Elodea nuttallii. Hg concentrations in subcellular fractions were 2 to 6 times higher in animals fed with cell wall than intracellular compartments. At the higher concentrations tested, the proportion of Hg in metal-sensitive fraction increased from 30.0 ± 6.1 to 41.0 ± 5.7% for individuals fed with intracellular compartment, while biologically detoxified metal fraction increased from 30.0 ± 6.1 to 50.0 ± 2.8% when fed with cell wall compartment. Data suggested that several thresholds of proteomic response are triggered by increased bioaccumulation in each subcellular fraction in correlation with Hg exclusively bound to the metal-sensitive fraction, while the increase of biologically detoxified metal likely had a cost for fitness. Proteomics analysis supported that the different binding sites and speciation in shoots subsequently resulted in different fate and cellular toxicity pathways to consumers. Our data confirmed that Hg bound in cell walls of plants can be assimilated by G. fossarum, which is consistent with its feeding strategy, hence pointing cell walls as a significant source for Hg transfers and toxicity in primary consumers. The high accumulation of Hg in macrophytes makes them a risk for food web transfer in shallow ecosystems. The present results allowed gaining new insights into the effects and uptake mechanisms of MeHg in aquatic primary consumers.

Co-expression network analysis identifies novel molecular pathways associated with cadmium and pyriproxyfen testicular toxicity in Gammarus fossarum

Omics approaches are continuously providing new clues on the mechanisms of action of contaminants in species of environmental relevance, contributing to the emergence of molecular ecotoxicology. Co-expression network approaches represent a suitable methodological framework for studying the rich content of omics datasets. This study aimed to find evidence of key pathways and proteins related to the testicular toxicity in the sentinel crustacean species Gammarus fossarum exposed to endocrine disruptors using a weighted protein co-expression network analysis. From a shotgun proteomics dataset of male gonads of G. fossarum organisms exposed to cadmium (Cd), pyriproxyfen (Pyr) and methoxyfenozide (Met) in laboratory conditions, four distinct modules were identified as significantly correlated to contaminants' exposure. Protein set enrichment analysis identified modules involved in cytoskeleton organization and oxidative stress response associated with the Cd exposure. The module associated with Pyr exposure was associated with endoplasmic reticulum stress (ER) response, and the module correlated with Met exposure was characterized by a significant proportion of amphipod-restricted proteins whose functions are still not characterized. Our results show that co-expression networks are efficient and adapted tools to identify new potential mode of actions from environmental sentinel species, such as G. fossarum, using a proteogenomic approach, even without an annotated genome.

Proteomics reveals sex-specific heat shock response of Baikal amphipod Eulimnogammarus cyaneus

The ancient Lake Baikal is the largest source of liquid freshwater on Earth and home to a unique fauna. Several hundred mostly cold-adapted endemic amphipod species inhabit Baikal, an ecosystem that is already being influenced by global change. In this study, we characterized the core proteome and heat stress-induced changes in a temperature-tolerant endemic amphipod, Eulimnogammarus cyaneus, using a proteogenomic approach (PRIDE dataset PXD013237) to unravel the molecular mechanisms of the observed adverse effects. As males were previously found to be much more tolerant to thermal stress, we placed special emphasis on differences between the sexes. For both sexes, we observed adaption of energy metabolism, cytoskeleton, lipid, and carbohydrate metabolism upon heat stress. In contrast, significant differences were determined in the molecular chaperone response. Females from the control conditions possessed significantly higher levels of heat shock proteins (HSP70, HSPb1, Hsc70-3), which, in contrast to males, were not further increased in response to heat stress. The inability of females to further increase heat shock protein synthesis in response to temperature stress may be due to sex-specific processes, such as egg production, requiring a large proportion of the available energy. As ovigerous females synthesize generally higher amounts of protein, they also need higher levels of molecular chaperones for the folding of these new proteins. Thus, the higher sensitivity of females to heat shock may be due to the lack of molecular chaperone molecules to counteract the heat-induced protein denaturation.

Can Proteomics Be Considered as a Valuable Tool to Assess the Toxicity of Nanoparticles in Marine Bivalves?

Exposure to nanoparticles (NPs) has been identified as a major concern for marine ecosystems. Because of their peculiar physico-chemical features, NPs are accumulated in marine organisms, which suffer a variety of adverse effects. In particular, bivalve mollusks represent a unique target for NPs, mainly because they are suspension-feeders with highly developed processes for cellular internalization of nano- and micrometric particles. Several studies have demonstrated that the uptake and the accumulation of NPs can induce sub-lethal effects towards marine bivalves. However, to understand the real risk of NP exposures the application of the so-called “omics” techniques (e.g., proteomics, genomics, metabolomics, lipidomics) has been suggested. In particular, proteomics has been used to study the effects of NPs and their mechanism(s) of action in marine bivalves, but to date its application is still limited. The present review aims at summarizing the state of the art concerning the application of proteomics as a tool to investigate the effects of nanoparticles on the proteome of marine bivalves, and to critically discuss the advantages and limitations of proteomics in this field of research. Relying on results obtained by studies that applied proteomics on bivalve tissues, proteomics application needs to be considered cautiously as a promising and valuable tool to shed light on toxicity and mechanism(s) of action of NPs. Although on one hand, the analysis of the current literature demonstrated undeniable strengths, potentiality and reliability of proteomics, on the other hand a number of limitations suggest that some gaps of knowledge need to be bridged, and methodological and technical improvements are necessary before proteomics can be readily and routinely applied to nanotoxicology studies.

Transgenerational inheritance of chemical-induced signature: A case study with simvastatin

The hypothesis that exposure to certain environmental chemicals during early life stages may disrupt reproduction across multiple non-exposed generations has significant implications for understanding disease etiology and adverse outcomes. We demonstrate here reproductive multi and transgenerational effects, at environmentally relevant levels, of one of the most prescribed human pharmaceuticals, simvastatin, in a keystone species, the amphipod Gammarus locusta. The transgenerational findings has major implications for hazard and risk assessment of pharmaceuticals and other contaminants of emerging concern given that transgenerational effects of environmental chemicals are not addressed in current hazard and risk assessment schemes. Considering that the mevalonate synthesis, one of the key metabolic pathways targeted by simvastatin, is highly conserved among metazoans, these results may also shed light on the potential transgenerational effects of simvastatin on other animals, including humans.

Continental-scale patterns of hyper-cryptic diversity within the freshwater model taxon Gammarus fossarum (Crustacea, Amphipoda)

Traditional morphological diagnoses of taxonomic status remain widely used while an increasing number of studies show that one morphospecies might hide cryptic diversity, i.e. lineages with unexpectedly high molecular divergence. This hidden diversity can reach even tens of lineages, i.e. hyper cryptic diversity. Even well-studied model-organisms may exhibit overlooked cryptic diversity. Such is the case of the freshwater crustacean amphipod model taxon Gammarus fossarum. It is extensively used in both applied and basic types of research, including biodiversity assessments, ecotoxicology and evolutionary ecology. Based on COI barcodes of 4926 individuals from 498 sampling sites in 19 European countries, the present paper shows (1) hyper cryptic diversity, ranging from 84 to 152 Molecular Operational Taxonomic Units, (2) ancient diversification starting already 26 Mya in the Oligocene, and (3) high level of lineage syntopy. Even if hyper cryptic diversity was already documented in G. fossarum, the present study increases its extent fourfold, providing a first continental-scale insight into its geographical distribution and establishes several diversification hotspots, notably south-eastern and central Europe. The challenges of recording hyper cryptic diversity in the future are also discussed.

High-multiplexed monitoring of protein biomarkers in the sentinel Gammarus fossarum by targeted scout-MRM assay, a new vision for ecotoxicoproteomics

Ecotoxicoproteomics employs mass spectrometry-based approaches centered on proteins of sentinel organisms to assess for instance, chemical toxicity in fresh water. In this study, we combined proteogenomics experiments and a novel targeted proteomics approach free from retention time scheduling called Scout-MRM. This methodology will enable the measurement of simultaneously changes in the relative abundance of multiple proteins involved in key physiological processes and potentially impacted by contaminants in the freshwater sentinel Gammarus fossarum. The development and validation of the assay were performed to target 157 protein biomarkers of this non-model organism. We carefully chose and validated the transitions to monitor using conventional parameters (linearity, repeatability, LOD, LOQ). Finally, the potential of the methodology is illustrated by measuring 277-peptide-plex assay (831 transitions) in sentinel animals exposed in natura to different agricultural sites potentially exposed to pesticide contamination. Multivariate data analyses highlighted the modulation of several key proteins involved in feeding and molting. This multiplex-targeted proteomics assay paves the way for the discovery and the use of a large panel of novel protein biomarkers in emergent ecotoxicological models for environmental monitoring in the future. BIOLOGICAL SIGNIFICANCE: The study contributed to the development of Scout-MRM for the high-throughput quantitation of a large panel of proteins in the Gammarus fossarum freshwater sentinel. Increasing the number of markers in ecotoxicoproteomics is of most interest to assess the impact of pollutants in freshwater organisms. The development and validation of the assay enabled the monitoring of a large panel of reporter peptides of exposed gammarids. To illustrate the applicability of the methodology, animals from different agricultural sites were analysed. The application of the assay highlighted the modulation of some biomarker proteins involved in key physiological pathways, such as molting, feeding and general stress response. Increasing multiplexing capabilities and field test will provide the development of diagnostic protein biomarkers for emergent ecotoxicological models in future environmental biomonitoring programs.

Combining proteogenomics and metaproteomics for deep taxonomic and functional characterization of microbiomes from a non-sequenced host

Metaproteomics of gut microbiomes from animal hosts lacking a reference genome is challenging. Here we describe a strategy combining high-resolution metaproteomics and host RNA sequencing (RNA-seq) with generalist database searching to survey the digestive tract of Gammarus fossarum, a small crustacean used as a sentinel species in ecotoxicology. This approach provides a deep insight into the full range of biomasses and metabolic activities of the holobiont components, and differentiates between the intestine and hepatopancreatic caecum.

Proteogenomics-Guided Evaluation of RNA-Seq Assembly and Protein Database Construction for Emergent Model Organisms

Proteogenomics is gaining momentum as, today, genomics, transcriptomics, and proteomics can be readily performed on any new species. This approach allows key alterations to molecular pathways to be identified when comparing conditions. For animals and plants, RNA-seq-informed proteomics is the most popular means of interpreting tandem mass spectrometry spectra acquired for species for which the genome has not yet been sequenced. It relies on high-performance de novo RNA-seq assembly and optimized translation strategies. Here, several pre-treatments for Illumina RNA-seq reads before assembly are explored to translate the resulting contigs into useful polypeptide sequences. Experimental transcriptomics and proteomics datasets acquired for individual Gammarus fossarum freshwater crustaceans are used, the most relevant procedure is defined by the ratio of MS/MS spectra assigned to peptide sequences. Removing reads with a mean quality score of less than 17-which represents a single probable nucleotide error on 150-bp reads-prior to assembly, increases the proteomics outcome. The best translation using Transdecoder is achieved with a minimal open reading frame length of 50 amino acids and systematic selection of ORFs longer than 900 nucleotides. Using these parameters, transcriptome assembly and translation informed by proteomics pave the way to further improvements in proteogenomics.

Complete transcriptome assembly and annotation of a critically important amphipod species in freshwater ecotoxicological risk assessment: Gammarus fossarum

Because of their crucial role in ecotoxicological risk assessment, amphipods (Crustacea) are commonly employed as model species in a wide range of studies. However, despite their ecological importance, their genome has not yet been completely annotated and molecular mechanisms underlying key pathways, such as the serotonin pathway, in development of ecotoxicological biomarkers of exposure to neuroactive pharmaceuticals are still poorly understood. Furthermore, genetic similarities and discrepancies with other model arthropods (e.g., Drosophila melanogaster) have not been completely clarified. In this report, we present a new transcriptome assembly of Gammarus fossarum, an important amphipod species, widespread in Central Europe. RNA-Seq with Illumina HiSeq technology was used to analyse samples extracted from total internal tissues. We used the Trinity and Trinotate software suites for transcriptome assembly and annotation, respectively. The quality of this assembly and the affiliated targeted homology searches greatly enrich the molecular knowledge on this species. Because of the lack of publicly available molecular information on the serotonin pathway, we also highlighted sequence homologies and divergences of the genes encoding the serotonin pathway components of the well-annotated arthropod D. melanogaster, and Crustacea with the corresponding genes of our assembly. An inferior number of hits was found when running a BLAST analysis of both D. melanogaster and Crustacea mRNA sequences encoding serotonin receptors available in GenBank against the total assembly, compared to other serotonin pathway components. A lack of information on important components for serotonin biosynthesis and vesicle endocytosis (i.e., tryptophan hydroxylase and vesicular monoamine transporter) in Crustacea was also brought to light. Our results will provide an extensive transcriptional resource for this important species in ecotoxicological risk assessment and highlight the need for a more detailed categorization of neuronal pathways components in invertebrates.

De novo transcriptomes of 14 gammarid individuals for proteogenomic analysis of seven taxonomic groups

Gammarids are amphipods found worldwide distributed in fresh and marine waters. They play an important role in aquatic ecosystems and are well established sentinel species in ecotoxicology. In this study, we sequenced the transcriptomes of a male individual and a female individual for seven different taxonomic groups belonging to the two genera Gammarus and Echinogammarus: Gammarus fossarum A, G. fossarum B, G. fossarum C, Gammarus wautieri, Gammarus pulex, Echinogammarus berilloni, and Echinogammarus marinus. These taxa were chosen to explore the molecular diversity of transcribed genes of genotyped individuals from these groups. Transcriptomes were de novo assembled and annotated. High-quality assembly was confirmed by BUSCO comparison against the Arthropod dataset. The 14 RNA-Seq-derived protein sequence databases proposed here will be a significant resource for proteogenomics studies of these ecotoxicologically relevant non-model organisms. These transcriptomes represent reliable reference sequences for whole-transcriptome and proteome studies on other gammarids, for primer design to clone specific genes or monitor their specific expression, and for analyses of molecular differences between gammarid species.

Measurement(s)	transcription profiling assay
Technology Type(s)	RNA sequencing
Factor Type(s)	sex • species
Sample Characteristic - Organism	Gammarus • Echinogammarus
Sample Characteristic - Environment	habitat

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.9777905

Comparative proteomics in the wild: Accounting for intrapopulation variability improves describing proteome response in a Gammarus pulex field population exposed to cadmium

High-throughput proteomics can be performed on animal sentinels for discovering key molecular biomarkers signing the physiological response and adaptation of organisms. Ecotoxicoproteomics is today amenable by means of proteogenomics to small arthropods such as Gammarids which are well known sentinels of aquatic environments. Here, we analysed two regional Gammarus pulex populations to characterize the potential proteome divergence induced in one site by natural bioavailable mono-metallic contamination (cadmium) compared to a non-contaminated site. Two RNAseq-derived protein sequence databases were established previously on male and female individuals sampled from the reference site. Here, individual proteomes were acquired on 10 male and 10 female paired organisms sampled from each site. Proteins involved in protein lipidation, carbohydrate metabolism, proteolysis, innate immunity, oxidative stress response and lipid transport were found more abundant in animals exposed to cadmium, while hemocyanins were found in lower abundance. The intrapopulation proteome variability of long-term exposed G. pulex was inflated relatively to the non-contaminated population. These results show that, while remaining a challenge for such organisms with not yet sequenced genomes, taking into account intrapopulation variability is important to better define the molecular players induced by toxic stress in a comparative field proteomics approach.

Discovery of Novel Crustacean and Cephalopod Flaviviruses: Insights into the Evolution and Circulation of Flaviviruses between Marine Invertebrate and Vertebrate Hosts

Most described flaviviruses (family Flaviviridae) are disease-causing pathogens of vertebrates maintained in zoonotic cycles between mosquitoes or ticks and vertebrate hosts. Poor sampling of flaviviruses outside vector-borne flaviviruses such as Zika virus and dengue virus has presented a narrow understanding of flavivirus diversity and evolution. In this study, we discovered three crustacean flaviviruses (Gammarus chevreuxi flavivirus, Gammarus pulex flavivirus, and Crangon crangon flavivirus) and two cephalopod flaviviruses (Southern Pygmy squid flavivirus and Firefly squid flavivirus). Bayesian and maximum likelihood phylogenetic methods demonstrate that crustacean flaviviruses form a well-supported clade and share a more closely related ancestor with terrestrial vector-borne flaviviruses than with classical insect-specific flaviviruses. In addition, we identify variants of Wenzhou shark flavivirus in multiple gazami crab (Portunus trituberculatus) populations, with active replication supported by evidence of an active RNA interference response. This suggests that Wenzhou shark flavivirus moves horizontally between sharks and gazami crabs in ocean ecosystems. Analyses of the mono- and dinucleotide composition of marine flaviviruses compared to that of flaviviruses with known host status suggest that some marine flaviviruses share a nucleotide bias similar to that of vector-borne flaviviruses. Furthermore, we identify crustacean flavivirus endogenous viral elements that are closely related to elements of terrestrial vector-borne flaviviruses. Taken together, these data provide evidence of flaviviruses circulating between marine vertebrates and invertebrates, expand our understanding of flavivirus host range, and offer potential insights into the evolution and emergence of terrestrial vector-borne flaviviruses.IMPORTANCE Some flaviviruses are known to cause disease in vertebrates and are typically transmitted by blood-feeding arthropods such as ticks and mosquitoes. While an ever-increasing number of insect-specific flaviviruses have been described, we have a narrow understanding of flavivirus incidence and evolution. To expand this understanding, we discovered a number of novel flaviviruses that infect a range of crustaceans and cephalopod hosts. Phylogenetic analyses of these novel marine flaviviruses suggest that crustacean flaviviruses share a close ancestor to all terrestrial vector-borne flaviviruses, and squid flaviviruses are the most divergent of all known flaviviruses to date. Additionally, our results indicate horizontal transmission of a marine flavivirus between crabs and sharks. Taken together, these data suggest that flaviviruses move horizontally between invertebrates and vertebrates in ocean ecosystems. This study demonstrates that flavivirus invertebrate-vertebrate host associations have arisen in flaviviruses at least twice and may potentially provide insights into the emergence or origin of terrestrial vector-borne flaviviruses.

Co-expression network analysis identifies gonad- and embryo-associated protein modules in the sentinel species Gammarus fossarum

Next generation sequencing and mass spectrometry technologies have recently expanded the availability of whole transcriptomes and proteomes beyond classical model organisms in molecular biology, even in absence of an annotated genome. However, the fragmented nature of transcriptomic and proteomic data reduces the ability to interpret the data, notably in non-model organisms. Network-based approaches may help extracting important biological information from -omics datasets. The reproductive cycle of the freshwater crustacean Gammarus fossarum.provides an excellent case study to test the relevance of a network analysis in non-model organisms. Here, we illustrated how the use of a co-expression network analysis (based on Weighted Gene Co-expression Network Analysis algorithm, WGCNA) allowed identifying protein modules whose expression profiles described germ cell maturation and embryonic development in the freshwater crustacean Gammarus fossarum. Proteome datasets included testes, ovaries or embryos samples at different maturation or developmental stages, respectively. We identified an embryonic module correlated with mid-developmental stages corresponding to the organogenesis and it was characterized by enrichment in proteins involved in RNA editing and splicing. An ovarian module was enriched in vitellogenin-like proteins and clottable proteins, confirming the diversity of proteins belonging to the large lipid transfer family involved in oocytes maturations in this freshwater amphipod. Moreover, our results found evidence of a fine-tuned regulation between energy production by glycolysis and actin-myosin-dependent events in G. fossarum spermatogenesis. This study illustrates the importance of applying systems biology approaches to emergent animal models to improve the understanding of the molecular mechanisms regulating important physiological events with ecological relevance.

The immune system of the freshwater zebra mussel, Dreissena polymorpha, decrypted by proteogenomics of hemocytes and plasma compartments

The immune system of bivalves is of great interest since it reflects the health status of these organisms during stressful conditions. While immune molecular responses are well documented for marine bivalves, few information is available for continental bivalves such as the zebra mussel, Dreissena polymorpha. A proteogenomic approach was conducted on both hemocytes and plasma to identified immune proteins of this non-model species. Combining transcriptomic sequences with mass spectrometry data acquired on proteins is a relevant strategy since 3020 proteins were identified, representing the largest protein inventory for this sentinel organism. Functional annotation and gene ontology (GO) analysis performed on the identified proteins described the main molecular players of hemocytes and plasma in immunity. GO analysis highlights the complementary immune functions of these two compartments in the management of micro-organisms. Functional annotation revealed new mechanisms in the immune defence of the zebra mussel. Proteins rarely observed in the hemolymph of bivalves were pinpointed such as natterin-like and thaumatin-like proteins. Furthermore, the high abundance of complement-related proteins observed in plasma suggested a strong implication of the complement system in the immune defence of D. polymorpha. This work brings a better understanding of the molecular mechanisms involved in zebra mussel immunity. SIGNIFICANCE: Although the molecular mechanisms of marine bivalves are widely investigated, little information is known for continental bivalves. Moreover, few proteomic studies described the complementarity of both hemolymphatic compartments (cellular and plasmatic) in the immune defence of invertebrates. The recent proteogenomics concept made it possible to discover proteins in non-model organisms. Here, we propose a proteogenomic strategy with the zebra mussel, a key sentinel species for biomonitoring of freshwater, to identify and describe the molecular actors involved in the immune system in both hemocytes and plasma compartments. More widely, this study provided new insight into bivalve immunity.

Identification, expression, and endocrine-disruption of three ecdysone-responsive genes in the sentinel species Gammarus fossarum

Taking advantage of a large transcriptomic dataset recently obtained in the sentinel crustacean amphipod Gammarus fossarum, we developed an approach based on sequence similarity and phylogenetic reconstruction to identify key players involved in the endocrine regulation of G. fossarum. Our work identified three genes of interest: the nuclear receptors RXR and E75, and the regulator broad-complex (BR). Their involvement in the regulation of molting and reproduction, along with their sensitivity to chemical contamination were experimentally assessed by studying gene expression during the female reproductive cycle, and after laboratory exposure to model endocrine disrupting compounds (EDCs): pyriproxyfen, tebufenozide and piperonyl butoxide. RXR expression suggested a role of this gene in ecdysis and post-molting processes. E75 presented two expression peaks that suggested a role in vitellogenesis, and molting. BR expression showed no variation during molting/reproductive cycle. After exposure to the three EDCs, a strong inhibition of the inter-molt E75 peak was observed with tebufenozide, and an induction of RXR after exposure to pyriproxyfen and piperonyl butoxide. These results confirm the implication of RXR and E75 in hormonal regulation of female reproductive cycles in G. fossarum and their sensitivity towards EDCs opens the possibility of using them as specific endocrine disruption biomarkers.

Prediction of a peptidome for the ecotoxicological model Hyalella azteca (Crustacea; Amphipoda) using a de novo assembled transcriptome

Due to its sensitivity to many environmental and anthropogenic stressors, including a wide range of chemical compounds, Hyalella azteca, a freshwater amphipod, has emerged as one of the most commonly used invertebrates for ecotoxicological assessment.Peptidergic signaling systems are key components in the control of organism-environment interactions, and there is a growing literature suggesting that they are targets of a number of aquatic toxicants.Interestingly, and despite its model species status in the field of ecotoxicology, little is known about the peptide hormones of H. azteca.Here, a transcriptome was produced for this species using the de novo assembler Trinity and mined for sequences encoding putative peptide precursors; the transcriptome was assembled from 460,291,636 raw reads and consists of 133,486 unique transcripts.Seventy-six sequences encoding peptide pre/preprohormones were identified from this transcriptome, allowing for the prediction of 202 distinct peptides, which included members of the allatostatin A, allatostatin B, allatostatin C, allatotropin, bursicon, CCHamide, corazonin, crustacean cardioactive peptide, crustacean hyperglycemic hormone/molt-inhibiting hormone, ecdysis-triggering hormone, eclosion hormone, elevenin, FMRFamide-like peptide, glycoprotein hormone, GSEFLamide, inotocin, leucokinin, myosuppressin, neuropeptide F, orcokinin, orcomyotropin, pigment dispersing hormone, proctolin, pyrokinin, red pigment concentrating hormone, RYamide, short neuropeptide F, SIFamide, sulfakinin, tachykinin-related peptide and trissin families.These peptides expand the known peptidome for H. azteca approximately nine-fold, forming a strong foundation for future studies of peptidergic control, including disruption by aquatic toxicants, in this important ecotoxicological model.

Evidence of low dose effects of the antidepressant fluoxetine and the fungicide prochloraz on the behavior of the keystone freshwater invertebrate Gammarus pulex

In recent years, behavior-related endpoints have been proposed as rapid and reliable ecotoxicological tools for risk assessment. In particular, the use of detritivores to test the toxicity of pollutants through feeding is currently becoming a well-known method. Experiments combining feeding with other behavioral endpoints can provide relevant information about direct and indirect toxicological effects of chemicals. We carried out a feeding experiment with the shredder Gammarus pulex in order to detect indirect (through leaf conditioning) and direct effects (through water exposure) of two pollutants at environmentally relevant concentrations: the fungicide prochloraz (6 μg/L) and the antidepressant fluoxetine (100 ng/L). Prochloraz inhibited fungal growth on leaves, but it did not affect either the microbial breakdown rates or the C:N ratio of the leaves. Individuals of G. pulex that were fed with treated leaves presented lower consumption rates, not only those fed with prochloraz-treated leaves, but also those fed with fluoxetine-treated leaves, and those fed with the mixture-treated leaves. Mixed-effects models revealed that the swimming velocity of the amphipods after the experiment was modulated by the exposure to fluoxetine, and also by the exposure to prochloraz. We demonstrate that both the antidepressant and the fungicide may cause significant sublethal effects at low concentrations. The combination of behavioral endpoints together with the application of mixed models provided a useful tool for early detection of the effects of toxicity mixtures in freshwater ecosystems.

Ecotoxico-Proteomics for Aquatic Environmental Monitoring: First in Situ Application of a New Proteomics-Based Multibiomarker Assay Using Caged Amphipods

As a proof of principle, a selected reaction monitoring (SRM) mass spectrometry-based methodology was applied to the simultaneous quantification of dozens of protein biomarkers in caged amphipods (Gammarus fossarum). We evaluated the suitability of the methodology to assess complex field contaminations through its application in the framework of a regional river monitoring network. Thanks to the high throughput acquisition of biomarker levels in G. fossarum exposed in four reference and 13 contaminated sites, we analyzed the individual responses of 38 peptides reporting for 25 proteins of interest in 170 organisms. Responses obtained in contaminated sites included inductions of vitellogenin-like proteins in male organisms, inductions of Na⁺K⁺/ATPases, and strong inhibitions of molt-related proteins such as chitinase and JHE-carboxylesterase. Proteins from detoxification and immunity processes were also found modulated in abundance. Summarizing, the results presented here show that the SRM strategy developed for multibiomarker measurement paves a very promising way to define multiple indicators of the health status of sentinel organisms for environmental hazard assessment.

Assessing the relevance of a multiplexed methodology for proteomic biomarker measurement in the invertebrate species Gammarus fossarum: A physiological and ecotoxicological study

Recently, a protein sequence database was built specifically for the sentinel non-model species Gammarus fossarum using a proteogenomics approach. A quantitative multiplexed targeted proteomics assay (using Selected Reaction Monitoring mass spectrometry) was then developed for a fast and simultaneous quantification of dozens of biomarker peptides specific of this freshwater sentinel crustacean species. In order to assess the relevance of this breakthrough methodology in ecotoxicology, the response patterns of a panel of 26 peptides reporting for 20 proteins from the Gammarus fossarum proteome with putative key functional roles (homeostasis, osmoregulation, nutrition, reproduction, molting,…) were recorded through male and female reproductive cycles and after exposure to environmental concentrations of cadmium and lead in laboratory-controlled conditions. Based on these results, we validated the implication of annotated vtg-like peptides in the oogenesis process, and the implication of Na⁺/K⁺ ATPase proteins in the molt cycle of organisms. Upon metal (cadmium and lead) contamination, peptides belonging to proteins annotated as involved in antioxidant and detoxification functions, immunity and molting were significantly down-regulated. Overall, this multiplex assay allowed gaining relevant insights upon disruption of different main functions in the sentinel species Gammarus fossarum. This breakthrough methodology in ecotoxicology offers a valid and high throughput alternative to currently used protocols, paving the way for future practical applications of proteogenomics-derived protein biomarkers in chemical risk assessment and environmental monitoring.

Characterisation of the circulating acellular proteome of healthy sheep using LC-MS/MS-based proteomics analysis of serum

BACKGROUND

Unlike humans, there is currently no publicly available reference mass spectrometry-based circulating acellular proteome data for sheep, limiting the analysis and interpretation of a range of physiological changes and disease states. The objective of this study was to develop a robust and comprehensive method to characterise the circulating acellular proteome in ovine serum.

METHODS

Serum samples from healthy sheep were subjected to shotgun proteomic analysis using nano liquid chromatography nano electrospray ionisation tandem mass spectrometry (nanoLC-nanoESI-MS/MS) on a quadrupole time-of-flight instrument (TripleTOF® 5600+, SCIEX). Proteins were identified using ProteinPilot™ (SCIEX) and Mascot (Matrix Science) software based on a minimum of two unmodified highly scoring unique peptides per protein at a false discovery rate (FDR) of 1% software by searching a subset of the Universal Protein Resource Knowledgebase (UniProtKB) database (http://www.uniprot.org). PeptideShaker (CompOmics, VIB-UGent) searches were used to validate protein identifications from ProteinPilot™ and Mascot.

RESULTS

ProteinPilot™ and Mascot identified 245 and 379 protein groups (IDs), respectively, and PeptideShaker validated 133 protein IDs from the entire dataset. Since Mascot software is considered the industry standard and identified the most proteins, these were analysed using the Protein ANalysis THrough Evolutionary Relationships (PANTHER) classification tool revealing the association of 349 genes with 127 protein pathway hits. These data are available via ProteomeXchange with identifier PXD004989.

CONCLUSIONS

These results demonstrated for the first time the feasibility of characterising the ovine circulating acellular proteome using nanoLC-nanoESI-MS/MS. This peptide spectral data contributes to a protein library that can be used to identify a wide range of proteins in ovine serum.

Molecular adaptation in the world's deepest-living animal: Insights from transcriptome sequencing of the hadal amphipod Hirondellea gigas

The Challenger Deep in the Mariana Trench is the deepest point in the oceans of our planet. Understanding how animals adapt to this harsh environment characterized by high hydrostatic pressure, food limitation, dark and cold is of great scientific interest. Of the animals dwelling in the Challenger Deep, amphipods have been captured using baited traps. In this study, we sequenced the transcriptome of the amphipod Hirondellea gigas collected at a depth of 10,929 m from the East Pond of the Challenger Deep. Assembly of these sequences resulted in 133,041 contigs and 22,046 translated proteins. Functional annotation of these contigs was made using the go and kegg databases. Comparison of these translated proteins with those of four shallow-water amphipods revealed 10,731 gene families, of which 5659 were single-copy orthologs. Base substitution analysis on these single-copy orthologs showed that 62 genes are positively selected in H. gigas, including genes related to β-alanine biosynthesis, energy metabolism and genetic information processing. For multiple-copy orthologous genes, gene family expansion analysis revealed that cold-inducible proteins (i.e., transcription factors II A and transcription elongation factor 1) as well as zinc finger domains are expanded in H. gigas. Overall, our results indicate that genetic adaptation to the hadal environment by H. gigas may be mediated by both gene family expansion and amino acid substitutions of specific proteins.

Integrating transcriptomic and proteomic data for accurate assembly and annotation of genomes

Complementing genome sequence with deep transcriptome and proteome data could enable more accurate assembly and annotation of newly sequenced genomes. Here, we provide a proof-of-concept of an integrated approach for analysis of the genome and proteome of Anopheles stephensi, which is one of the most important vectors of the malaria parasite. To achieve broad coverage of genes, we carried out transcriptome sequencing and deep proteome profiling of multiple anatomically distinct sites. Based on transcriptomic data alone, we identified and corrected 535 events of incomplete genome assembly involving 1196 scaffolds and 868 protein-coding gene models. This proteogenomic approach enabled us to add 365 genes that were missed during genome annotation and identify 917 gene correction events through discovery of 151 novel exons, 297 protein extensions, 231 exon extensions, 192 novel protein start sites, 19 novel translational frames, 28 events of joining of exons, and 76 events of joining of adjacent genes as a single gene. Incorporation of proteomic evidence allowed us to change the designation of more than 87 predicted “noncoding RNAs” to conventional mRNAs coded by protein-coding genes. Importantly, extension of the newly corrected genome assemblies and gene models to 15 other newly assembled Anopheline genomes led to the discovery of a large number of apparent discrepancies in assembly and annotation of these genomes. Our data provide a framework for how future genome sequencing efforts should incorporate transcriptomic and proteomic analysis in combination with simultaneous manual curation to achieve near complete assembly and accurate annotation of genomes.

GAPP: A Proteogenomic Software for Genome Annotation and Global Profiling of Post-translational Modifications in Prokaryotes

Although the number of sequenced prokaryotic genomes is growing rapidly, experimentally verified annotation of prokaryotic genome remains patchy and challenging. To facilitate genome annotation efforts for prokaryotes, we developed an open source software called GAPP for genome annotation and global profiling of post-translational modifications (PTMs) in prokaryotes. With a single command, it provides a standard workflow to validate and refine predicted genetic models and discover diverse PTM events. We demonstrated the utility of GAPP using proteomic data from Helicobacter pylori, one of the major human pathogens that is responsible for many gastric diseases. Our results confirmed 84.9% of the existing predicted H. pylori proteins, identified 20 novel protein coding genes, and corrected four existing gene models with regard to translation initiation sites. In particular, GAPP revealed a large repertoire of PTMs using the same proteomic data and provided a rich resource that can be used to examine the functions of reversible modifications in this human pathogen. This software is a powerful tool for genome annotation and global discovery of PTMs and is applicable to any sequenced prokaryotic organism; we expect that it will become an integral part of ongoing genome annotation efforts for prokaryotes. GAPP is freely available at https://sourceforge.net/projects/gappproteogenomic/.

Ovary and embryo proteogenomic dataset revealing diversity of vitellogenins in the crustacean Gammarus fossarum

Ovaries and embryos from sexually mature Gammarus fossarum were sampled at different stages of the reproductive cycle. The soluble proteome was extracted for five biological replicates and samples were subjected to trypsin digestion. The resulting peptides were analyzed by high resolution tandem mass spectrometry with a LTQ-Orbitrap XL instrument. The MS/MS spectra were assigned with a previously described RNAseq-derived G. fossarum database. The proteins highlighted by proteogenomics were monitored and their abundance kinetics over the different stages revealed a large panel of vitellogenins. Criteria were i) accumulation during oogenesis, ii) decrease during embryogenesis, iii) classified as female-specific, and iv) sequence similarity and phylogenetic analysis. The data accompanying the manuscript describing the database searches and comparative analysis (“High-throughput proteome dynamics for discovery of key proteins in sentinel species: unsuspected vitellogenins diversity in the crustacean Gammarus fossarum” by Trapp et al. [1]) have been deposited to the ProteomeXchange via the PRIDE repository with identifiers PRIDE: PXD001002.

Tissue-specific Proteogenomic Analysis of Plutella xylostella Larval Midgut Using a Multialgorithm Pipeline

The diamondback moth, Plutella xylostella (L.), is the major cosmopolitan pest of brassica and other cruciferous crops. Its larval midgut is a dynamic tissue that interfaces with a wide variety of toxicological and physiological processes. The draft sequence of the P. xylostella genome was recently released, but its annotation remains challenging because of the low sequence coverage of this branch of life and the poor description of exon/intron splicing rules for these insects. Peptide sequencing by computational assignment of tandem mass spectra to genome sequence information provides an experimental independent approach for confirming or refuting protein predictions, a concept that has been termed proteogenomics. In this study, we carried out an in-depth proteogenomic analysis to complement genome annotation of P. xylostella larval midgut based on shotgun HPLC-ESI-MS/MS data by means of a multialgorithm pipeline. A total of 876,341 tandem mass spectra were searched against the predicted P. xylostella protein sequences and a whole-genome six-frame translation database. Based on a data set comprising 2694 novel genome search specific peptides, we discovered 439 novel protein-coding genes and corrected 128 existing gene models. To get the most accurate data to seed further insect genome annotation, more than half of the novel protein-coding genes, i.e. 235 over 439, were further validated after RT-PCR amplification and sequencing of the corresponding transcripts. Furthermore, we validated 53 novel alternative splicings. Finally, a total of 6764 proteins were identified, resulting in one of the most comprehensive proteogenomic study of a nonmodel animal. As the first tissue-specific proteogenomics analysis of P. xylostella, this study provides the fundamental basis for high-throughput proteomics and functional genomics approaches aimed at deciphering the molecular mechanisms of resistance and controlling this pest.

Proteogenomics for understanding oncology: recent advances and future prospects

The concept of proteogenomics has emerged rapidly as a valuable approach to integrate mass spectrometry-derived proteomic data with genomic and transcriptomic data. It is used to harness the full potential of the former dataset in the discovery of potential biomarkers, therapeutic targets and novel proteins associated with various biological processes including diseases. Proteogenomic strategies have been successfully utilized to identify novel genes and redefine annotation of existing gene models in various genomes. In recent years, this approach has been extended to the field of cancer biology to unravel complexities in the tumor genomes and proteomes. Standard proteomics workflows employing translated cancer genomes and transcriptomes can potentially identify peptides from mutant proteins, splice variants and fusion proteins in the tumor proteome, which in addition to the currently available biomarker panels can serve as potential diagnostic and prognostic biomarkers, besides having therapeutic utility. This review focuses on the role of proteogenomics to understand cancer biology.

Data for comparative proteomics of ovaries from five non-model, crustacean amphipods

Ovaries were taken from five sexually mature amphipods: Gammarus fossarum, Gammarus pulex, Gammarus roeseli, Hyallela azteca and Parhyale hawaiensis. The soluble proteome extracted from individual pair of ovaries from five biological replicates was trypsin digested and the resulting peptides were analyzed by high resolution tandem mass spectrometry. The spectra were assigned with four protein sequence databases with different specificities: a RNAseq-derived G. fossarum database; a RNAseq-derived P. hawaiensis database; both originating from ovaries transcriptome; the Daphnia pulex database derived from whole-genome sequencing and the NCBInr database. The best interpretation was obtained for most animals with the specific RNA-seq protein database previously established by means of RNAseq carried out on G. fossarum. Proteins identified in the five amphipod species allow defining the core-proteome of female reproductive tissues of the Senticaudata suborder. The data accompanying the manuscript describing the database searches and comparative analysis Trapp et al., 2015 [1] have been deposited to the ProteomeXchange with identifiers PXD002253 (G. fossarum), PXD002254 (G. pulex), PXD002255 (G. roeseli), PXD002256 (H. Azteca), and PXD002257 (P. hawaiensis).

Proteogenomic insights into the core-proteome of female reproductive tissues from crustacean amphipods

As a result of the poor genome sequence coverage of crustacean amphipods, characterization of their evolutionary biology relies mostly on phenotypic traits. Here, we analyzed the proteome of ovaries from five amphipods, all from the Senticaudata suborder, with the objective to obtain insights into the core-proteome of female reproductive systems. These amphipods were from either the Gammarida infraorder: Gammarus fossarum, Gammarus pulex, Gammarus roeseli, or the Talitrida infraorder: Parhyale hawaiensis and Hyalella azteca. Ovaries from animals sampled at the end of their reproductive cycle were dissected. Their whole protein contents were extracted and their proteomes were recorded by high-throughput nanoLC-MS/MS with a high-resolution mass spectrometer. We interpreted tandem mass spectrometry data with the protein sequence resource from G. fossarum and P. hawaiensis, both recently established by RNA sequencing. The large molecular biodiversity within amphipods was assessed by the ratio of MS/MS spectra assigned for each sample, which tends to diverge rapidly along the taxonomic level considered. The core-proteome was defined as the proteins conserved along all samples, thus detectable by the homology-based proteomic assignment procedure. This specific subproteome may be further enriched in the future with the analysis of new species and update of the protein sequence resource.