Transcriptome analysis of the Yesso scallop, Patinopecten yessoensis gills in response to water temperature fluctuations

Warming and polymetallic stress induce proteomic and physiological shifts in the neurotoxic Alexandrium pacificum as possible response to global changes

Harmful Algal Blooms involving the dinoflagellate Alexandrium pacificum continue to increase in ecosystems suffering the climate warming and anthropogenic pressure. Changes in the total proteome and physiological traits of the Mediterranean A. pacificum SG C10-3 strain were measured in response to increasing temperature (24 °C, 27 °C, 30 °C) and trace metal contamination (Cu2+, Pb2+, Zn2+, Cd2+). Warming reduced the cell densities and maximal growth rate (μmax), but the strain persisted at 30 °C with more large cells. The polymetallic stress increased cell sizes, reduced cell growth at 24 °C-27 °C and it increased this at 30 °C. Toxin profiles showed a predominance of GTX4 (32-38 %), then C2 (11-34 %) or GTX6 (18-24 %) among the total Paralytic Shellfish Toxins, however these were modified under warming, showing increased contents in GTX1 (among the most toxic), GTX5, C1 and NeoSTX, while dc-NeoSTX and STX (among the most toxic) only appeared at 30 °C. Under polymetallic contamination, warming also increased contents in GTX5 and NeoSTX. In contrast, polymetallic stress, or warming had harmful effects on C2 contents. Proteins were more quantitatively produced by A. pacificum SG C10-3 under warming in accordance with the high levels of up-regulated proteins found in the total proteome in this condition. Polymetallic stress, only or combined with warming, led to low proteomic modifications (1 % or 4 %), whereas warming induced strong 52 % modified proteomic response, mainly based on up-regulated proteins involved in photosynthesis (light harvesting complex protein), carbohydrate metabolism (arylsulfatase) and translation (ribosomal proteins), and with the lesser down-regulated proteins principally associated with the lipid metabolism (type I polyketide synthase). Our results show that warming triggers a strong up-regulated A. pacificum SG C10-3 proteomic response, which, coupled to modified cell sizes and toxin profiles, could help it to withstand stress conditions. This could presage the success of A. pacificum in anthropized ecosystems submitted to global warming in which this dinoflagellate also might be more toxic.

Genome-Wide Identification, Characterization, and Transcriptional Profile of the HECT E3 Ubiquitin Ligase Gene Family in the Hard-Shelled Mussel Mytilus coruscus Gould

The homologous E6-AP carboxy-terminal structural domain (HECT) contained in E3 ubiquitin ligases (E3s) is a key factor in protein degradation and maintenance of cellular homeostasis in animals. However, the functional roles and evolutionary aspects of the HECT gene family in bivalve mussels remain unclear and warrant further investigation. In this study, we identified 22 HECT genes within the genome of Mytilus coruscus Gould, all containing a conserved HECT structural domain derived from dispersed repeats, distributed unevenly across 11 chromosomes. Phylogenetic analysis classified M. coruscus HECT genes into six major classes, with amino acid sequences within the same evolutionary clade displaying similar conserved motifs. Homology analysis with HECT genes of four bivalve species revealed that M. coruscus and Mytilus galloprovincialis possessed the largest number of homologous gene pairs, showing a significant correlation between the two in the evolution of the HECT gene family. Homology analysis with HECT genes of four bivalve species revealed that M. coruscus and M. galloprovincialis possessed the largest number of homologous gene pairs, showing a significant correlation between the two in the evolution of the HECT gene family. M. coruscus exhibited pronounced and specific expression in gills and blood tissues. Notably, Mco_UPL3 gene expression was significantly upregulated after 12 h of acute heat stress (33 °C) and 24 h of Vibrio injection (0.4 OD). Gene ontology analysis of the HECT genes in M. coruscus revealed that it is primarily enriched in protein modification and degradation functions. This suggests that HECT genes may play a key role in protein degradation and immunomodulation in M. coruscus. These findings offer valuable insights for the breeding of stress-tolerant traits in M. coruscus. In summary, our data shed light on the potential functions of HECT E3 ligases in response to heat stress and Vibrio infection, providing practical guidance for enhancing resilience through breeding in M. coruscus.

