Mutation in promoter region of a serine protease inhibitor confers Perkinsus marinus resistance in the eastern oyster (Crassostrea virginica)

Genome-Wide Analysis of Family I84 Protease Inhibitor Genes in Three Bivalves Reveals Important Information About the Protein Family's Evolution

Family I84 serine protease inhibitors are believed to be mollusk specific proteins involved in host defense. The molecular evolution of the family, however, remains to be understood. In this study, the genes of Family I84 protease inhibitors in 3 bivalves, Crassostrea gigas, Crassostrea virginica and Tegillarca granosa, were analyzed at the genomic level. A total of 66 Family I84 genes (22 in C. gigas, 28 in C. virginica and 16 in T. granosa) were identified from the 3 species. They distributed unevenly in the genomes involving 4 chromosomes in C. gigas and 5 chromosomes in C. virginica and T. granosa and some genes were tandemly duplicated. Most genes had 3 exons with 12 genes having 4 exons and 1 gene having 2 exons. All genes but 1 from C. gigas and 1 from T. granosa encoded peptides with a signal sequence at the N-terminus, and the properties of the predicted mature molecules were similar. Four conserved motifs were identified in the 66 amino acid sequences. Collinear analysis revealed higher collinearity between the 2 oyster species in general genes and in Family I84 genes. Phylogenetic analysis of the 66 genes with those previously reported from 3 other bivalves and 1 gastropod showed that Family I84 protease inhibitor genes from the same species tended to be grouped together in terminal branches of the constructed Maximum likelihood tree, but most internal nodes were poorly supported by the bootstrap values. In addition, differences in expression patterns between the genes of a same species were observed in the developmental stages and tissues of C. gigas and T. granosa. Moreover, the co-expression of genes within Family I84 and Family I84 genes with non-Family I84 were also detected in C. gigas and T. granosa. These results suggested that Family I84 protease inhibitor genes evolved by active duplications and structural and functional diversifications after the speciation of related mollusks, and the diversified protease inhibitor family was likely multifunctional.

Development and Evaluation of High-Density SNP Arrays for the Eastern Oyster Crassostrea virginica

The eastern oyster Crassostrea virginica is a major aquaculture species for the USA. The sustainable development of eastern oyster aquaculture depends upon the continued improvement of cultured stocks through advanced breeding technologies. The Eastern Oyster Breeding Consortium (EOBC) was formed to advance the genetics and breeding of the eastern oyster. To facilitate efficient genotyping needed for genomic studies and selection, the consortium developed two single-nucleotide polymorphism (SNP) arrays for the eastern oyster: one screening array with 566K SNPs and one breeders' array with 66K SNPs. The 566K screening array was developed based on whole-genome resequencing data from 292 oysters from Atlantic and Gulf of Mexico populations; it contains 566,262 SNPs including 47K from protein-coding genes with a marker conversion rate of 48.34%. The 66K array was developed using best-performing SNPs from the screening array, which contained 65,893 oyster SNPs including 22,984 genic markers with a calling rate of 99.34%, a concordance rate of 99.81%, and a much-improved marker conversion rate of 92.04%. Null alleles attributable to large indels were found in 13.1% of the SNPs, suggesting that copy number variation is pervasive. Both arrays provided easy identification and separation of selected stocks from wild progenitor populations. The arrays contain 31 mitochondrial SNPs that allowed unambiguous identification of Gulf mitochondrial genotypes in some Atlantic populations. The arrays also contain 756 probes from 13 oyster and human pathogens for possible detection. Our results show that marker conversion rate is low in high polymorphism species and that the two-step process of array development can greatly improve array performance. The two arrays will advance genomic research and accelerate genetic improvement of the eastern oyster by delineating genetic architecture of production traits and enabling genomic selection. The arrays also may be used to monitor pedigree and inbreeding, identify selected stocks and their introgression into wild populations, and assess the success of oyster restoration.

What Can Genetics Do for the Control of Infectious Diseases in Aquaculture?

Simple Summary

Infectious diseases place an economic burden on aquaculture and a limitation to its growth. This state-of-the-art review describes the application of genetics and genomics as novel tools to control infectious disease in aquaculture.

Abstract

Infectious diseases place an economic burden on aquaculture and a limitation to its growth. An innovative approach to mitigate their impact on production is breeding for disease resistance: selection for domestication, family-based selection, marker-assisted selection, and more recently, genomic selection. Advances in genetics and genomics approaches to the control of infectious diseases are key to increasing aquaculture efficiency, profitability, and sustainability and to reducing its environmental footprint. Interaction and co-evolution between a host and pathogen can, however, turn breeding to boost infectious disease resistance into a potential driver of pathogenic change. Parallel molecular characterization of the pathogen and its virulence and antimicrobial resistance genes is therefore essential to understand pathogen evolution over time in response to host immunity, and to apply appropriate mitigation strategies.

Pooled resequencing of larvae and adults reveals genomic variations associated with Ostreid herpesvirus 1 resistance in the Pacific oyster Crassostrea gigas

