High-throughput transcriptome sequencing of the cold seep mussel Bathymodiolus platifrons

Potential biomarkers of metal toxicity in deep-sea invertebrates - A critical review of the omics data

Deep-sea mining (DSM) activities are expected to release potentially toxic metal mixtures through the generation of sediment plumes to the marine environment. This may disrupt the normal functioning of biological mechanisms, adversely affecting deep-sea invertebrate organisms. It is thus essential to understand the ecotoxicological effects from these toxic elements in deep-sea organisms and the omics approaches applied to ecotoxicology are seen as promising tools. Here, we provide an overview of the principal biological modifications identified in deep-sea invertebrates when exposed to metals and critically evaluate the current knowledge and discuss which potential biomarkers may be useful after metal exposure. Most of the 50 omics studies on deep-sea invertebrates revised are comparative transcriptomes (n = 41). Forty-three potential biomarker candidates are highlighted from immune system, 46 from cellular metabolism and 29 from oxidative stress. The processes mostly affected by metal toxicity in deep-sea invertebrates are related to innate immune defense; sulfur, chitin, and catabolic metabolism; antioxidation; and detoxification. We acknowledge the current limitations and future perspectives for their uses and emphasize the need to invest in further ecotoxicological studies using the omics approaches.

Deciphering deep-sea chemosynthetic symbiosis by single-nucleus RNA-sequencing

Bathymodioline mussels dominate deep-sea methane seep and hydrothermal vent habitats and obtain nutrients and energy primarily through chemosynthetic endosymbiotic bacteria in the bacteriocytes of their gill. However, the molecular mechanisms that orchestrate mussel host-symbiont interactions remain unclear. Here, we constructed a comprehensive cell atlas of the gill in the mussel Gigantidas platifrons from the South China Sea methane seeps (1100 m depth) using single-nucleus RNA-sequencing (snRNA-seq) and whole-mount in situ hybridisation. We identified 13 types of cells, including three previously unknown ones, and uncovered unknown tissue heterogeneity. Every cell type has a designated function in supporting the gill's structure and function, creating an optimal environment for chemosynthesis, and effectively acquiring nutrients from the endosymbiotic bacteria. Analysis of snRNA-seq of in situ transplanted mussels clearly showed the shifts in cell state in response to environmental oscillations. Our findings provide insight into the principles of host-symbiont interaction and the bivalves' environmental adaption mechanisms.

To live or die: "Fine-tuning" adaptation revealed by systemic analyses in symbiotic bathymodiolin mussels from diverse deep-sea extreme ecosystems

Hydrothermal vents (HVs) and cold seeps (CSs) are typical deep-sea extreme ecosystems with their own geochemical characteristics to supply the unique living conditions for local communities. Once HVs or CSs stop emission, the dramatic environmental change would pose survival risks to deep-sea organisms. Up to now, limited knowledge has been available to understand the biological responses and adaptive strategy to the extreme environments and their transition from active to extinct stage, mainly due to the technical difficulties and lack of representative organisms. In this study, bathymodiolin mussels, the dominant and successful species surviving in diverse deep-sea extreme ecosystems, were collected from active and extinct HVs (Southwest Indian Ocean) or CSs (South China Sea) via two individual cruises. The transcriptomic analysis and determination of multiple biological indexes in stress defense and metabolic systems were conducted in both gills and digestive glands of mussels, together with the metagenomic analysis of symbionts in mussels. The results revealed the ecosystem- and tissue-specific transcriptional regulation in mussels, addressing the autologous adaptations in antioxidant defense, energy utilization and key compounds (i.e. sulfur) metabolism. In detail, the successful antioxidant defense contributed to conquering the oxidative stress induced during the unavoidable metabolism of xenobiotics commonly existing in the extreme ecosystems; changes in metabolic rate functioned to handle toxic matters in different surroundings; upregulated gene expression of sulfide:quinone oxidoreductase indicated an active sulfide detoxification in mussels from HVs and active stage of HVs & CSs. Coordinately, a heterologous adaptation, characterized by the functional compensation between symbionts and mussels in energy utilization, sulfur and carbon metabolism, was also evidenced by the bacterial metagenomic analysis. Taken together, a new insight was proposed that symbiotic bathymodiolin mussels would develop a "finetuning" strategy combining the autologous and heterologous regulations to fulfill the efficient and effective adaptations for successful survival.

Chromosome-level genome assembly of the deep-sea snail Phymorhynchus buccinoides provides insights into the adaptation to the cold seep habitat

BACKGROUND

The deep-sea snail Phymorhynchus buccinoides belongs to the genus Phymorhynchus (Neogastropoda: Raphitomidae), and it is a dominant specie in the cold seep habitat. As the environment of the cold seep is characterized by darkness, hypoxia and high concentrations of toxic substances such as hydrogen sulfide (H2S), exploration of the diverse fauna living around cold seeps will help to uncover the adaptive mechanisms to this unique habitat. In the present study, a chromosome-level genome of P. buccinoides was constructed and a series of genomic and transcriptomic analyses were conducted to explore its molecular adaptation mechanisms to the cold seep environments.

RESULTS

The assembled genome size of the P. buccinoides was approximately 2.1 Gb, which is larger than most of the reported snail genomes, possibly due to the high proportion of repetitive elements. About 92.0% of the assembled base pairs of contigs were anchored to 34 pseudo-chromosomes with a scaffold N50 size of 60.0 Mb. Compared with relative specie in the shallow water, the glutamate regulative and related genes were expanded in P. buccinoides, which contributes to the acclimation to hypoxia and coldness. Besides, the relatively high mRNA expression levels of the olfactory/chemosensory genes in osphradium indicate that P. buccinoides might have evolved a highly developed and sensitive olfactory organ for its orientation and predation. Moreover, the genome and transcriptome analyses demonstrate that P. buccinoides has evolved a sulfite-tolerance mechanism by performing H2S detoxification. Many genes involved in H2S detoxification were highly expressed in ctenidium and hepatopancreas, suggesting that these tissues might be critical for H2S detoxification and sulfite tolerance.

CONCLUSIONS

In summary, our report of this chromosome-level deep-sea snail genome provides a comprehensive genomic basis for the understanding of the adaptation strategy of P. buccinoides to the extreme environment at the deep-sea cold seeps.

Capacity of deep-sea corals to obtain nutrition from cold seeps aligned with microbiome reorganization

