Network analysis of oyster transcriptome revealed a cascade of cellular responses during recovery after heat shock

Transcriptomic Response of the Atlantic Surfclam (Spisula solidissima) to Acute Heat Stress

There is clear evidence that the oceans are warming due to anthropogenic climate change, and the northeastern coast of USA contains some of the fastest warming areas. This warming is projected to continue with serious biological and social ramifications for fisheries and aquaculture. One species particularly vulnerable to warming is the Atlantic surfclam (Spisula solidissima). The surfclam is a critically important species, linking marine food webs and supporting a productive, lucrative, and sustainable fishery. The surfclam is also emerging as an attractive candidate for aquaculture diversification, but the warming of shallow coastal farms threatens the expansion of surfclam aquaculture. Little is known about the adaptive potential of surfclams to cope with ocean warming. In this study, the surfclam transcriptome under heat stress was examined. Two groups of surfclams were subjected to heat stress to assess how artificial selection may alter gene expression. One group of clams had been selected for greater heat tolerance (HS) and the other was composed of random control clams (RC). After a 6-h exposure to 16 or 29 °C, gill transcriptome expression profiles of the four temperature/group combinations were determined by RNA sequencing and compared. When surfclams experienced heat stress, they exhibited upregulation of heat shock proteins (HSPs), inhibitors of apoptosis (IAPs), and other stress-response related genes. RC clams differentially expressed 1.7 times more genes than HS clams, yet HS clams had a stronger response of key stress response genes, including HSPs, IAPs, and genes involved with mitigating oxidative stress. The findings imply that the HS clams have a more effective response to heat stress after undergoing the initial selection event due to genetic differences created by the selection, epigenetic memory of the first heat shock, or both. This work provides insights into how surfclams adapt to heat stress and should inform future breeding programs that attempt to breed surfclam for greater heat tolerance, and ultimately bring greater resiliency to shellfish farms.

Investigation of Gene Networks in Three Components of Immune System Provides Novel Insights into Immune Response Mechanisms against Edwardsiella tarda Infection in Paralichthys olivaceus

As a quintessential marine teleost, Paralichthys olivaceus demonstrates vulnerability to a range of pathogens. Long-term infection with Edwardsiella tarda significantly inhibits fish growth and even induces death. Gills, blood, and kidneys, pivotal components of the immune system in teleosts, elicit vital regulatory roles in immune response processes including immune cell differentiation, diseased cell clearance, and other immunity-related mechanisms. This study entailed infecting P. olivaceus with E. tarda for 48 h and examining transcriptome data from the three components at 0, 8, and 48 h post-infection employing weighted gene co-expression network analysis (WGCNA) and protein-protein interaction (PPI) network analysis. Network analyses revealed a series of immune response processes after infection and identified multiple key modules and key, core, and hub genes including xpo1, src, tlr13, stat1, and mefv. By innovatively amalgamating WGCNA and PPI network methodologies, our investigation facilitated an in-depth examination of immune response mechanisms within three significant P. olivaceus components post-E. tarda infection. Our results provided valuable genetic resources for understanding immunity in P. olivaceus immune-related components and assisted us in further exploring the molecular mechanisms of E. tarda infection in teleosts.

Transcriptome Analysis of Reciprocal Hybrids Between Crassostrea gigas and C. angulata Reveals the Potential Mechanisms Underlying Thermo-Resistant Heterosis

Heterosis, also known as hybrid vigor, is widely used in aquaculture, but the molecular causes for this phenomenon remain obscure. Here, we conducted a transcriptome analysis to unveil the gene expression patterns and molecular bases underlying thermo-resistant heterosis in Crassostrea gigas ♀ × Crassostrea angulata ♂ (GA) and C. angulata ♀ × C. gigas ♂ (AG). About 505 million clean reads were obtained, and 38,210 genes were identified, of which 3779 genes were differentially expressed between the reciprocal hybrids and purebreds. The global gene expression levels were toward the C. gigas genome in the reciprocal hybrids. In GA and AG, 95.69% and 92.00% of the differentially expressed genes (DEGs) exhibited a non-additive expression pattern, respectively. We observed all gene expression modes, including additive, partial dominance, high and low dominance, and under- and over-dominance. Of these, 77.52% and 50.00% of the DEGs exhibited under- or over-dominance in GA and AG, respectively. The over-dominance DEGs common to reciprocal hybrids were significantly enriched in protein folding, protein refolding, and intrinsic apoptotic signaling pathway, while the under-dominance DEGs were significantly enriched in cell cycle. As possible candidate genes for thermo-resistant heterosis, GRP78, major egg antigen, BAG, Hsp70, and Hsp27 were over-dominantly expressed, while MCM6 and ANAPC4 were under-dominantly expressed. This study extends our understanding of the thermo-resistant heterosis in oysters.