The Ostreid herpesvirus 1 (OsHV-1) is a lethal pathogen of the Pacific oyster (Crassostrea gigas), an important aquaculture species. To understand the genetic architecture of the defense against the pathogen, we studied genomic variations associated with herpesvirus-caused mortalities by pooled whole-genome resequencing of before and after-mortality larval samples as well as dead and surviving adults from a viral challenge. Analysis of the resequencing data identified 5,271 SNPs and 1,883 genomic regions covering 3,111 genes in larvae, and 18,692 SNPs and 28,314 regions covering 4,863 genes in adults that were significantly associated with herpesvirus-caused mortalities. Only 1,653 of the implicated genes were shared by larvae and adults, suggesting that the antiviral response or resistance in larvae and adults involves different sets of genes or differentiated members of expanded gene families. Combined analyses with previous transcriptomic data from challenge experiments revealed that transcription of many mortality-associated genes was also significantly upregulated by herpesvirus infection confirming their importance in antiviral response. Key immune response genes especially those encoding antiviral receptors such as TLRs and RLRs displayed strong association between variation in regulatory region and herpesvirus-caused mortality, suggesting they may confer resistance through transcriptional modulation. These results point to previously undescribed genetic mechanisms for disease resistance at different developmental stages and provide candidate polymorphisms and genes that are valuable for understanding antiviral immune responses and breeding for herpesvirus resistance.

Genome Wide Identification and Expression Profiling Indicate Expansion of Family I84 Protease Inhibitor via Gene Tandem Duplication and Divergence in Razor Clam Sinonovacula constricta

Family I84 protease inhibitors represent a novel family in the MEROPS peptidase database and are likely unique for molluscan host defense. Two Family I84 members, scSI-1 and scSI-2, were reported from the razor clam Sinonovacula constricta in a previous research. In the present study, 12 additional genes, named scSI-3 to scSI-14, were identified via genome wide sequence analyses. Among them, 10 genes were predicted to have a signal sequence, but one (scSI-7) was not. Besides, one sequence (scSI-14) was likely to encode a prematurely terminated peptide. The predicted mature peptides shared characteristics including 12 conserved cysteine residues, isoelectric points of 4.98 to 6.11, and molecular weights of 7.1 to 9.3 kDa with previously reported family members. Four motifs were characterized in 13 predicted mature peptides (with exception of scSI-14), which shared two to four conserved cysteine residues, are possibly to form two functional domain comprised 6 cysteine residues, respectively. At genomic level, all the 14 razor clam Family I84 genes were organized into 3 exons and 2 introns; 13 of them clustered in 3 regions of 100 kb on 3 separate chromosomes, suggesting tandem duplications of related genes. The promoter region of all the 14 genes was predicted to share some transcription factor binding sites, in particular those responsive to pathological and physiological stimuli, but no shared motifs were identified. Analyses also revealed differences in expression patterns among the genes. One gene in a tandem duplicated gene pairs usually showed a higher expression level than the other whereas non-tandem duplicated genes exhibited a higher degree of correlation in expression level. In addition, 8 of the 14 genes demonstrated higher level of expression in Vibrio tolerant clams than in non-tolerant clams following challenges with Vibrio parahaemolyticus. These results generated important information about the evolution of Family I84 protease inhibitors in S. constricta.

Perkinsus marinus suppresses in vitro eastern oyster apoptosis via IAP-dependent and caspase-independent pathways involving TNFR, NF-kB, and oxidative pathway crosstalk

The protozoan parasite Perkinsus marinus causes Dermo disease in eastern oysters, Crassostrea virginica, and can suppress apoptosis of infected hemocytes using incompletely understood mechanisms. This study challenged hemocytes in vitro with P. marinus for 1 h in the presence or absence of caspase inhibitor Z-VAD-FMK or Inhibitor of Apoptosis protein (IAP) inhibitor GDC-0152. Hemocytes exposure to P. marinus significantly reduced granulocyte apoptosis, and pre-incubation with Z-VAD-FMK did not affect P. marinus-induced apoptosis suppression. Hemocyte pre-incubation with GDC-0152 prior to P. marinus challenge further reduced apoptosis of granulocytes with engulfed parasite, but not mitochondrial permeabilization. This suggests P. marinus-induced apoptosis suppression may be caspase-independent, affect an IAP-involved pathway, and occur downstream of mitochondrial permeabilization. P. marinus challenge stimulated hemocyte differential expression of oxidation-reduction, TNFR, and NF-kB pathways. WGCNA analysis of P. marinus expression in response to hemocyte exposure revealed correlated protease, kinase, and hydrolase expression that could contribute to P. marinus-induced apoptosis suppression.

Transcriptomic Response to Perkinsus marinus in Two Crassostrea Oysters Reveals Evolutionary Dynamics of Host-Parasite Interactions

Infectious disease outbreaks are causing widespread declines of marine invertebrates including corals, sea stars, shrimps, and molluscs. Dermo is a lethal infectious disease of the eastern oyster Crassostrea virginica caused by the protist Perkinsus marinus. The Pacific oyster Crassostrea gigas is resistant to Dermo due to differences in the host-parasite interaction that is not well understood. We compared transcriptomic responses to P. marinus challenge in the two oysters at early and late infection stages. Dynamic and orchestrated regulation of large sets of innate immune response genes were observed in both species with remarkably similar patterns for most orthologs, although responses in C. virginica were stronger, suggesting strong or over-reacting immune response could be a cause of host mortality. Between the two species, several key immune response gene families differed in their expansion, sequence variation and/or transcriptional response to P. marinus, reflecting evolutionary divergence in host-parasite interaction. Of note, significant upregulation of inhibitors of apoptosis (IAPs) was observed in resistant C. gigas but not in susceptible C. virginica, suggesting upregulation of IAPs is an active defense mechanism, not a passive response orchestrated by P. marinus. Compared with C. gigas, C. virginica exhibited greater expansion of toll-like receptors (TLRs) and positive selection in P. marinus responsive TLRs. The C1q domain containing proteins (C1qDCs) with the galactose-binding lectin domain that is involved in P. marinus recognition, were only present and significantly upregulated in C. virginica. These results point to previously undescribed differences in host defense genes between the two oyster species that may account for the difference in susceptibility, providing an expanded portrait of the evolutionary dynamics of host-parasite interaction in lophotrochozoans that lack adaptive immunity. Our findings suggest that C. virginica and P. marinus have a history of coevolution and the recent outbreaks may be due to increased virulence of the parasite.