Cold seeps in the deep sea harbor various animals that have adapted to utilize seepage chemicals with the aid of chemosynthetic microbes that serve as primary producers. Corals are among the animals that live near seep habitats and yet, there is a lack of evidence that corals gain benefits and/or incur costs from cold seeps. Here, we focused on Callogorgia delta and Paramuricea sp. type B3 that live near and far from visual signs of currently active seepage at five sites in the deep Gulf of Mexico. We tested whether these corals rely on chemosynthetically-derived food in seep habitats and how the proximity to cold seeps may influence; (i) coral colony traits (i.e., health status, growth rate, regrowth after sampling, and branch loss) and associated epifauna, (ii) associated microbiome, and (iii) host transcriptomes. Stable isotope data showed that many coral colonies utilized chemosynthetically derived food, but the feeding strategy differed by coral species. The microbiome composition of C. delta, unlike Paramuricea sp., varied significantly between seep and non-seep colonies and both coral species were associated with various sulfur-oxidizing bacteria (SUP05). Interestingly, the relative abundances of SUP05 varied among seep and non-seep colonies and were strongly correlated with carbon and nitrogen stable isotope values. In contrast, the proximity to cold seeps did not have a measurable effect on gene expression, colony traits, or associated epifauna in coral species. Our work provides the first evidence that some corals may gain benefits from living near cold seeps with apparently limited costs to the colonies. Cold seeps provide not only hard substrate but also food to cold-water corals. Furthermore, restructuring of the microbiome communities (particularly SUP05) is likely the key adaptive process to aid corals in utilizing seepage-derived carbon. This highlights that those deep-sea corals may upregulate particular microbial symbiont communities to cope with environmental gradients.

Transcriptome analysis of gene expression profiling from the deep sea in situ to the laboratory for the cold seep mussel Gigantidas haimaensis

Background

The deep-sea mussel Gigantidas haimaensis is a representative species from the Haima cold seep ecosystem in the South China Sea that establishes endosymbiosis with chemotrophic bacteria. During long-term evolution, G. haimaensis has adapted well to the local environment of cold seeps. Until now, adaptive mechanisms responding to environmental stresses have remained poorly understood.

Results

In this study, transcriptomic analysis was performed for muscle tissue of G. haimaensis in the in situ environment (MH) and laboratory environment for 0 h (M0), 3 h (M3) and 9 h (M9), and 187,368 transcript sequences and 22,924 annotated differentially expressed genes (DEGs) were generated. Based on Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, these DEGs were enriched with a broad spectrum of biological processes and pathways, including those associated with antioxidants, apoptosis, chaperones, immunity and metabolism. Among these significantly enriched pathways, protein processing in the endoplasmic reticulum and metabolism were the most affected metabolic pathways. These results may imply that G. haimaensis struggles to support the life response to environmental change by changing gene expression profiles.

Conclusion

The present study provides a better understanding of the biological responses and survival strategies of the mussel G. haimaensis from deep sea in situ to the laboratory environment.

Molecular characterization of peptidoglycan recognition proteins from Mytilus coruscus

Mytilus shows great immune resistance to various bacteria from the living waters, indicating a complex immune recognition mechanism against various microbes. Peptidoglycan recognition proteins (PGRPs) play an important role in the defense against invading microbes via the recognition of the immunogenic substance peptidoglycan (PGN). Therefore, eight PGRPs were identified from the gill transcriptome of Mytilus coruscus. The sequence features, expression pattern in various organs and larval development stages, and microbes induced expression profiles of these Mytilus PGRPs were determined. Our data revealed the constitutive expression of PGRPs in various organs with relative higher expression level in immune-related organs. The expression of PGRPs is developmentally regulated, and most PGRPs are undetectable in larvae stages. The expression level of most PGRPs was significantly increased with in vivo microbial challenges, showing strong response to Gram-positive strain in gill and digestive gland, strong response to Gram-negative strain in hemocytes, and relative weaker response to fungus in the three tested organs. In addition, the function analysis of the representative recombinant expressed PGRP (rMcPGRP-2) confirmed the antimicrobial and agglutination activities, showing the immune-related importance of PGRP in Mytilus. Our work suggests that Mytilus PGRPs can act as pattern recognition receptors to recognize the invading microorganisms and the antimicrobial effectors during the innate immune response of Mytilus.

Adaptation and evolution of the sea anemone Alvinactis sp. to deep-sea hydrothermal vents: A comparison using transcriptomes

Sea anemones are diverse and ecologically successful members of Anthozoa. They are often found in intertidal and shallow waters, although a few of them inhabit harsher living conditions, such as deep-sea hydrothermal vents. Here, we sequenced the transcriptome of the vent sea anemone Alvinactis sp., which was collected from Edmond vent along the central Indian Ocean ridge at a depth of 3275 m, to explore the molecular mechanisms related to adaptation to vents. Compared with another deep-sea anemone (Paraphelliactis xishaensis) and five shallow water sea anemones, a total of 117 positively selected genes and 46 significantly expanded gene families were found in Alvinactis sp. specifically that may be related to its vent-specific aspect of adaptation. In addition, 127 positively selected genes and 23 significantly expanded gene families that were found in both Alvinactis sp. and P. xishaensis. Among these, vent-specific adaptations of Alvinactis sp. may involve genetic alterations in peroxisome, ubiquitin-mediated protein degradation, oxidative phosphorylation, and endocytosis, and its deep-sea adaptation may involve changes in genetic information processing. Differentially expressed genes between Alvinactis sp. and the deep-sea anemone P. xishaensis were enriched in a variety of pathways related to adaptation, such as energy metabolism, genetic information processing, endocytosis, and peroxisomes. Overall, we provided the first transcriptome of sea anemones that inhabit vents, which enriches our knowledge of deep-sea hydrothermal vent adaptation and the diversity of sea anemones.

Chromosome-Level Genome Assembly of the Butter Clam Saxidomus purpuratus

