Transcriptome analysis reveals a rich gene set related to innate immunity in the Eastern oyster (Crassostrea virginica)

Proteomic profile of Ostrea edulis haemolymph in response to bonamiosis and identification of candidate proteins as resistance markers

European flat oyster Ostrea edulis populations have suffered extensive mortalities caused by bonamiosis. The protozoan parasite Bonamia ostreae is largely responsible for this disease in Europe, while its congener B. exitiosa has been detected more recently in various European countries. Both of these intracellular parasites are able to survive and proliferate within haemocytes, the main cellular effectors of the immune system in molluscs. Two-dimensional electrophoresis was used to compare the haemolymph protein profile between Bonamia spp.-infected and non-infected oysters within 3 different stocks, a Galician stock of oysters selected for resistance against bonamiosis, a non-selected Galician stock and a selected Irish stock. Thirty-four proteins with a presumably relevant role in the oyster-Bonamia spp. interaction were identified; they were involved in major metabolic pathways, such as energy production, respiratory chain, oxidative stress, signal transduction, transcription, translation, protein degradation and cell defence. Furthermore, the haemolymph proteomic profiles of the non-infected oysters of the 2 Galician stocks were compared. As a result, 7 proteins representative of the non-infected Galician oysters selected for resistance against bonamiosis were identified; these 7 proteins could be considered as candidate markers of resistance to bonamiosis, which should be further assessed.

De novo assembly, gene annotation, and marker development using Illumina paired-end transcriptome sequencing in the Crassadoma gigantea

Crassadoma gigantea is an important commercial marine bivalve species in Baja California and Mexico. In this study, we have applied RNA-Seq technology to profile the transcriptome of the C. gigantea for the first time. A total of 80,832,518 raw reads were produced from a Illumina HiSeq4000 platform, and 77,306,198 (95.64%) clean reads were generated after trimming the adaptor sequences. The transcriptome assembled into 158,855 transcripts with an N50 size of 1995 bp and an average size of 1008 bp. A number of DNA repair related genes, such as MSH3, EGF, TGF, IGF, FGF, encoding different groups of growth factors were found in the transcriptome data of C. gigantean. In addition, immune related genes Toll-like receptor (TLR) including TLR1, TLR2, TLR3, TLR4, TLR5, TLR6, TLR7, TLR8, and TLR9 was also observed in C. gigantean. A set of 12 polymorphic microsatellite loci were firstly developed and characterized in C. gigantea. The results show that the number of alleles and expected heterozygosity ranged from 3 to 9 and from 0.254 to 0.820, respectively. The average polymorphic information content was 0.790. These microsatellite loci will facilitate future studies of population structure and conservation genetics in this species.

The oyster immunity

Oysters, the common name for a number of different bivalve molluscs, are the worldwide aquaculture species and also play vital roles in the function of ecosystem. As invertebrate, oysters have evolved an integrated, highly complex innate immune system to recognize and eliminate various invaders via an array of orchestrated immune reactions, such as immune recognition, signal transduction, synthesis of antimicrobial peptides, as well as encapsulation and phagocytosis of the circulating haemocytes. The hematopoietic tissue, hematopoiesis, and the circulating haemocytes have been preliminary characterized, and the detailed annotation of the Pacific oyster Crassostrea gigas genome has revealed massive expansion and functional divergence of innate immune genes in this animal. Moreover, immune priming and maternal immune transfer are reported in oysters, suggesting the adaptability of invertebrate immunity. Apoptosis and autophagy are proved to be important immune mechanisms in oysters. This review will summarize the research progresses of immune system and the immunomodulation mechanisms of the primitive catecholaminergic, cholinergic, neuropeptides, GABAergic and nitric oxidase system, which possibly make oysters ideal model for studying the origin and evolution of immune system and the neuroendocrine-immune regulatory network in lower invertebrates.

Hepatic transcriptome analysis of juvenile GIFT tilapia (Oreochromis niloticus), fed diets supplemented with different concentrations of resveratrol

The GIFT (Genetically Improved Farmed Tilapia) tilapia, Oreochromis niloticus, is cultured widely for the production of freshwater fish in China. Streptococcosis, which is related to pathogenic infections, occurs frequently in juvenile and adult female GIFT individuals. Resveratrol (RES) has been used in feed to control these infections in freshwater tilapia. To address the effects of RES on tilapia, we used high-throughput RNA sequencing technology (RNA-Seq, HiSeq. 2500) to explore the global transcriptomic response and specific involvement of hepatic mRNA of juvenile O. niloticus fed with diets containing different concentrations of (0, 0.025, 0.05, and 0.1g/kg) RES. A total of > 24,513,018 clean reads were generated and then assembled into 23,244 unigenes. The unigenes were annotated by comparing them against non-redundant protein sequence (Nr), non-redundant nucleotide (Nt), Swiss-Prot, Pfam, Gene Ontology database (GO), Clusters of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, and 12,578 unigenes were annotated to the GO database. A total of 1444 (0.025g/kg RES), 1526 (0.05g/kg RES), and 3135 (0.1g/kg RES) genes were detected as significant differentially expressed genes (DEGs), when compared with the controls. A total of 6 (0.025 vs 0.05g/kg RES), 19 (0.025 vs 0.1g/kg RES), and 124 (0.05 vs 0.1g/kg RES) genes were detected as significant DEGs. Six genes, including dnah7x1, sox4, fam46a, hsp90a, ddit4, and nmrk2, were associated with an immune response. These findings provide information on the innate immunity of GIFT and might contribute to the development of strategies for the effective management of diseases and long-term sustainability of O. niloticus culture.

Transcriptome analysis of Mytilus coruscus hemocytes in response to Vibrio alginnolyficus infection

As an economically important bivalve, the Mytilus coruscus is cultured widely in the eastern coast of China. In recent years, this bivalve has been seriously affected by the pathogenic infections. To elucidate the host defense mechanisms of M. coruscus against pathogenic challenge, the hemocyte transcriptomes of M. coruscus before and after Vibrio alginnolyficus infection were analyzed using the deep-sequencing platform Illumina/HiSeq-2500, meanwhile the differentially expressed genes (DEGs) were investigated. In total, 130,031,083 clean reads were obtained and then assembled into 63,942 unigenes with an average length of 810 bp and an N50 of 1056 bp. Unigenes were annotated by comparing against nr, Swiss-Prot, KEGG, COG, KOG, GO, and Pfam databases, and 27,345 unigenes (42.77%) were annotated in at least one database. After bacterial challenge, 1270 and 265 genes were identified as remarkably up- or down-regulated, respectively, amongst 1154 were associated with 122 pathways, including classical immune-related pathways, such as 'Toll-like receptor signaling', 'the complement cascades', 'MAPK signaling pathway', 'Apoptosis' and 'Wnt signaling pathway'. Besides, nine genes which were differently-expressed immuno-related were confirmed by using quantitative real-time PCR. These findings would provide new insights on the M. coruscus innate immunity, based on which, some novel strategies for management of diseases and long-term sustainability of M. coruscus culture could be developed.

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.

Integrated application of transcriptomics and metabolomics provides insights into glycogen content regulation in the Pacific oyster Crassostrea gigas

BACKGROUND

The Pacific oyster Crassostrea gigas is an important marine fishery resource, which contains high levels of glycogen that contributes to the flavor and the quality of the oyster. However, little is known about the molecular and chemical mechanisms underlying glycogen content differences in Pacific oysters. Using a homogeneous cultured Pacific oyster family, we explored these regulatory networks at the level of the metabolome and the transcriptome.

RESULTS