Protective role of trehalose in the Namib desert ant, Ocymyrmex robustior

Over recent decades, increasing attention has been paid to how low-molecular-weight molecules affect thermal tolerance in animals. Although the disaccharide sugar trehalose is known to serve as a thermal protectant in unicellular organisms, nothing is known about its potential role in insects. In this study, we investigated the effect of trehalose on heat tolerance in the Namib desert ant, Ocymyrmex robustior, one of the most thermotolerant animals found in terrestrial ecosystems. First, we tested whether a trehalose-supplemented diet increased worker survival following exposure to heat stress. Second, we assessed the degree of protein damage by comparing protein aggregation levels for trehalose-supplemented workers and control workers. Third, we compared the expression levels of three genes involved in trehalose metabolism. We found that trehalose supplementation significantly enhanced worker heat tolerance, increased metabolic levels of trehalose and reduced protein aggregation under conditions of heat stress. Expression levels of the three genes varied in a manner that was consistent with the maintenance of trehalose in the hemolymph and tissues under conditions of heat stress. Altogether, these results suggest that increased trehalose concentration may help protect Namib desert ant individuals against heat stress. More generally, they highlight the role played by sugar metabolites in boosting tolerance in extremophiles.

Dynamic transcriptome analysis reveals the gene network of gonadal development from the early history life stages in dwarf surfclam Mulinia lateralis

BACKGROUND

Gonadal development is driven by a complex genetic cascade in vertebrates. However, related information remains limited in molluscs owing to the long generation time and the difficulty in maintaining whole life cycle in the lab. The dwarf surfclam Mulinia lateralis is considered an ideal bivalve model due to the short generation time and ease to breed in the lab.

RESULTS

To gain a comprehensive understanding of gonadal development in M. lateralis, we conducted a combined morphological and molecular analysis on the gonads of 30 to 60 dpf. Morphological analysis showed that gonad formation and sex differentiation occur at 35 and 40-45 dpf, respectively; then the gonads go through gametogenic cycle. Gene co-expression network analysis on 40 transcriptomes of 35-60 dpf gonads identifies seven gonadal development-related modules, including two gonad-forming modules (M6, M7), three sex-specific modules (M14, M12, M11), and two sexually shared modules (M15, M13). The modules participate in different biological processes, such as cell communication, glycan biosynthesis, cell cycle, and ribosome biogenesis. Several hub transcription factors including SOX2, FOXZ, HSFY, FOXL2 and HES1 are identified. The expression of top hub genes from sex-specific modules suggests molecular sex differentiation (35 dpf) occurs earlier than morphological sex differentiation (40-45 dpf).

CONCLUSION

This study provides a deep insight into the molecular basis of gonad formation, sex differentiation and gametogenesis in M. lateralis, which will contribute to a comprehensive understanding of the reproductive regulation network in molluscs.

Identification of stress-related genes by co-expression network analysis based on the improved turbot genome

Turbot (Scophthalmus maximus), commercially important flatfish species, is widely cultivated in Europe and China. With the continuous expansion of the intensive breeding scale, turbot is exposed to various stresses, which greatly impedes the healthy development of turbot industry. Here, we present an improved high-quality chromosome-scale genome assembly of turbot using a combination of PacBio long-read and Illumina short-read sequencing technologies. The genome assembly spans 538.22 Mb comprising 27 contigs with a contig N50 size of 25.76 Mb. Annotation of the genome assembly identified 104.45 Mb repetitive sequences, 22,442 protein-coding genes and 3,345 ncRNAs. Moreover, a total of 345 stress responsive candidate genes were identified by gene co-expression network analysis based on 14 published stress-related RNA-seq datasets consisting of 165 samples. Significantly improved genome assembly and stress-related candidate gene pool will provide valuable resources for further research on turbot functional genome and stress response mechanism, as well as theoretical support for the development of molecular breeding technology for resistant turbot varieties.