Family I84 protease inhibitors likely constitute a Mollusca-specific protein family functioning in host defense

Protease inhibitors are proteins or small polypeptides functioning in numerous biological processes in all organisms. The I84 family of protease inhibitors in the MEROPS database represents a novel protease inhibitor family that has been reported in 2 bivalves, Crassostrea virginica and Sinonovacula constricta, and is believed to play a role in host defense. In the present study, 7 new members of Family I84 were identified in 2 bivalves, Meretrix meretrix and Mytilus galloprovincialis, and 1 gastropod, Haliotis discus hannai, at the mRNA level via cDNA cloning. The expression patterns of the newly identified genes varied in response to salinity stresses and pathogen-associated molecular pattern stimulations, suggesting their involvement in the host defense of related species. Additionally, analyses of sequence data in public databases did not reveal any Family I84 protease inhibitor molecules in non-molluscan animals. The results indicated that Family I84 protease inhibitors are likely mollusk specific, constituting a unique host defense mechanism in molluscan species.

Vital Role of Glutamate Dehydrogenase Gene in Ammonia Detoxification and the Association Between its SNPs and Ammonia Tolerance in Sinonovacula constricta

Increasing evidence has revealed accumulated ammonia will cause adverse effects on the growth, reproduction, and survival of aquatic animals. As a marine benthic mollusk, the razor clam Sinonovacula constricta shows better growth and survival under high ammonia nitrogen environment. However, little is known about its adaptation mechanisms to high ammonia stress in an integrated mariculture system. In this study, we analyzed the association between the polymorphism of glutamate dehydrogenase gene (GDH), a key gene involved in ammonia nitrogen detoxification, and ammonia tolerance. The results showed that 26 and 22 single-nucleotide polymorphisms (SNPs) of GDH in S. constricta (denoted as Sc-GDH) were identified from two geographical populations, respectively. Among them, two SNPs (c.323T > C and c.620C > T) exhibited a significant and strong association with ammonia tolerance, suggesting that Sc-GDH gene could serve as a potential genetic marker for molecular marker–assisted selection to increase survival rate and production of S. constricta. To observe the histological morphology and explore the histocellular localization of Sc-GDH, by paraffin section and hematoxylin–eosin staining, the gills were divided into gill filament (contains columnar and flattened cells) and gill cilia, whereas hepatopancreas was made up of individual hepatocytes. The results of immunohistochemistry indicated that the columnar cells of gill filaments and the endothelial cells of hepatocytes were the major sites for Sc-GDH secretion. Under ammonia stress (180 mg/L), the expression levels of Sc-GDH were extremely significantly downregulated at 24, 48, 72, and 96 h (P < 0.01) after RNA interference. Thus, we can speculate that Sc-GDH gene may play an important role in the defense process against ammonia stress. Overall, these findings laid a foundation for further research on the adaptive mechanisms to ammonia–nitrogen tolerance for S. constricta.

Regulation Between HSF1 Isoforms and HSPs Contributes to the Variation in Thermal Tolerance Between Two Oyster Congeners

Heat shock transcription factor 1 (HSF1) plays an important role in regulating heat shock, which can activate heat shock proteins (HSPs). HSPs can protect organisms from thermal stress. Oysters in the intertidal zone can tolerate thermal stress. The Pacific oyster (Crassostrea gigas gigas) and Fujian oyster (C. gigas angulata)—allopatric subspecies with distinct thermal tolerances—make good study specimens for analyzing and comparing thermal stress regulation. We cloned and compared HSF1 isoforms, which is highly expressed under heat shock conditions in the two subspecies. The results revealed that two isoforms (HSF1a and HSF1d) respond to heat shock in both Pacific and Fujian oysters, and different heat shock conditions led to various combinations of isoforms. Subcellular localization showed that isoforms gathered in the nucleus when exposed to heat shock. The co-immunoprecipitation revealed that HSF1d can be a dimer. In addition, we selected HSPs that are expressed under the heat shock response, according to the RNA-seq and proteomic analyses. For the HSPs, we analyzed the coding part and the promoter sequences. The result showed that the domains of HSPs are conserved in two subspecies, but the promoters are significantly different. The Dual-Luciferase assay showed that the induced expression isoform HSF1d had the highest activity in C. gigas gigas, while the constitutively-expressed HSF1a was most active in C. gigas angulata. In addition, variation in the level of HSP promoters appeared to be correlated with gene expression. We argue that this gene is regulated based on the different expression levels between the two subspecies’ responses to heat shock. In summary, various stress conditions can yield different HSF1 isoforms and respond to heat shock in both oyster subspecies. Differences in how the isoforms and promoter are activated may contribute to their differential expressions. Overall, the results comparing C. gigas gigas and C. gigas angulata suggest that these isoforms have a regulatory relationship under heat shock, providing valuable information on the thermal tolerance mechanism in these commercially important oyster species.

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.

A functional polymorphism in the promoter of RhoB is associated with susceptibility to Vibrio anguillarum in turbot (Scophthalmus maximus)