Oysters with the highest and lowest natural glycogen content were selected for differential transcriptome and metabolome analysis. We identified 1888 differentially-expressed genes, seventy-five differentially-abundant metabolites, which are part of twenty-seven signaling pathways that were enriched using an integrated analysis of the interaction between the differentially-expressed genes and the differentially-abundant metabolites. Based on these results, we found that a high expression of carnitine O-palmitoyltransferase 2 (CPT2), indicative of increased fatty acid degradation, is associated with a lower glycogen content. Together, a high level of expression of phosphoenolpyruvate carboxykinase (PEPCK), and high levels of glucogenic amino acids likely underlie the increased glycogen production in high-glycogen oysters. In addition, the higher levels of the glycolytic enzymes hexokinase (HK) and pyruvate kinase (PK), as well as of the TCA cycle enzymes malate dehydrogenase (MDH) and pyruvate carboxylase (PYC), imply that there is a concomitant up-regulation of energy metabolism in high-glycogen oysters. High-glycogen oysters also appeared to have an increased ability to cope with stress, since the levels of the antioxidant glutathione peroxidase enzyme 5 (GPX5) gene were also increased.

CONCLUSION

Our results suggest that amino acids and free fatty acids are closely related to glycogen content in oysters. In addition, oysters with a high glycogen content have a greater energy production capacity and a greater ability to cope with stress. These findings will not only provide insights into the molecular mechanisms underlying oyster quality, but also promote research into the molecular breeding of oysters.

FunctionAnnotator, a versatile and efficient web tool for non-model organism annotation

Along with the constant improvement in high-throughput sequencing technology, an increasing number of transcriptome sequencing projects are carried out in organisms without decoded genome information and even on environmental biological samples. To study the biological functions of novel transcripts, the very first task is to identify their potential functions. We present a web-based annotation tool, FunctionAnnotator, which offers comprehensive annotations, including GO term assignment, enzyme annotation, domain/motif identification and predictions for subcellular localization. To accelerate the annotation process, we have optimized the computation processes and used parallel computing for all annotation steps. Moreover, FunctionAnnotator is designed to be versatile, and it generates a variety of useful outputs for facilitating other analyses. Here, we demonstrate how FunctionAnnotator can be helpful in annotating non-model organisms. We further illustrate that FunctionAnnotator can estimate the taxonomic composition of environmental samples and assist in the identification of novel proteins by combining RNA-Seq data with proteomics technology. In summary, FunctionAnnotator can efficiently annotate transcriptomes and greatly benefits studies focusing on non-model organisms or metatranscriptomes. FunctionAnnotator, a comprehensive annotation web-service tool, is freely available online at: http://fa.cgu.edu.tw/. This new web-based annotator will shed light on field studies involving organisms without a reference genome.

High throughput sequencing of RNA transcriptomes in Ruditapes philippinarum identifies genes involved in osmotic stress response

Ruditapes philippinarum, is an economically important marine bivalve species. The ability to cope with low salinity stress is quite important for the survival of aquatic species under natural conditions. In this study, the transcriptional response of the Manila clam to low salinity stress was characterized using RNA sequencing. The transcriptomes of a low salinity-treatment group (FRp1, FRp2), which survived under low salinity stress, and control group (SRp1, SRp2), which was not subjected to low salinity stress, were sequenced with the Illumina HiSeq platform. A total of 196,578 unigenes were generated. GO and KEGG analyses revealed that signal transduction, immune response, cellular component organization or biogenesis, and energy production processes were the most highly enriched pathways among the genes that were differentially expressed under low salinity stress. All these pathways could be assigned to the following biological functions in the low salinity tolerant Manila clam: signal response to low salinity stress, antioxidant response, intracellular free amino acid transport and metabolism, energy production and conversion, cell signaling pathways, and regulation of ionic content and cell volume. In summary, this is the first study using high-throughput sequencing to identify gene targets that could explain osmotic regulation mechanisms under salinity stress in R. philippinarum.

A Genome-Wide Association Study Identifies the Genomic Region Associated with Shell Color in Yesso Scallop, Patinopecten yessoensis

The shell color polymorphism widely exists in economic shellfish, which not only results in a better visual perception but also shows great value as an economic trait for breeding. Small numbers of reddish-orange shell Yesso scallops, Patinopecten yessoensis, were found in cultured populations compared to the brown majority. In this study, a genome-wide association study was conducted to understand the genetic basis of shell color. Sixty-six 2b-RAD libraries with equal numbers of reddish-orange and brown shell individuals were constructed and sequenced using the Illumina HiSeq 2000 platform. A total of 322,332,684 high-quality reads were obtained, and the average sequencing depth was 18.4×. One genomic region on chromosome 11 that included 239 single-nucleotide polymorphisms (SNPs) was identified as significantly associated with shell color. After verification by high-resolution melting in another population, two SNPs were selected as specific loci for reddish-orange shell color. These two SNPs could be used to improve the selective breeding progress of true-breeding strains with complete reddish-orange scallops. In addition, within the significantly associated genomic region, candidate genes were identified using marker sequences to search the draft genome of Yesso scallop. Three genes (LDLR, FRIS, and FRIY) with known functions in carotenoid metabolism were identified. Further study using high-performance liquid chromatography proved that the relative level of carotenoids in the reddish-orange shells was 40 times higher than that in the brown shells. These results suggested that the accumulation of carotenoids contributes to the formation of reddish-orange shells.

Effects of Florida Red Tides on histone variant expression and DNA methylation in the Eastern oyster Crassostrea virginica

Massive algal proliferations known as Harmful Algal Blooms (HABs) represent one of the most important threats to coastal areas. Among them, the so-called Florida Red Tides (FRTs, caused by blooms of the dinoflagellate Karenia brevis and associated brevetoxins) are particularly detrimental in the southeastern U.S., causing high mortality rates and annual losses in excess of $40 million. The ability of marine organisms to cope with environmental stressors (including those produced during HABs) is influenced by genetic and epigenetic mechanisms, the latter resulting in phenotypic changes caused by heritable modifications in gene expression, without involving changes in the genetic (DNA) sequence. Yet, studies examining cause-effect relationships between environmental stressors, specific epigenetic mechanisms and subsequent responses are still lacking. The present work contributes to increase this knowledge by investigating the effects of Florida Red Tides on two types of mechanisms participating in the epigenetic memory of Eastern oysters: histone variants and DNA methylation. For that purpose, a HAB simulation was conducted in laboratory conditions, exposing oysters to increasing concentrations of K. brevis. The obtained results revealed, for the first time, the existence of H2A.X, H2A.Z and macroH2A genes in this organism, encoding histone variants potentially involved in the maintenance of genome integrity during responses to the genotoxic effect of brevetoxins. Additionally, an increase in H2A.X phosphorylation (γH2A.X, a marker of DNA damage) and a decrease in global DNA methylation were observed as the HAB simulation progressed. Overall, the present work provides a basis to better understand how epigenetic mechanisms participate in responses to environmental stress in marine invertebrates, opening new avenues to incorporate environmental epigenetics approaches into management and conservation programs.

Transcriptome profiles of the clam Meretrix petechialis hepatopancreas in response to Vibrio infection

Microbial diseases have received much attention due to their enormous destruction of aquaculture, and Vibrio parahaemolyticus is one of the main pathogens that cause bacterial disease in the clam Meretrix petechialis. To better understand the molecular mechanisms of the immune response to Vibrio in M. petechialis, RNA-Seq was applied to explore global expression changes of hepatopancreas from this clam after Vibrio challenge. There were 199,318,966 clean reads obtained by Illumina sequencing, which were further assembled into 214,577 transcripts, and then 147,255 unigenes with an N50 of 1393 bp were identified. Gene ontology (GO) analysis revealed 21 biological process subcategories, 15 cellular component subcategories and 12 molecular function subcategories. A total of 8358 unigenes were mapped onto 267 biological signaling pathways by KEGG, among which there were 16 pathways related to the immune system. In total, 206 differentially expressed genes (DEGs) were identified, including 113 up-regulated unigenes and 93 down-regulated unigenes. In these DEGs, 96 DEGs were annotated in at least one database, accounting for 46.60% of all significant DEGs. To validate the transcriptome dataset, 15 DEGs were selected for real-time qPCR confirmation and the results showed that expression patterns of 13 genes (86.7%) agreed well with the RNA-Seq analysis. Fourteen of the 206 DEGs were annotated to be immune-related genes, and we examined the expression patterns of four immune-related DEGs using clams post immersion challenge. This study enriched the M. petechialis transcriptome database and provided insight into the immune response of M. petechialis against Vibrio infection.

