Development of a Pacific oyster (Crassostrea gigas) 31,918-feature microarray: identification of reference genes and tissue-enriched expression patterns

m6A Profile Dynamics Indicates Regulation of Oyster Development by m6A-RNA Epitranscriptomes

The N6-methylation of RNA adenosines (N6-methyladenosine, m6A) is an important regulator of gene expression with critical implications in vertebrate and insect development. However, the developmental significance of epitranscriptomes in lophotrochozoan organisms remains unknown. Using methylated RNA immunoprecipitation sequencing (MeRIP-seq), we generated transcriptome-wide m6A-RNA methylomes covering the entire development of the oyster from oocytes to juveniles. Oyster RNA classes display specific m6A signatures, with messenger RNAs (mRNAs) and long non-coding RNAs (lncRNAs) exhibiting distinct profiles and being highly methylated compared to transposable element (TE) transcripts. Epitranscriptomes are dynamic and correspond to the chronological steps of development (cleavage, gastrulation, organogenesis, and metamorphosis), with minimal mRNA and lncRNA methylation at the morula stage followed by a global increase. mRNA m6A levels are correlated with transcript levels, and shifts in methylation profiles correspond to expression kinetics. Differentially methylated transcripts cluster according to embryo-larval stages and bear the corresponding developmental functions (cell division, signal transduction, morphogenesis, and cell differentiation). The m6A level of TE transcripts is also regulated and peaks during the gastrulation. We demonstrate that m6A-RNA methylomes are dynamic and associated with gene expression regulation during oyster development. The putative epitranscriptome implication in the cleavage, maternal-to-zygotic transition, and cell differentiation in a lophotrochozoan model brings new insights into the control and evolution of developmental processes.

Isolation of multidrug-resistant Escherichia coli and Salmonella spp. from sulfonamide-treated diarrheic calves

Background and Aim

The bovine industry is threatened by one of the most serious and deadly enteric diseases, calf diarrhea, particularly in developing nations like Bangladesh. In this context, bacterial resistance to antimicrobial drugs and its detrimental consequences have become a critical public health issue that is difficult to address globally. This study aimed to isolate and identify Escherichia coli and Salmonella spp. with their antibiogram and antibiotic resistance gene detection from sulfonamide-treated diarrheic calves.

Materials and Methods

Twelve diarrheic calves suffering from calf diarrhea in a dairy farm were selected and a total of 36 fecal samples were aseptically collected directly from rectum before, during, and at the end of treatment for each calf to determine the total viable count, total E. coli count and total Salmonella count. A polymerase chain reaction was used for the specific detection of E. coli and Salmonella genus targeting fliC and invA genes, respectively. Antibiotic sensitivity test of the isolated E. coli and Salmonella spp. were performed by the disk diffusion method for eight commonly used antibiotics.

Results

A total of 36 E. coli (100%) and 12 Salmonella spp. (33%) were isolated from the samples and were confirmed by polymerase chain reaction. Total viable count was found to be ranged from 35 × 10⁷ to 99 × 10¹⁰ colony-forming unit (CFU)/g fecal sample before starting sulfonamide treatment, 34 × 10⁵ to 25 × 10¹⁰ CFU/g during treatment with sulfonamide, and 48 × 10³ to 69 × 10¹⁰ CFU/g immediately after completion of sulfonamide treatment. Total E. coli count was found to be ranged from 4 × 10⁴ to 36 × 10¹⁰ CFU/g, 24 × 10⁴ to 23 × 10⁸ CFU/g, and 13 × 10⁴ to 85 × 10¹⁰ CFU/g, whereas total Salmonella count was found to be ranged from 16 × 10⁶ to 18.5 × 10¹¹ CFU/g, 15 × 10⁴ to 44 × 10⁷ CFU/g, and 13.2 × 10⁵ to 21 × 10¹⁰ CFU/g fecal sample before starting sulfonamide treatment, during treatment with sulfonamide immediately after completion of sulfonamide treatment, respectively. The in vitro antibiotic sensitivity test showed that all the E. coli and Salmonella spp. isolated from diarrheic calves (100%) contained multidrug-resistant (MDR) phenotypes. Escherichia coli isolates were found 100% resistant to amoxicillin (AMX), cefuroxime, cephalexin (CN), erythromycin (ERY), and tetracycline (TET); whereas 94.4%, 86.1%, and 77.8% isolates were resistant to doxycycline (DOX), moxifloxacin (MOF), and gentamycin (GEN), respectively. In case of Salmonella isolates, all were found 100% resistant to AMX, CN, and ERY; whereas 91.7% of resistance was observed for DOX, MOF, cefuroxime, GEN, and TET. Based on the molecular screening of the antibiotic resistance genes, tetA gene was present in 83.3% of the isolated E. coli and 75% of the isolated Salmonella strains, whereas 83.3% E. coli and 79.2% Salmonella isolates contained blaTEM gene.

Conclusion

These findings suggest that MDR E. coli and Salmonella spp. might be responsible for calf scouring, which is challenging to treat with antibiotics or sulfonamide drugs alone. Therefore, it is important to check the antibiotic sensitivity pattern to select a suitable antibiotic for the treatment of calf scoring. A suitable antibiotic or combination of an antibiotic and sulfonamide could be effective against E. coli and Salmonella spp. responsible for calf scouring.

Gonadal transcriptomes associated with sex phenotypes provide potential male and female candidate genes of sex determination or early differentiation in Crassostrea gigas, a sequential hermaphrodite mollusc

Background

In the animal kingdom, mollusca is an important phylum of the Lophotrochozoa. However, few studies have investigated the molecular cascade of sex determination/early gonadal differentiation within this phylum. The oyster Crassostrea gigas is a sequential irregular hermaphrodite mollusc of economic, physiological and phylogenetic importance. Although some studies identified genes of its sex-determining/−differentiating pathway, this particular topic remains to be further deepened, in particular with regard to the expression patterns. Indeed, these patterns need to cover the entire period of sex lability and have to be associated to future sex phenotypes, usually impossible to establish in this sequential hermaphrodite. This is why we performed a gonadal RNA-Seq analysis of diploid male and female oysters that have not changed sex for 4 years, sampled during the entire time-window of sex determination/early sex differentiation (stages 0 and 3 of the gametogenetic cycle). This individual long-term monitoring gave us the opportunity to explain the molecular expression patterns in the light of the most statistically likely future sex of each oyster.

Results

The differential gene expression analysis of gonadal transcriptomes revealed that 9723 genes were differentially expressed between gametogenetic stages, and 141 between sexes (98 and 43 genes highly expressed in females and males, respectively). Eighty-four genes were both stage- and sex-specific, 57 of them being highly expressed at the time of sex determination/early sex differentiation. These 4 novel genes including Trophoblast glycoprotein-like, Protein PML-like, Protein singed-like and PREDICTED: paramyosin, while being supported by RT-qPCR, displayed sexually dimorphic gene expression patterns.

Conclusions

This gonadal transcriptome analysis, the first one associated with sex phenotypes in C. gigas, revealed 57 genes highly expressed in stage 0 or 3 of gametogenesis and which could be linked to the future sex of the individuals. While further study will be needed to suggest a role for these factors, some could certainly be original potential actors involved in sex determination/early sex differentiation, like paramyosin and could be used to predict the future sex of oysters.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07838-1.