Herein, we provide the first whole-genome sequence of the purple butter clam (Saxidomus purpuratus), an economically important bivalve shellfish. Specifically, we sequenced and de novo assembled the genome of Sa. purpuratus based on PromethION long reads and Hi-C data. The 978-Mb genome of Sa. purpuratus comprises 19 chromosomes with 36,591 predicted protein-coding genes. The N50 length of Sa. purpuratus genome is 52 Mb, showing the highest continuous assembly among bivalve genomes. The Benchmarking by Universal Single-Copy Orthologs assessment indicated that 95.07% of complete metazoan universal single-copy orthologs (n = 954) were present in the assembly. Approximately 51% of Sa. purpuratus genome comprises repetitive sequences. Based on the high-quality Sa. purpuratus genome, we resolved half of the immune-associated genes, namely, scavenger receptor (SR) proteins, which are collinear to those in the closely related Cyclina sinensis genome. This finding suggested a high degree of conservation among immune-associated genes. Twenty-two (19%) SR proteins are tandemly duplicated in Sa. purpuratus genome, suggesting putative convergence evolution. Overall, Sa. purpuratus genome provides a new resource for the discovery of economically important traits and immune-response genes.

Life in the Dark: Phylogenetic and Physiological Diversity of Chemosynthetic Symbioses

Possibly the last discovery of a previously unknown major ecosystem on Earth was made just over half a century ago, when researchers found teaming communities of animals flourishing two and a half kilometers below the ocean surface at hydrothermal vents. We now know that these highly productive ecosystems are based on nutritional symbioses between chemosynthetic bacteria and eukaryotes and that these chemosymbioses are ubiquitous in both deep-sea and shallow-water environments. The symbionts are primary producers that gain energy from the oxidation of reduced compounds, such as sulfide and methane, to fix carbon dioxide or methane into biomass to feed their hosts. This review outlines how the symbiotic partners have adapted to living together. We first focus on the phylogenetic and metabolic diversity of these symbioses and then highlight selected research directions that could advance our understanding of the processes that shaped the evolutionary and ecological success of these associations. Expected final online publication date for the Annual Review of Microbiology, Volume 75 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Biochemical and metabolic responses of the deep-sea mussel Bathymodiolus platifrons to cadmium and copper exposure

Greater interest in commercial deep-sea mining has been accompanied by mounting environmental concerns, including metal contamination resulting from mining activities. However, little is known about the toxic effects of metal exposure on deep-sea life. Given its ability to accumulate metals from the surrounding environment, its wide distribution at both vents and seeps, and its high abundance, the deep-sea mussel Bathymodiolus platifrons could serve as an ideal model to investigate the toxicological responses of deep-sea organisms to metal exposure. Here, we evaluated metal accumulation, traditional metal-related biomarkers, namely acid phosphatase (ACP), alkaline phosphatase (AKP), superoxide dismutase, catalase, reduced glutathione, metallothioneins, and malondialdehyde, as well as metabolic profiles in the gills of B. platifrons after a 7-day exposure to copper (100 μg/L), cadmium (500 μg/L), or copper-plus-cadmium treatments (100 μg/L Cu and 500 μg/L Cd). Metal exposure concentrations selected in this study can be found in deep-sea hydrothermal environments. Metal exposure resulted in significant metal accumulation in the gills of the mussel, indicating that B. platifrons has promise for use as an indicator of deep-sea metal pollution levels. Traditional biomarkers (AKP, ACP, and measured antioxidants) revealed cellular injury and oxidative stress in mussels following metal exposure. Metabolic responses in the three treatment groups indicated that metal exposure perturbed osmoregulation, energy metabolism, and nucleotide metabolism in mussels, in a response marked by differentially altered levels of amino acids, hypotaurine, betaine, succinate, glucose 6-phosphate, fructose 6-phosphate, guanosine, guanosine 5'-monophosphate, and inosine. Nevertheless, several uniquely altered metabolites were found in each treatment exposure group, suggesting dissimilar modes of toxicity between the two metal types. In the Cd-exposed group, the monosaccharide D-allose, which is involved in suppressing mitochondrial ROS production, was downregulated, a response consistent with oxidative stress in Cd-exposed B. platifrons. In the Cu-exposed group, the detected alterations in dopamine, dopamine-related, and serotonin-related metabolites together suggest disturbed neurotransmission in Cu-exposed B. platifrons. In the Cu-plus-Cd group, we detected a decline in fatty acid levels, implying that exposure to both metals jointly exerted a negative influence on the physiological functioning of the mussel. To the best of our knowledge, this is the first study to investigate changes in metabolite profiles in Bathymodiolus mussels exposed to metal. The findings reported here advance our understanding of the adverse impact of metal exposure on deep-sea life and can inform deep-sea mining assessments through the use of multiple biomarkers.

A Toll-like receptor identified in Gigantidas platifrons and its potential role in the immune recognition of endosymbiotic methane oxidation bacteria

Symbiosis with chemosynthetic bacteria is an important ecological strategy for the deep-sea megafaunas including mollusks, tubeworms and crustacean to obtain nutrients in hydrothermal vents and cold seeps. How the megafaunas recognize symbionts and establish the symbiosis has attracted much attention. Bathymodiolinae mussels are endemic species in both hydrothermal vents and cold seeps while the immune recognition mechanism underlying the symbiosis is not well understood due to the nonculturable symbionts. In previous study, a lipopolysaccharide (LPS) pull-down assay was conducted in Gigantidas platifrons to screen the pattern recognition receptors potentially involved in the recognition of symbiotic methane-oxidizing bacteria (MOB). Consequently, a total of 208 proteins including GpTLR13 were identified. Here the molecular structure, expression pattern and immune function of GpTLR13 were further analyzed. It was found that GpTLR13 could bind intensively with the lipid A structure of LPS through surface plasmon resonance analysis. The expression alternations of GpTLR13 transcripts during a 28-day of symbiont-depletion assay were investigated by real-time qPCR. As a result, a robust decrease of GpTLR13 transcripts was observed accompanying with the loss of symbionts, implying its participation in symbiosis. In addition, GpTLR13 transcripts were found expressed exclusively in the bacteriocytes of gills of G. platifrons by in situ hybridization. It was therefore speculated that GpTLR13 may be involved in the immune recognition of symbiotic methane-oxidizing bacteria by specifically recognizing the lipid A structure of LPS. However, the interaction between GpTLR13 and symbiotic MOB was failed to be addressed due to the nonculturable symbionts. Nevertheless, the present result has provided with a promising candidate as well as a new approach for the identification of symbiont-related genes in Bathymodiolinae mussels.

Molecular analyses of the gill symbiosis of the bathymodiolin mussel Gigantidas platifrons