Aquaculture genomics, genetics and breeding in the United States: current status, challenges, and priorities for future research

Advancing the production efficiency and profitability of aquaculture is dependent upon the ability to utilize a diverse array of genetic resources. The ultimate goals of aquaculture genomics, genetics and breeding research are to enhance aquaculture production efficiency, sustainability, product quality, and profitability in support of the commercial sector and for the benefit of consumers. In order to achieve these goals, it is important to understand the genomic structure and organization of aquaculture species, and their genomic and phenomic variations, as well as the genetic basis of traits and their interrelationships. In addition, it is also important to understand the mechanisms of regulation and evolutionary conservation at the levels of genome, transcriptome, proteome, epigenome, and systems biology. With genomic information and information between the genomes and phenomes, technologies for marker/causal mutation-assisted selection, genome selection, and genome editing can be developed for applications in aquaculture. A set of genomic tools and resources must be made available including reference genome sequences and their annotations (including coding and non-coding regulatory elements), genome-wide polymorphic markers, efficient genotyping platforms, high-density and high-resolution linkage maps, and transcriptome resources including non-coding transcripts. Genomic and genetic control of important performance and production traits, such as disease resistance, feed conversion efficiency, growth rate, processing yield, behaviour, reproductive characteristics, and tolerance to environmental stressors like low dissolved oxygen, high or low water temperature and salinity, must be understood. QTL need to be identified, validated across strains, lines and populations, and their mechanisms of control understood. Causal gene(s) need to be identified. Genetic and epigenetic regulation of important aquaculture traits need to be determined, and technologies for marker-assisted selection, causal gene/mutation-assisted selection, genome selection, and genome editing using CRISPR and other technologies must be developed, demonstrated with applicability, and application to aquaculture industries.Major progress has been made in aquaculture genomics for dozens of fish and shellfish species including the development of genetic linkage maps, physical maps, microarrays, single nucleotide polymorphism (SNP) arrays, transcriptome databases and various stages of genome reference sequences. This paper provides a general review of the current status, challenges and future research needs of aquaculture genomics, genetics, and breeding, with a focus on major aquaculture species in the United States: catfish, rainbow trout, Atlantic salmon, tilapia, striped bass, oysters, and shrimp. While the overall research priorities and the practical goals are similar across various aquaculture species, the current status in each species should dictate the next priority areas within the species. This paper is an output of the USDA Workshop for Aquaculture Genomics, Genetics, and Breeding held in late March 2016 in Auburn, Alabama, with participants from all parts of the United States.

De novo transcriptome analysis and microsatellite marker development for population genetic study of a serious insect pest, Rhopalosiphum padi (L.) (Hemiptera: Aphididae)

The bird cherry-oat aphid, Rhopalosiphum padi (L.), is one of the most abundant aphid pests of cereals and has a global distribution. Next-generation sequencing (NGS) is a rapid and efficient method for developing molecular markers. However, transcriptomic and genomic resources of R. padi have not been investigated. In this study, we used transcriptome information obtained by RNA-Seq to develop polymorphic microsatellites for investigating population genetics in this species. The transcriptome of R. padi was sequenced on an Illumina HiSeq 2000 platform. A total of 114.4 million raw reads with a GC content of 40.03% was generated. The raw reads were cleaned and assembled into 29,467 unigenes with an N50 length of 1,580 bp. Using several public databases, 82.47% of these unigenes were annotated. Of the annotated unigenes, 8,022 were assigned to COG pathways, 9,895 were assigned to GO pathways, and 14,586 were mapped to 257 KEGG pathways. A total of 7,936 potential microsatellites were identified in 5,564 unigenes, 60 of which were selected randomly and amplified using specific primer pairs. Fourteen loci were found to be polymorphic in the four R. padi populations. The transcriptomic data presented herein will facilitate gene discovery, gene analyses, and development of molecular markers for future studies of R. padi and other closely related aphid species.

Transcriptomic analysis of Crassostrea sikamea × Crassostrea angulata hybrids in response to low salinity stress

Hybrid oysters often show heterosis in growth rate, weight, survival and adaptability to extremes of salinity. Oysters have also been used as model organisms to study the evolution of host-defense system. To gain comprehensive knowledge about various physiological processes in hybrid oysters under low salinity stress, we performed transcriptomic analysis of gill tissue of Crassostrea sikamea ♀ × Crassostrea angulata♂ hybrid using the deep-sequencing platform Illumina HiSeq. We exploited the high-throughput technique to delineate differentially expressed genes (DEGs) in oysters maintained in hypotonic conditions. A total of 199,391 high quality unigenes, with average length of 644 bp, were generated. Of these 35 and 31 genes showed up- and down-regulation, respectively. Functional categorization and pathway analysis of these DEGs revealed enrichment for immune mechanism, apoptosis, energy metabolism and osmoregulation under low salinity stress. The expression patterns of 41 DEGs in hybrids and their parental species were further analyzed by quantitative real-time PCR (qRT-PCR). This study will serve as a platform for subsequent gene expression analysis regarding environmental stress. Our findings will also provide valuable information about gene expression to better understand the immune mechanism, apoptosis, energy metabolism and osmoregulation in hybrid oysters under low salinity stress.

Infectious diseases of marine molluscs and host responses as revealed by genomic tools

More and more infectious diseases affect marine molluscs. Some diseases have impacted commercial species including MSX and Dermo of the eastern oyster, QPX of hard clams, withering syndrome of abalone and ostreid herpesvirus 1 (OsHV-1) infections of many molluscs. Although the exact transmission mechanisms are not well understood, human activities and associated environmental changes often correlate with increased disease prevalence. For instance, hatcheries and large-scale aquaculture create high host densities, which, along with increasing ocean temperature, might have contributed to OsHV-1 epizootics in scallops and oysters. A key to understanding linkages between the environment and disease is to understand how the environment affects the host immune system. Although we might be tempted to downplay the role of immunity in invertebrates, recent advances in genomics have provided insights into host and parasite genomes and revealed surprisingly sophisticated innate immune systems in molluscs. All major innate immune pathways are found in molluscs with many immune receptors, regulators and effectors expanded. The expanded gene families provide great diversity and complexity in innate immune response, which may be key to mollusc's defence against diverse pathogens in the absence of adaptive immunity. Further advances in host and parasite genomics should improve our understanding of genetic variation in parasite virulence and host disease resistance.

Moving from Histopathology to Molecular Tools in the Diagnosis of Molluscs Diseases of Concern under EU Legislation

One of the main factors limiting molluscs production is the presence of pathogens and diseases. Disease agent transfer via transfers of live molluscs has been a major cause of disease outbreaks and epizootics. Because of that, the European Union has adopted several decisions and directives, the last in 2006 (2006/88/EC) to control movements of marine organisms over the European countries. Once the disease is established in a determined area its eradication is a complicated task because life cycle of pathogens are not completely known and only a good and early diagnosis of the disease could be the most appropriate way to deal with it. Besides, molluscs do not have an adaptive immune response and vaccination strategies are not possible. Molluscs listed diseases under EU legislation are mainly protozoan parasites, that's why histological techniques are recognized for their diagnosis. However, molecular techniques are being increasingly used primarily as confirmatory techniques of the presence of the pathogens but also in disease monitoring programs. Research perspectives are mainly focussed in the optimization, of the already described techniques to gain in sensitivity and sensibility and in the development of new molecular biology techniques (quantitative real time PCRs), that are faster and easier to apply and that allow a positive diagnosis even in early stages of infection. However, molecular tools detect DNA sequences of the pathogen which does not imply that pathogen is viable in the cell host and the infection is established. Consequently, it needs to be validated against other techniques, such as histology or in situ hybridization, so that its reliability can be determined.

Sequencing and de novo assembly of visceral mass transcriptome of the critically endangered land snail Satsuma myomphala: Annotation and SSR discovery