Vitellogenin gene expression in marine mussels exposed to ethinylestradiol: No induction at the transcriptional level

Vitellogenin (Vtg), a large multidomain protein precursor of egg-yolk proteins, is used as an endocrine disruption biomarker in fish, and in the last decades, its use has been extended to invertebrates like mollusks. However, it remains unclear whether invertebrate endocrine system produces Vtg in response to estrogens, like it occurs in oviparous vertebrates. In a previous study, no evidence of induction of Vtg expression at protein level was found in gonads of the marine mussel Mytilus galloprovincialis after exposure to the estrogenic chemical 17α-ethinylestradiol (EE2). In the present follow-up study, it was investigated whether there is any effect of EE2 on Vtg abundance at transcriptional level in M. galloprovincialis gonads. To this aim, RT-qPCR analysis targeting three different domains of Vtg transcript was performed on gonads of mussels that were exposed either 4 or 24 days to 100 ng/L EE2. In addition, several reference genes were analysed and a selection of these for potential use in further RT-qPCR analyses on mussel male and female gonads is provided. Results showed higher expression in females than in males for the three analysed Vtg domains, and no evidence of Vtg mRNA induction due to EE2 either in females or males. The present results, together with those obtained from previous analysis at protein level, support that Vtg is not an adequate biomarker for xenoestrogenicity in marine mussels. Additionally, nucleotide sequences of Vtg transcripts of three closely-related species from Mytilus edulis complex (M. galloprovincialis, M. edulis and M. trossulus) are provided and compared with Vtg sequences from other mollusk species to assess the level of conservation and evolutionary relationships among species.

De novo transcriptome analysis of the mussel Perna viridis after exposure to the toxic dinoflagellate Prorocentrum lima

Diarrheic shellfish poisoning (DSP) toxins are produced by harmful microalgae and accumulate in bivalve mollusks, causing various toxicity. These toxic effects appear to abate with increasing DSP concentration and longer exposure time, however, the underlying mechanisms remain unclear. To explore the underlying molecular mechanisms, de novo transcriptome analysis of the digestive gland of Perna viridis was performed after Prorocentrum lima exposure. RNA-seq analysis showed that 1886 and 237 genes were up- and down-regulated, respectively after 6 h exposure to P. lima, while 265 genes were up-regulated and 217 genes were down-regulated after 96 h compared to the control. These differentially expressed genes mainly involved in Nrf2 signing pathways, immune stress, apoptosis and cytoskeleton, etc. Combined with qPCR results, we speculated that the mussel P. viridis might mainly rely on glutathione S-transferase (GST) and ABC transporters to counteract DSP toxins during short-term exposure. However, longer exposure of P. lima could activate the Nrf2 signaling pathway and inhibitors of apoptosis protein (IAP), which in turn reduced the damage of DSP toxins to the mussel. DSP toxins could induce cytoskeleton destabilization and had some negative impact on the immune system of bivalves. Collectively, our findings uncovered the crucial molecular mechanisms and the regulatory metabolic nodes that underpin the defense mechanism of bivalves against DSP toxins and also advanced our current understanding of bivalve defense mechanisms.

Variation in global transcriptomic response to Perkinsus marinus infection among eastern oyster families highlights potential mechanisms of disease resistance

Dermo disease, caused by the protozoan parasite Perkinsus marinus, negatively impacts wild and cultured Eastern oyster populations, yet our knowledge of the mechanistic bases for parasite pathogenicity and the Eastern oyster's response to it is limited. To better understand host responses to the parasite and identify molecular mechanisms underlying disease-resistance phenotypes, we experimentally challenged two families exhibiting divergent Dermo-resistance phenotypes with the parasite, generated global expression profiles using RNAseq and identified differentially expressed transcripts between control and challenged oysters from each family at multiple time points post-parasite injection. The susceptible and resistant families exhibited strikingly different transcriptomic responses to the parasite over a 28-day time period. The resistant family exhibited a strong, focused, early response to P. marinus infection, where many significantly upregulated transcripts were associated with the biological processes "regulation of proteolysis" and "oxidation-reduction process." P. marinus virulence factors are mainly comprised of proteases that facilitate parasite invasion and weaken host humoral defenses, thus host upregulation of transcripts associated with negative regulation of proteolysis is consistent with a Dermo-resistant phenotype. In contrast, the susceptible family mounted a very weak, disorganized, initial response to the parasite. Few transcripts were differentially expressed between control and injected oysters, and no functional enrichment was detected among them. At the final 28 d time point 2450 differentially expressed transcripts were identified and were associated with either "G-protein coupled receptor activity" (upregulated) or "microtubule-based process" (downregulated). A handful of protease inhibitors were differentially expressed between control and injected susceptible oysters, but this function was not enriched in the susceptible data set. The differential expression patterns observed in this study provide valuable insight into the functional basis of Dermo resistance and suggest that the timing of expression is just as important as the transcripts being expressed.

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.

Histone Methylation Participates in Gene Expression Control during the Early Development of the Pacific Oyster Crassostrea gigas

Histone methylation patterns are important epigenetic regulators of mammalian development, notably through stem cell identity maintenance by chromatin remodeling and transcriptional control of pluripotency genes. But, the implications of histone marks are poorly understood in distant groups outside vertebrates and ecdysozoan models. However, the development of the Pacific oyster Crassostrea gigas is under the strong epigenetic influence of DNA methylation, and Jumonji histone-demethylase orthologues are highly expressed during C. gigas early life. This suggests a physiological relevance of histone methylation regulation in oyster development, raising the question of functional conservation of this epigenetic pathway in lophotrochozoan. Quantification of histone methylation using fluorescent ELISAs during oyster early life indicated significant variations in monomethyl histone H3 lysine 4 (H3K4me), an overall decrease in H3K9 mono- and tri-methylations, and in H3K36 methylations, respectively, whereas no significant modification could be detected in H3K27 methylation. Early in vivo treatment with the JmjC-specific inhibitor Methylstat induced hypermethylation of all the examined histone H3 lysines and developmental alterations as revealed by scanning electronic microscopy. Using microarrays, we identified 376 genes that were differentially expressed under methylstat treatment, which expression patterns could discriminate between samples as indicated by principal component analysis. Furthermore, Gene Ontology revealed that these genes were related to processes potentially important for embryonic stages such as binding, cell differentiation and development. These results suggest an important physiological significance of histone methylation in the oyster embryonic and larval life, providing, to our knowledge, the first insights into epigenetic regulation by histone methylation in lophotrochozoan development.

Systematic identification and validation of the reference genes from 60 RNA-Seq libraries in the scallop Mizuhopecten yessoensis

Background

Reverse transcription quantitative PCR (RT-qPCR) is widely used for gene expression analysis in various organisms. Its accuracy largely relies on the stability of reference genes, making reference gene selection a vital step in RT-qPCR experiments. However, previous studies in mollusks only focused on the reference genes widely used in vertebrates.

Results

