Proteogenomic reconstruction of organ-specific metabolic networks in an environmental sentinel species, the amphipod Gammarus fossarum.

Metabolic pathways are affected by the impacts of environmental contaminants underlying a large variability of toxic effects across different species. However, the systematic reconstruction of metabolic pathways remains limited in environmental sentinel species due to the lack of available genomic data in many taxa of animal diversity. In this study we used a multi-omics approach to reconstruct the most comprehensive map of metabolic pathways for a crustacean model in biomonitoring, the amphipod Gammarus fossarum in order to improve the knowledge of the metabolism of this sentinel species. We revisited the assembly of RNA-seq data by de novo approaches to reduce RNA contaminants and transcript redundancy. We also acquired extensive mass spectrometry shotgun proteomic data on several organs from a reference population of G. fossarum males and females to identify organ-specific metabolic profiles. The G. fossarum metabolic pathway reconstruction (available through the metabolic database GamfoCyc) was performed by adapting the genomic tool CycADS and we identified 377 pathways representing 7630 annotated enzymes, 2610 enzymatic reactions and the expression of 858 enzymes was experimentally validated by proteomics. To our knowledge, our analysis provides for the first time a systematic metabolic pathway reconstruction and the proteome profiles of these pathways at the organ level in this sentinel species. As an example, we show an elevated abundance in enzymes involved in ATP biosynthesis and fatty acid beta-oxidation indicative of the high-energy requirement of the gills, or the key anabolic and detoxification role of the hepatopancreatic caeca, as exemplified by the specific expression of the retinoid biosynthetic pathways and glutathione synthesis. In conclusion, the multi-omics data integration performed in this study provides new resources to investigate metabolic processes in crustacean amphipods and their role in mediating the effects of environmental contaminant exposures in sentinel species. SYNOPSIS: This study provide the first evidence that it is possible to combine multiple omics data to exhaustively describe the metabolic network of a model species in ecotoxicology, Gammarus fossarum, for which a reference genome is not yet available.

Organ-oriented proteogenomics functional atlas of three aquatic invertebrate sentinel species.

Proteogenomic methodologies have enabled the identification of protein sequences in wild species without annotated genomes, shedding light on molecular mechanisms affected by pollution. However, proteomic resources for sentinel species are limited, and organ-level investigations are necessary to expand our understanding of their molecular biology. This study presents proteomic resources obtained from proteogenomic analyses of key organs (hepatopancreas, gills, hemolymph) from three established aquatic sentinel invertebrate species of interest in ecotoxicological/ecological research and environmental monitoring: Gammarus fossarum, Dreissena polymorpha, and Palaemon serratus. Proteogenomic analyses identified thousands of proteins for each species, with over 90% of them being annotated to putative function. Functional analysis validated the relevance of the proteomic atlases by revealing similarities in functional annotation of catalogues of proteins across analogous organs in the three species, while deep contrasts between functional profiles are delimited across different organs in the same organism. These organ-level proteomic atlases are crucial for future research on these sentinel animals, aiding in the evaluation of aquatic environmental risks and providing a valuable resource for ecotoxicological studies.

Male Differentiation in the Marine Copepod Oithona nana Reveals the Development of a New Nervous Ganglion and Lin12-Notch-Repeat Protein-Associated Proteolysis.

Copepods are among the most numerous animals, and they play an essential role in the marine trophic web and biogeochemical cycles. The genus Oithona is described as having the highest density of copepods. The Oithona male paradox describes the activity states of males, which are obliged to alternate between immobile and mobile phases for ambush feeding and mate searching, respectively, while the female is less mobile and feeds less. To characterize the molecular basis of this sexual dimorphism, we combined immunofluorescence, genomics, transcriptomics, and protein-protein interaction approaches and revealed the presence of a male-specific nervous ganglion. Transcriptomic analysis showed male-specific enrichment for nervous system development-related transcripts. Twenty-seven Lin12-Notch Repeat domain-containing protein coding genes (LDPGs) of the 75 LDPGs identified in the genome were specifically expressed in males. Furthermore, some LDPGs coded for proteins with predicted proteolytic activity, and proteases-associated transcripts showed a male-specific enrichment. Using yeast double-hybrid assays, we constructed a protein-protein interaction network involving two LDPs with proteases, extracellular matrix proteins, and neurogenesis-related proteins. We also hypothesized possible roles of the LDPGs in the development of the lateral ganglia through helping in extracellular matrix lysis, neurites growth guidance, and synapses genesis.

Investigating population-scale allelic differential expression in wild populations of Oithona similis (Cyclopoida, Claus, 1866).