Satsuma myomphala is critically endangered through loss of natural habitats, predation by natural enemies, and indiscriminate collection. It is a protected species in Korea but lacks genomic resources for an understanding of varied functional processes attributable to evolutionary success under natural habitats. For assessing the genetic information of S. myomphala, we performed for the first time, de novo transcriptome sequencing and functional annotation of expressed sequences using Illumina Next-Generation Sequencing (NGS) platform and bioinformatics analysis. We identified 103,774 unigenes of which 37,959, 12,890, and 17,699 were annotated in the PANM (Protostome DB), Unigene, and COG (Clusters of Orthologous Groups) databases, respectively. In addition, 14,451 unigenes were predicted under Gene Ontology functional categories, with 4581 assigned to a single category. Furthermore, 3369 sequences with 646 having Enzyme Commission (EC) numbers were mapped to 122 pathways in the Kyoto Encyclopedia of Genes and Genomes Pathway database. The prominent protein domains included the Zinc finger (C2H2-like), Reverse Transcriptase, Thioredoxin-like fold, and RNA recognition motif domain. Many unigenes with homology to immunity, defense, and reproduction-related genes were screened in the transcriptome. We also detected 3120 putative simple sequence repeats (SSRs) encompassing dinucleotide to hexanucleotide repeat motifs from >1kb unigene sequences. A list of PCR primers of SSR loci have been identified to study the genetic polymorphisms. The transcriptome data represents a valuable resource for further investigations on the species genome structure and biology. The unigenes information and microsatellites would provide an indispensable tool for conservation of the species in natural and adaptive environments.

Transcriptomic evaluation of the American oyster, Crassostrea virginica, deployed during the Deepwater Horizon oil spill: Evidence of an active hydrocarbon response pathway

Estuarine organisms were impacted by the Deepwater Horizon oil spill which released ∼5 million barrels of crude oil into the Gulf of Mexico in the spring and summer of 2010. Crassostrea virginica, the American oyster, is a keystone species in these coastal estuaries and is routinely used for environmental monitoring purposes. However, very little is known about their cellular and molecular responses to hydrocarbon exposure. In response to the spill, a monitoring program was initiated by deploying hatchery-reared oysters at three sites along the Alabama and Mississippi coast (Grand Bay, MS, Fort Morgan, AL, and Orange Beach, AL). Oysters were deployed for 2-month periods at five different time points from May 2010 to May 2011. Gill and digestive gland tissues were harvested for gene expression analysis and determination of aliphatic and polycyclic aromatic hydrocarbon (PAH) concentrations. To facilitate identification of stress response genes that may be involved in the hydrocarbon response, a nearly complete transcriptome was assembled using Roche 454 and Illumina high-throughput sequencing from RNA samples obtained from the gill and digestive gland tissues of deployed oysters. This effort resulted in the assembly and annotation of 27,227 transcripts comprised of a large assortment of stress response genes, including members of the aryl hydrocarbon receptor (AHR) pathway, Phase I and II biotransformation enzymes, antioxidant enzymes and xenobiotic transporters. From this assembly several potential biomarkers of hydrocarbon exposure were chosen for expression profiling, including the AHR, two cytochrome P450 1A genes (CYP1A-like 1 and CYP1A-like 2), Cu/Zn superoxide dismutase (CuZnSOD), glutathione S-transferase theta (GST theta) and multidrug resistance protein 3 (MRP3). Higher expression levels of GST theta and MRP3 were observed in gill tissues from all three sites during the summer to early fall 2010 deployments. Linear regression analysis indicated a statistically significant relationship between total PAH levels in digestive gland tissue samples with CYP1A-like 2, CuZnSOD, GST theta and MRP3 induction. These observations provide evidence of a potentially conserved AHR pathway in invertebrates and yield new insight into the development of novel biomarkers for use in environmental monitoring activities.

De Novo Transcriptome Analysis of Differential Functional Gene Expression in Largemouth Bass (Micropterus salmoides) after Challenge with Nocardia seriolae

Largemouth bass (Micropterus salmoides) are common hosts of an epizootic bacterial infection by Nocardia seriolae. We conducted transcriptome profiling of M. salmoides to understand the host immune response to N. seriolae infection, using the Illumina sequencing platform. De novo assembly of paired-end reads yielded 47,881 unigenes, the total length, average length, N50, and GC content of which were 49,734,288, 1038, 1983 bp, and 45.94%, respectively. Annotation was performed by comparison against non-redundant protein sequence (NR), non-redundant nucleotide (NT), Swiss-Prot, Clusters of Orthologous Groups (COG), Kyoto Encyclopaedia of Genes and Genomes (KEGG), Gene Ontology (GO), and Interpro databases, yielding 28,964 (NR: 60.49%), 36,686 (NT: 76.62%), 24,830 (Swissprot: 51.86%), 8913 (COG: 18.61%), 20,329 (KEGG: 42.46%), 835 (GO: 1.74%), and 22,194 (Interpro: 46.35%) unigenes. Additionally, 8913 unigenes were classified into 25 Clusters of Orthologous Groups (KOGs) categories, and 20,329 unigenes were assigned to 244 specific signalling pathways. RNA-Seq by Expectation Maximization (RSEM) and PossionDis were used to determine significantly differentially expressed genes (False Discovery Rate (FDR) < 0.05) and we found that 1384 were upregulated genes and 1542 were downregulated genes, and further confirmed their regulations using reverse transcription quantitative PCR (RT-qPCR). Altogether, these results provide information on immune mechanisms induced during bacterial infection in largemouth bass, which may facilitate the prevention of nocardiosis.

Transcriptome sequencing and de novo characterization of Korean endemic land snail, Koreanohadra kurodana for functional transcripts and SSR markers