As an isoform of Rho family GTPases, RhoB plays a pivotal role in cytoskeletal organization, cell proliferation, apoptosis and immune response. However, the regulatory mechanisms of RhoB expression in aquatic animals are still unknown. In the present study, we first construct Vibrio anguillarum infection model in S. maximus, including susceptible and resistant individuals. Then the temporal expression of RhoB was detected after V. anguillarum challenge using qRT-PCR and found that RhoB transcripts were significantly induced in the liver, gill and blood despite of differential expression levels and responsive time points. In addition, the mRNA levels of RhoB in resistant individuals were significantly higher than in susceptible ones. The length of 2083 bp sequences of RhoB promoter was cloned and characterized. Moreover, DNA methylation of the RhoB promoter was measured by bisulfite sequencing (BSP) and hypo-methylated was detected in the CpG islands. Three SNPs (-1590, -1575 and -1449) and two haplotypes in the promoter region of RhoB were identified to be associated with V. anguillarum resistance in turbot by association analysis in group 17-R and 17-S. Deletion analysis indicated that these SNPs could negatively mediate the activity of RhoB promoter. Site-directed mutagenesis and qRT-PCR of individuals with different genotypes demonstrated that -1575 T/A polymorphism affected promoter activity. Further study showed that this mutation altered the binding site of the transcription factor CREB. Co-transfection of SmCREB and RhoB promoter was performed in HEK293T cells which confirmed the -1575 allelic differences on transcriptional activity, with the susceptibility allele showing reduced activity. Taken together, our findings implicate that losing of binding of CREB to SmRhoB promoter due to -1575T/A polymorphisms enhances SmRhoB expression in resistant turbot, which provide insights into the effect of SmRhoB expression in response to V. anguillarum infection.

Signatures of selection for bonamiosis resistance in European flat oyster (Ostrea edulis): New genomic tools for breeding programs and management of natural resources

The European flat oyster (Ostrea edulis) is a highly appreciated mollusk with an important aquaculture production throughout the 20th century, in addition to playing an important role on coastal ecosystems. Overexploitation of natural beds, habitat degradation, introduction of non-native species, and epidemic outbreaks have severely affected this important resource, particularly, the protozoan parasite Bonamia ostreae, which is the main concern affecting its production and conservation. In order to identify genomic regions and markers potentially associated with bonamiosis resistance, six oyster beds distributed throughout the European Atlantic coast were sampled. Three of them have been exposed to this parasite since the early 1980s and showed some degree of innate resistance (long-term affected group, LTA), while the other three were free of B. ostreae at least until sampling date (naïve group, NV). A total of 14,065 SNPs were analyzed, including 37 markers from candidate genes and 14,028 from a medium-density SNP array. Gene diversity was similar between LTA and NV groups suggesting no genetic erosion due to long-term exposure to the parasite, and three population clusters were detected using the whole dataset. Tests for divergent selection between NV and LTA groups detected the presence of a very consistent set of 22 markers, located within a putative single genomic region, which suggests the presence of a major quantitative trait locus associated with B. ostreae resistance. Moreover, 324 outlier loci associated with factors other than bonamiosis were identified allowing fully discrimination of all the oyster beds. A practical tool which included the 84 highest discriminative markers for tracing O. edulis populations was developed and tested with empirical data. Results reported herein could assist the production of stocks with improved resistance to bonamiosis and facilitate the management of oyster beds for recovery production and ecosystem services provided by this species.

Pathogen proteases and host protease inhibitors in molluscan infectious diseases

The development of infectious diseases represents an outcome of dynamic interactions between the disease-producing agent's pathogenicity and the host's self-defense mechanism. Proteases secreted by pathogenic microorganisms and protease inhibitors produced by host species play an important role in the process. This review aimed at summarizing major findings in research on pathogen proteases and host protease inhibitors that had been proposed to be related to the development of mollusk diseases. Metalloproteases and serine proteases respectively belonging to Family M4 and Family S8 of the MEROPS system are among the most studied proteases that may function as virulence factors in mollusk pathogens. On the other hand, a mollusk-specific family (Family I84) of novel serine protease inhibitors and homologues of the tissue inhibitor of metalloprotease have been studied for their potential in the molluscan host defense. In addition, research at the genomic and transcriptomic levels showed that more proteases of pathogens and protease inhibitor of hosts are likely involved in mollusk disease processes. Therefore, the pathological significance of interactions between pathogen proteases and host protease inhibitors in the development of molluscan infectious diseases deserves more research efforts.

Testing plasma subtilisin inhibitory activity as a selective marker for dermo resistance in eastern oysters

Recent findings have suggested that eastern oyster plasma possesses inhibitors of the protease subtilisin, which play a role in the host defense against Perkinsus marinus, a protist parasite causing dermo. A study was conducted to determine whether plasma subtilisin inhibitory activity (PSIA) could be used as a selective marker in breeding programs for dermo resistance. Eastern oysters Crassostrea virginica from 2 wild Louisiana populations shown to differ in dermo resistance were collected and their PSIA was measured. Three groups of oysters were established to spawn from each population. One group was composed of randomly sampled oysters (i.e. unselected) and the other 2 groups were composed of oysters with the highest or lowest PSIA. After spawning, progenies were deployed in October 2014 in a dermo endemic area and sampled quarterly for 2 yr to measure their mortality, growth, P. marinus infection intensity, condition index, PSIA, and the gene expression of 3 subtilisin inhibitors (cvSI-1, cvSI-2, and cvSI-3). Oyster cumulative mortalities of the progenies of all groups increased both years from April to October, concomitant with increasing P. marinus infection intensities. Mortalities and P. marinus infection intensities differed markedly between the 2 populations, but differences between the unselected and selected groups of each population were limited. Measurements of PSIA and cvSI-1, cvSI-2, and cvSI-3 gene expressions between the progenies of all groups showed few differences. CvSI-1 gene expression in surviving oysters of the most susceptible population was increased at the end of the study, adding additional support to the potential role of cvSI-1 defense against P. marinus.