Gills de novo assembly reveals oxidative stress, unfolded protein, and immune response on red cusk-eel (Genypterus chilensis) under thermal stress

The heat waves on the South Pacific coast could lead to thermal stress in native fish. The red cusk-eel (Genypterus chilensis) is relevant for Chilean artisanal fisheries and aquaculture diversification. This study examined the effect of high-temperature stress in the gills of G. chilensis in control (14 °C) and high-temperature stress (19 °C) conditions. High-temperature stress induces a significant increase in gills cortisol levels. Additionally, oxidative damage was observed in gills (protein carbonylation and lipoperoxidation). RNA-seq data was used to build the first transcriptome assembly of gills in this species (23,656 annotated transcripts). A total of 1138 down-regulated and 1531 up-regulated transcripts were observed in response to high-temperature stress in gills. The enrichment analysis showed immune response and replication enriched processes (on down-regulated transcripts), and processes related to the folding of proteins, endoplasmic reticulum, and transporter activity (on up-regulated transcripts). The present study showed how gills could be affected by high-temperature stress.

Comparative transcriptomic analysis revealed changes in multiple signaling pathways involved in protein degradation in the digestive gland of Mytilus coruscus during high-temperatures

As a result of global warming, the Mytilus coruscus living attached in the intertidal zone experience extreme and fluctuating changes in temperature, and extreme temperature changes are causing mass mortality of intertidal species. This study explores the transcriptional response of M. coruscus at different temperatures (18 °C, 26 °C, and 33 °C) and different times (0, 12, and 24 h) of action by analyzing the potential temperature of the intertidal zone. In response to high temperatures, several signaling pathways in M. coruscus, ribosome, endocytosis, endoplasmic reticulum stress, protein degradation, and lysosomes, interact to counter the adverse effects of high temperatures on protein homeostasis. Increased expression of key genes, including heat shock proteins (Hsp70, Hsp20, and Hsp110), Lysosome-associated membrane glycoprotein (LAMP), endoplasmic reticulum chaperone (BiP), and baculoviral IAP repeat-containing protein 7 (BIRC7), may further mitigate the effects of heat stress and delay mortality in M. coruscus. These results reveal changes in multiple signaling pathways involved in protein degradation during high-temperature stress, which will contribute to our overall understanding of the molecular mechanisms underlying the response of M. coruscus to high-temperature stress.

Tissue-Specific and Time-Dependent Expressions of PC4s in Bay Scallop (Argopecten irradians irradians) Reveal Function Allocation in Thermal Response

Transcriptional coactivator p15 (PC4) encodes a structurally conserved but functionally diverse protein that plays crucial roles in RNAP-II-mediated transcription, DNA replication and damage repair. Although structures and functions of PC4 have been reported in most vertebrates and some invertebrates, the PC4 genes were less systematically identified and characterized in the bay scallop Argopecten irradians irradians. In this study, five PC4 genes (AiPC4s) were successfully identified in bay scallops via whole-genome scanning through in silico analysis. Protein structure and phylogenetic analyses of AiPC4s were conducted to determine the identities and evolutionary relationships of these genes. Expression levels of AiPC4s were assessed in embryos/larvae at all developmental stages, in healthy adult tissues and in different tissues (mantles, gills, hemocytes and hearts) being processed under 32 °C stress with different time durations (0 h, 6 h, 12 h, 24 h, 3 d, 6 d and 10 d). Spatiotemporal expression profiles of AiPC4s suggested the functional roles of the genes in embryos/larvae at all developmental stages and in healthy adult tissues in bay scallop. Expression regulations (up- and down-) of AiPC4s under high-temperature stress displayed both tissue-specific and time-dependent patterns with function allocations, revealing that AiPC4s performed differentiated functions in response to thermal stress. This work provides clues of molecular function allocation of PC4 in scallops in response to thermal stress and helps in illustrating how marine bivalves resist elevated seawater temperature.