The Korean endemic land snail Koreanohadra kurodana (Gastropoda: Bradybaenidae) found in humid areas of broadleaf forests and shrubs have been considered vulnerable as the number of individuals are declining in recent years. The species is poorly characterized at the genomic level that limits the understanding of functions at the molecular and genetics level. In the present study, we performed de novo transcriptome sequencing to produce a comprehensive transcript dataset of visceral mass tissue of K. kurodana by the Illumina paired-end sequencing technology. Over 234 million quality reads were assembled to a total of 315,924 contigs and 191,071 unigenes, with an average and N50 length of 585.6 and 715 bp and 678 and 927 bp, respectively. Overall, 36.32 % of the unigenes found matches to known protein/nucleotide sequences in the public databases. The direction of the unigenes to functional categories was determined using COG, GO, KEGG, and InterProScan protein domain search. The GO analysis search resulted in 22,967 unigenes (12.02 %) being categorized into 40 functional groups. The KEGG annotation revealed that metabolism pathway genes were enriched. The most prominent protein motifs include the zinc finger, ribonuclease H, reverse transcriptase, and ankyrin repeat domains. The simple sequence repeats (SSRs) identified from >1 kb length of unigenes show a dominancy of dinucleotide repeat motifs followed with tri- and tetranucleotide motifs. A number of unigenes were putatively assessed to belong to adaptation and defense mechanisms including heat shock proteins 70, Toll-like receptor 4, AMP-activated protein kinase, aquaporin-2, etc. Our data provide a rich source for the identification and functional characterization of new genes and candidate polymorphic SSR markers in K. kurodana. The availability of transcriptome information ( http://bioinfo.sch.ac.kr/submission/ ) would promote the utilization of the resources for phylogenetics study and genetic diversity assessment.

De novo assembly and annotation of the European abalone Haliotis tuberculata transcriptome

The European abalone Haliotis tuberculata is a delicacy and consequently a commercially valuable gastropod species. Aquaculture production and wild populations are subjected to multiple climate-associated stressors and anthropogenic pressures, including rising sea-surface temperatures, ocean acidification and an emerging pathogenic Vibrio infection. Transcript expression data provides a valuable resource for understanding abalone responses to variation in the biotic and abiotic environment. To generate an extensive transcriptome, we performed next-generation sequencing of RNA on larvae exposed to temperature and pH variation and on haemolymph of adults from two wild populations after experimental infection with Vibrio harveyi. We obtained more than 1.5 billion raw paired-end reads, which were assembled into 328,519 contigs. Filtration and clustering produced a transcriptome of 41,099 transcripts, of which 10,626 (25.85%) were annotated with Blast hits, and 7380 of these were annotated with Gene Ontology (GO) terms in Blast2Go. A differential expression analysis comparing all samples from the two life stages identified 5690 and 10,759 transcripts with significantly higher expression in larvae and adult haemolymph respectively. This is the greatest sequencing effort yet in the Haliotis genus, and provides the first high-throughput transcriptomic resource for H. tuberculata.

Combined exposure to pyrene and fluoranthene and their molecular effects on the Sydney rock oyster, Saccostrea glomerata

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitously detected in the water column, associated with particulate matter or in the tissue of marine organisms such as molluscs. PAH exposure and their resultant bioaccumulation in molluscs can cause a range of serious physiological effects in the affected animals. To examine the molecular response of these xenobiotics in bivalves, Sydney rock oysters (Saccostrea glomerata) were exposed to pyrene and fluoranthene for seven days. Chemical analysis of the soft-tissue of PAH stressed S. glomerata confirmed that pyrene and fluoranthene could be bioaccumulated by these oysters. RNA-Seq analysis of PAH-exposed S. glomerata showed a total of 765 transcripts differentially expressed between control and PAH-stressed oysters. Closer examination of the transcripts revealed a range genes encoding enzymes involved in PAH detoxification (e.g. cytochrome P450), innate immune responses (e.g. pathogen recognition, phagocytosis) and protein synthesis. Overall, pyrene and fluoranthene exposure appears to have resulted in a suppression of pathogen recognition and some protein synthesis processes, whereas transcripts of genes encoding proteins involved in clearance of cell debris and some transcripts of genes involved in PAH detoxification were induced in response to the stressors. Pyrene and fluoranthene exposure thus invoked a complex molecular response in S. glomerata, with results suggesting that oysters focus on removing the stressors from their system and dealing with the downstream effects of PAH exposure, potentially at the exclusion of other, less immediate concerns (e.g. protection from infection).

Transcriptomic responses to low temperature stress in the Manila clam, Ruditapes philippinarum

The Manila clam, Ruditapes philippinarum, is an economically important shellfish in marine aquaculture, with a broad thermal tolerance. The ability to cope with cold stress is quite important for the survival of aquatic species under natural conditions. A cold-tolerant clam that can survive the winter at temperatures below 0 °C might extend our understanding of the mechanisms underlying the response to cold stress. In this study, the transcriptional response of the Manila clam to cold stress (-1 °C) was characterized using RNA sequencing. The transcriptomes of a cold-treatment (O) group of clams, which survived under cold stress, and the control group (OC2), which was not subjected to cold stress, were sequenced with the Illumina HiSeq platform. In all, 148,593 unigenes were generated. Compared with the unigene expression profile of the control group, 1760 unigenes were up regulated and 2147 unigenes were down regulated in the O group. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses revealed that signal transduction, mitochondrial metabolism, cellular component organization or biogenesis, and energy production processes were the most highly enriched pathways among the genes that were differentially expressed under cold stress. All these pathways could be assigned to the following biological functions in the cold-tolerant Manila clam: signal response to cold stress, antioxidant response, cell proliferation, and energy production.

Transcriptome Analysis of the Sydney Rock Oyster, Saccostrea glomerata: Insights into Molluscan Immunity

BACKGROUND

Oysters have important ecological functions in their natural environment, acting as global carbon sinks and improving water quality by removing excess nutrients from the water column. During their life-time oysters are exposed to a variety of pathogens that can cause severe mortality in a range of oyster species. Environmental stressors encountered in their habitat can increase the susceptibility of oysters to these pathogens and in general have been shown to impact on oyster immunity, making immune parameters expressed in these marine animals an important research topic.

RESULTS

Paired-end Illumina high throughput sequencing of six S. glomerata tissues exposed to different environmental stressors resulted in a total of 484,121,702 paired-end reads. When reads and assembled transcripts were compared to the C. gigas genome, an overall low level of similarity at the nucleotide level, but a relatively high similarity at the protein level was observed. Examination of the tissue expression pattern showed that some transcripts coding for cathepsins, heat shock proteins and antioxidant proteins were exclusively expressed in the haemolymph of S. glomerata, suggesting a role in innate immunity. Furthermore, analysis of the S. glomerata ORFs showed a wide range of genes potentially involved in innate immunity, from pattern recognition receptors, components of the Toll-like signalling and apoptosis pathways to a complex antioxidant defence mechanism.

CONCLUSIONS

This is the first large scale RNA-Seq study carried out in S. glomerata, showing the complex network of innate immune components that exist in this species. The results confirmed that many of the innate immune system components observed in mammals are also conserved in oysters; however, some, such as the TLR adaptors MAL, TRIF and TRAM are either missing or have been modified significantly. The components identified in this study could help explain the oysters' natural resilience against pathogenic microorganisms encountered in their natural environment.

Multi-species protein similarity clustering reveals novel expanded immune gene families in the eastern oyster Crassostrea virginica

Comparative genomics research in non-model species has highlighted how invertebrate hosts possess complex diversified repertoires of immune molecules. The levels of diversification in particular immune gene families appear to differ between invertebrate lineages and even between species within lineages, reflecting differences not only in evolutionary histories, but also in life histories, environmental niches, and pathogen exposures. The goal of this research was to identify immune-related gene families experiencing high levels of diversification in eastern oysters, Crassostrea virginica. Families containing 1) transcripts differentially expressed in eastern oysters in response to bacterial challenge and 2) a larger number of transcripts compared to other species included those coding for the C1q and C-type lectin domain containing proteins (C1qDC and CTLDC), GTPase of the immune-associated proteins (GIMAP), scavenger receptors (SR), fibrinogen-C domain containing proteins (also known as FREPs), dopamine beta-hydrolase (DBH), interferon-inducible 44 (IFI44), serine protease inhibitors, apextrin, and dermatopontin. Phylogenetic analysis of two of the families significantly expanded in bivalves, IFI44 and GIMAP, showed a patchy distribution within both protostomes and deuterostomes, suggesting multiple independent losses and lineage-specific expansions. Increased availability of genomic information for a broader range of non-model species broadly distributed through vertebrate and invertebrate phyla will likely lead to improved knowledge on mechanisms of immune-gene diversification.

De Novo Transcriptome Analysis of Two Seahorse Species (Hippocampus erectus and H. mohnikei) and the Development of Molecular Markers for Population Genetics

Seahorse conservation has been performed utilizing various strategies for many decades, and the deeper understanding of genomic information is necessary to more efficiently protect the germplasm resources of seahorse species. However, little genetic information about seahorses currently exists in the public databases. In this study, high-throughput RNA sequencing for two seahorse species, Hippocampus erectus and H. mohnikei, was carried out, and de novo assembly generated 37,506 unigenes for H. erectus and 36,113 unigenes for H. mohnikei. Among them, 17,338 (46.23%) unigenes for H. erectus and 17,900 (49.57%) for H. mohnikei were successfully annotated based on the information available from the public databases. Through comparing the unigenes of two seahorse species, 7,802 candidate orthologous genes were identified and 5,268 genes among them could be annotated. In addition, gene ontology analysis of two species was similarly performed on biological processes, cellular components, and molecular functions. Twenty-four and twenty-one unigenes in H. erectus and H. mohnikei were annotated in the biosynthesis of unsaturated fatty acids pathways, and both seahorses lacked the Δ12 and Δ15 desaturases. Total of 8,992 and 9,116 SSR loci were obtained from H. erectus and H. mohnikei unigenes, respectively. Dozens of SSR were developed and then applied to assess the population genetic diversity, as well as cross-amplified in a related species, H. trimaculatus. The HO and HE values of the tested populations for H. erectus, H. mohnikei, and H. trimaculatus were medium. These resources would facilitate the conservation of the species through a better understanding of the genomics and comparative genome analysis within the Hippocampus genus.