In this study, we conducted the transcriptome-wide identification of reference genes in the bivalve mollusk Mizuhopecten yessoensis based on 60 transcriptomes covering early development, adult tissues and gonadal development. A total of 964, 1210 and 2097 candidate reference genes were identified, respectively, resulting in a core set of 568 genes. Functional enrichment analysis showed that these genes are significantly overrepresented in Gene Ontology (GO) terms or Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways related to ribosomes, energy production, etc. Six genes (RS23, EF1A, NDUS4, SELR1, EIF3F, and OLA1) were selected from the candidate genes for RT-qPCR validation, together with 6 commonly used reference genes (ACT, CYTC, HEL, EF1B, GAPDH and RPL16). Stability analyses using geNorm, NormFinder and the comparative delta-Ct method revealed that the new candidate reference genes are more stable than the traditionally used genes, and ACT and CYTC are not recommended under either of the three circumstances. There was a significant correlation between the Ct of RT-qPCR and the log₂(TPM) of RNA-Seq data (Ct = − 0.94 log₂(TPM) + 29.67, R² = 0.73), making it easy to estimate the Ct values from transcriptome data prior to RT-qPCR experiments.

Conclusion

Our study represents the first transcriptome-wide identification of reference genes for early development, adult tissues, and gonadal development in the Yesso scallop and will benefit gene expression studies in other bivalve mollusks.

Electronic supplementary material

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

Molecular evolution and functional characterisation of insulin related peptides in molluscs: Contributions of Crassostrea gigas genomic and transcriptomic-wide screening

Insulin Related Peptides (IRPs) belong to the insulin superfamily and possess a typical structure with two chains, B and A, linked by disulphide bonds. As the sequence conservation is usually low between members, IRPs are classified according to the number and position of their disulphide bonds. In molluscan species, the first IRPs identified, named Molluscan Insulin-related Peptides (MIPs), exhibit four disulphide bonds. The genomic and transcriptomic data screening in the Pacific oyster Crassostrea gigas (Mollusc, Bivalvia) allowed us to identify six IRP sequences belonging to three structural groups. Cg-MIP1 to 4 have the typical structure of MIPs with four disulphide bonds. Cg-ILP has three disulphide bonds like vertebrate Insulin-Like Peptides (ILPs). The last one, Cg-MILP7 has a significant homology with Drosophila ILP7 (DILP7) associated with two additional cysteines allowing the formation of a fourth disulphide bond. The phylogenetic analysis points out that ILPs may be the most ancestral form. Moreover, it appears that ILP7 orthologs are probably anterior to lophotrochozoa and ecdysozoa segregation. In order to investigate the diversity of physiological functions of the oyster IRPs, we combine in silico expression data, qPCR measurements and in situ hybridization. The Cg-ilp transcript, mainly detected in the digestive gland and in the gonadal area, is potentially involved in the control of digestion and gametogenesis. The expression of Cg-mip4 is mainly associated with the larval development. The Cg-mip transcript shared by the Cg-MIP1, 2 and 3, is mainly expressed in visceral ganglia but its expression was also observed in the gonads of mature males. This pattern suggested the key roles of IRPs in the control of sexual reproduction in molluscan species.

De novo male gonad transcriptome draft for the marine mussel Perumytilus purpuratus with a focus on its reproductive-related proteins

Perumytilus purpuratus is a marine mussel considered a bioengineer species with a broad distribution in the Pacific and Atlantic coast of South America. Studies have shown two geographically and genetically differentiated subpopulations at molecular level and in sperm morphological traits. To open avenues for molecular research on P. purpuratus, a global de novo transcriptome from gonadal tissue of mature males was sequenced using the Illumina platform. From a total of 126.38 million reads, 37,765 transcripts were successfully annotated. BUSCO analysis determined a level of 89% completeness for the assembled transcriptome. The functional gene ontology (GO) annotation indicated that, in terms of abundance, the transcripts related with molecular function were the most represented, followed by those related with biological process and cellular components. Additionally, a subset of GO annotations generated using the "sperm" term resulted in a total of 1,294 sequences where the biological process category was the more represented, with transcripts strongly associated to sperm-processes required for fertilization, and with processes where the sperm-egg interaction could be implicated. Our work will contribute to the evolutionary understanding of the molecular mechanisms related to tissue-specific functions. This work reports the first male gonad transcriptome for the mussel P. purpuratus, generating a useful transcriptomic resource for this species and other closely related mytilids.

Morphological and molecular criteria allow the identification of putative germ stem cells in a lophotrochozoan, the Pacific oyster Crassostrea gigas

While our knowledge of bivalve gametogenesis recently progressed, data on early stages of gametogenesis remain to be developed, especially when dealing with germinal stem cells (GSC) and their niche in these organisms. Here, we wish to develop a strategy to identify putative GSC in Pacific oyster Crassostrea gigas based on morphological criteria combined with vasa marker expression. A histological quantitative approach, based on stereology, allowed us to identify two types of early germ cells in the germinal epithelium, one presenting round nuclei and the other irregular ones. Both early germ cell types present slightly condensed chromatin in nucleus, are vasa-positive and the Oyvlg (oyster vasa-like gene) expression in these cells is recorded throughout the whole gametogenesis process. The microenvironment of an early germ cell in oyster includes an associated somatic cell presenting an immunolabeling for BMP2/4 and a close myoid cell. In agreement with the GSC characteristics in other species, we postulate that putative germ stem cells in C. gigas correspond to the early germ cell type with irregular nucleus shape; those early germ cells with a round nucleus may consist in progenitors.

fam20C participates in the shell formation in the pearl oyster, Pinctada fucata

Kinase-family with sequence similarity 20, member C (Fam20C) is a protein kinase, which can phosphorylate biomineralization related proteins in vertebrate animals. However, the function of Fam20C in invertebrate animals especially the role in biomineralization is still unknown. Herein, we cloned the cDNA of fam20C from the pearl oyster, Pinctada fucata. It is showed that the expression of fam20C in the mantle edge was much higher than other tissues. In situ hybridization showed that fam20C was expressed mostly in the outer epithelial cells of the middle fold, indicating it may play important roles in the shell formation. Besides, fam20C expression increased greatly in the D-shape stage of pearl oyster development, when the shell was first formed. During the shell repair process, the expression level of fam20C increased 1.5 times at 6 h after shell notching. Knockdown of fam20C in vivo by RNA interference resulted in abnormally stacking of calcium carbonate crystals at the edges of nacre tablets, showing direct evidence that fam20C participates in the shell formation. This study provides an insight into the role of kinase protein in the shell formation in mollusk and broaden our understanding of biomineralization mechanism.

Transcriptional response after exposure to domoic acid-producing Pseudo-nitzschia in the digestive gland of the mussel Mytilus galloprovincialis