Long-term affected flat oyster (Ostrea edulis) haemocytes show differential gene expression profiles from naïve oysters in response to Bonamia ostreae

European flat oyster (Ostrea edulis) production has suffered a severe decline due to bonamiosis. The responsible parasite enters in oyster haemocytes, causing an acute inflammatory response frequently leading to death. We used an immune-enriched oligo-microarray to understand the haemocyte response to Bonamia ostreae by comparing expression profiles between naïve (NS) and long-term affected (AS) populations along a time series (1 d, 30 d, 90 d). AS showed a much higher response just after challenge, which might be indicative of selection for resistance. No regulated genes were detected at 30 d in both populations while a notable reactivation was observed at 90 d, suggesting parasite latency during infection. Genes related to extracellular matrix and protease inhibitors, up-regulated in AS, and those related to histones, down-regulated in NS, might play an important role along the infection. Twenty-four candidate genes related to resistance should be further validated for selection programs aimed to control bonamiosis.

Two I84 family protease inhibitors from Chinese razor clams Sinonovacula constricta expressed in response to environmental challenges

Protease inhibitors play critical roles in numerous biological processes including host defense in all multicellular organisms. Eighty three evolutionary families of protease inhibitors are currently accommodated in the MEROPS database and the I84 family currently consists of 3 novel serine protease inhibitors from the eastern oyster Crassostrea virginica. In this study, we identified 2 new I84 family members from the Chinese razor clam Sinonovacula constricta, scSI-1 and scSI-2, using cDNA cloning and sequencing. The scSI-1 cDNA consisted of 494 bp with a 282 bp ORF encoding a 93-amino acid polypeptide that was predicted to have a 19-amino acid signal peptide and a 74-residue mature protein with a calculated molecular mass of 8248.5 Da. The scSI-2 cDNA was 490 bp long with a 273 bp ORF encoding a 90-amino acid polypeptide that was predicted to have an 18-amino acid signal peptide and a 72-residue nature protein with a calculated molecular mass of 7528.4 Da. ScSI-1 and scSI-2 shared high sequence similarity with the 3 known members of I84 family and both expressed primarily in the clam digestive glands. Protease inhibitory activity in the clam plasma also exhibited the signature kinetic characteristics of the I84 members from the oyster. In addition, levels of scSI-1 and scSI-2 gene expression in digestive glands and the protease inhibitory activity in plasma elevated significantly in clams challenged by bacterial injections and Vibrio harveyi was more effective than Staphylococcus epidermidis in inducing the gene expression and plasma protease inhibitory activity. Moreover, drastic changes of salinity and temperature also caused significant changes in the gene expression and plasma activity. These results indicated that scSI-1 and scSI-2 represented 2 new members of the I84 family and they likely play a role in clam host defense against infections and in reactions against physiochemical stressors.

The polymorphism of chicken-type lysozyme gene in Japanese flounder (Paralichthys olivaceus) and its association with resistance/susceptibility to Listonella anguillarum

Lysozyme is a crucially spread hydrolase in organisms that can defend against bacterial infection in innate immunity. In this study, we successfully sequenced the coding region of chicken-type lysozyme gene (PoLysC) in Paralichthys olivaceus and identified nine single nucleotide polymorphisms (SNPs). We then amplified the 2500 bp promoter region of lysozyme and identified the eight sites of polymorphisms. All SNPs were genotyped between susceptible and resistance groups after Listonella anguillarum challenge. One of these SNP sites in the codon of PoLysC was genotyped and determined to be a significant marker by analyzing its distribution in the susceptible and resistant groups. As a nonsynonymous mutation, the frequency of 140G/C genotype in the resistant group was higher (67.74%) than that in the susceptible group (32.26%). The linkage between SNP140 and polymorphisms in the promoter region was also studied. Results revealed that the frequency of haplotype CC_-536/CC_-1200/GG₁₄₀ in the resistance group was significantly higher than that in the susceptible group. The quantitative expression of lysozyme gene in the resistant group was also higher than that in the susceptible group. This finding indicated that the linkage between polymorphism -536 and -1200 sites in promoter and SNP140 in codon sequence was associated with the resistance of P. olivaceus to L. anguillarum. All these results suggest that the mutations in promoter and coding region were related to changes in PoLysC for resisting L. anguillarum. The haplotype CC_-536/CC_-1200/GG₁₄₀ was a potential marker and can thus be applied to selective breeding for the disease resistance of P. olivaceus.

Identification of cvSI-3 and evidence for the wide distribution and active evolution of the I84 family of protease inhibitors in mollusks

Protease inhibitors are an extremely diverse group of proteins that control the proteolytic activities of proteases and play a crucial role in biological processes including host defenses. The I84 family of protease inhibitors in the MEROPS database currently consists of cvSI-1 and cvSI-2, two novel serine protease inhibitors purified and characterized from the eastern oyster Crassostrea virginica plasma and believed to play a role in host defense and disease resistance. In the present study, a third member of I84 family, named cvSI-3, was identified from C. virginica by cDNA cloning and sequencing. The full cvSI-3 cDNA was composed of 342 bp including a 255 bp open reading frame (ORF) that encodes an 84-amino acid peptide. The mature cvSI-3 molecule was predicted to have 68 amino acid residues after removal of a 16-amino acid signal peptide, with a calculated molecular mass of 7724.5 Da and a theoretical isoelectric point (pI) of 6.28. CvSI-3 amino acid sequence shared 41% identity with cvSI-2 and 37% identity with cvSI-1, which included 12 conserved cysteines. Quantitative real-time PCR determined that cvSI-3 gene expressed primarily in oyster digestive glands. Real-time PCR also detected that cvSI-1, cvSI-2 and cvSI-3 expression levels in digestive glands varied significantly, with cvSI-2 showing the highest expression level and cvSI-3 the lowest. Additionally, a significant correlation was detected between cvSI-2 and cvSI-3 mRNAs levels. Searches into sequence databases using cvSI-1, cvSI-2 and cvSI-3 as queries retrieved ESTs suggesting the possible existence of at least 9 more I84 family members in eastern oysters and of I84 family protease inhibitors in various bivalve and gastropod species. Moreover, orthologs of all C. virginica I84 family members or potential member genes were found to be present in the C. gigas genome, and their distributions among species provided important information about the evolution of the I84 family of protease inhibitors. It appears that the I84 family of protease inhibitors is widely distributed and actively evolving in the Phylum Mollusca.

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.

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.

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.