Microplastics influence physiological processes, growth and reproduction in the Manila clam, Ruditapes philippinarum

Microplastics (<5 mm) are widely distributed in marine environments and pose a serious threat to bivalves. Here, the ingestion and accumulation of polystyrene microplastics (PS microplastics, diameters 5 and 10 μm) by the Manila clam, Ruditapes philippinarum, and their impacts on physiological processes, growth and reproduction were studied. The results showed that both PS microplastics were ingested by the Manila clam and accumulated in their gills, hepatopancreases and intestines. Furthermore, the accumulation of 5 and 10 μm PS microplastics significantly increased the rates of respiration and excretion while significantly decreasing feeding and absorption efficiency (AE), leading to a dramatically reduced amount of energy available for growth (SfG) and ultimately led to slower growth. The dynamic energy budget (DEB) model predicts that PS microplastic exposure for 200 days would cause lower shell/flesh growth rates and reproductive potentiality. Transcriptomic profiles support these results, as carbon and protein metabolism and oxytocin and insulin-related signaling pathways were significantly altered in clams in response to PS microplastics. This study provides evidence that microplastics strongly affect the physiological activities, energy allocation, growth and reproduction of filter-feeding bivalves.

Integrated transcriptomics and metabolomics analyses reveal benzo[a]pyrene enhances the toxicity of mercury to the Manila clam, Ruditapes philippinarum

Mercury (Hg²⁺) and benzo[a]pyrene (BaP) are ubiquitous and persistent pollutants with multiple toxicities in bivalve molluscs. Here, the toxicological responses in the gills of Manila clams, Ruditapes philippinarum, to Hg²⁺ (10 μg L^-1), BaP (3 μg L^-1), and their mixture were analysed using transcriptomics and metabolomics approaches. Comparisons of the transcriptomes and metabolomes of Hg²⁺-and/or BaP-treated clams with control animals revealed the involvement of the detoxification metabolism, immune defence, energy-related pathways, and osmotic regulation in the stress response of R. philippinarum. Exposure to Hg²⁺ alone primarily enhanced the detoxification and energy metabolic pathways by significantly increasing the expression of genes associated with heat-shock proteins and oxidative phosphorylation. However, co-exposure to Hg²⁺ and BaP caused greater immunotoxicity and disrupted detoxification metabolism, the TCA cycle, glycolysis, and ATP generation. The expression levels of cytochrome P450 1A1 (CYP1A1), multidrug resistance-associated protein 1 (MRP1), and myosin (MYO), and the activity of electron transport system (ETS) in gills were detected, supporting the underlying toxic mechanisms of Hg²⁺ and BaP. We suggest that the presence of BaP enhances the toxicity of Hg²⁺ by 1) hampering the detoxification of Hg²⁺, 2) increasing the immunotoxicity of Hg²⁺, and 3) constraining energy availability for clams.

Transcriptomic analyses provide insights into the adaptive responses to heat stress in the ark shells, Scapharca subcrenata

The ark shell, Scapharca subcrenata, is susceptible to high temperature which may lead to mass mortality in hot summers. Herein, we conducted the transcriptomic analyses of haemocytes in ark shells under thermal stress, to reveal the underlying molecular mechanisms of heat resistance in these animals. The results showed that a total of 7773, 11,500 and 13,046 unigenes were expressed differentially at 12, 24 and 48 h post thermal stress, respectively. The expression levels of key DEGs as revealed by RNA-seq were confirmed by quantitative real-time PCR. GO and KEGG enrichment analyses showed that the DEGs were mainly associated with apoptosis, NF-kappa B signaling pathway, TNF signaling pathway and RIG-I-like receptor signaling pathway. Among the DEGs, 40 were candidate heat stress response-related genes and 169 were identified to be involved in antioxidant defense, cell detoxification, protein metabolism and endoplasmic reticulum stress responses. It seemed that ark shells may adapt to short term thermal stress through regulation of protein metabolism, DNA replication and anti-apoptotic system. However, if the stress sustains, it may cause irreparable injury gradually in the animals due to oxygen limitation and metabolic dysregulation. Noteworthily, the expression of DEGs involved in protein biosynthesis and proteolysis was significantly elevated in ark shells under heat stress. These findings may provide preliminary insights into the molecular response of ark shells to acute thermal stress and lay the groundwork for marker-assisted selection of heat-resistant strains in S. subcrenata.