Bivalve molluscs are filter feeding species that can accumulate biotoxins in their body tissues during harmful algal blooms. Amnesic Shellfish Poisoning (ASP) is caused by species of the diatom genus Pseudo-nitzschia, which produces the toxin domoic acid. The Mytilus galloprovincialis digestive gland transcriptome was de novo assembled based on the sequencing of 12 cDNA libraries, six obtained from control mussels and six from mussels naturally exposed to domoic acid-producing diatom Pseudo-nitzschia australis. After de novo assembly 94,727 transcripts were obtained, with an average length of 1015 bp and a N50 length of 761 bp. The assembled transcripts were clustered (homology > 90%) into 69,294 unigenes. Differential gene expression analysis was performed (DESeq2 algorithm) in the digestive gland following exposure to the toxic algae. A total of 1158 differentially expressed unigenes (absolute fold change > 1.5 and p-value < 0.05) were detected: 686 up-regulated and 472 down-regulated. Several membrane transporters belonging to the family of the SLC (solute carriers) were over-expressed in exposed mussels. Functional enrichment was performed using Pfam annotations obtained from the genes differentially expressed, 37 Pfam families were found to be significantly (FDR adjusted p-value < 0.1) enriched. Some of these families (sulfotransferases, aldo/keto reductases, carboxylesterases, C1q domain and fibrinogen C-terminal globular domain) could be putatively involved in detoxification processes, in the response against of the oxidative stress and in immunological processes. Protein network analysis with STRING algorithm found alteration of the Notch signaling pathway under the action of domoic acid-producing Pseudo-nitzschia. In conclusion, this study provides a high quality reference transcriptome of M. galloprovincialis digestive gland and identifies potential genes involved in the response to domoic acid.

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.

Selection of housekeeping genes as internal controls for quantitative RT-PCR analysis of the veined rapa whelk (Rapana venosa)

Background

The veined rapa whelk Rapana venosa is an important commercial shellfish in China and quantitative real-time PCR (qRT-PCR) has become the standard method to study gene expression in R. venosa. For accurate and reliable gene expression results, qRT-PCR assays require housekeeping genes as internal controls, which display highly uniform expression in different tissues or stages of development. However, to date no studies have validated housekeeping genes in R. venosa for use as internal controls for qRT-PCR.

Methods

In this study, we selected the following 13 candidate genes for suitability as internal controls: elongation factor-1α (EF-1α), α-actin (ACT), cytochrome c oxidase subunit 1 (COX1), nicotinamide adenine dinucleotide dehydrogenase (ubiquinone) 1α subcomplex subunit 7 (NDUFA7), 60S ribosomal protein L5 (RL5), 60S ribosomal protein L28 (RL28), glyceraldehyde 3-phosphate dehydrogenase (GAPDH), β-tubulin (TUBB), 40S ribosomal protein S25 (RS25), 40S ribosomal protein S8 (RS8), ubiquitin-conjugating enzyme E2 (UBE2), histone H3 (HH3), and peptidyl-prolyl cis-trans isomerase A (PPIA). We measured the expression levels of these 13 candidate internal controls in eight different tissues and twelve larvae developmental stages by qRT-PCR. Further analysis of the expression stability of the tested genes was performed using GeNorm and RefFinder algorithms.

Results

Of the 13 candidate genes tested, we found that EF-1α was the most stable internal control gene in almost all adult tissue samples investigated with RL5 and RL28 as secondary choices. For the normalization of a single specific tissue, we suggested that EF-1α and NDUFA7 are the best combination in gonad, as well as COX1 and RL28 for intestine, EF-1α and RL5 for kidney, EF-1α and COX1 for gill, EF-1α and RL28 for Leiblein and mantle, EF-1α, RL5, and NDUFA7 for liver, GAPDH, PPIA, and RL28 for hemocyte. From a developmental perspective, we found that RL28 was the most stable gene in all developmental stages measured, and COX1 and RL5 were appropriate secondary choices. For the specific developmental stage, we recommended the following combination for normalization, PPIA, RS25, and RL28 for stage 1, RL5 and RL28 for stage 2 and 5, RL28 and NDUFA7 for stage 3, and PPIA and TUBB for stage 4.

Discussion

Our results are instrumental for the selection of appropriately validated housekeeping genes for use as internal controls for gene expression studies in adult tissues or larval development of R. venosa in the future.

Selection and evaluation of new reference genes for RT-qPCR analysis in Epinephelus akaara based on transcriptome data

Groupers are an economically important fish species in world fishery markets. Because many studies using RT-qPCR have addressed gene expression in groupers, appropriate reference genes are required to obtain reliable and accurate results. In this study, the most suitable reference genes were identified from eleven candidate genes of one of the most valuable species, Epinephelus akaara, in a range of different experimental conditions. Using the software packages geNorm, NormFinder, BestKeeper and refFinder, three traditionally used reference genes, β-actin (β-ACT), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), and beta-2-microglobulin (B2M), were identified as not suitable for E. akaara gene expression studies, whereas two newly identified reference genes, conserved oligomeric Golgi complex subunit 5 (Cog5) and brefeldin a-inhibited guanine nucleotide-exchange protein 1 (ARFGEF1), could be universally applied under all the tested conditions. These data provide the foundation for more precise results in RT-qPCR studies of gene expression in E. akaara and other Epinephelus species.

Construction of an Ostrea edulis database from genomic and expressed sequence tags (ESTs) obtained from Bonamia ostreae infected haemocytes: Development of an immune-enriched oligo-microarray

The flat oyster, Ostrea edulis, is one of the main farmed oysters, not only in Europe but also in the United States and Canada. Bonamiosis due to the parasite Bonamia ostreae has been associated with high mortality episodes in this species. This parasite is an intracellular protozoan that infects haemocytes, the main cells involved in oyster defence. Due to the economical and ecological importance of flat oyster, genomic data are badly needed for genetic improvement of the species, but they are still very scarce. The objective of this study is to develop a sequence database, OedulisDB, with new genomic and transcriptomic resources, providing new data and convenient tools to improve our knowledge of the oyster's immune mechanisms. Transcriptomic and genomic sequences were obtained using 454 pyrosequencing and compiled into an O. edulis database, OedulisDB, consisting of two sets of 10,318 and 7159 unique sequences that represent the oyster's genome (WG) and de novo haemocyte transcriptome (HT), respectively. The flat oyster transcriptome was obtained from two strains (naïve and tolerant) challenged with B. ostreae, and from their corresponding non-challenged controls. Approximately 78.5% of 5619 HT unique sequences were successfully annotated by Blast search using public databases. A total of 984 sequences were identified as being related to immune response and several key immune genes were identified for the first time in flat oyster. Additionally, transcriptome information was used to design and validate the first oligo-microarray in flat oyster enriched with immune sequences from haemocytes. Our transcriptomic and genomic sequencing and subsequent annotation have largely increased the scarce resources available for this economically important species and have enabled us to develop an OedulisDB database and accompanying tools for gene expression analysis. This study represents the first attempt to characterize in depth the O. edulis haemocyte transcriptome in response to B. ostreae through massively sequencing and has aided to improve our knowledge of the immune mechanisms of flat oyster. The validated oligo-microarray and the establishment of a reference transcriptome will be useful for large-scale gene expression studies in this species.

Assessment of housekeeping genes as internal references in quantitative expression analysis during early development of oyster