Measurement(s)	whole genome sequencing
Technology Type(s)	PacBio long-read and Illumina short-read sequencing technologies

Single-molecule long-read (SMRT) transcriptome sequencing of Mercenaria mercenaria reveals a powerful anti-apoptotic system critical for air exposure endurance

Mercenaria mercenaria is an economically important clam species and exhibits an outstanding resistance to multiple environmental stressors. However, our understanding of their stress adaptability is limited due to a lack of genomic information, such as transcriptome resources. In this study, single-molecule long-read (SMRT) mRNA sequencing was performed to obtain the full-length gill transcriptome reference sequences of M. mercenaria under air exposure stress. In all, 14.5 G subreads were obtained and assembled into 64,603 unigenes, among which 50,613 were successfully annotated. Additionally, 56,295 SSRs, 1457 transcription factors, and 5924 lncRNAs were identified in M. mercenaria transcriptome. Furthermore, numerous apoptosis-related transcripts were identified according to Swiss-Prot annotation and their numbers were counted. We also found that most apoptosis-related transcripts exhibited typical domains of a certain protein family through conserved domain prediction. Additionally, eight typical sequences related to apoptosis pathway were detected by RT-PCR, with the aim to show the sequential variation of gene expression levels under air exposure. These results implied that the complicated apoptosis system, especially the powerful anti-apoptotic system was critical for M. mercenaria to endure air exposure.

The Members of the Highly Diverse Crassostrea gigas Integrin Family Cooperate for the Generation of Various Immune Responses

Studies on invertebrate immune receptors can provide insights into characteristics specific to innate immune system. Here, eight α and three β integrins are identified from an invertebrate, the Pacific oyster Crassostrea gigas, and their possible immune functions are studied. Oyster α/β integrins exhibit a higher degree of sequence and structural variability than the members from Homo sapiens and Drosophila melanogaster. The analysis reveals that oyster RGD- and laminin-binding receptor homologs are present in the phylogenetic tree of α integrins, but the other six oyster α integrins mainly form a species-specific branch; meanwhile, oyster β integrins are clustered with insect β integrins but distinct from a member from the mollusk Biomphalaria glabrata. Although phylogenetically lacking the important α integrin branches of LDV-binding, PS3-type, and αI-containing integrins, oyster integrins can bind to most ECM ligands, including RGDCP, LDVCP, GFOGERCP, and laminin protein in a distinct binding pattern. Besides, oyster integrins are distributed in different hemocyte subpopulations, while only specific integrins are selectively involved in hemocyte phagocytosis, migration, and encapsulation, and some of them participate in more than one immune response in a sophisticated pattern. Especially, oyster β integrins are arranged in the core to mediate complex immune responses, unlike the counterparts in humans that mainly depend on αI-containing integrins to incite immune reactions. This study represents the first comprehensive attempt to reveal the structural and evolutionary features of the integrin family and their involvement in cellular immune responses in the non-model invertebrate C. gigas and sheds light on the characteristics specific to the innate immune system in the integrin family.

Identification of HSF1 Target Genes Involved in Thermal Stress in the Pacific Oyster Crassostrea gigas by ChIP-seq

The Pacific oyster Crassostrea gigas, a commercially important species inhabiting the intertidal zone, facing enormous temperature fluctuations. Therefore, it is important to identify candidate genes and key regulatory relationships associated with thermal tolerance, which can aid the molecular breeding of oysters. Heat shock transcription factor 1 (HSF1) plays an important role in the thermal stress resistance. However, the regulatory relationship between the expansion of heat shock protein (HSP) HSP 70 and HSF1 is not yet clear in C. gigas. In this study, we analyzed genes regulated by HSF1 in response to heat shock by chromatin immunoprecipitation followed by sequencing (ChIP-seq), determined the expression patterns of target genes by qRT-PCR, and validated the regulatory relationship between one HSP70 and HSF1. We found 916 peaks corresponding to HSF1 binding sites, and these peaks were annotated to the nearest genes. In Gene Ontology analysis, HSF1 target genes were related to signal transduction, energy production, and response to biotic stimulus. Four HSP70 genes, two HSP40 genes, and one small HSP gene exhibited binding to HSF1. One HSP70 with a binding site in the promoter region was validated to be regulated by HSF1 under heat shock. These results provide a basis for future studies aimed at determining the mechanisms underlying thermal tolerance and provide insights into gene regulation in the Pacific oyster.