Clam focal and systemic immune responses to QPX infection revealed by RNA-seq technology

Background

The hard clam Mercenaria mercenaria is an important seafood species widely exploited along the eastern coasts of the United States and play a crucial role in coastal ecology and economy. Severe hard clam mortalities have been associated with the protistan parasite QPX (Quahog Parasite Unknown). QPX infection establishes in pallial organs with the lesions typically characterized as nodules, which represent inflammatory masses formed by hemocyte infiltration and encapsulation of parasites. QPX infection is known to induce host changes on both the whole-organism level and at specific lesion areas, which imply systemic and focal defense responses, respectively. However, little is known about the molecular mechanisms underlying these alterations.

Results

RNA-seq was performed using Illumina Hiseq 2000 (641 Million 100 bp reads) to characterize M. mercenaria focal and systemic immune responses to QPX. Transcripts were assembled and the expression levels were compared between nodule and healthy tissues from infected clams, and between these and tissues from healthy clams. De novo assembly reconstructed a consensus transcriptome of 62,980 sequences that was functionally-annotated. A total of 3,131 transcripts were identified as differentially expressed in different tissues. Results allowed the identification of host immune factors implicated in the systemic and focal responses against QPX and unraveled the pathways involved in parasite neutralization. Among transcripts significantly modulated upon host-pathogen interactions, those involved in non-self recognition, signal transduction and defense response were over-represented. Alterations in pathways regulating hemocyte focal adhesion, migration and apoptosis were also demonstrated.

Conclusions

Our study is the first attempt to thoroughly characterize M. mercenaria transcriptome and identify molecular features associated with QPX infection. It is also one of the first studies contrasting focal and systemic responses to infections in invertebrates using high-throughput sequencing. Results identified the molecular signatures of clam systemic and focal defense responses, to collectively mediate immune processes such as hemocyte recruitment and local inflammation. These investigations improve our understanding of bivalve immunity and provide molecular targets for probing the biological bases of clam resistance towards QPX.

Electronic supplementary material

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

Molecular Basis for Adaptation of Oysters to Stressful Marine Intertidal Environments

Oysters that occupy estuarine and intertidal habitats have well-developed stress tolerance mechanisms to tolerate harsh and dynamically changing environments. In this review, we summarize common pathways and genomic features in oyster that are responsive to environmental stressors such as temperature, salinity, hypoxia, air exposure, pathogens, and anthropogenic pollutions. We first introduce the key genes involved in several pathways, which constitute the molecular basis for adaptation to stress. We use genome analysis to highlight the strong cellular homeostasis system, a unique adaptive characteristic of oysters. Next, we provide a global view of features of the oyster genome that contribute to stress adaptation, including oyster-specific gene expansion, highly inducible expression, and functional divergence. Finally, we review the consequences of interactions between oysters and the environment from ecological and evolutionary perspectives by discussing mass mortality and adaptive divergence among populations and related species of the genus Crassostrea. We conclude with prospects for future study.

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.

Three novel clade B serine protease inhibitors from disk abalone, Haliotis discus discus: Molecular perspectives and responses to immune challenges and tissue injury

Serine protease inhibitors (SERPINs) control cellular protease activity in order to maintain cellular homeostasis. The immune and inflammatory responses of invertebrate clade B SERPINs have not been widely reported. In the present study, three proteins with high similarity to clade B SERPINs, referred to as AbSERPIN-1, AbSERPIN-2 and AbSERPIN-3, were identified from disk abalone (Haliotis discus discus). While AbSERPIN-1 (399 aa) was of a typical size for this protein class, AbSERPIN-2 (506 aa) and AbSERPIN-3 (532 aa) were relatively larger. Bioinformatic analysis revealed the characteristic SERPIN domain in each AbSERPIN. In addition, the N-terminal region of both AbSERPIN-2 and AbSERPIN-3 contained a predicted low complexity region (LCR) and a signal peptide, suggesting that these proteins are secretory proteins and are, thus, novel peptides. Tertiary structural models of the AbSERPINs highlighted their structural and functional conservation. Ubiquitous expression of AbSERPIN transcripts was evaluated by quantitative real time PCR (qPCR) analysis in seven tissue types. AbSERPIN-1, AbSERPIN-2, and AbSERPIN-3 transcript levels were highest in mantle, hemocytes, and muscles, respectively. Temporal expression analysis revealed that AbSERPINs were significantly (P < 0.05) elevated in hemocytes during the early/middle stages following the injection of a bacterial pathogen (Vibrio parahaemolyticus or Listeria monocytogenes) or an immuno-stimulant (lipopolysaccharide). Moreover, mantle tissue injury led to significant changes in the temporal expression of AbSERPIN mRNA. Specifically, transcription of AbSERPIN-1 and AbSERPIN-3 was considerably up-regulated, while expression of AbSERPIN-2 was suppressed. These results suggest a potential role of AbSERPINs in response to pathogen invasion and tissue injury in disk abalone.

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.