Identification of Immune-Related Genes and Development of SSR/SNP Markers from the Spleen Transcriptome of Schizothorax prenanti

Schizothorax prenanti (S. prenanti) is mainly distributed in the upstream regions of the Yangtze River and its tributaries in China. This species is indigenous and commercially important. However, in recent years, wild populations and aquacultures have faced the serious challenges of germplasm variation loss and an increased susceptibility to a range of pathogens. Currently, the genetics and immune mechanisms of S. prenanti are unknown, partly due to a lack of genome and transcriptome information. Here, we sought to identify genes related to immune functions and to identify molecular markers to study the function of these genes and for trait mapping. To this end, the transcriptome from spleen tissues of S. prenanti was analyzed and sequenced. Using paired-end reads from the Illumina Hiseq2500 platform, 48,517 transcripts were isolated from the spleen transcriptome. These transcripts could be clustered into 37,785 unigenes with an N50 length of 2,539 bp. The majority of the unigenes (35,653, 94.4%) were successfully annotated using non-redundant nucleotide sequence analysis (nt), and the non-redundant protein (nr), Swiss-Prot, Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. KEGG pathway assignment identified more than 500 immune-related genes. Furthermore, 7,545 putative simple sequence repeats (SSRs), 857,535 single nucleotide polymorphisms (SNPs), and 53,481 insertion/deletion (InDels) were detected from the transcriptome. This is the first reported high-throughput transcriptome analysis of S. prenanti, and it provides valuable genetic resources for the investigation of immune mechanisms, conservation of germplasm, and molecular marker-assisted breeding of S. prenanti.

Sequencing, De Novo Assembly, and Annotation of the Transcriptome of the Endangered Freshwater Pearl Bivalve, Cristaria plicata, Provides Novel Insights into Functional Genes and Marker Discovery

BACKGROUND

The freshwater mussel Cristaria plicata (Bivalvia: Eulamellibranchia: Unionidae), is an economically important species in molluscan aquaculture due to its use in pearl farming. The species have been listed as endangered in South Korea due to the loss of natural habitats caused by anthropogenic activities. The decreasing population and a lack of genomic information on the species is concerning for environmentalists and conservationists. In this study, we conducted a de novo transcriptome sequencing and annotation analysis of C. plicata using Illumina HiSeq 2500 next-generation sequencing (NGS) technology, the Trinity assembler, and bioinformatics databases to prepare a sustainable resource for the identification of candidate genes involved in immunity, defense, and reproduction.

RESULTS

The C. plicata transcriptome analysis included a total of 286,152,584 raw reads and 281,322,837 clean reads. The de novo assembly identified a total of 453,931 contigs and 374,794 non-redundant unigenes with average lengths of 731.2 and 737.1 bp, respectively. Furthermore, 100% coverage of C. plicata mitochondrial genes within two unigenes supported the quality of the assembler. In total, 84,274 unigenes showed homology to entries in at least one database, and 23,246 unigenes were allocated to one or more Gene Ontology (GO) terms. The most prominent GO biological process, cellular component, and molecular function categories (level 2) were cellular process, membrane, and binding, respectively. A total of 4,776 unigenes were mapped to 123 biological pathways in the KEGG database. Based on the GO terms and KEGG annotation, the unigenes were suggested to be involved in immunity, stress responses, sex-determination, and reproduction. A total of 17,251 cDNA simple sequence repeats (cSSRs) were identified from 61,141 unigenes (size of >1 kb) with the most abundant being dinucleotide repeats.

CONCLUSIONS

This dataset represents the first transcriptome analysis of the endangered mollusc, C. plicata. The transcriptome provides a comprehensive sequence resource for the conservation of genetic information in this species and enrichment of the genetic database. The development of molecular markers will assist in the genetic improvement of C. plicata.

Alterations of the immune transcriptome in resistant and susceptible hard clams (Mercenaria mercenaria) in response to Quahog Parasite Unknown (QPX) and temperature

Quahog Parasite Unknown (QPX) is a fatal protistan parasite that causes severe losses in the hard clam (Mercenaria mercenaria) fisheries along the northeastern coast of the US. Field and laboratory studies of QPX disease have demonstrated a major role for water temperature and M. mercenaria genetic origin in disease development. Infections are more likely to occur at cold temperatures, with clam stocks originating from southern states being more susceptible than clams from northern origin where disease is enzootic. Even though the influence of temperature on QPX infection have been examined in susceptible and resistant M. mercenaria at physiological and cellular scales, the underlying molecular mechanisms associated with host-pathogen interactions remain largely unknown. This study was carried out to explore the molecular changes in M. mercenaria in response to temperature and QPX infection on the transcriptomic level, and also to compare molecular responses between susceptible and resistant clam stocks. A M. mercenaria oligoarray (15 K Agilent) platform was produced based on our previously generated transcriptomic data and was used to compare gene expression profiles in naive and QPX-infected susceptible (Florida stock) and resistant (Massachusetts) clams maintained at temperatures favoring disease development (13 °C) or clam healing (21 °C). In addition, transcriptomic changes reflecting focal (the site of infection, mantle) and systemic (circulating hemocytes) responses were also assessed using the oligoarray platform. Results revealed significant regulation of multiple biological pathways by temperature and QPX infection, mainly associated with immune recognition, microbial killing, protein synthesis, oxidative protection and metabolism. Alterations were widely systemic with most changes in gene expression revealed in hemocytes, highlighting the role of circulating hemocytes as the first line of defense against pathogenic stress. A large number of complement-related recognition molecules with fibrinogen or C1q domains were shown to be specially induced following QPX challenge, and the expression of these molecules was significantly higher in resistant clams as compared to susceptible ones. These highly variable immune proteins may be potent candidate molecular markers for future study of M. mercenaria resistance against QPX. Beyond the specific case of clam response to QPX, this study also provides insights into the primitive complement-like system in the hard clam.

Identification of functional genes involved in Cd(2+) response of Chinese surf clam (Mactra chinensis) through transcriptome sequencing

The Chinese surf clam Mactra chinensis is an economically important bivalve species in the coastal waters of Liaoning and Shandong Province, China. In this study, we carried out transcriptome sequencing to develop molecular resources for M. chinensis and conducted an acute test of Cd(2+) stimulation through quantitative real-time PCR (qRT-PCR) to analyze the relative expression of six functional genes. A total of 100,839 transcripts and 56,712 unigenes were obtained from 39.9 million filtered reads and 21,305 unigenes were annotated by hitting against NCBI database. According to the results of qRT-PCR, heat shock protein 22 (Hsp22) and cytochrome P450 (CYP450(2C31)) were inhibited in the low concentration, and induced in the high concentration of Cd(2+); thioredoxin peroxidase (TPx-A) was at normal level in low concentration, but induced in high concentration of Cd(2+); glutathione peroxidase A (GPA), glutathione peroxidase 1 (GPA1) and Mn superoxide dismutase gene (MnSOD) were down-regulated when exposed to any treatment groups. Expression levels of the six functional genes following Cd(2+) exposure indicated that these genes were linked to environmental stress. Moreover, the present work enriched the molecule genetic data of M. chinensis.

Analysis of synonymous codon usage patterns in sixty-four different bivalve species

Synonymous codon usage bias (CUB) is a defined as the non-random usage of codons encoding the same amino acid across different genomes. This phenomenon is common to all organisms and the real weight of the many factors involved in its shaping still remains to be fully determined. So far, relatively little attention has been put in the analysis of CUB in bivalve mollusks due to the limited genomic data available. Taking advantage of the massive sequence data generated from next generation sequencing projects, we explored codon preferences in 64 different species pertaining to the six major evolutionary lineages in Bivalvia. We detected remarkable differences across species, which are only partially dependent on phylogeny. While the intensity of CUB is mild in most organisms, a heterogeneous group of species (including Arcida and Mytilida, among the others) display higher bias and a strong preference for AT-ending codons. We show that the relative strength and direction of mutational bias, selection for translational efficiency and for translational accuracy contribute to the establishment of synonymous codon usage in bivalves. Although many aspects underlying bivalve CUB still remain obscure, we provide for the first time an overview of this phenomenon in this large, commercially and environmentally important, class of marine invertebrates.