Transcriptomic Profiling Provides Insights into Inbreeding Depression in Yesso Scallop Patinopecten yessoensis

Inbreeding often causes a decline in biological fitness, known as inbreeding depression. In genetics study, inbreeding coefficient f gives the proportion by which the heterozygosity of an individual is reduced by inbreeding. With the development of high-throughput sequencing, researchers were able to perform deep approaches to investigate which genes are affected by inbreeding and reveal some molecular underpinnings of inbreeding depression. As one commercially important species, Yesso scallop Patinopecten yessoensis confront the same dilemma of inbreeding depression. To examine how inbreeding affects gene expression, we compared the transcriptome of two experimentally selfing families with inbreeding coefficient f reached 0.5 as well as one natural population (f ≈ 0) of P. yessoensis. A total of 24 RNA-Seq libraries were constructed using scallop adductor muscle, and eventually 676.56 M (96.85%) HQ reads were acquired. Based on differential gene analysis, we were able to identify nine common differentially expressed genes (DEGs) across the top-ranked 30 DEGs in both selfing families in comparation with the natural population. Remarkable, through weighted gene co-expression network analysis (WGCNA), five common DEGs were found enriched in the most significant inbreeding related functional module M14 (FDR = 1.64E-156), including SREBP1, G3BP2, SBK1, KIAA1161, and AATs-Glupro. These five genes showed significantly higher expression in self-bred progeny. Suggested by the genetic functional analysis, up-regulated SREBP1, G3BP2, and KIAA1161 may suggest a perturbing lipid metabolism, a severe inframammary reaction or immune response, and a stress-responsive behavior. Besides, the significant higher SBK1 and AATs-Glupro may reflect the abnormal cellular physiological situation. Together, these genetic aberrant transcriptomic performances may contribute to inbreeding depression in P. yessoensis, deteriorating the stress tolerance and survival phenotype in self-bred progeny. Our results would lay a foundation for further comprehensive understanding of bivalve inbreeding depression, which may potentially benefit the genetic breeding for scallop aquaculture.

RNAi based transcriptome suggests genes potentially regulated by HSF1 in the Pacific oyster Crassostrea gigas under thermal stress

Background

The Pacific oyster Crassostrea gigas is an important fishery resource that is sensitive to temperature fluctuations. Thus, it has evolved a protection mechanism against heat stress by increasing the expression of the gene coding for heat shock protein (HSP) 70 under elevated temperatures. In other animals, heat shock response is a transcriptional response driven by the heat shock transcription factor 1 (HSF1) and thermal stress can trigger HSP70 expression to protect the organism via HSF1. However, the regulatory relationship between HSF1 and HSP remains unclear in Pacific oyster. Therefore, in the present study, we examined the transcriptomic response of several to thermal stress following HSF1 interference.

Results

We identified 150 genes responsive to heat shock including seven HSP genes, six of which belonging to the group of 17 HSP genes enriched in response to heat shock, according to weighted gene co-expression network analysis (WGCNA). The other gene was enriched in the module correlated with HSF1 interference. In addition, we found 48 and 47 genes that were upregulated and downregulated by HSF1 in response to heat shock, respectively. In the upregulated genes, we identified one HSP70 potentially regulated by HSF1 in response to heat shock. Furthermore, based on differentially expressed genes and WGCNA analyses, we found that the hypoxia signaling pathway was enriched under heat shock conditions. Five genes were then selected to detect dynamic changes through time. The results suggested that gene expression was correlated with HSF1 expression. The regulation of HSP70 by HSF1 was preliminarily confirmed by binding site predictions and by a dual luciferase assay.

Conclusions