Impact of Ringer's Solution Challenge Stress to Immunostimulatory Experiment, Insights From Japanese Flounder

Ringer's or phosphate buffer saline (PBS) solution buffer usually was used as dilution butter in intraperitoneal injection. Stress could activate immune response, inflammatory response and glycogen metabolic process. The impact of solution buffer as a stressor to immune system was ignored in immunostimulatory experiment. In this report, we tested the hypothesis that the innate immune response and glycogen metabolic process were altered when it were challenged with Ringer's in Japanese flounder (Paralichthys olivaceus). RNA-seq was performed after challenge with Ringer's at 8 h and 48 h. The data revealed that the expression profiles of blood, gill, and kidney were significantly changed. Differentially expressed genes (DEGs) were identified, and energy metabolic and immune-related genes were up-regulated or down-regulated obviously. GO and KEGG analyses showed that DEGs were mainly enriched in innate immune terms and pathways. Weighted gene co-expression networks analysis (WGCNA) also indicated the highest association module with stress. A total of 16 genes were detected in the gray module, which were immune-related and metabolic-related genes. These results provided fundamental information on intraperitoneal injection with solution buffer. It offered useful clues to further explore the functional mechanism of stress and immunity.

Host Defense Effectors Expressed by Hemocytes Shape the Bacterial Microbiota From the Scallop Hemolymph

The interaction between host immune response and the associated microbiota has recently become a fundamental aspect of vertebrate and invertebrate animal health. This interaction allows the specific association of microbial communities, which participate in a variety of processes in the host including protection against pathogens. Marine aquatic invertebrates such as scallops are also colonized by diverse microbial communities. Scallops remain healthy most of the time, and in general, only a few species are fatally affected on adult stage by viral and bacterial pathogens. Still, high mortalities at larval stages are widely reported and they are associated with pathogenic Vibrio. Thus, to give new insights into the interaction between scallop immune response and its associated microbiota, we assessed the involvement of two host antimicrobial effectors in shaping the abundances of bacterial communities present in the scallop Argopecten purpuratus hemolymph. To do this, we first characterized the microbiota composition in the hemolymph from non-stimulated scallops, finding both common and distinct bacterial communities dominated by the Proteobacteria, Spirochaetes and Bacteroidetes phyla. Next, we identified dynamic shifts of certain bacterial communities in the scallop hemolymph along immune response progression, where host antimicrobial effectors were expressed at basal level and early induced after a bacterial challenge. Finally, the transcript silencing of the antimicrobial peptide big defensin ApBD1 and the bactericidal/permeability-increasing protein ApLBP/BPI1 by RNA interference led to an imbalance of target bacterial groups from scallop hemolymph. Specifically, a significant increase in the class Gammaproteobacteria and the proliferation of Vibrio spp. was observed in scallops silenced for each antimicrobial. Overall, our results strongly suggest that scallop antimicrobial peptides and proteins are implicated in the maintenance of microbial homeostasis and are key molecules in orchestrating host-microbiota interactions. This new evidence depicts the delicate balance that exists between the immune response of A. purpuratus and the hemolymph microbiota.

Structure and functional analysis reveal an important regulated role of arginine kinase in Patinopecten yessoensis under low pH stress