De novo assembly and characterization of foot transcriptome and microsatellite marker development for Paphia textile

Paphia textile is an important, aquaculture bivalve clam species distributed mainly in China, Philippines, and Malaysia. Recent studies of P. textile have focused mainly on artificial breeding and nutrition analysis, and the transcriptome and genome of P. textile have rarely been reported. In this work, the transcriptome of P. textile foot tissue was sequenced on an Illumina HiSeq™ 2000 platform. A total of 20,219,795 reads were generated, resulting in 4.08 Gb of raw data. The raw reads were cleaned and assembled into 54,852 unigenes with an N50 length of 829 bp. Of these unigenes, 38.92% were successfully annotated based on their matches to sequences in seven public databases. Among the annotated unigenes, 14,571 were assigned Gene Ontology terms, 5448 were classified to Clusters of Orthologous Groups categories, and 6738 were mapped to 228 pathways in the Kyoto Encyclopedia of Genes and Genomes database. For functional marker development, 5605 candidate simple sequence repeats were identified in the transcriptome and 80 primer pairs were selected randomly and amplified in a wild population of P. textile. A total of 36 loci that exhibited obvious repeat length polymorphisms were detected. The transcriptomic data and microsatellite markers will provide valuable resources for future functional gene analyses, genetic map construction, and quantitative trait loci mapping in P. textile.

Integration of Transcriptomic and Proteomic Approaches Provides a Core Set of Genes for Understanding of Scallop Attachment

Attachment is an essential physiological process in life histories of many marine organisms. Using a combination of transcriptomic and proteomic approach, scallop byssal proteins (Sbps) and their associated regulatory network genes were investigated for the first time. We built the first scallop foot transcriptome library, and 75 foot-specific genes were identified. Through integration of transcriptomic-proteomic approach, seven unique Sbps were identified. Of them, three showed significant amino acid sequence homology to known proteins. In contrast, the rest did not show significant protein matches, indicating they are possibly novel proteins. Our transcriptomic and proteomic analyses also suggest that post-translational modification may be one of the significant features for Sbps as well. Taken together, our study provides the first multidimensional collection of a core set of genes that may be potentially involved in scallop byssal attachment.

IL-17 signaling components in bivalves: Comparative sequence analysis and involvement in the immune responses

The recent discovery of soluble immune-regulatory molecules in invertebrates takes advantage of the rapid growth of next generation sequencing datasets. Following protein domain searches in the transcriptomes of 31 bivalve spp. and in few available mollusk genomes, we retrieved 59 domains uniquely identifying interleukin 17 (IL-17) and 96 SEFIR domains typical of IL-17 receptors and CIKS/ACT1 proteins acting downstream in the IL-17 signaling pathway. Compared to the Chordata IL-17 family members, we confirm a separate clustering of the bivalve domain sequences and a consistent conservation pattern of amino acid residues. Analysis performed at transcript and genome level allowed us to propose an updated view of the components outlining the IL-17 signaling pathway in Mytilus galloprovincialis and Crassostrea gigas (in both species, homology modeling reduced the variety of IL-17 domains to only two 3D structures). Digital expression analysis indicated more heterogeneous expression levels for the mussel and oyster IL-17 ligands than for IL-17 receptors and CIKS/CIKSL proteins. Besides, new qPCR analyses confirmed such gene expression trends in hemocytes and gills of mussels challenged with heat-killed bacteria. These results uphold the involvement of an ancient IL-17 signaling pathway in the bivalve immune responses and, likewise in humans, suggest the possibility of distinctive modulatory roles of individual IL-17s/IL-17 receptors. Overall, the common evidence of pro-inflammatory cytokines and inter-related intracellular signaling pathways in bivalves definitely adds complexity to the invertebrate immunity.

Expression Divergence of Duplicate Genes in the Protein Kinase Superfamily in Pacific Oyster

Gene duplication has been proposed to serve as the engine of evolutionary innovation. It is well recognized that eukaryotic genomes contain a large number of duplicated genes that evolve new functions or expression patterns. However, in mollusks, the evolutionary mechanisms underlying the divergence and the functional maintenance of duplicate genes remain little understood. In the present study, we performed a comprehensive analysis of duplicate genes in the protein kinase superfamily using whole genome and transcriptome data for the Pacific oyster. A total of 64 duplicated gene pairs were identified based on a phylogenetic approach and the reciprocal best BLAST method. By analyzing gene expression from RNA-seq data from 69 different developmental and stimuli-induced conditions (nine tissues, 38 developmental stages, eight dry treatments, seven heat treatments, and seven salty treatments), we found that expression patterns were significantly correlated for a number of duplicate gene pairs, suggesting the conservation of regulatory mechanisms following divergence. Our analysis also identified a subset of duplicate gene pairs with very high expression divergence, indicating that these gene pairs may have been subjected to transcriptional subfunctionalization or neofunctionalization after the initial duplication events. Further analysis revealed a significant correlation between expression and sequence divergence (as revealed by synonymous or nonsynonymous substitution rates) under certain conditions. Taken together, these results provide evidence for duplicate gene sequence and expression divergence in the Pacific oyster, accompanying its adaptation to harsh environments. Our results provide new insights into the evolution of duplicate genes and their expression levels in the Pacific oyster.

Transcriptome analysis of the pearl oyster (Pinctada fucata) hemocytes in response to Vibrio alginolyticus infection

The pearl oyster Pinctada fucata is cultured widely for production of marine pearls in China, while mass mortalities, likely related to pathogenic infections, have occurred frequently in juvenile, mother and operated oysters. To address this issue, understanding host defense mechanisms of P. fucata against pathogenic challenge is extremely important. In the present study, a comparative analysis of hemocyte transcriptomes of P. fucata before and after Vibrio alginolyticus infection was conducted using the Illumina/Hiseq-2000 RNA-Seq technology. A total of 56,345,139 clean reads were generated and then assembled into 74,007 unigenes with an average length of 680 bp and an N50 of 1197 bp. Unigenes were annotated by comparing against non-redundant protein sequence (nr), non-redundant nucleotide (nt), Swiss-Prot, Pfam, Gene Ontology database (GO), Clusters of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, and 29,615 unigenes (40.01%) were annotated in at least one database. There were 636 genes (518 up-regulated and 118 down-regulated) that were significantly differentially expressed after bacterial challenge, and among which 369 were associated with 122 pathways, including classical immune-related pathways, such as 'MAPK signaling pathway', 'Chemokine signaling pathway', 'Apoptosis' and 'Wnt signaling pathway'. These findings provide information on the pearl oyster innate immunity and may contribute to developing strategies for management of diseases and long-term sustainability of P. fucata culture.

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.

Developing tools for the study of molluscan immunity: The sequencing of the genome of the eastern oyster, Crassostrea virginica

The eastern oyster, Crassostrea virginica, provides important ecological and economical services, making it the target of restoration projects and supporting a significant fishery/aquaculture industry with landings valued at more than $100 million in 2012 in the United States of America. Due to the impact of infectious diseases on wild, restored, and cultured populations, the eastern oyster has been the focus of studies on host-pathogen interactions and immunity, as well as the target of selective breeding efforts for disease resistant oyster lines. Despite these efforts, relatively little is known about the genetic basis of resistance to diseases or environmental stress, not only in eastern oyster, but also in other molluscan species of commercial interest worldwide. In order to develop tools and resources to assist in the elucidation of the genomic basis of traits of commercial, biological, and ecological interest in oysters, a team of genome and bioinformatics experts, in collaboration with the oyster research community, is sequencing, assembling, and annotating the first reference genome for the eastern oyster and producing an exhaustive transcriptome from a variety of oyster developmental stages and tissues in response to a diverse set of environmentally-relevant stimuli. These transcriptomes and reference genome for the eastern oyster, added to the already available genome and transcriptomes for the Pacific oyster (Crassostrea gigas) and other bivalve species, will be an essential resource for the discovery of candidate genes and markers associated with traits of commercial, biological, and ecologic importance in bivalve molluscs, including those related to host-pathogen interactions and immunity.