The early development of mollusks exhibits important characteristics from the developmental and evolutionary perspective. With the increasing number of genome-wide studies, accurate analyses of quantitative gene expression during development are impeded by the lack of validated reference genes. To improve the situation, in this study, we analyzed the expression stability of seven candidate housekeeping genes during early development of the Pacific oyster Crassostrea gigas: actin, glyceraldehyde-3-phosphate dehydrogenase (gapdh), α subunit of elongation factor 1 (elf1α), adp-ribosylation factor 1 (arf1), heterogeneous nuclear ribonucleoprotein q, ubiquitin-conjugating enzyme e2d2 and ribosomal protein s18. We focused on 11 stages from oocyte to D-veliger, which include crucial developmental processes such as axis determination, gastrulation and shell formation. Gene expression stabilities were assessed with the three commonly used programs geNorm, NormFinder and BestKeeper. Although the results obtained with the three programs varied to some extent, in general, arf1, elf1α and gapdh were highly ranked and actin was poorly ranked. This analysis also indicated that multiple genes should be used for normalization, and we concluded that arf1-elf1α-gapdh should be used as internal references. The findings of this study will help researchers to obtain accurate results in future quantitative gene expression analysis of development in bivalve mollusks.

Quantitative analysis of oyster larval proteome provides new insights into the effects of multiple climate change stressors

The metamorphosis of planktonic larvae of the Pacific oyster (Crassostrea gigas) underpins their complex life-history strategy by switching on the molecular machinery required for sessile life and building calcite shells. Metamorphosis becomes a survival bottleneck, which will be pressured by different anthropogenically induced climate change-related variables. Therefore, it is important to understand how metamorphosing larvae interact with emerging climate change stressors. To predict how larvae might be affected in a future ocean, we examined changes in the proteome of metamorphosing larvae under multiple stressors: decreased pH (pH 7.4), increased temperature (30 °C), and reduced salinity (15 psu). Quantitative protein expression profiling using iTRAQ-LC-MS/MS identified more than 1300 proteins. Decreased pH had a negative effect on metamorphosis by down-regulating several proteins involved in energy production, metabolism, and protein synthesis. However, warming switched on these down-regulated pathways at pH 7.4. Under multiple stressors, cell signaling, energy production, growth, and developmental pathways were up-regulated, although metamorphosis was still reduced. Despite the lack of lethal effects, significant physiological responses to both individual and interacting climate change related stressors were observed at proteome level. The metamorphosing larvae of the C. gigas population in the Yellow Sea appear to have adequate phenotypic plasticity at the proteome level to survive in future coastal oceans, but with developmental and physiological costs.

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.

Additive transcriptomic variation associated with reproductive traits suggest local adaptation in a recently settled population of the Pacific oyster, Crassostrea gigas

Background

Originating from Northeast Asia, the Pacific oyster Crassostrea gigas has been introduced into a large number of countries for aquaculture purpose. Following introduction, the Pacific oyster has turned into an invasive species in an increasing number of coastal areas, notably recently in Northern Europe.

Methods

To explore potential adaptation of reproductive traits in populations with different histories, we set up a common garden experiment based on the comparison of progenies from two populations of Pacific oyster sampled in France and Denmark and their hybrids. Sex ratio, condition index and microarray gene expression in gonads, were analyzed in each progeny (n = 60).

Results

A female-biased sex-ratio and a higher condition index were observed in the Danish progeny, possibly reflecting an evolutionary reproductive strategy to increase the potential success of natural recruitment in recently settled population. Using multifarious statistical approaches and accounting for sex differences we identified several transcripts differentially expressed between the Danish and French progenies, for which additive genetic basis is suspected (showing intermediate expression levels in hybrids, and therefore additivity). Candidate transcripts included mRNA coding for sperm quality and insulin metabolism, known to be implicated in coordinated control and success of reproduction.

Conclusions

Observed differences suggest that adaptation of invasive populations might have occurred during expansion acting on reproductive traits, and in particular on a female-biased sex-ratio, gamete quality and fertility.

Electronic supplementary material

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

A Label-Free Proteomic Analysis on Competent Larvae and Juveniles of the Pacific Oyster Crassostrea gigas

Current understandings on the molecular mechanisms underlying bivalve metamorphosis are still fragmentary, and a comprehensive description is required. In this study, using a large-scale label-free proteomic approach, we described and compared the proteomes of competent larvae (CL) and juveniles (JU) of the Pacific oyster, Crassostrea gigas. A total of 788 proteins were identified: 392 in the CL proteome and 636 in the JU proteome. Gene Ontology analysis of the proteome from each sample revealed active metabolic processes in both stages. Further quantitative analyses revealed 117 proteins that were differentially expressed between the two samples. These proteins were divided into eight groups: cytoskeleton and cell adhesion, protein synthesis and degradation, immunity and stress response, development of particular tissues, signal regulation, metabolism and energy supply, transport, and other proteins. A certification experiment using real-time PCR assay confirmed 20 of 30 examined genes exhibited the same trends at the mRNA and protein levels. The differentially expressed proteins may play roles in tissue remodeling, signal transduction, and organ development during and after metamorphosis. Novel roles were proposed for some differentially expressed proteins, such as chymotrypsin. The results of this work provide an overview of metamorphosis and post-metamorphosis development of C. gigas at the protein level. Future studies on the functions of the differentially expressed proteins will help to obtain a more in-depth understanding of bivalve metamorphosis.

Applications of biological tools or biomarkers in aquatic biota: A case study of the Tamar estuary, South West England

Biological systems are the ultimate recipients of pollutant-induced damage. Consequently, our traditional reliance on analytical tools is not enough to assess ecosystem health. Biological responses or biomarkers are therefore also considered to be important tools for environmental hazard and risk assessments. Due to historical mining, other anthropogenic activities, and its conservational importance (e.g. NATURA sites, SACs), the Tamar estuary in South West England is an ideal environment in which to examine applications of such biological tools. This review presents a thorough and critical evaluation of the different biological tools used in the Tamar estuary thus far, while also discussing future perspectives for biomarker studies from a global perspective. In particular, we focus on the challenges which hinder applications of biological tools from being more readily incorporated into regulatory frameworks, with the aim of enabling both policymakers and primary stakeholders to maximise the environmental relevance and regulatory usefulness of such tools.

The use of -omic tools in the study of disease processes in marine bivalve mollusks

Our understanding of disease processes and host-pathogen interactions in model species has benefited greatly from the application of medium and high-throughput genomic, metagenomic, epigenomic, transcriptomic, and proteomic analyses. The rate at which new, low-cost, high-throughput -omic technologies are being developed has also led to an expansion in the number of studies aimed at gaining a better understanding of disease processes in bivalves. This review provides a catalogue of the genetic and -omic tools available for bivalve species and examples of how -omics has contributed to the advancement of marine bivalve disease research, with a special focus in the areas of immunity, bivalve-pathogen interactions, mechanisms of disease resistance and pathogen virulence, and disease diagnosis. The analysis of bivalve genomes and transcriptomes has revealed that many immune and stress-related gene families are expanded in the bivalve taxa examined thus far. In addition, the analysis of proteomes confirms that responses to infection are influenced by epigenetic, post-transcriptional, and post-translational modifications. The few studies performed in bivalves show that epigenetic modifications are non-random, suggesting a role for epigenetics in regulating the interactions between bivalves and their environments. Despite the progress -omic tools have enabled in the field of marine bivalve disease processes, there is much more work to be done. To date, only three bivalve genomes have been sequenced completely, with assembly status at different levels of completion. Transcriptome datasets are relatively easy and inexpensive to generate, but their interpretation will benefit greatly from high quality genome assemblies and improved data analysis pipelines. Finally, metagenomic, epigenomic, proteomic, and metabolomic studies focused on bivalve disease processes are currently limited but their expansion should be facilitated as more transcriptome datasets and complete genome sequences become available for marine bivalve species.