Although the deep-sea bathymodiolin mussels have been intensively studied as a model of animal-bacteria symbiosis, it remains challenging to assess the host-symbiont interactions due to the complexity of the symbiotic tissue-the gill. Using cold-seep mussel Gigantidas platifrons as a model, we isolated the symbiont harboring bacteriocytes and profiled the transcriptomes of the three major parts of the symbiosis-the gill, the bacteriocyte, and the symbiont. This breakdown of the complex symbiotic tissue allowed us to characterize the host-symbiont interactions further. Our data showed that the gill's non-symbiotic parts play crucial roles in maintaining and protecting the symbiosis; the bacteriocytes supply the symbiont with metabolites, control symbiont population, and shelter the symbiont from phage infection; the symbiont dedicates to the methane oxidation and energy production. This study demonstrates that the bathymodiolin symbiosis interacts at the tissue, cellular, and molecular level, maintaining high efficiency and harmonic chemosynthetic micro niche.

The complete mitochondrial genome of a cold seep gastropod Phymorhynchus buccinoides (Neogastropoda: Conoidea: Raphitomidae)

Phymorhynchus is a genus of deep-sea snails that are most distributed in hydrothermal vent or cold seep environments. In this study, we presented the complete mitochondrial genome of P. buccinoides, a cold seep snail from the South China Sea. It is the first mitochondrial genome of a cold seep member of the superfamily Conoidea. The mitochondrial genome is 15,764 bp in length, and contains 13 protein-coding genes (PCGs), 2 rRNA genes, and 22 tRNA genes. These genes are encoded on the positive strand, except for 8 tRNA genes that are encoded on the negative strand. The start codon ATG and 3 types of stop codons, TAA, TAG and the truncated termination codon T, are used in the 13 PCGs. All 13 PCGs in the 26 species of Conoidea share the same gene order, while several tRNA genes have been translocated. Phylogenetic analysis revealed that P. buccinoides clustered with Typhlosyrinx sp., Eubela sp., and Phymorhynchus sp., forming the Raphitomidae clade, with high support values. Positive selection analysis showed that a residue located in atp6 (18 S) was identified as the positively selected site with high posterior probabilities, suggesting potential adaption to the cold seep environment. Overall, our data will provide a useful resource on the evolutionary adaptation of cold seep snails for future studies.

Metal adaptation strategies of deep-sea Bathymodiolus mussels from a cold seep and three hydrothermal vents in the West Pacific

Deep-sea Bathymodiolus mussels are ubiquitous in most cold seeps and hydrothermal fields, where they have adapted to various toxic environments including high metal exposure. However, there is scarce knowledge of metal accumulation and metal-related biomarkers in B. mussels. Here, we present data for metal concentrations (Ag, Cd, Cr, Cu, Fe, Mn, Pb, and Zn) and metal related biomarkers (superoxide dismutase-SOD, catalase-CAT, glutathione peroxidase-GPX, glutathione-GSH, metallothioneins-MTs, and lipid peroxidation-LPO) in different tissues of B. mussels from four different deep-sea geochemical settings: one cold seep and three vent fields in the West Pacific Ocean. Results showed that mussel gills generally exhibited higher metal enrichment than the mantle. Mussels from hydrothermal vents usually had higher metal concentrations (Fe, Cr, Cd, and Pb) than those from cold seep, which could be related to their higher contents in fluids or sediments. However, despite quite different metals loads among the geochemical environment settings, Mn, Zn, and Cu concentrations varied over a smaller range across the sampling sites, implying biological regulation by deep-sea mussels for these elements. Several statistically significant correlations were observed between SOD, CAT, GSH, MTs, and metal levels in analyzed tissues. Although the vent ecosystem is harsher than the cold seep ecosystem, according to our results their mussels' biomarker levels were not so different. This finding suggests that some adaptive or compensatory mechanisms may occur in chronically polluted deep-sea mussels. Principal component analysis allowed for distinguishing different deep-sea settings, indicating that B. mussels are robust indicators of their living environments. We also compared our results with those reported for coastal mussels. To our best knowledge, this is the first integrated study to report metal accumulation and metal-related biomarkers in the deep-sea B. mussels from the West Pacific.

Insights into the strategy of micro-environmental adaptation: Transcriptomic analysis of two alvinocaridid shrimps at a hydrothermal vent

Diffusing fluid at a deep-sea hydrothermal vent creates rapid, acute physico-chemical gradients that correlate strongly with the distribution of the vent fauna. Two alvinocaridid shrimps, Alvinocaris longirostris and Shinkaicaris leurokolos occupy distinct microhabitats around these vents and exhibit different thermal preferences. S. leurokolos inhabits the central area closer to the active chimney, while A. longirostris inhabits the peripheral area. In this study, we screened candidate genes that might be involved in niche separation and microhabitat adaptation through comparative transcriptomics. The results showed that among the top 20% of overexpressed genes, gene families related to protein synthesis and structural components were much more abundant in S. leurokolos compared to A. longirostris. Moreover, 15 out of 25 genes involved in cellular carbohydrate metabolism were related to trehalose biosynthesis, versus 1 out of 5 in A. longirostris. Trehalose, a non-reducing disaccharide, is a multifunctional molecule and has been proven to act as a protectant responsible for thermotolerance in Saccharomyces cerevisiae. Putative positively selected genes involved in chitin metabolism and the immune system (lectin, serine protease and antimicrobial peptide) were enriched in S. leurokolos. In particular, one collagen and two serine proteases were found to have experienced strong positive selection. In addition, sulfotransferase-related genes were both overexpressed and positively selected in S. leurokolos. Finally, genes related to structural proteins, immune proteins and protectants were overexpressed or positively selected. These characteristics could represent adaptations of S. leurokolos to its microhabitat, which need to be confirmed by more evidence, such as data from large samples and different development stages of these alvinocaridid shrimps.

Identification and characterization of endosymbiosis-related immune genes in deep-sea mussels Gigantidas platifrons