Our results revealed that the expression of HSP70 and HSP20 was initially triggered after 2 h of heat shock, and one of the HSP70 genes was potentially regulated by HSF1. From these results, it is evident that not all heat-inducible genes were triggered simultaneously in response to heat shock stress. Overall, the results revealed a possible HSF1–HSP regulatory relationship in Pacific oyster, providing valuable information on the mechanisms of thermal tolerance in this commercially important oyster.

Electronic supplementary material

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

Transcriptomic responses to extreme low salinity among locally adapted populations of Olympia oyster (Ostrea lurida)

The Olympia oyster (Ostrea lurida) is a foundation species inhabiting estuaries along the North American west coast. In California estuaries, O. lurida is adapted to local salinity regimes and populations differ in low salinity tolerance. In this study, oysters from three California populations were reared for two generations in a laboratory common garden and subsequently exposed to low salinity seawater. Comparative transcriptomics was then used to understand species-level responses to hyposmotic stress and population-level mechanisms underlying divergent salinity tolerances. Gene expression patterns indicate Olympia oysters are sensitive to hyposmotic stress: All populations respond to low salinity by up-regulating transcripts indicative of protein unfolding, DNA damage and cell cycle arrest after sub-lethal exposure. Among O. lurida populations, transcriptomic profiles differed constitutively and in response to low salinity. Despite two generations in common-garden conditions, transcripts encoding apoptosis modulators were constitutively expressed at significantly different levels in the most tolerant population. Expression of cell death regulators may facilitate cell fate decisions when salinity declines. Following low salinity exposure, oysters from the more tolerant population expressed a small number of mRNAs at significantly higher levels than less tolerant populations. Proteins encoded by these transcripts regulate ciliary activity within the mantle cavity and may function to prolong valve closure and reduce mortality in low salinity seawater. Collectively, gene expression patterns suggest sub-lethal impacts of hyposmotic stress in Olympia oysters are considerable and that even oysters with greater low salinity tolerance may be vulnerable to future freshwater flooding events.

The transcriptional response of the Pacific oyster Crassostrea gigas against acute heat stress

The Pacific oyster, Crassostrea gigas, has evolved sophisticated mechanisms to adapt the changing ambient conditions, and protect themselves from stress-induced injuries. In the present study, the expression profiles of mRNA transcripts in the haemocytes of oysters under heat stress were examined to reveal the possible mechanism of heat stress response. There were 23,315, 23,904, 23,123 and 23,672 transcripts identified in the haemocytes of oysters cultured at 25 °C for 0, 6, 12, and 24 h (designed as B, H6, H12, H24), respectively. And 22,330 differentially expressed transcripts (DTs) were yielded in the pairwise comparisons between the above four samples, which corresponded to 8074 genes. There were 9, 12 and 22 Gene Ontology (GO) terms identified in the DT pairwise comparison groups of H6_B, H12_H6 and H24_H12, respectively, and the richest GO terms in biological process category were cellular catabolic process, translational initiation and apoptotic process, respectively. There were 108, 102 and 102 KEGG pathways successfully retrieved from DTs comparison groups DTH6_B, DTH12_H6 and DTH24_H12, respectively, among which 93 pathways were shared by all three comparison groups, and most of them were related to metabolism of protein, carbohydrate and fat. The expression patterns of 12 representative heat stress response-relevant genes detected by quantitative real-time PCR (qRT-PCR) were similar to those obtained from transcriptome analysis. By flow cytometric analysis, the apoptosis rate of haemocytes increased significantly after oysters were treated at 25 °C for 24 h and recovered at 4 °C for 12 h (p < 0.05) and 36 h (p < 0.01), and it also increased significantly when the heat treatment lasted to 60 h (p < 0.01). The present results indicated that, when oysters encountered short term heat stress, the expression of genes related to energy metabolism, as well as unfolded protein response (UPR) and anti-apoptotic system, were firstly regulated to maintain basic life activities, and then a large number of genes involved in stabilizing protein conformation and facilitating further protein refolding were activated to repair the stress injury. However, the stress injury gradually became irreparable with the stress persisting, and apoptosis was activated when the heat treatment prolonged to 24 h. The information was useful to better understand the molecular mechanism of heat stress response and develop strategies for the improvement of oyster survival rate during summer high-temperature period.