An updated molecular basis for mussel immunity

Non-self recognition with the consequent tolerance or immune reaction is a crucial process to succeed as living organisms. At the same time the interactions between host species and their microbiome, including potential pathogens and parasites, significantly contribute to animal life diversity. Marine filter-feeding bivalves, mussels in particular, can survive also in heavily anthropized coastal waters despite being constantly surrounded by microorganisms. Based on the first outline of the Mytilus galloprovincialis immunome dated 2011, the continuously growing transcript data and the recent release of a draft mussel genome, we explored the available sequence data and scientific literature to reinforce our knowledge on the main gene-encoded elements of the mussel immune responses, from the pathogen recognition to its clearance. We carefully investigated molecules specialized in the sensing and targeting of potential aggressors, expected to show greater molecular diversification, and outlined, whenever relevant, the interconnected cascades of the intracellular signal transduction. Aiming to explore the diversity of extracellular, membrane-bound and intracellular pattern recognition receptors in mussel, we updated a highly complex immune system, comprising molecules which are described here in detail for the first time (e.g. NOD-like receptors) or which had only been partially characterized in bivalves (e.g. RIG-like receptors). Overall, our comparative sequence analysis supported the identification of over 70 novel full-length immunity-related transcripts in M. galloprovincialis. Nevertheless, the multiplicity of gene functions relevant to immunity, the involvement of part of them in other vital processes, and also the lack of a refined mussel genome make this work still not-exhaustive and support the development of more specific studies.

The complexity of apoptotic cell death in mollusks: An update

Apoptosis is a type of programmed cell death that produces changes in cell morphology and in biochemical intracellular processes without inflammatory reactions. The components of the apoptotic pathways are conserved throughout evolution. Caspases are key molecules involved in the transduction of the death signal and are responsible for many of the biochemical and morphological changes associated with apoptosis. Nowadays, It is known that caspases are activated through two major apoptotic pathways (the extrinsic or death receptor pathway and the intrinsic or mitochondrial pathway), but there are also evidences of at least other alternative pathway (the perforin/granzyme pathway). Apoptosis in mollusks seems to be similar in complexity to apoptosis in vertebrates but also has unique features maybe related to their recurrent exposure to environmental changes, pollutants, pathogens and also related to the sedentary nature of some stages in the life cycle of mollusks bivalves and gastropods. As in other animals, apoptotic process is involved in the maintenance of tissue homeostasis and also constitutes an important immune response that can be triggered by a variety of stimuli, including cytokines, hormones, toxic insults, viruses, and protozoan parasites. The main goal of this work is to present the current knowledge of the molecular mechanisms of apoptosis in mollusks and to highlight those steps that need further study.

Transcriptome analysis reveals strong and complex antiviral response in a mollusc

Viruses are highly abundant in the oceans, and how filter-feeding molluscs without adaptive immunity defend themselves against viruses is not well understood. We studied the response of a mollusc Crassostrea gigas to Ostreid herpesvirus 1 µVar (OsHV-1μVar) infections using transcriptome sequencing. OsHV-1μVar can replicate extremely rapidly after challenge of C. gigas as evidenced by explosive viral transcription and DNA synthesis, which peaked at 24 and 48 h post-inoculation, respectively, accompanied by heavy oyster mortalities. At 120 h post-injection, however, viral gene transcription and DNA load, and oyster mortality, were greatly reduced indicating an end of active infections and effective control of viral replication in surviving oysters. Transcriptome analysis of the host revealed strong and complex responses involving the activation of all major innate immune pathways that are equipped with expanded and often novel receptors and adaptors. Novel Toll-like receptor (TLR) and MyD88-like genes lacking essential domains were highly up-regulated in the oyster, possibly interfering with TLR signal transduction. RIG-1/MDA5 receptors for viral RNA, interferon-regulatory factors, tissue necrosis factors and interleukin-17 were highly activated and likely central to the oyster's antiviral response. Genes related to anti-apoptosis, oxidation, RNA and protein destruction were also highly up-regulated, while genes related to anti-oxidation were down-regulated. The oxidative burst induced by the up-regulation of oxidases and severe down-regulation of anti-oxidant genes may be important for the destruction of viral components, but may also exacerbate oyster mortality. This study provides unprecedented insights into antiviral response in a mollusc. The mobilization and complex regulation of expanded innate immune-gene families highlights the oyster genome's adaptation to a virus-rich marine environment.

Transcriptomics Analysis of Crassostrea hongkongensis for the Discovery of Reproduction-Related Genes

BACKGROUND

The reproductive mechanisms of mollusk species have been interesting targets in biological research because of the diverse reproductive strategies observed in this phylum. These species have also been studied for the development of fishery technologies in molluscan aquaculture. Although the molecular mechanisms underlying the reproductive process have been well studied in animal models, the relevant information from mollusks remains limited, particularly in species of great commercial interest. Crassostrea hongkongensis is the dominant oyster species that is distributed along the coast of the South China Sea and little genomic information on this species is available. Currently, high-throughput sequencing techniques have been widely used for investigating the basis of physiological processes and facilitating the establishment of adequate genetic selection programs.

RESULTS

The C.hongkongensis transcriptome included a total of 1,595,855 reads, which were generated by 454 sequencing and were assembled into 41,472 contigs using de novo methods. Contigs were clustered into 33,920 isotigs and further grouped into 22,829 isogroups. Approximately 77.6% of the isogroups were successfully annotated by the Nr database. More than 1,910 genes were identified as being related to reproduction. Some key genes involved in germline development, sex determination and differentiation were identified for the first time in C.hongkongensis (nanos, piwi, ATRX, FoxL2, β-catenin, etc.). Gene expression analysis indicated that vasa, nanos, piwi, ATRX, FoxL2, β-catenin and SRD5A1 were highly or specifically expressed in C.hongkongensis gonads. Additionally, 94,056 single nucleotide polymorphisms (SNPs) and 1,699 simple sequence repeats (SSRs) were compiled.

CONCLUSIONS

Our study significantly increased C.hongkongensis genomic information based on transcriptomics analysis. The group of reproduction-related genes identified in the present study constitutes a new tool for research on bivalve reproduction processes. The large group of molecular markers discovered in this study will be useful for population screening and marker assisted selection programs in C.hongkongensis aquaculture.

Repertoire and evolution of TNF superfamily in Crassostrea gigas: implications for expansion and diversification of this superfamily in Mollusca

Tumor necrosis factor superfamily (TNFSF) members represent a group of cytokines participating in diverse immunological, pathological and developmental pathways. However, compared with deuterostomia and cnidaia, the composition and evolution of TNF homologous in protostomia are still not well understood. In the present study, a total of 81 TNF superfamily (TNFSF) genes from 15 mollusk species, including 23 TNFSF genes from Crassostrea gigas, were surveyed by genome-wide bioinformatics analysis. The phylogenetic analysis showed that 14 out of 23 C. gigas TNFSF genes in five clades exhibited orthologous relationships with Pinctada fucata TNFSF genes. Moreover, there were 15 C. gigas TNFSF genes located in oyster-specific clusters, which were contributed by small-scaled tandem and/or segmental duplication events in oyster. By comparing the sequences of duplicated TNFSF pairs, exon loss and variant in exon/intron length were revealed as the major modes of divergence in gene structure. Most of the duplicated C. gigas TNFSF pairs were evolved under purifying selection with consistent tissue expression patterns, implying functional constraint shaped diversification. This study demonstrated the expansion and early divergence of TNF superfamily in C. gigas, which provides potential insight into revealing the evolution and function of this superfamily in mollusk.