Deep-sea mussels of the subfamily Bathymodiolinae are common and numerically dominant species widely distributed in cold seeps and hydrothermal vents. During long-time evolution, deep-sea mussels have evolved to be well adapted to the local environment of cold seeps and hydrothermal vents by various ways, especially by establishing endosymbiosis with chemotrophic bacteria. However, biological processes underlying the establishment and maintenance of symbiosis between host mussels and symbionts are largely unclear. In the present study, Gigantidas platifrons genes possibly involved in the symbiosis with methane oxidation symbionts were identified and characterized by Lipopolysaccharide (LPS) pull-down and in situ hybridization. Five immune related proteins including Toll-like receptor 2 (TLR2), integrin, vacuolar sorting protein (VSP), matrix metalloproteinase 1 (MMP1), and leucine-rich repeat (LRR-1) were identified by LPS pull-down assay. These five proteins were all conserved in either molecular sequences or functional domains and known to be key molecules in host immune recognition, phagocytosis, and lysosome-mediated digestion. Furthermore, in situ hybridization of LRR-1, TLR2 and VSP genes was conducted to investigate their expression patterns in gill tissues of G. platifrons. Consequently, LRR-1, TLR2, and VSP genes were found expressed exclusively in the bacteriocytes of G. platifrons. Therefore, it was suggested that TLR2, integrin, VSP, MMP1, and LRR-1 might be crucial molecules in the symbiosis between G. platifrons and methane oxidation bacteria by participating in symbiosis-related immune processes.

Identification of cells expressing two peptidoglycan recognition proteins in the gill of the vent mussel, Bathymodiolus septemdierum

In symbiotic systems in which symbionts are transmitted horizontally, hosts must accept symbionts from the environment while defending themselves against invading pathogenic microorganisms. How they distinguish pathogens from symbionts and how the latter evade host immune defences are not clearly understood. Recognition of foreign materials is one of the most critical steps in stimulating immune responses, and pattern recognition receptors (PRRs) play vital roles in this process. In this study, we focused on a group of highly conserved PRRs, peptidoglycan recognition proteins (PGRPs), in the deep-sea mussel, Bathymodiolus septemdierum, which harbours chemosynthetic bacteria in their gill epithelial cells. We isolated B. septemdierum PGRP genes BsPGRP-S and BsPGRP-L, which encode a short- and a long-type PGRP, respectively. The short-type PGRP has a signal peptide and was expressed in the asymbiotic goblet mucous cells in the gill epithelium, whereas the long-type PGRP was predicted to include a transmembrane domain and was expressed in gill bacteriocytes. Based on these findings, we hypothesize that the secreted and transmembrane PGRPs are engaged in host defence against pathogenic bacteria and/or in the regulation of symbiosis via different cellular localizations and mechanisms.

Stress Adapted Mollusca and Nematoda Exhibit Convergently Expanded Hsp70 and AIG1 Gene Families

We recently sequenced the genome of the first subterrestrial metazoan, the nematode Halicephalobus mephisto. A central finding was a dramatic expansion of genes encoding avrRpt2 induced gene (AIG1), and 70 kDa heat shock (Hsp70) domains. While the role of Hsp70 in thermotolerance is well established, the contribution of AIG1 is much more poorly characterized, though in plants some members of this family are heat-induced. Hypothesizing that this dual domain expansion may constitute a general biosignature of thermal stress adaptation, here we examine a number of genomes, finding that expansion of both AIG1 and Hsp70 is common in bivalves. Phylogenetic analysis reveals that the bivalve-specific Hsp70 protein expansion groups with H. mephisto sequences. Our identification of the same gene expansions in bivalves and a nematode implies that this biosignature may be a general stress adaptation strategy for protostomes, particularly those organisms that cannot escape their stressful environments. We hypothesize that the two families play largely complementary mechanistic roles, with Hsp70 directly refolding heat-denatured proteins while AIG1 promotes cellular and organismal survival by inhibiting apoptosis.

Transcriptomic analysis reveals insights into deep-sea adaptations of the dominant species, Shinkaia crosnieri (Crustacea: Decapoda: Anomura), inhabiting both hydrothermal vents and cold seeps

Background

Hydrothermal vents and cold seeps are typical deep-sea chemosynthetically-driven ecosystems that allow high abundance of specialized macro-benthos. To gather knowledge about the genetic basis of adaptation to these extreme environments, species shared between different habitats, especially for the dominant species, are of particular interest. The galatheid squat lobster, Shinkaia crosnieri Baba and Williams, 1998, is one of the few dominant species inhabiting both deep-sea hydrothermal vents and cold seeps. In this study, we performed transcriptome analyses of S. crosnieri collected from the Iheya North hydrothermal vent (HV) and a cold seep in the South China Sea (CS) to provide insights into how this species has evolved to thrive in different deep-sea chemosynthetic ecosystems.

Results

We analyzed 5347 orthologs between HV and CS to identify genes under positive selection through the maximum likelihood approach. A total of 82 genes were identified to be positively selected and covered diverse functional categories, potentially indicating their importance for S. crosnieri to cope with environmental heterogeneity between deep-sea vents and seeps. Among 39,806 annotated unigenes, a large number of differentially expressed genes (DEGs) were identified between HV and CS, including 339 and 206 genes significantly up-regulated in HV and CS, respectively. Most of the DEGs associated with stress response and immunity were up-regulated in HV, possibly allowing S. crosnieri to increase its capability to manage more environmental stresses in the hydrothermal vents.

Conclusions

We provide the first comprehensive transcriptomic resource for the deep-sea squat lobster, S. crosnieri, inhabiting both hydrothermal vents and cold seeps. A number of stress response and immune-related genes were positively selected and/or differentially expressed, potentially indicating their important roles for S. crosnieri to thrive in both deep-sea vents and cold seeps. Our results indicated that genetic adaptation of S. crosnieri to different deep-sea chemosynthetic environments might be mediated by adaptive evolution of functional genes related to stress response and immunity, and alterations in their gene expression that lead to different stress resistance. However, further work is required to test these proposed hypotheses. All results can constitute important baseline data for further studies towards elucidating the adaptive mechanisms in deep-sea crustaceans.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5753-7) contains supplementary material, which is available to authorized users.

High rates of apoptosis visualized in the symbiont-bearing gills of deep-sea Bathymodiolus mussels