Conservation and divergence of mitochondrial apoptosis pathway in the Pacific oyster, Crassostrea gigas

Apoptosis is considered a crucial part of the host defense system in oysters according to previous reports; however, the exact process by which this occurs remains unclear. Besides, mitochondrial apoptosis is the primary method of apoptosis in vertebrate cells, but has been poorly studied in invertebrates and is quite controversial. In this study, we investigated the molecular mechanism of mitochondrial apoptosis in the Pacific oyster Crassostrea gigas. Notably, we show that most key elements involved in the vertebrate mitochondrial apoptosis pathway – including mitochondrial outer membrane permeabilization, cytochrome c release, and caspase activation – are also present in C. gigas. In contrast, the lack of Bcl-2 homology 3-only subfamily members and apoptotic protease activating factor-1 (APAF-1) protein revealed evolutionary diversity from other phyla. Our results support that mitochondrial apoptosis in animals predates the emergence of vertebrates, but suggest that an unexpectedly diverse mitochondrial apoptosis pathway may exist in invertebrates. In addition, our work provided new clues for an improved understanding of how bivalve acclimate themselves to an inconstant environment.

Gene Co-Expression Network Analysis Reveals the Correlation Patterns Among Genes in Euryhaline Adaptation of Crassostrea gigas

The Pacific oyster Crassostrea gigas is a dominant aquaculture species in many intertidal zones throughout the Pacific and Atlantic Oceans and can tolerate a wide range of salinity. Studying the gene expression profiles of oyster gills had found differentially expressed genes (DEGs) involved in salinity tolerance. A systematic study of cellular response to salinity stress may provide insights into the mechanism of acquired salinity tolerance. Here, weighted gene co-expression network analysis (WGCNA) was carried out using RNA-seq data from gill transcriptome in response to different salinity. A total of 25,463 genes were parsed into 22 gene modules, of which 5 gene modules were identified as salinity-related modules. Brown module was the only one significantly correlated with salinity and free amino acids (FAAs) contents, which was associated with cellular metabolism, biosynthesis of amino acids, oxidation reduction, electron transport, nitrogen compound metabolism, and others. The enriched pathways in brown module were mainly about FAAs metabolism. The other four modules were significantly correlated with certain FAAs, and were over-represented in certain salinity. These results indicated that C. gigas triggered different FAAs in different salinity stress. This study represents the first RNA-seq gene network analysis in oysters responding to different salinity stresses. These results provide a systems-level framework to help understand the complexity of cellular process in response to osmotic stress and show the function and regulated genes of different FAAs at the molecular level.

Expression Characterization of Stress Genes Under High and Low Temperature Stresses in the Pacific Oyster, Crassostrea gigas

As a characteristic sessile inhabitant of the intertidal zone, the Pacific oyster Crassostrea gigas occupies one of the most physically stressful environments on earth. With high exposure to terrestrial conditions, oysters must tolerate broad fluctuations in temperature range. However, oysters' cellular and molecular responses to temperature stresses have not been fully characterized. Here, we analyzed oyster transcriptome data under high and low temperatures. We also identified over 30 key temperature stress-responsive candidate genes, which encoded stress proteins such as heat shock proteins and apoptosis-associated proteins. The expression characterization of these genes under short-term cold and hot environments (5 and 35 °C) and long-term cold environments (5 °C) was detected by quantitative real-time PCR. Most of these genes reached expression peaks during the recovery stage after 24 h of heat stress, and these genes were greatly induced around day 3 in long-term cold stress while responded little to short-term cold stress. In addition, in the second heat stress after 2 days of recovery, oysters showed milder expression in these genes and a lower mortality rate, which indicated the existence of plasticity in the oyster's response to heat stress. We confirmed that homeostatic flexibility and anti-apoptosis might be crucial centers of temperature stress responses in oysters. Furthermore, we analyzed stress gene families in 11 different species and found that the linage-specific expansion of stress genes might be implicated in adaptive evolution. These results indicated that both plasticity and evolution played an important role in the stress response adaptation of oysters.

GigaTON: an extensive publicly searchable database providing a new reference transcriptome in the pacific oyster Crassostrea gigas

BACKGROUND