Arginine kinase (AK), an important member of the phosphokinase family, is involved in temporal and spatial adenosine triphosphate (ATP) buffering systems. AK plays an important role in physiological function and metabolic regulations, in particular tissues with high and fluctuating energy demands. In present study, four AK genes were firstly identified from Yesso scallop (Patinopecten yessoensis) genome, respectively named PyAK1-4. PyAKs have highly conserved structures with a six-exon/five-exon structure, except for PyAK3. PyAK3 contains an unusual two-domain structure and a "bridge intron" between the two domains, which may originate from gene duplication and subsequent fusion. Phylogenetic analysis showed that all PyAKs belonged to an AK supercluster together with other AK proteins from Mollusca, Platyhelminthes, Arthropoda, and Nematode. A transcriptome database demonstrated that PyAK3 and PyAK4 were the main functional executors with high expression level during larval development and in adult tissues, while PyAK1 and PyAK2 were expressed at a low level. Furthermore, both PyAK2 and PyAK3 showed notably high expression in the male gonad, and PyAK4 was broadly expressed in almost all tissues with the highest level in striated muscle, indicating a tissue-specific expression pattern of PyAKs. In addition, quantitative real-time PCR results demonstrated that the expression of PyAK2, PyAK3 and PyAK4 were significantly upregulated in response to pH stress, especially in an extremely acidifying condition (pH 6.5), revealing the possible involvement of PyAKs in energetic homeostasis during environmental changes. Collectively, a comprehensive analysis of PyAKs was conducted in P. yessoensis. The diversity of PyAKs and their specific expression patterns promote a better understanding of energy metabolism in the growth, development and environmental response of P. yessoensis.

How will farmed populations of freshwater fish deal with the extreme climate scenario in 2100? Transcriptional responses of Colossoma macropomum from two Brazilian climate regions

Tambaqui (Colossoma macropomum Cuvier, 1818) is an endemic fish of the Amazon and Orinoco basins, and it is the most economically important native species in Brazil being raised in five climatically distinct regions. In the face of current global warming, environmental variations in farm ponds represent additional challenges that may drive new adaptive regional genetic variations among broodstocks of tambaqui. In an experimental context based on the high-emission scenario of the 5th Intergovernmental Panel on Climate Change (IPCC) report, we used two farmed tambaqui populations to test this hypothesis. RNA-seq transcriptome analysis was performed in the liver of juvenile tambaqui from northern (Balbina Experimental Station, Balbina, AM) and southeastern (Brumado Fish Farming, Mogi Mirim, SP) Brazilian regions kept for 30 days in artificial environmental rooms mimicking the current and extreme climate scenarios. Three Illumina MiSeq runs produced close to 120 million 500 bp paired-end reads; 191,139 contigs were assembled with N50 = 1595. 355 genes were differentially expressed for both populations in response to the extreme scenario. After enrichment analysis, each population presented a core set of genes to cope with climate change. Northern fish induced genes related to the cellular response to stress, activation of MAPK activity, response to unfolded protein, protein metabolism and cellular response to DNA damage stimuli. Genes biologically involved in regulating cell proliferation, protein stabilisation and protein ubiquitination for degradation through the ubiquitin-proteasome system were downregulated. Genes associated with biological processes, including the cellular response to stress, MAPK cascade activation, homeostatic processes and positive regulation of immune responses were upregulated in southeastern fish. The downregulated genes were related to cytoskeleton organisation, energy metabolism, and the regulation of transcription and biological rhythms. Our findings reveal the signatures of promising candidate genes involved in the regional plasticity of each population of tambaqui in dealing with upcoming climate changes.

Immunity-related genes and signaling pathways under hypoxic stresses in Haliotis diversicolor: a transcriptome analysis