Differential expression of serine protease inhibitors 1 and 2 in Crassostrea corteziensis and C. virginica infected with Perkinsus marinus

Proliferation of Perkinsus marinus (Dermo) in vitro is inhibited by the action of 2 serine protease inhibitors belonging to the I-84 family. We compared the levels of expression of serine protease inhibitors 1 and 2 (SPI-1 and SPI-2) in 2 oyster species (Crassostrea virginica and C. corteziensis) inoculated with the parasite P. marinus. C. virginica is well known to be susceptible to this parasite, whereas C. corteziensis is apparently more tolerant. Oysters were inoculated with trophozoites (1 × 106 trophozoites oyster-1) of P. marinus while control oysters were injected with saline solution. Oysters were maintained in a closed water system for 2 wk. The oysters were then sacrificed and parasite burden, histological damage, and gene expression were evaluated. The results showed that the challenged oysters presented a significant increase in parasite burden, which generated histological damage in digestive gland and gills. Quantitative PCR detected significant differences in SPI-1 and SPI-2 expression levels in the 2 oyster species, with C. corteziensis showing higher expression levels than C. virginica as a response to P. marinus inoculation. Our results provide valuable information for the understanding of the defense response in C. corteziensis and a possible explanation for its tolerance to the parasite.

Infectious diseases affect marine fisheries and aquaculture economics

Seafood is a growing part of the economy, but its economic value is diminished by marine diseases. Infectious diseases are common in the ocean, and here we tabulate 67 examples that can reduce commercial species' growth and survivorship or decrease seafood quality. These impacts seem most problematic in the stressful and crowded conditions of aquaculture, which increasingly dominates seafood production as wild fishery production plateaus. For instance, marine diseases of farmed oysters, shrimp, abalone, and various fishes, particularly Atlantic salmon, cost billions of dollars each year. In comparison, it is often difficult to accurately estimate disease impacts on wild populations, especially those of pelagic and subtidal species. Farmed species often receive infectious diseases from wild species and can, in turn, export infectious agents to wild species. However, the impact of disease export on wild fisheries is controversial because there are few quantitative data demonstrating that wild species near farms suffer more from infectious diseases than those in other areas. The movement of exotic infectious agents to new areas continues to be the greatest concern.

Polymorphisms of heat shock protein 90 (Hsp90) in the sea cucumber Apostichopus japonicus and their association with heat-resistance

Heat shock protein 90 (Hsp90) functions as a molecular chaperone and plays an important role in the resistance of organisms to stress, particularly heat-stress. In our study, 12 exons and 11 introns of hsp90 were identified in the sea cucumber Apostichopus japonicus. Twenty-two single nucleotide polymorphisms (SNPs), including three non-synonymous mutations, were detected in the exons. Susceptible and resistant individuals were distinguished using a high-temperature (32 °C) challenge experiment. Three blocks with high linkage disequilibrium were detected among these SNPs. Five of the twenty-two SNPs were shown to be significantly associated with susceptibility/resistance to high temperature by correlation analysis (chi-square test, P < 0.05). To confirm the importance of these five SNPs, a heat-resistance strain (HRS) was selected through three generations. Using the common population as the control group, it was shown that the distributions of genotypes and alleles of SNP e10-1 and e11-6 were significantly different between the two groups (P < 0.05). SNP e10-1 was trimorphic, with three alleles (A, C and T) and five genotypes (AA, CC, AT, CT and AC). The allele frequency of SNP e2-3 was also significantly associated with this trait (P < 0.05). This is the first demonstration of SNPs related to heat-resistance in A. japonicus and supports the use of SNP markers in the selective breeding of sea cucumbers.

Identification of MMV Malaria Box inhibitors of Perkinsus marinus using an ATP-based bioluminescence assay

"Dermo" disease caused by the protozoan parasite Perkinsus marinus (Perkinsozoa) is one of the main obstacles to the restoration of oyster populations in the USA. Perkinsus spp. are also a concern worldwide because there are limited approaches to intervention against the disease. Based on the phylogenetic affinity between the Perkinsozoa and Apicomplexa, we exposed Perkinsus trophozoites to the Medicines for Malaria Venture Malaria Box, an open access compound library comprised of 200 drug-like and 200 probe-like compounds that are highly active against the erythrocyte stage of Plasmodium falciparum. Using a final concentration of 20 µM, we found that 4 days after exposure 46% of the compounds were active against P. marinus trophozoites. Six compounds with IC50 in the µM range were used to compare the degree of susceptibility in vitro of eight P. marinus strains from the USA and five Perkinsus species from around the world. The three compounds, MMV666021, MMV665807 and MMV666102, displayed a uniform effect across Perkinsus strains and species. Both Perkinsus marinus isolates and Perkinsus spp. presented different patterns of response to the panel of compounds tested, supporting the concept of strain/species variability. Here, we expanded the range of compounds available for inhibiting Perkinsus proliferation in vitro and characterized Perkinsus phenotypes based on their resistance to six compounds. We also discuss the implications of these findings in the context of oyster management. The Perkinsus system offers the potential for investigating the mechanism of action of the compounds of interest.

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