The Pacific oyster, Crassostrea gigas, is one of the most important aquaculture shellfish resources worldwide. Important efforts have been undertaken towards a better knowledge of its genome and transcriptome, which makes now C. gigas becoming a model organism among lophotrochozoans, the under-described sister clade of ecdysozoans within protostomes. These massive sequencing efforts offer the opportunity to assemble gene expression data and make such resource accessible and exploitable for the scientific community. Therefore, we undertook this assembly into an up-to-date publicly available transcriptome database: the GigaTON (Gigas TranscriptOme pipeliNe) database.

DESCRIPTION

We assembled 2204 million sequences obtained from 114 publicly available RNA-seq libraries that were realized using all embryo-larval development stages, adult organs, different environmental stressors including heavy metals, temperature, salinity and exposure to air, which were mostly performed as part of the Crassostrea gigas genome project. This data was analyzed in silico and resulted into 56621 newly assembled contigs that were deposited into a publicly available database, the GigaTON database. This database also provides powerful and user-friendly request tools to browse and retrieve information about annotation, expression level, UTRs, splice and polymorphism, and gene ontology associated to all the contigs into each, and between all libraries.

CONCLUSIONS

The GigaTON database provides a convenient, potent and versatile interface to browse, retrieve, confront and compare massive transcriptomic information in an extensive range of conditions, tissues and developmental stages in Crassostrea gigas. To our knowledge, the GigaTON database constitutes the most extensive transcriptomic database to date in marine invertebrates, thereby a new reference transcriptome in the oyster, a highly valuable resource to physiologists and evolutionary biologists.

Immune and stress responses in oysters with insights on adaptation

Oysters are representative bivalve molluscs that are widely distributed in world oceans. As successful colonizers of estuaries and intertidal zones, oysters are remarkably resilient against harsh environmental conditions including wide fluctuations in temperature and salinity as well as prolonged air exposure. Oysters have no adaptive immunity but can thrive in microbe-rich estuaries as filter-feeders. These unique adaptations make oysters interesting models to study the evolution of host-defense systems. Recent advances in genomic studies including sequencing of the oyster genome have provided insights into oyster's immune and stress responses underlying their amazing resilience. Studies show that the oyster genomes are highly polymorphic and complex, which may be key to their resilience. The oyster genome has a large gene repertoire that is enriched for immune and stress response genes. Thousands of genes are involved in oyster's immune and stress responses, through complex interactions, with many gene families expanded showing high sequence, structural and functional diversity. The high diversity of immune receptors and effectors may provide oysters with enhanced specificity in immune recognition and response to cope with diverse pathogens in the absence of adaptive immunity. Some members of expanded immune gene families have diverged to function at different temperatures and salinities or assumed new roles in abiotic stress response. Most canonical innate immunity pathways are conserved in oysters and supported by a large number of diverse and often novel genes. The great diversity in immune and stress response genes exhibited by expanded gene families as well as high sequence and structural polymorphisms may be central to oyster's adaptation to highly stressful and widely changing environments.

The genome-wide identification of mitogen-activated protein kinase kinase (MKK) genes in Yesso scallop Patinopecten yessoensis and their expression responses to bacteria challenges

Mitogen-activated protein kinase kinases (MKK) are the essential components of the evolutionarily conserved MAPK signaling cascade, which regulates a variety of cellular activities and innate immune responses. Although MKK genes have been extensively studied in various vertebrate and invertebrate species, they have not been systematically characterized in bivalves. In this study, we identified and characterized five MKK genes (PyMKK1/2, PyMKK4, PyMKK5, PyMKK3/6 and PyMKK7) in the Yesso scallop (Patinopecten yessoensis). Phylogenetic and protein structural analyses were conducted to determine their identities and evolutionary relationships. To gain insights into the possible roles of MKK genes during scallop innate immune responses, quantitative real-time PCR (qRT-PCR) was used to investigate their expression profiles during different developmental stages in samples taken from healthy adult tissues and hemocytes after Micrococcus luteus and Vibrio anguillarum bacterial infections. The Yesso scallop MKKs (PyMKKs) were found to have highly conserved structural features compared to the MKK genes from other invertebrate species. Using qRT-PCR analysis, three distinct expression patterns were detected among the PyMKKs over the course of ten different developmental stages. In adult scallops, the majority of the PyMKKs were highly expressed in mantle, gill, muscle and hemocytes. The differential expression patterns of the five PyMKKs after M. luteus (Gram-positive) and V. anguillarum (Gram-negative) bacterial infections suggested their possible involvement in the innate immune response and provide the foundation and resource for the further study on innate immune response of MAPK signal pathway in mollusk.