An immune-enriched oligo-microarray analysis of gene expression in Manila clam (Venerupis philippinarum) haemocytes after a Perkinsus olseni challenge

Parasites of the genus Perkinsus cause high mortality and economic losses in bivalves commonly produced in global aquaculture. Although the immune responses of oysters and clams naturally infected with Perkinsus marinus or Perkinsus olseni have been extensively studied, there is not much information on host response at the early stages of infection. In this study, we analysed how P. olseni influences the gene expression profiles of haemocytes from the Manila clam (Venerupis philippinarum) using temporal experimental infections and an immune-enriched microarray. We identified an early phase of infection that was characterised by no mortality and by the increased expression of genes associated with pathogen recognition, production of nitrogen radicals and antimicrobial activity. Cellular processes such as inhibition of serine proteases and proliferation were also involved in this early response. This phase was followed by an intermediate stage, when the pathogen was most likely multiplying and infecting new areas of the body, and animals began to die. In this stage, many genes related to cell movement were over-expressed. Thirty days after infection metabolic pathway genes were the most affected. Apoptosis appears to be important during pathogenesis. Our results provide novel observations of the broader innate immune response triggered by P. olseni at different infection stages.

A family of variable immunoglobulin and lectin domain containing molecules in the snail Biomphalaria glabrata

Technical limitations have hindered comprehensive studies of highly variable immune response molecules that are thought to have evolved due to pathogen-mediated selection such as fibrinogen-related proteins (FREPs) from Biomphalaria glabrata. FREPs combine upstream immunoglobulin superfamily (IgSF) domains with a C-terminal fibrinogen-related domain (FreD) and participate in reactions against trematode parasites. From RNAseq data we assembled a de novo reference transcriptome of B. glabrata to investigate the diversity of FREP transcripts. This study increased over two fold the number of bonafide FREP subfamilies and revealed important sequence diversity within FREP12 subfamily. We also report the discovery of related molecules that feature one or two IgSF domains associated with different C-terminal lectin domains, named C-type lectin-related proteins (CREPs) and Galectin-related protein (GREP). Together, the highly similar FREPs, CREPs and GREP were designated VIgL (Variable Immunoglobulin and Lectin domain containing molecules).

Transcriptomic profiling of gametogenesis in triploid Pacific Oysters Crassostrea gigas: towards an understanding of partial sterility associated with triploidy

BACKGROUND

Triploidy can occur in many animal species but is often lethal. Among invertebrates, amphibians and fishes, triploids are viable although often sterile or infertile. Most triploids of the Pacific oyster Crassostrea gigas are almost sterile (named "3nβ") yet a low but significant proportion show an advanced gametogenesis (named "3nα"). These oysters thus constitute an interesting model to study the effect of triploidy on germ cell development. We used microarrays to compare the gonad transcriptomes of diploid 2n and the abovementioned triploid 3nβ and 3nα male and female oysters throughout gametogenesis.

RESULTS

All triploids displayed an upregulation of genes related to DNA repair and apoptosis and a downregulation of genes associated with cell division. The comparison of 3nα and 3nβ transcriptomes with 2n revealed the likely involvement of a cell cycle checkpoint during mitosis in the successful but delayed development of gonads in 3nα individuals. In contrast, a disruption of sex differentiation mechanisms may explain the sterility of 3nβ individuals with 3nβ females expressing male-specific genes and 3nβ males expressing female-specific genes.

CONCLUSIONS

The disruption of sex differentiation and mitosis may be responsible for the impaired gametogenesis of triploid Pacific oysters. The function of the numerous candidate genes identified in our study should now be studied in detail in order to elucidate their role in sex determination, mitosis/meiosis control, pachytene cell cycle checkpoint, and the control of DNA repair/apoptosis.

Selection of reference genes as internal controls for gene expression in tissues of red abalone Haliotis rufescens (Mollusca, Vetigastropoda; Swainson, 1822)

The red abalone Haliotis rufescens is one of the most important species for aquaculture in Baja California, México, and despite this, few gene expression studies have been done in tissues such as gill, head and gonad. For this purpose, reverse transcription and quantitative real time PCR (RT-qPCR) is a powerful tool for gene expression evaluation. For a reliable analysis, however, it is necessary to select and validate housekeeping genes that allow proper transcription quantification. Stability of nine housekeeping genes (ACTB, BGLU, TUBB, CY, GAPDH, HPRTI, RPL5, SDHA and UBC) was evaluated in different tissues of red abalone (gill, head and gonad/digestive gland). Four-fold serial dilutions of cDNA (from 25 ngμL(-1) to 0.39 ngμL(-1)) were used to prepare the standard curve, and it showed gene efficiencies between 0.95 and 0.99, with R(2)=0.99. geNorm and NormFinder analysis showed that RPL5 and CY were the most stable genes considering all tissues, whereas in gill HPRTI and BGLU were most stable. In gonad/digestive gland, RPL5 and TUBB were the most stable genes with geNorm, while SDHA and HPRTI were the best using NormFinder. Similarly, in head the best genes were RPL5 and UBC with geNorm, and GAPDH and CY with NormFinder. The technical variability analysis with RPL5 and abalone gonad/digestive gland tissue indicated a high repeatability with a variation coefficient within groups ≤ 0.56% and between groups ≤ 1.89%. These results will help us for further research in reproduction, thermoregulation and endocrinology in red abalone.

Dual transcriptomics of virus-host interactions: comparing two Pacific oyster families presenting contrasted susceptibility to ostreid herpesvirus 1

BACKGROUND

Massive mortality outbreaks affecting Pacific oyster (Crassostrea gigas) spat in various countries have been associated with the detection of a herpesvirus called ostreid herpesvirus type 1 (OsHV-1). However, few studies have been performed to understand and follow viral gene expression, as it has been done in vertebrate herpesviruses. In this work, experimental infection trials of C. gigas spat with OsHV-1 were conducted in order to test the susceptibility of several bi-parental oyster families to this virus and to analyze host-pathogen interactions using in vivo transcriptomic approaches.

RESULTS

The divergent response of these oyster families in terms of mortality confirmed that susceptibility to OsHV-1 infection has a significant genetic component. Two families with contrasted survival rates were selected. A total of 39 viral genes and five host genes were monitored by real-time PCR. Initial results provided information on (i) the virus cycle of OsHV-1 based on the kinetics of viral DNA replication and transcription and (ii) host defense mechanisms against the virus.

CONCLUSIONS

In the two selected families, the detected amounts of viral DNA and RNA were significantly different. This result suggests that Pacific oysters are genetically diverse in terms of their susceptibility to OsHV-1 infection. This contrasted susceptibility was associated with dissimilar host gene expression profiles. Moreover, the present study showed a positive correlation between viral DNA amounts and the level of expression of selected oyster genes.

Gonad transcriptome analysis of pearl oyster Pinctada margaritifera: identification of potential sex differentiation and sex determining genes

BACKGROUND