As a benthic filter-feeder of estuaries, the eastern oyster, Crassostrea virginica, faces tremendous exposure to microbial pathogens. How eastern oysters without adaptive immunity survive in pathogen-rich environments is of fundamental interest, but studies on its immune system are hindered by the lack of genomic resources. We sequenced the transcriptome of an adult oyster with short Illumina reads and assembled 66,229 contigs with a N50 length of 1,503 bp. The assembly covered 89.4 % of published ESTs and 97.9 % of mitochondrial genes demonstrating its quality. A set of 39,978 contigs and unigenes (>300 bp) were identified and annotated by searching public databases. Analysis of the gene set yielded a diverse set of 657 genes related to innate immunity, including many pertaining to pattern recognition, effectors, signal transduction, cytokines, and apoptosis. Gene families encoding C1q domain containing proteins, CTLD, IAPs, Ig_I-set, and TRAFs expanded in C. virginica and Crassostrea gigas. Many key genes of the apoptosis system including IAP, BAX, BAC-2, caspase, FADD, and TNFR were identified, suggesting C. virginica possess advanced apoptosis and apoptosis-regulating systems. Our results show that short Illumina reads can produce transcriptomes of highly polymorphic genomes with coverage and integrity comparable to that from longer 454 reads. The expansion and high diversity in gene families related to innate immunity, point to a complex defense system in the lophotrochozoan C. virginica, probably in adaptation to a pathogen-rich environment.

Immunological responses of the mangrove oysters Crassostrea gasar naturally infected by Perkinsus sp. in the Mamanguape Estuary, Paraíba state (Northeastern, Brazil)

Perkinsus genus includes protozoan parasites of marine mollusks, especially bivalves. In the last four years, this parasite has been detected in mangrove oysters Crassostrea rhizophorae and Crassostrea gasar from the Northeastern region of Brazil. Hemocytes are the key cells of the oyster immune system, being responsible for a variety of cellular and humoral reactions, such as phagocytosis, encapsulation and the release of several effector molecules that control the invasion and proliferation of microorganisms. In Brazil, there is little information on perkinsosis and none on the immune responses of native oysters' species against Perkinsus spp. The objective of this study was to determine the effects of natural infection by Perkinsus sp. on the immunological parameters of mangrove oysters C. gasar cultured in the Mamanguape River Estuary (Paraíba, Brazil). Adults oysters (N = 40/month) were sampled in December 2011, March, May, August and October 2012. Gills were removed and used to determine the presence and intensity of the Perkinsus sp. infection, according to a scale of four levels (1-4), using the Ray's fluid thioglycollate medium assay. Immunological parameters were measured in hemolymph samples by flow cytometry, including: total hemocyte count (THC), differential hemocyte count (DHC), cell mortality, phagocytic capacity, and production of Reactive Oxygen Species (ROS). The plasma was used to determine the hemagglutination activity. The results showed the occurrence of Perkinsus sp. with the highest mean prevalence (93.3%) seen so far in oyster populations in Brazil. Despite that, no oyster mortality was associated. In contrast, we observed an increase in hemocyte mortality and a suppression of two of the main defense mechanisms, phagocytosis and ROS production in infected oysters. The increase in the percentage of blast-like cells on the hemolymph, and the increase in THC in oysters heavily infected (at the maximum intensity, 4) suggest an induction of hemocytes proliferation. The immunological parameters varied over the studied months, which may be attributed to the dynamics of infection by Perkinsus sp. The results of the present study demonstrate that Perkinsus sp. has a deleterious effect on C. gasar immune system, mainly in high intensities, which likely renders oysters more susceptible to other pathogens and diseases.

Polymorphisms of anti-lipopolysaccharide factors in the swimming crab Portunus trituberculatus and their association with resistance/susceptibility to Vibrio alginolyticus

Anti-lipopolysaccharide factor (ALF) is an important antimicrobial peptide (AMP) that can bind and neutralize major component of Gram-negative bacteria cell wall, lipopolysaccharide (LPS). Seven isoforms of anti-lipopolysaccharide factors (PtALF1-7) were previously identified from the swimming crab Portunus trituberculatus in our laboratory. Here, polymorphisms of PtALF1-7 were detected and their association with resistance/susceptibility to Vibrio alginolyticus (a main Gram-negative bacteria causing high mortality in P. trituberculatus) were investigated. We identified 127, 96, 103, 53 and 158 single nucleotide polymorphisms (SNPs) in genomic fragments of PtALF1-3, PtALF4, PtALF5, PtALF6 and PtALF7, respectively. Among them, totally sixteen SNPs were significantly associated with resistance/susceptibility to V. alginolyticus (P < 0.05). Of these sixteen SNPs, most were located in introns and noncoding exons, while two synonymous SNPs and one nonsynonymous SNP were in coding exons. Additionally, simple sequence repeats (SSRs) were only identified in introns and noncoding exons of PtALF4, PtALF5 and PtALF7. Although no significant difference of allele frequencies was found, these SSRs had different polymorphic alleles according to the repeat number between susceptible and resistant stocks. After further confirmation, polymorphisms investigated here might be applied as potential molecular markers for future selection of resistant strains to diseases caused by Gram-negative bacteria.

A Review of Current State of Knowledge Concerning Perkinsus marinus Effects on Crassostrea virginica (Gmelin) (the Eastern Oyster)

The eastern oyster, Crassostrea virginica (Gmelin), is both an important component of our estuaries and an important farmed food animal along the east and south coasts of the United States. Its populations have been significantly diminished in the wild due to decades of overfishing beginning in the 1890s. Unfortunately, in 1950, a new disease in eastern oysters caused by the protistan agent, Perkinsus marinus, was identified. The disease, resulting from infection with this protozoan, leads to high mortality of both wild and cultured eastern oysters. Current restoration efforts are hampered by the disease, as is the aquaculture of this economically important food. The parasite infects hemocytes and causes hemolytic anemia and general degeneration of the tissues, leading to death. Ongoing research efforts are attempting to develop oysters resistant to the disease. Transport regulations exist in may states. Infection with P. marinus is listed as a reportable disease by the World Health Organization.