Next generation sequence analysis of the transcriptome of Sydney rock oysters (Saccostrea glomerata) exposed to a range of environmental stressors

Sydney rock oysters (Saccostrea glomerata) were exposed to environmental stressors at contaminated field sites or in a controlled laboratory setting. RNA seq transcriptome data were generated for the gill and digestive gland using Roche's 454 pyrosequencing technology. 28,685 contigs were de novo assembled which encoded 11,671 different protein products. The data will act as a reference for future studies in ecology, immunology and environmental toxicology.

Sequencing-based gene network analysis provides a core set of gene resource for understanding thermal adaptation in Zhikong scallop Chlamys farreri

Marine organisms are commonly exposed to variable environmental conditions, and many of them are under threat from increased sea temperatures caused by global climate change. Generating transcriptomic resources under different stress conditions are crucial for understanding molecular mechanisms underlying thermal adaptation. In this study, we conducted transcriptome-wide gene expression profiling of the scallop Chlamys farreri challenged by acute and chronic heat stress. Of the 13 953 unique tags, more than 850 were significantly differentially expressed at each time point after acute heat stress, which was more than the number of tags differentially expressed (320-350) under chronic heat stress. To obtain a systemic view of gene expression alterations during thermal stress, a weighted gene coexpression network was constructed. Six modules were identified as acute heat stress-responsive modules. Among them, four modules involved in apoptosis regulation, mRNA binding, mitochondrial envelope formation and oxidation reduction were downregulated. The remaining two modules were upregulated. One was enriched with chaperone and the other with microsatellite sequences, whose coexpression may originate from a transcription factor binding site. These results indicated that C. farreri triggered several cellular processes to acclimate to elevated temperature. No modules responded to chronic heat stress, suggesting that the scallops might have acclimated to elevated temperature within 3 days. This study represents the first sequencing-based gene network analysis in a nonmodel aquatic species and provides valuable gene resources for the study of thermal adaptation, which should assist in the development of heat-tolerant scallop lines for aquaculture.

Cyclic-AMP mediated regulation of ABCB mRNA expression in mussel haemocytes

Background

The multixenobiotic resistance system (MXR) allows aquatic organisms to cope with their habitat despite high pollution levels by over-expressing membrane and intracellular transporters, including the P-glycoprotein (Pgp). In mammals transcription of the ABCB1 gene encoding Pgp is under cAMP/PKA-mediated regulation; whether this is true in mollusks is not fully clarified.

Methodology/Principal Findings

cAMP/PKA regulation and ABCB mRNA expression were assessed in haemocytes from Mediterranean mussels (Mytilus galloprovincialis) exposed in vivo for 1 week to 0.3 ng/L fluoxetine (FX) alone or in combination with 0.3 ng/L propranolol (PROP). FX significantly decreased cAMP levels and PKA activity, and induced ABCB mRNA down-regulation. FX effects were abolished in the presence of PROP. In vitro experiments using haemocytes treated with physiological agonists (noradrenaline and serotonin) and pharmacological modulators (PROP, forskolin, dbcAMP, and H89) of the cAMP/PKA system were performed to obtain clear evidence about the involvement of the signaling pathway in the transcriptional regulation of ABCB. Serotonin (5-HT) decreased cAMP levels, PKA activity and ABCB mRNA expression but increased the mRNA levels for a putative 5-HT₁ receptor. Interestingly, 5-HT₁ was also over-expressed after in vivo exposures to FX. 5-HT effects were counteracted by PROP. Forskolin and dbcAMP increased PKA activity as well as ABCB mRNA expression; the latter effect was abolished in the presence of the PKA inhibitor H89.

Conclusions

This study provides the first direct evidence for the cAMP/PKA-mediated regulation of ABCB transcription in mussels.