Acclimation allowed by variation in gene or allele expression in natural populations is increasingly understood as a decisive mechanism, as much as adaptation, for species evolution. However, for small eukaryotic organisms, as species from zooplankton, classical methods face numerous challenges. Here, we propose the concept of allelic differential expression at the population-scale (psADE) to investigate the variation in allele expression in natural populations. We developed a novel approach to detect psADE based on metagenomic and metatranscriptomic data from environmental samples. This approach was applied on the widespread marine copepod, Oithona similis, by combining samples collected during the Tara Oceans expedition (2009-2013) and de novo transcriptome assemblies. Among a total of 25,768 single nucleotide variants (SNVs) of O. similis, 572 (2.2%) were affected by psADE in at least one population (FDR < 0.05). The distribution of SNVs under psADE in different populations is significantly shaped by population genomic differentiation (Pearson r = 0.87, p = 5.6 × 10-30), supporting a partial genetic control of psADE. Moreover, a significant amount of SNVs (0.6%) were under both selection and psADE (p < .05), supporting the hypothesis that natural selection and psADE tends to impact common loci. Population-scale allelic differential expression offers new insights into the gene regulation control in populations and its link with natural selection.

Discovering millions of plankton genomic markers from the Atlantic Ocean and the Mediterranean Sea.

Comparison of the molecular diversity in all plankton populations present in geographically distant water columns may allow for a holistic view of the connectivity, isolation and adaptation of organisms in the marine environment. In this context, a large-scale detection and analysis of genomic variants directly in metagenomic data appeared as a powerful strategy for the identification of genetic structures and genes under natural selection in plankton. Here, we used discosnp++, a reference-free variant caller, to produce genetic variants from large-scale metagenomic data and assessed its accuracy on the copepod Oithona nana in terms of variant calling, allele frequency estimation and population genomic statistics by comparing it to the state-of-the-art method. discosnp ++ produces variants leading to similar conclusions regarding the genetic structure and identification of loci under natural selection. discosnp++ was then applied to 120 metagenomic samples from four size fractions, including prokaryotes, protists and zooplankton sampled from 39 tara Oceans sampling stations located in the Atlantic Ocean and the Mediterranean Sea to produce a new set of marine genomic markers containing more than 19 million of variants. This new genomic resource can be used by the community to relocate these markers on their plankton genomes or transcriptomes of interest. This resource will be updated with new marine expeditions and the increase of metagenomic data (availability: http://bioinformatique.rennes.inria.fr/taravariants/).

Chitin distribution in the Oithona digestive and reproductive systems revealed by fluorescence microscopy.

Among copepods, which are the most abundant animals on Earth, the genus Oithona is described as one of the most numerous and plays a major role in the marine food chain and biogeochemical cycles, particularly through the excretion of chitin-coated fecal pellets. Despite the morphology of several Oithona species is well known, knowledge of its internal anatomy and chitin distribution is still limited. To answer this problem, Oithona nana and O. similis individuals were stained by Wheat Germ Agglutinin-Fluorescein IsoThioCyanate (WGA-FITC) and DiAmidino-2-PhenylIndole (DAPI) for fluorescence microscopy observations. The image analyses allowed a new description of the organization and chitin content of the digestive and reproductive systems of Oithona male and female. Chitin microfibrils were found all along the digestive system from the stomach to the hindgut with a higher concentration at the peritrophic membrane of the anterior midgut. Several midgut shrinkages were observed and proposed to be involved in faecal pellet shaping and motion. Amorphous chitin structures were also found to be a major component of the ducts and seminal vesicles and receptacles. The rapid staining protocol we proposed allowed a new insight into the Oithona internal anatomy and highlighted the role of chitin in the digestion and reproduction. This method could be applied to a wide range of copepods in order to perform comparative anatomy analyses.

New insights into global biogeography, population structure and natural selection from the genome of the epipelagic copepod Oithona.

In the epipelagic ocean, the genus Oithona is considered as one of the most abundant and widespread copepods and plays an important role in the trophic food web. Despite its ecological importance, little is known about Oithona and cyclopoid copepods genomics. Therefore, we sequenced, assembled and annotated the genome of Oithona nana. The comparative genomic analysis integrating available copepod genomes highlighted the expansions of genes related to stress response, cell differentiation and development, including genes coding Lin12-Notch-repeat (LNR) domain proteins. The Oithona biogeography based on 28S sequences and metagenomic reads from the Tara Oceans expedition showed the presence of O. nana mostly in the Mediterranean Sea (MS) and confirmed the amphitropical distribution of Oithona similis. The population genomics analyses of O. nana in the Northern MS, integrating the Tara Oceans metagenomic data and the O. nana genome, led to the identification of genetic structure between populations from the MS basins. Furthermore, 20 loci were found to be under positive selection including four missense and eight synonymous variants, harbouring soft or hard selective sweep patterns. One of the missense variants was localized in the LNR domain of the coding region of a male-specific gene. The variation in the B-allele frequency with respect to the MS circulation pattern showed the presence of genomic clines between O. nana and another undefined Oithona species possibly imported through Atlantic waters. This study provides new approaches and results in zooplankton population genomics through the integration of metagenomic and oceanographic data.