Black pearl farming is based on culture of the blacklip pearl oyster Pinctada margaritifera (Mollusca, lophotrochozoa), a protandrous hermaphrodite species. At first maturation, all individuals are males. The female sex appears progressively from two years old, which represents a limitation for broodstock conditioning for aquaculture production. In marine mollusks displaying hermaphroditic features, data on sexual determinism and differentiation, including the molecular sex determining cascade, are scarce. To increase genomic resources and identify the molecular mechanisms whereby gene expression may act in the sexual dimorphism of P. margaritifera, we performed gonad transcriptome analysis.

RESULTS

The gonad transcriptome of P. margaritifera was sequenced from several gonadic samples of males and females at different development stages, using a Next-Generation-Sequencing method and RNAseq technology. After Illumina sequencing, assembly and annotation, we obtained 70,147 contigs of which 62.2% shared homologies with existing protein sequences, and 9% showed functional annotation with Gene Ontology terms. Differential expression analysis identified 1,993 differentially expressed contigs between the different categories of gonads. Clustering methods of samples revealed that the sex explained most of the variation in gonad gene expression. K-means clustering of differentially expressed contigs showed 815 and 574 contigs were more expressed in male and female gonads, respectively. The analysis of these contigs revealed the presence of known specific genes coding for proteins involved in sex determinism and/or differentiation, such as dmrt and fem-1 like for males, or foxl2 and vitellogenin for females. The specific gene expression profiles of pmarg-fem1-like, pmarg-dmrt and pmarg-foxl2 in different reproductive stages (undetermined, sexual inversion and regression) suggest that these three genes are potentially involved in the sperm-oocyte switch in P. margaritifera.

CONCLUSIONS

The study provides a new transcriptomic tool to study reproduction in hermaphroditic marine mollusks. It identifies sex differentiation and potential sex determining genes in P. margaritifera, a protandrous hermaphrodite species.

Gene expression profile analysis of Manila clam (Ruditapes philippinarum) hemocytes after a Vibrio alginolyticus challenge using an immune-enriched oligo-microarray

Background

The Manila clam (Ruditapes philippinarum) is a cultured bivalve with worldwide commercial importance, and diseases cause high economic losses. For this reason, interest in the immune genes in this species has recently increased. The present work describes the construction of the first R. philippinarum microarray containing immune-related hemocyte sequences and its application to study the gene transcription profiles of hemocytes from clams infected with V. alginolyticus through a time course.

Results

The complete set of sequences from R. philippinarum available in the public databases and the hemocyte sequences enriched in immune transcripts were assembled successfully. A total of 12,156 annotated sequences were used to construct the 8 ×15 k oligo-microarray. The microarray experiments yielded a total of 579 differentially expressed transcripts. Using the gene expression results, the associated Gene Ontology terms and the enrichment analysis, we found different response mechanisms throughout the experiment. Genes related to signaling, transcription and apoptosis, such as IL-17D, NF-κB or calmodulin, were typically expressed as early as 3 hours post-challenge (hpc), while characteristic immune genes, such as PGRPs, FREPs and defense proteins appeared later at 8 hpc. This immune-triggering response could have affected a high number of processes that seemed to be activated 24 hpc to overcome the Vibrio challenge, including the expression of many cytoskeleton molecules, which is indicative of the active movement of hemocytes. In fact functional studies showed an increment in apoptosis, necrosis or cell migration after the infection. Finally, 72 hpc, activity returned to normal levels, and more than 50% of the genes were downregulated in a negative feedback of all of the previously active processes.

Conclusions

Using a new version of the R. philippinarum oligo-microarray, a putative timing for the response against a Vibrio infection was established. The key point to overcome the challenge seemed to be 8 hours after the challenge, when we detected immune functions that could lead to the destruction of the pathogen and the activation of a wide variety of processes related to homeostasis and defense. These results highlight the importance of a fast response in bivalves and the effectiveness of their innate immune system.

Reproduction-related genes in the pearl oyster genome

Molluscan reproduction has been a target of biological research because of the various reproductive strategies that have evolved in this phylum. It has also been studied for the development of fisheries technologies, particularly aquaculture. Although fundamental processes of reproduction in other phyla, such as vertebrates and arthropods, have been well studied, information on the molecular mechanisms of molluscan reproduction remains limited. The recently released draft genome of the pearl oyster Pinctada fucata provides a novel and powerful platform for obtaining structural information on the genes and proteins involved in bivalve reproduction. In the present study, we analyzed the pearl oyster draft genome to screen reproduction-related genes. Analysis was mainly conducted for genes reported from other molluscs for encoding orthologs of reproduction-related proteins in other phyla. The gene search in the P. fucata gene models (version 1.1) and genome assembly (version 1.0) were performed using Genome Browser and BLAST software. The obtained gene models were then BLASTP searched against a public database to confirm the best-hit sequences. As a result, more than 40 gene models were identified with high accuracy to encode reproduction-related genes reported for P. fucata and other molluscs. These include vasa, nanos, doublesex- and mab-3-related transcription factor, 5-hydroxytryptamine (5-HT) receptors, vitellogenin, estrogen receptor, and others. The set of reproduction-related genes of P. fucata identified in the present study constitute a new tool for research on bivalve reproduction at the molecular level.

Bivalve omics: state of the art and potential applications for the biomonitoring of harmful marine compounds

The extraordinary progress experienced by sequencing technologies and bioinformatics has made the development of omic studies virtually ubiquitous in all fields of life sciences nowadays. However, scientific attention has been quite unevenly distributed throughout the different branches of the tree of life, leaving molluscs, one of the most diverse animal groups, relatively unexplored and without representation within the narrow collection of well established model organisms. Within this Phylum, bivalve molluscs play a fundamental role in the functioning of the marine ecosystem, constitute very valuable commercial resources in aquaculture, and have been widely used as sentinel organisms in the biomonitoring of marine pollution. Yet, it has only been very recently that this complex group of organisms became a preferential subject for omic studies, posing new challenges for their integrative characterization. The present contribution aims to give a detailed insight into the state of the art of the omic studies and functional information analysis of bivalve molluscs, providing a timely perspective on the available data resources and on the current and prospective applications for the biomonitoring of harmful marine compounds.

Physiological change under OsHV-1 contamination in Pacific oyster Crassostrea gigas through massive mortality events on fields

BACKGROUND

Massive mortalities have been observed in France since 2008 on spat and juvenile Pacific oysters, Crassostrea gigas. A herpes virus called OsHV-1, easily detectable by PCR, has been implicated in the mortalities as demonstrated by the results of numerous field studies linking mortality with OsHV-1 prevalence. Moreover, experimental infections using viral particles have documented the pathogenicity of OsHV-1 but the physiological responses of host to pathogen are not well known.

RESULTS

The aim of this study was to understand mechanisms brought into play against the virus during infection in the field. A microarray assay has been developed for a major part of the oyster genome and used for studying the host transcriptome across mortality on field. Spat with and without detectable OsHV-1 infection presenting or not mortality respectively were compared by microarray during mortality episodes. In this study, a number of genes are regulated in the response to pathogen infection on field and seems to argue to an implication of the virus in the observed mortality. The result allowed establishment of a hypothetic scheme of the host cell's infection by, and response to, the pathogen.