Symbiosis between Bathymodiolus and Gammaproteobacteria allows these deep-sea mussels to live in toxic environments such as hydrothermal vents and cold seeps. The quantity of endosymbionts within the gill-bacteriocytes appears to vary according to the hosts environment; however, the mechanisms of endosymbiont population size regulation remain obscure. We investigated the possibility of a control of endosymbiont density by apoptosis, a programmed cell death, in three mussel species. Fluorometric TUNEL and active Caspase-3-targeting antibodies were used to visualize and quantify apoptotic cells in mussel gills. To control for potential artefacts due to depressurization upon specimen recovery from the deep-sea, the apoptotic rates between mussels recovered unpressurised, versus mussels recovered in a pressure-maintaining device, were compared in two species from hydrothermal vents on the Mid-Atlantic Ridge: Bathymodiolus azoricus and B. puteoserpentis. Results show that pressurized recovery had no significant effect on the apoptotic rate in the gill filaments. Apoptotic levels were highest in the ciliated zone and in the circulating hemocytes, compared to the bacteriocyte zone. Apoptotic gill-cells in B. aff. boomerang from cold seeps off the Gulf of Guinea show similar distribution patterns. Deep-sea symbiotic mussels have much higher rates of apoptosis in their gills than the coastal mussel Mytilus edulis, which lacks chemolithoautotrophic symbionts. We discuss how apoptosis might be one of the mechanisms that contribute to the adaptation of deep-sea mussels to toxic environments and/or to symbiosis.

Two Novel Short Peptidoglycan Recognition Proteins (PGRPs) From the Deep Sea Vesicomyidae Clam Archivesica packardana: Identification, Recombinant Expression and Bioactivity

Vesicomyidae clams are common species living in cold seeps, which incorporates symbiotic bacteria into their body maintaining endosymbiosis relationship. As members of pattern recognition receptor (PRR) family, peptidoglycan recognition proteins (PGRPs) recognize pathogen associated molecular patterns and play an important role in innate immunity. In present study, two short PGRPs (ApPGRP-1 and -2) were first identified from Vesicomyidae clam Archivesica packardana. Sequences analysis showed that they have both conserved Zn²⁺ binding sites (H-H-C) and amidase catalytic sites (H-Y-H-T-C), and phylogenetic tree indicated that they clustered with short PGRPs of other molluscs. PGN assay showed that ApPGRPs could bind Lys-type PGN from Staphylococcus aureus and Dap-type PGN from Bacillus subtilis, and revealed amidase activity with selective zinc ion dependence. rApPGRP-1 and -2 (recombinant ApPGRP-1 and -2) could bind six bacteria with a broad spectrum and had both zinc-dependent and -independent bactericidal activity. ApPGRPs had the complete functions of effectors and partial functions of receptors from PGRPs. Further analyses showed that ApPGRPs from A. packardana might be involved in the endosymbiosis relationship between the host clam and endosymbiotic bacteria as a regulator. The results of these experiments suggested that ApPGRPs were involved in cold seep clams’ immune response. This study provides basic information for further research on the immune mechanisms of deep sea organisms.

Trace metals in oysters: molecular and cellular mechanisms and ecotoxicological impacts

Oysters are important benthic bivalves in coastal and estuarine environments. They are widely farmed due to their rapid growth and taste; they are also widely applied in environmental monitoring of coastal pollution due to their accumulation of contaminants. Most importantly, oysters are among the few marine organisms that are considered to be hyper-accumulators of many toxic metals, such as cadmium, copper and zinc. As such, there is a tremendous call to study the interactions between metals and oysters, especially due to the increasing metal pollution in many coastal and estuarine waters. Over the past decades, many studies have focused on metal accumulation in oysters as well as the ecotoxicological effects of metals on oysters. In this review, we summarize the recent progress in our understanding of the molecular and cellular mechanisms of metal accumulation, sequestration and toxicity in oysters. Applications of modern technologies such as omics and nanoscale imaging have added significantly to our knowledge of metal biology in oysters. Variations between different metals also demonstrate the diversity of the interactions between oysters and metals. Despite this recent progress, however, there is a need for further study of the molecular mechanisms of metal uptake and toxicity as well as the joint effects of metal mixtures on oyster populations. Oysters have higher numbers of stress responsive genes than most animals, which may have been induced by gene duplication during the evolution of their intertidal environmental adaptations. The divergent expression of stress responsive genes may explain the different tolerances for metals among different species. These fundamental studies may eventually provide promising solutions for reducing toxic metal concentrations in oysters for safe consumption by humans. To conclude, the complexity of life history and metal chemistry of oysters coupled with emerging pollution and application of modern techniques represents an important and exciting research area in modern ecotoxicology.

pH Stress-Induced Cooperation between Rhodococcus ruber YYL and Bacillus cereus MLY1 in Biodegradation of Tetrahydrofuran

Microbial consortia consisting of cooperational strains exhibit biodegradation performance superior to that of single microbial strains and improved remediation efficiency by relieving the environmental stress. Tetrahydrofuran (THF), a universal solvent widely used in chemical and pharmaceutical synthesis, significantly affects the environment. As a refractory pollutant, THF can be degraded by some microbial strains under suitable conditions. There are often a variety of stresses, especially pH stress, that inhibit the THF-degradation efficiency of microbial consortia. Therefore, it is necessary to study the molecular mechanisms of microbial cooperational degradation of THF. In this study, under conditions of low pH (initial pH = 7.0) stress, a synergistic promotion of the THF degradation capability of the strain Rhodococcus ruber YYL was found in the presence of a non-THF degrading strain Bacillus cereus MLY1. Metatranscriptome analysis revealed that the low pH stress induced the strain YYL to up-regulate the genes involved in anti-oxidation, mutation, steroid and bile acid metabolism, and translation, while simultaneously down-regulating the genes involved in ATP production. In the co-culture system, strain MLY1 provides fatty acids, ATP, and amino acids for strain YYL in response to low pH stress during THF degradation. In return, YYL shares the metabolic intermediates of THF with MLY1 as carbon sources. This study provides the preliminary mechanism to understand how microbial consortia improve the degradation efficiency of refractory furan pollutants under environmental stress conditions.

Comparative Transcriptomic Analysis Reveals Candidate Genes and Pathways Involved in Larval Settlement of the Barnacle Megabalanus volcano