Due to increased temperatures and aquaculture density, thermal and hypoxia stresses have become serious problems for the aquaculture of abalone Haliotis diversicolor. Stresses lead to immunosuppression, which can cause severe negative impacts on aquaculture farms. To study the mechanism of immunosuppression after hypoxia stress and bacterial challenge, transcriptomes of H. diversicolor hemocytes involved in immunity were profiled. A total of 307,395,572 clean reads were generated and assembled into 99,774 unigenes. KEGG analysis indicated that 225 unigenes with immunologic function were mapped into immune-related pathways. Expression of 41 unigenes measured by quantitative real-time PCR (qRT-PCR) showed consistent results with that of transcriptome analysis. When exposure challenge of Vibrio parahaemolyticus, it is indicated that the PI3K-AKT, MAPK, NF-κB and P53 signal pathways were involved in the hypoxia-induced immunosuppression of H. diversicolor. Furthermore, when the AKT gene (HdAKT) was inhibited by double-stranded RNA (dsRNA), expression levels of HdAKT was lower than the blank and control group in hemocytes at 4 h, 12 h and 24 h (p < 0.05).

De novo assembly, characterization of tissue-specific transcriptomes and identification of immune related genes from the scallop Argopecten purpuratus

The scallop Argopecten purpuratus is one of the most economically important cultured mollusks on the coasts from Chile and Peru but its production has declined, in part, due to the emergence of mass mortality events of unknown origin. Driven by this scenario, increasing progress has been made in recent years in the comprehension of immune response mechanisms in this species. However, it is still not entirely understood how different mucosal interfaces participate and cooperate with the immune competent cells, the hemocytes, in the immune defense. Thus, in this work we aimed to characterize the transcriptome of three tissues with immune relevance from A. purpuratus by next-generation sequencing and de novo transcriptome assembly. For this, 18 cDNA libraries were constructed from digestive gland, gills and hemocytes tissues of scallops from different immune conditions and sequenced by the Illumina HiSeq4000 platform. A total of 967.964.884 raw reads were obtained and 967.432.652 clean reads were generated. The clean reads were de novo assembled into 46.601 high quality contigs and 32.299 (69.31%) contigs were subsequently annotated. In addition, three de novo specific assemblies were performed from clean reads obtained from each tissue cDNA libraries for their comparison. Gene ontology (GO) and KEGG analyses revealed that annotated sequences from digestive gland, gills and hemocytes could be classified into both general and specific subcategory terms and known biological pathways, respectively, according to the tissue nature. Finally, several immune related candidate genes were identified, and the differential expression of tissue-specific genes was established, suggesting they could display specific roles in the host defense. The data presented in this study provide the first insight into the tissue specific transcriptome profiles of A. purpuratus, which should be considered for further research on the interplay between the hemocytes and mucosal immune responses.

Transcriptome analysis and discovery of genes involved in immune pathways in Solen strictus (Gould, 1861) under Vibrio anguillarum

The Solen strictus (Gould, 1861) has been recognized as an important marine economic bivalve in Eastern and Southeast Asia. To gain a better understanding of the S. strictus immune system and its related genes in response to bacterial infections, we performed a comparative gene transcription analysis from S. strictus with Vibrio anguillarum through RNA-Seq technology, meanwhile the differentially expressed genes (DEGs) were investigated. After assembly, a total of 195,774 transcripts with an average length of 996 bp were obtained. Total 153,038 unigenes were annotated in the nr, Swiss-Prot, KEGG, COG, KOG, GO and Pfam databases, and 56,597 unigenes (36.98%) were annotated in at least one database. After bacterial challenge, there were 1588 significant differentially expressed genes (DEGs) between the challenged and control groups, including 999 up-regulated and 589 down-regulated genes. All the DEGs were classified into three gene ontology categories, and allocated to 225 KEGG pathways. Immune-related genes were detected from immune system pathways among the top 20 enriched pathways, such as Toll-like receptor signaling, RIG-I-like receptor signaling and NOD-like receptor signaling pathway. In addition, 56,079 potential simple sequence repeats (SSRs) and 1,031,521 candidate single nucleotide polymorphisms (SNPs) were detected and identified in the S. strictus transcriptome. Results of the present study will provide valuable theoretical resources for future genetic and genomic research on S. strictus. The research results will be helpful for improving the efficiency and quality of artificial breeding, establishing genetic linkage map, and enhancing health management for this species.