CONCLUSIONS

This response shows a "sensu stricto" innate immunity through genic regulation of the virus OsHV-1 life cycle, but also others biological processes resulting to complex interactions between host and pathogens in general.

Shotgun proteomics as a viable approach for biological discovery in the Pacific oyster

Shotgun proteomics offers an efficient means to characterize proteins in a complex mixture, particularly when sufficient genomic resources are available. In order to assess the practical application of shotgun proteomics in the Pacific oyster (Crassostrea gigas), liquid chromatography coupled with tandem mass spectrometry was used to characterize the gill proteome. Using information from the recently published Pacific oyster genome, 1043 proteins were identified. Biological samples (n = 4) and corresponding technical replicates (three) were similar in both specific proteins identified and expression, as determined by normalized spectral abundance factor. A majority of the proteins identified (703) were present in all biological samples. Functional analysis of the protein repertoire illustrates that these proteins represent a wide range of biological processes, supporting the dynamic function of the gill. These insights are important for understanding environmental influences on the oyster, because the gill tissue acts as the interface between the oyster and its environment. In silico analysis indicated that this sequencing effort identified a large proportion of the complete gill proteome. Together, these data demonstrate that shotgun sequencing is a viable approach for biological discovery and will play an important role in future studies of oyster physiology.

Rapid mitochondrial adjustments in response to short-term hypoxia and re-oxygenation in the Pacific oyster Crassostrea gigas

As oxygen concentrations in marine coastal habitats can fluctuate rapidly and drastically, sessile marine organisms like the oyster Crassostrea gigas can experience marked and rapid oxygen variations. In this study, we investigated the responses of oyster gill mitochondria to short-term hypoxia (3 h and 12 h, at 1.7 mg O2 L-1) and subsequent re-oxygenation. Mitochondrial respiratory rates (state 3 and 4 stimulated by glutamate), phosphorylation efficiency (RCR ratio and ADP/O) were measured. Cytochrome c oxidase activity (CCO) and cytochrome concentrations (a, b, c1, c) were measured to investigate the rearrangements of respiratory chain subunits. The potential implication of an alternative oxidase (AOX) was investigated using an inhibitor of the respiratory chain (antimycin A) and through gene expression analysis in gills and digestive gland. Results indicate a down regulation of mitochondrial capacity, with 60% inhibition of respiratory rates after 12 h of hypoxia. RCR ratio remained stable, while ADP/O increased after 12 h of hypoxia and 1 h of re-oxygenation, suggesting increased phosphorylation efficiency. CCO showed a fast and remarkable increase of its catalytic activity only after 3 h of hypoxia. AOX mRNA levels showed similar patterns in gill and digestive gland, and were up regulated after 12 and 24 h of hypoxia and during re-oxygenation. Results suggest a set of controls in regulating mitochondrial functions in response to oxygen fluctuations and demonstrate the fast and extreme plasticity of oyster mitochondria in response to oxygen variations.

Gametogenesis in the Pacific oyster Crassostrea gigas: a microarrays-based analysis identifies sex and stage specific genes

Background

The Pacific oyster Crassostrea gigas (Mollusca, Lophotrochozoa) is an alternative and irregular protandrous hermaphrodite: most individuals mature first as males and then change sex several times. Little is known about genetic and phenotypic basis of sex differentiation in oysters, and little more about the molecular pathways regulating reproduction. We have recently developed and validated a microarray containing 31,918 oligomers (Dheilly et al., 2011) representing the oyster transcriptome. The application of this microarray to the study of mollusk gametogenesis should provide a better understanding of the key factors involved in sex differentiation and the regulation of oyster reproduction.

Methodology/Principal Findings

Gene expression was studied in gonads of oysters cultured over a yearly reproductive cycle. Principal component analysis and hierarchical clustering showed a significant divergence in gene expression patterns of males and females coinciding with the start of gonial mitosis. ANOVA analysis of the data revealed 2,482 genes differentially expressed during the course of males and/or females gametogenesis. The expression of 434 genes could be localized in either germ cells or somatic cells of the gonad by comparing the transcriptome of female gonads to the transcriptome of stripped oocytes and somatic tissues. Analysis of the annotated genes revealed conserved molecular mechanisms between mollusks and mammals: genes involved in chromatin condensation, DNA replication and repair, mitosis and meiosis regulation, transcription, translation and apoptosis were expressed in both male and female gonads. Most interestingly, early expressed male-specific genes included bindin and a dpy-30 homolog and female-specific genes included foxL2, nanos homolog 3, a pancreatic lipase related protein, cd63 and vitellogenin. Further functional analyses are now required in order to investigate their role in sex differentiation in oysters.

Conclusions/Significance

This study allowed us to identify potential markers of early sex differentiation in the oyster C. gigas, an alternative hermaphrodite mollusk. We also provided new highly valuable information on genes specifically expressed by mature spermatozoids and mature oocytes.

Identification of a tubulin-α gene specifically expressed in testis and adductor muscle during stable reference gene selection in the hermaphrodite gonad of the lion's paw scallop Nodipecten subnodosus

For non-model species, as many used for aquaculture, with minimal or no genomic information, relative quantification of gene expression studies requires preliminary research including the isolation of potential reference genes and the identification of those stably expressed under the biological conditions of interest. Here we report on the isolation of five partial gene sequences from gonad tissue cDNA in the functional hermaphrodite scallop Nodipecten subnodosus to be evaluated as reference genes: 18S-rRNA, riboprotein l8 (rp-l8), actin-β (act-β), elongation factor 1α (ef-1α) and alpha-tubulin-α (tub-α). We found that 18S-rRNA was stably expressed independently of the priming method used to reverse transcribe RNA to cDNA, oligo-dT or random hexamer. Stability analysis for the five putative reference genes with geNorm and NormFinder indicated that 18S together with rp-l8 were the most stable genes for normalization of gene expression during gonad development in both, male and female sexual regions of the hermaphrodite N. subnodosus. The least stable gene was tub-α, showing a biased expression profile between sexual regions of the gonad, therefore this gene was analyzed thereafter as a target gene together with vitellogenin (vit) and a DEAD-box RNA helicase (dbx) gene. Relative expression, estimated by normalization with the combination of 18S and rp-l8 as reference genes, indicated that as gonad development advanced two of the target genes were up-regulated, tub-α in the male region and vit in the female region. Whereas an increased expression was expected during development for vit for its known role in vitellogenesis, the increased expression of tub-α in the male sexual region was unexpected, and pointed toward this gene being a testis-specific α-tubulin isotype. Further analyses of gene expression among tissues indicated that tub-α is specifically and highly expressed in the male gonad, although expression in adductor muscle was also observed at significantly lower levels. The existence of testis specific α- and β-tubulins has been previously reported in other taxa, relating their function to sperm axoneme formation. Tissue-specific tubulin genes, particularly their promoters, have recently found an application as native promoters for transgene tissue-specific expression in research and reproductive control of insect plagues. The third target gene, a putative member of the DEAD-box RNA helicase family (dbx), showed no changes in expression during gonad development or between sexual regions, therefore it was chosen to discuss the different statistical inferences resulting from the arbitrary use of 'randomly chosen' reference genes when normalizing gene expression.