Megabalanus barnacle is one of the model organisms for marine biofouling research. However, further elucidation of molecular mechanisms underlying larval settlement has been hindered due to the lack of genomic information thus far. In the present study, cDNA libraries were constructed for cyprids, the key stage for larval settlement, and adults of Megabalanus volcano. After high-throughput sequencing and de novo assembly, 42,620 unigenes were obtained with a N50 value of 1532 bp. These unigenes were annotated by blasting against the NCBI non-redundant (nr), Swiss-Prot, Cluster of Orthologous Groups (COG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Finally, 19,522, 15,691, 14,459, and 10,914 unigenes were identified correspondingly. There were 22,158 differentially expressed genes (DEGs) identified between two stages. Compared with the cyprid stage, 8241 unigenes were down-regulated and 13,917 unigenes were up-regulated at the adult stage. The neuroactive ligand-receptor interaction pathway (ko04080) was significantly enriched by KEGG enrichment analysis of the DEGs, suggesting that it possibly involved in larval settlement. Potential functions of three conserved allatostatin neuropeptide-receptor pairs and two light-sensitive opsin proteins were further characterized, indicating that they might regulate attachment and metamorphosis at cyprid stage. These results provided a deeper insight into the molecular mechanisms underlying larval settlement of barnacles.

Bias and Correction in RNA-seq Data for Marine Species

RNA-seq is a recently developed approach widely used for transcriptome profiling in biological analyses that use next-generation sequencing technologies. Accurate estimation of gene expression levels is critical for answering biological questions. Here, we show that the commonly used measure of gene expression levels, fragments per kilobase of transcript per million mapped reads (FPKM), is biased in transcript length, GC content, and dinucleotide frequencies in the RNA-seq analysis of marine species. We used a generalized linear model to correct the observed biases of FPKM. We used RNA-seq data sets from eight species obtained by different sequencing methods to evaluate the correction methods. Our work contributes to the understanding of potential technical artifacts in RNA-seq experiments for marine species, and presents a means by which more accurate gene expression measures can be obtained.

Comparative immunogenomics of molluscs

Comparative immunology, studying both vertebrates and invertebrates, provided the earliest descriptions of phagocytosis as a general immune mechanism. However, the large scale of animal diversity challenges all-inclusive investigations and the field of immunology has developed by mostly emphasizing study of a few vertebrate species. In addressing the lack of comprehensive understanding of animal immunity, especially that of invertebrates, comparative immunology helps toward management of invertebrates that are food sources, agricultural pests, pathogens, or transmit diseases, and helps interpret the evolution of animal immunity. Initial studies showed that the Mollusca (second largest animal phylum), and invertebrates in general, possess innate defenses but lack the lymphocytic immune system that characterizes vertebrate immunology. Recognizing the reality of both common and taxon-specific immune features, and applying up-to-date cell and molecular research capabilities, in-depth studies of a select number of bivalve and gastropod species continue to reveal novel aspects of molluscan immunity. The genomics era heralded a new stage of comparative immunology; large-scale efforts yielded an initial set of full molluscan genome sequences that is available for analyses of full complements of immune genes and regulatory sequences. Next-generation sequencing (NGS), due to lower cost and effort required, allows individual researchers to generate large sequence datasets for growing numbers of molluscs. RNAseq provides expression profiles that enable discovery of immune genes and genome sequences reveal distribution and diversity of immune factors across molluscan phylogeny. Although computational de novo sequence assembly will benefit from continued development and automated annotation may require some experimental validation, NGS is a powerful tool for comparative immunology, especially increasing coverage of the extensive molluscan diversity. To date, immunogenomics revealed new levels of complexity of molluscan defense by indicating sequence heterogeneity in individual snails and bivalves, and members of expanded immune gene families are expressed differentially to generate pathogen-specific defense responses.

Complexities of gene expression patterns in natural populations of an extremophile fish (Poecilia mexicana, Poeciliidae)

Variation in gene expression can provide insights into organismal responses to environmental stress and physiological mechanisms mediating adaptation to habitats with contrasting environmental conditions. We performed an RNA-sequencing experiment to quantify gene expression patterns in fish adapted to habitats with different combinations of environmental stressors, including the presence of toxic hydrogen sulphide (H2 S) and the absence of light in caves. We specifically asked how gene expression varies among populations living in different habitats, whether population differences were consistent among organs, and whether there is evidence for shared expression responses in populations exposed to the same stressors. We analysed organ-specific transcriptome-wide data from four ecotypes of Poecilia mexicana (nonsulphidic surface, sulphidic surface, nonsulphidic cave and sulphidic cave). The majority of variation in gene expression was correlated with organ type, and the presence of specific environmental stressors elicited unique expression differences among organs. Shared patterns of gene expression between populations exposed to the same environmental stressors increased with levels of organismal organization (from transcript to gene to physiological pathway). In addition, shared patterns of gene expression were more common between populations from sulphidic than populations from cave habitats, potentially indicating that physiochemical stressors with clear biochemical consequences can constrain the diversity of adaptive solutions that mitigate their adverse effects. Overall, our analyses provided insights into transcriptional variation in a unique system, in which adaptation to H2 S and darkness coincide. Functional annotations of differentially expressed genes provide a springboard for investigating physiological mechanisms putatively underlying adaptation to extreme environments.

Characterization of the IRF2 proteins isolated from the deep-sea mussel Bathymodiolus platifrons and the shallow-water mussel Modiolus modiolus

Interferon regulatory factors (IRFs) are transcription factors that play important roles in immune defense, stress response, hematopoietic differentiation, and cell apoptosis. IRFs of invertebrate organisms and their functions remain largely unexplored. In the present study, for the first time new IRFs (BpIRF2 and MmIRF2) were identified in the deep-sea mussel Bathymodiolus platifrons and the shallow-water mussel Modiolus modiolus. The open reading frame of BpIRF2 and MmIRF2 encoded putative proteins of 354 and 348 amino acids, respectively. Comparison and phylogenetic analysis revealed that both IRF2 proteins were new identified invertebrate IRF molecular. As transcriptional factors, both BpIRF2 and MmIRF2 could activate the interferon-stimulated response element-containing promoter and BpIRF2 could interact with itself. Moreover, both BpIRF2 and MmIRF2 were localized to the cytoplasm and nucleus. Collectively, these results demonstrated that IRF2 proteins might be crucial in the innate immunity of deep-sea and shallow-water mussels.

Exploring the molecular basis of adaptive evolution in hydrothermal vent crab Austinograea alayseae by transcriptome analysis

Elucidating the genetic mechanisms of adaptation to the hydrothermal vent in organisms at genomic level is significant for understanding the adaptive evolution process in the extreme environment. We performed RNA-seq on four different tissues of a vent crab species, Austinograea alayseae, producing 725,461 unigenes and 134,489 annotated genes. Genes related to sensory, circadian rhythm, hormone, hypoxia stress, metal detoxification and immunity were identified. It was noted that in the degenerated eyestalk, transcription of phototransduction related genes which are important for retinal function was greatly reduced; three crucial neuropeptide hormones, one molt-inhibiting and two crustacean hyperglycemic hormone precursors were characterized with conserved domains; hypoxia-inducible factor 1 and two novel isoforms of metallothioneins in the vent crabs were discovered. An analysis of 6,932 orthologs among three crabs A. alayseae, Portunus trituberculutus and Eriocheir sinensis revealed 19 positive selected genes (PSGs). Most of the PSGs were involved in immune responses, such as crustins and anti-lipopolysaccharide factor, suggesting their function in the adaptation to environment. The characterization of the first vent crab transcriptome provides abundant resources for genetic and evolutionary studies of this species, and paves the way for further investigation of vent adaptation process in crabs.

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.

Comparative transcriptome analysis of Rimicaris sp. reveals novel molecular features associated with survival in deep-sea hydrothermal vent

Shrimp of the family Alvinocarididae are the predominant megafauna of deep-sea hydrothermal vents. However, genome information on this family is currently unavailable. In the present study, by employing Illumina sequencing, we performed the first de novo transcriptome analysis of the gills of the shrimp Rimicaris sp. from the hydrothermal vent in Desmos, Manus Basin. The analysis was conducted in a comparative manner with the shrimp taken directly from the vent (GR samples) and the shrimp that had been maintained for ten days under normal laboratory condition (mGR samples). Among the 128,938 unigenes identified, a large number of differentially expressed genes (DEGs) between the GR and mGR samples were detected, including 2365 and 1607 genes significantly upregulated and downregulated, respectively, in GR. The DEGs covered diverse functional categories. Most of the DEGs associated with immunity were downregulated in GR, while most of the DEGs associated with sulfur metabolism and detoxification were upregulated in GR. These results provide the first comprehensive transcriptomic resource for hydrothermal vent Rimicaris and revealed varied categories of genes likely involved in deep-sea survival.

Adaptation and evolution of deep-sea scale worms (Annelida: Polynoidae): insights from transcriptome comparison with a shallow-water species

Polynoid scale worms (Polynoidae, Annelida) invaded deep-sea chemosynthesis-based ecosystems approximately 60 million years ago, but little is known about their genetic adaptation to the extreme deep-sea environment. In this study, we reported the first two transcriptomes of deep-sea polynoids (Branchipolynoe pettiboneae, Lepidonotopodium sp.) and compared them with the transcriptome of a shallow-water polynoid (Harmothoe imbricata). We determined codon and amino acid usage, positive selected genes, highly expressed genes and putative duplicated genes. Transcriptome assembly produced 98,806 to 225,709 contigs in the three species. There were more positively charged amino acids (i.e., histidine and arginine) and less negatively charged amino acids (i.e., aspartic acid and glutamic acid) in the deep-sea species. There were 120 genes showing clear evidence of positive selection. Among the 10% most highly expressed genes, there were more hemoglobin genes with high expression levels in both deep-sea species. The duplicated genes related to DNA recombination and metabolism, and gene expression were only enriched in deep-sea species. Deep-sea scale worms adopted two strategies of adaptation to hypoxia in the chemosynthesis-based habitats (i.e., rapid evolution of tetra-domain hemoglobin in Branchipolynoe or high expression of single-domain hemoglobin in Lepidonotopodium sp.).

Adaptation to deep-sea chemosynthetic environments as revealed by mussel genomes

Hydrothermal vents and methane seeps are extreme deep-sea ecosystems that support dense populations of specialized macro-benthos such as mussels. But the lack of genome information hinders the understanding of the adaptation of these animals to such inhospitable environments. Here we report the genomes of a deep-sea vent/seep mussel (Bathymodiolus platifrons) and a shallow-water mussel (Modiolus philippinarum). Phylogenetic analysis shows that these mussel species diverged approximately 110.4 million years ago. Many gene families, especially those for stabilizing protein structures and removing toxic substances from cells, are highly expanded in B. platifrons, indicating adaptation to extreme environmental conditions. The innate immune system of B. platifrons is considerably more complex than that of other lophotrochozoan species, including M. philippinarum, with substantial expansion and high expression levels of gene families that are related to immune recognition, endocytosis and caspase-mediated apoptosis in the gill, revealing presumed genetic adaptation of the deep-sea mussel to the presence of its chemoautotrophic endosymbionts. A follow-up metaproteomic analysis of the gill of B. platifrons shows methanotrophy, assimilatory sulfate reduction and ammonia metabolic pathways in the symbionts, providing energy and nutrients, which allow the host to thrive. Our study of the genomic composition allowing symbiosis in extremophile molluscs gives wider insights into the mechanisms of symbiosis in other organisms such as deep-sea tubeworms and giant clams.

Linking genomic responses of gonads with reproductive impairment in marine medaka (Oryzias melastigma) exposed chronically to the chemopreventive and antifouling agent, 3,3'-diindolylmethane (DIM)

3,3'-Diindolylmethane (DIM) has been promoted as an effective chemopreventive and antifouling additive. However, the concurrent risks or side effects of DIM are not fully understood, especially on tissues responsive to estrogen. Therefore, this study employed marine medaka (Oryzias melastigma) as a test model to evaluate relative safety and explore mechanisms of toxic action of DIM on development and function of gonad after chronic (28days) aqueous exposure to relatively low doses (0μg/L or 8.5μg/L). Integration of comprehensive toxicogenomic analysis at the transcriptome and proteome levels with apical endpoints, such as production of eggs and swimming performance of larvae, elucidated the molecular linkage in gonad from bottom up along the reproductive adverse outcome pathway. A series of sequential changes at the transcript and protein levels were linked to lesser fecundity and viability of larvae exposed to DIM. Anomalous production of vitellogenin (VTG) and eggshell proteins in testis confirmed the estrogenic potency of DIM. In the ovary, although storage of VTG was greater, lesser expressions of cathepsin enzymes blocked cleavage and incorporation of VTG into oocytes as yolk, which acted together with lower eggshell proteins to inhibit maturation of primary oocyte and thus contributed to impairment of fecundity. Overall, this study demonstrated that exposure to DIM impaired reproductive fitness. Diverse molecular initiating changes in gonads were linked to apical endpoints that could be used in assessment of risks posed by DIM on gametogenesis. In combination with chemical stability and potent endocrine disruption, the results of this study can inform decisions about the use of DIM either as chemopreventive or antifouling agent.