Single nucleotide polymorphisms in i-type lysozyme gene and their correlation with vibrio-resistance and growth of clam Meretrix meretrix based on the selected resistance stocks

LysM-containing proteins function in the resistance of Litopenaeus vannamei against Vibrio parahaemolyticus infection

Lysin motif (LysM) is a functional domain that can bind to peptidoglycans, chitin and their derivatives. The LysM-containing proteins participate in multiple biological processes, such as the hydrolysis of bacterial cell walls and the perception of PAMPs in plants and high animals. In the present study, two genes encoding LysM-containing proteins, designated as LvLysM1 and LvLysM2, were identified in the Pacific white shrimp, Litopenaeus vannamei, and their functions during Vibrio infection were analyzed. The open-reading frame (ORF) of LvLysM1 was 795 bp, only encoding a LysM domain at the N-terminal region. The ORF of LvLysM2 was 834 bp, encoding a LysM domain at the central region and a transmembrane region at the C-terminal region. Both LvLysM1 and LvLysM2 were widely transcribed in all tested shrimp tissues. Enzyme-linked immunosorbent assay (ELISA) showed that the recombinant protein of LvLysM2 could bind to different bacterial polysaccharides, while LvLysM1 showed no direct binding activity. The transcripts of LvLysMs in gills increased significantly after infection with Vibrio parahaemolyticus. When LvLysM1 or LvLysM2 was knocked down by dsRNA, the mortality of shrimp was significantly increased after infection with Vibrio parahaemolyticus. Interestingly, some SNPs existed in these two genes were apparently correlated with the VpAHPND resistance of shrimp. These results suggested that LvLysM1 and LvLysM2 might contribute to the disease resistance of shrimp. The data provide new knowledge about the function of LysM-containing proteins in shrimp and potential genetic markers for disease resistance breeding.

Characterization of a novel activating protein-1 (AP-1) gene and the association of its single nucleotide polymorphisms with vibrio resistance in Tegillarca granosa

The blood clam Tegillarca granosa is a commercial marine bivalve of economic value, accounting for approximately 50% of clam production in China. In recent years, the yield of blood clams has been threatened by bacterial infections caused by marine Vibrio species that thrive under a rising sea temperature. The transcription factor activating protein-1 (AP-1) is emerging as an important player in the innate immunity of marine bivalves against viral or bacterial infections. In this study, the full-length cDNA of a novel T. granosa AP-1 (TgAP-1) was cloned for the first time. The 1591-bp cDNA encoded a protein of 292 amino acid residues with a calculated molecular weight of 32.8 kDa. The TgAP-1 protein contained an N-terminal Jun domain and a C-terminal basic region leucine zipper domain typically found in Jun proteins (a subfamily of AP-1 proteins). TgAP-1 was ubiquitously expressed in T. granosa, with the highest expression detected in the gill and foot, followed by the mantle, hemolymph, and hepatopancreas. Exposure to Vibrio harveyi induced TgAP-1 expression in gill tissues and the expression levels of TgAP-1 of resistant blood clams were always lower than that of control population whether Vibro infection or not. A total of 18 single nucleotide polymorphisms (SNPs) of TgAP-1 were detected in T. granosa. SNP-typing and haplotyping of resistant and susceptible populations revealed that six SNPs (AG type of TgSNP-1, GA type of TgSNP-2, TG type of TgSNP-4, CT type of TgSNP-7, AG type of TgSNP-11, and GA type of TgSNP-12) and four haplotypes (fHap2, fHap3, fHap6, and fHap7) were significantly associated with V. harveyi resistance. Risk assessment showed that fHap2 (CG) and fHap7 (GA) were associated with an increased resistance, while fHap3 (CT) and fHap6 (AG) were associated with an increased susceptibility. The results from this study supported a potential role of TgAp-1 in the anti-Vibro immunity of T. granosa. The discovery of the genetic molecular markers and haplotypes related to Vibrio resistance can provide guidance for selective breeding of T. granosa in the future.

Identification and Expression Characterization of the Smad3 Gene and SNPs Associated with Growth Traits in the Hard Clam (Meretrix meretrix)

It has been demonstrated that the sekelsky mothers against decapentaplegic homolog 3 (Smad3) plays an important role in the growth and development of vertebrates. However, little is known about the association between the Smad3 gene and the growth traits of mollusks. In this study, Smad3 from the hard clam Meretrix meretrix (Mm-Smad3) was cloned, characterized, and screened for growth-related single nucleotide polymorphisms (SNPs) in its exons. The full-length cDNA of Mm-Smad3 was 1938 bp, encoding a protein with 428 amino acid residues. The protein sequence included an MH1 (27–135 aa) and MH2 domain (233–404 aa). Promoter analysis showed that the promoter sequence of Mm-Smad3 was 2548 bp, and the binding sites of Pit-1a, Antp, Hb, and other transcription factors are related to the growth and development of hard clams. The phylogenetic tree was divided into two major clusters, including mollusks and vertebrate. The expression level of Mm-Smad3 was predominantly detected in the mantle and foot, while extremely less expression was observed in the digestive gland. The low expression level of Mm-Smad3 was detected at the stages of unfertilized mature eggs, fertilized eggs, four-cell embryos, blastula, gastrulae, trochophore, and D-shaped larvae, whereas an opposite trend was observed regarding the highest expression at the umbo larvae stage (p < 0.05). In the mantle repair experiment, the time-course expression profiles showed that compared to the expression level at 0 h, Mm-Smad3 significantly decreased at 6 h (p < 0.05) but increased at 12 and 48 h. Further, the association analysis identified 11 SNPs in the exons of Mm-Smad3, of which three loci (c.597 C > T, c.660 C > T, c.792 A > T) were significantly related to the growth traits of clam (p < 0.05). Overall, our findings indicated that Mm-Smad3 is a growth-related gene and the detected SNP sites provide growth-related markers for molecular marker-assisted breeding of this species.

Potential for Genetic Improvement of Resistance to Perkinsus olseni in the Manila Clam, Ruditapes philippinarum, Using DNA Parentage Assignment and Mass Spawning

The Manila clam Ruditapes philippinarum, a major cultured shellfish species, is threatened by infection with the microparasite Perkinsus olseni, whose prevalence increases with high water temperatures. Under the current trend of climate change, the already severe effects of this parasitic infection might rapidly increase the frequency of mass mortality events. Treating infectious diseases in bivalves is notoriously problematic, therefore selective breeding for resistance represents a key strategy for mitigating the negative impact of pathogens. A crucial step in initiating selective breeding is the estimation of genetic parameters for traits of interest, which relies on the ability to record parentage and accurate phenotypes in a large number of individuals. Here, to estimate the heritability of resistance against P. olseni, a field experiment mirroring conditions in industrial clam production was set up, a genomic tool was developed for parentage assignment, and parasite load was determined through quantitative PCR. A mixed-family cohort of potentially 1,479 clam families was produced in a hatchery by mass spawning of 53 dams and 57 sires. The progenies were seeded in a commercial clam production area in the Venice lagoon, Italy, where high prevalence of P. olseni had previously been reported. Growth and parasite load were monitored every month and, after 1 year, more than 1,000 individuals were collected for DNA samples and phenotype recording. A pooled sequencing approach was carried out using DNA samples from the hatchery broodstock and from a Venice lagoon clam population, providing candidate markers used to develop a 245-SNP panel. Parentage assignment for 246 F1 individuals showed sire and dam representation were high (75 and 85%, respectively), indicating a very limited risk of inbreeding. Moderate heritability (0.23 ± 0.11-0.35 ± 0.13) was estimated for growth traits (shell length, shell weight, total weight), while parasite load showed high heritability, estimated at 0.51 ± 0.20. No significant genetic correlations were found between growth-associated traits and parasite load. Overall, the preliminary results provided by this study show high potential for selecting clams resistant to parasite load. Breeding for resistance may help limit the negative effects of climate change on clam production, as the prevalence of the parasite is predicted to increase under a future scenario of higher temperatures. Finally, the limited genetic correlation between resistance and growth suggests that breeding programs could incorporate dual selection without negative interactions.

Identification of Single Nucleotide Polymorphisms Related to the Resistance Against Acute Hepatopancreatic Necrosis Disease in the Pacific White Shrimp Litopenaeus vannamei by Target Sequencing Approach

Acute hepatopancreatic necrosis disease (AHPND) is a major bacterial disease in Pacific white shrimp Litopenaeus vannamei farming, which is caused by Vibrio parahaemolyticus. AHPND has led to a significant reduction of shrimp output since its outbreak. Selective breeding of disease-resistant broodstock is regarded as a key strategy in solving the disease problem. Understanding the relationship between genetic variance and AHPND resistance is the basis for marker-assisted selection in shrimp. The purpose of this study was to identify single nucleotide polymorphisms (SNPs) associated with the resistance against AHPND in L. vannamei. In this work, two independent populations were used for V. parahaemolyticus challenge and the resistant or susceptible shrimp were evaluated according to the survival time after Vibrio infection. The above two populations were genotyped separately by a SNP panel designed based on the target sequencing platform using a pooling strategy. The SNP panel contained 508 amplicons from DNA fragments distributed evenly along the genome and some immune-related genes of L. vannamei. By analyzing the allele frequency in the resistant and susceptible groups, 30 SNPs were found to be significantly associated with the resistance of the shrimp against V. parahaemolyticus infection (false discovery rate corrected at P < 0.05). Three SNPs were further validated by individual genotyping in all samples of population 1. Our study illustrated that target sequencing and pooling sequencing were effective in identifying the markers associated with economic traits, and the SNPs identified in this study could be used as molecular markers for breeding disease-resistant shrimp.

Identification and immune-related analysis of SNPs in Litopenaeus vannamei Toll3 receptor

Tolls and Toll-like receptors (TLRs), as innate immune-recognition receptors that recognize molecular patterns associated with microbial pathogens, play a critical role in antimicrobial immune responses. Here, we report on single nucleotide polymorphisms (SNPs) of Litopenaeus vannamei Toll3 (LvToll3). Multiple sequence alignment of the L. vannamei Toll3 Leucine rich repeat C-terminal domain (LvToll3-LRR-CT) with other L. vannamei Tolls LRR-CT domains showed 39.23% - 43.96% homology at the nucleic acid level and 20.31% - 30.00% identity at the amino acid level. Analysis of different shrimp tissues by polymerase chain reaction denaturing gradient gel electrophoresis (PCR-DGGE) revealed that LvToll3-LRR-CT had genetic polymorphisms at both the genomic deoxyribonucleic acid (gDNA) and complementary deoxyribonucleic acid (cDNA) levels. Further, high-throughput sequencing analysis confirmed the presence of 8 non-synonymous SNP (nsSNP) and 1 nsSNPs with frequency greater than 1% at the gDNA level, while 13 nsSNPs and 2 nsSNPs with frequency greater than 1% at the cDNA level. In silico analysis revealed that the α-helix secondary structure and tertiary structure of LvToll3 changed when 3 SNPs (C2039T, T2041C, T2228C) were mutated. Interestingly, 2 novel bands on PCR-DGGE, which were identified as 2 nsSNPs (C2140A, T2186A) were observed following challenge with Streptococcus iniae but not with Vibrio parahaemolyticus or White spot syndrome virus (WSSV). Moreover, the secondary and tertiary structures of LvToll3 changed when the nsSNP T2186A was mutated. The present findings therefore provide novel insight into the molecular basis of shrimp innate immune response to pathogens through the generation of specific SNPs.

SNPs in the Toll1 receptor of Litopenaeus vannamei are associated with immune response

Tolls and Toll-like receptors (TLRs) are important regulators in the innate immune system and their genetic variations usually affect the host's susceptibility/resistance to pathogen infections. In this study, we report on the single nucleotide polymorphisms (SNPs) of Toll1 in Litopenaeus vannamei (LvToll1) and how this is associated with immune response. PCR-DGGE analysis revealed genetic polymorphisms in LvToll1 at both the genomic DNA (gDNA) and cDNA levels. Using high-throughput sequencing, 223 SNPs were identified at the gDNA level, of which 145 were non-synonymous SNP (nsSNP), with 3 nsSNPs having frequency over 1%. On the other hand, 60 SNPs were identified at the cDNA level including 38 nsSNPs and 4 nsSNPs with frequency over 1%. Upon challenging shrimps with Streptococcus iniae, Vibrio parahaemolyticus and white spot syndrome virus (WSSV), LvToll1 was shown to generate 6, 4 and 4 novel bands, respectively when analyzed with PCR-DGGE. Sequencing analysis of these bands showed that they contained 6, 4 and 2 nsSNPs, respectively. Moreover, the nsSNP C1526T was detected in S. iniae-resistant but not in susceptible shrimps. Most significantly, the C1526T mutation could shorten the α-helix of the LRR domain and was predicted to affect the function of LvToll1, indicating that SNP C1526T might be associated with shrimp's resistance to pathogen infections. In sum, our findings here reveal that the genetic polymorphisms of Toll receptor are linked with the immune response to pathogen infections in L. vannamei.

Identification of an MITF gene and its polymorphisms associated with the Vibrio resistance trait in the clam Meretrix petechialis

Microphthalmia-associated transcription factor (MITF) regulates the transcription of its target genes by binding to their promoters. In this study, an MITF gene, MpMITF was identified in the clam Meretrix petechialis. The full-length cDNA of MpMITF is 3564 bp with an ORF of 1365 bp. The deduced amino acid sequence consists of a conserved functional structure of bHLH-LZ, which could bind with E-box. The mRNA and protein expression levels of MpMITF were significantly up-regulated 6 h post-Vibrio injection. The mRNA expression of MpMITF increased on day 2 and peaked on day 10 post-Vibrio immersion. Furthermore, MpMITF expression was significantly up-regulated in most resistant families of clams (P < 0.05) but did not change significantly in most susceptive families of clams after the Vibrio immersion challenge. These results suggest that, in clams, MpMITF participates in the immune response against a Vibrio infection. Genotyping in two clam groups with different resistant levels to Vibrio parahaemolyticus (i.e., 11-R and 11-S), thirteen SNPs and five haplotypes were detected in the DNA sequence of MpMITF, of which five SNPs and two haplotypes were associated with Vibrio resistance. Four SNPs (SNP2, 5, 6 and 13) and one haplotype (Hap1) were further confirmed to be associated with Vibrio resistance in M. petechialis by association analysis in different clam families. This study deepens the understanding of MITF in marine bivalves and provides potential candidate markers for resistance selection in the clam M. petechialis.

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.

The polymorphisms of a MIF gene and their association with Vibrio resistance in the clam Meretrix meretrix

Macrophage migration inhibitory factor (MIF) is an important proinflammatory cytokine that mediates both innate and adaptive immune responses. In this study, a homolog of MIF was identified in the clam Meretrix meretrix. Ten SNPs in the DNA partial sequence of MmMIF were found to be significantly associated with Vibrio resistance (P < 0.05). Distinct expression patterns of MmMIF among different haplotypes were observed after Vibrio challenge. The results showed that haplotypes did not affect MmMIF expression in the negative control group, while the expression of MmMIF in clams with Hap1 and Hap1/Hap2 was significantly lower than that with Hap2 at 24 h in the PBS-injected group but significantly higher than that with Hap2 in the Vibrio-injected group. The results also indicate that Hap1 and Hap1/Hap2 can specifically respond to mechanical stimulation while Hap2 can specifically respond to Vibrio infection. The effect of a missense mutation was detected by site-directed mutagenesis using fusion expression of the protein, which showed that the SNP g.737 (I > V) has no effect on redox activity and tautomerase activity. These studies identified a potential marker that is enriched in Vibrio-resistant clams that can be used for the genetic breeding of Meretrix meretrix.

Identification of the MmeHairy gene and expression analysis affected by two SNPs in the 3'-untranslated region in the clam Meretrix meretrix

As a bHLH transcriptional repressor, Hairy-related proteins can bind to DNA sites in target gene promoters and negatively regulate gene transcription. In the present study, the full-length cDNA of Hairy was obtained from the clam Meretrix meretrix (MmeHairy), and two SNPs in the 3'-untranslated region (UTR) of this gene, SNP1066 and 1067, were identified and characterized. Multiple sequence alignment and phylogenetic analysis revealed that MmeHairy belongs to the Hairy protein subfamily. Analysis of tissue expression patterns showed that the mRNA of MmeHairy had the highest expression level in the hepatopancreas. The expression levels of MmeHairy were up-regulated in the hepatopancreas after Vibrio challenge. Genotyping and quantitative analysis showed that the mRNA levels of MmeHairy were significantly different among individual clams with different genotypes at SNP1066 and 1067 (P < 0.05), which indicated that these two SNP loci may affect the expression of MmeHairy and could be used as candidate markers for future selection in M. meretrix breeding programs.

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.

Adaptation to enemy shifts: rapid resistance evolution to local Vibrio spp. in invasive Pacific oysters

One hypothesis for the success of invasive species is reduced pathogen burden, resulting from a release from infections or high immunological fitness of invaders. Despite strong selection exerted on the host, the evolutionary response of invaders to newly acquired pathogens has rarely been considered. The two independent and genetically distinct invasions of the Pacific oyster Crassostrea gigas into the North Sea represent an ideal model system to study fast evolutionary responses of invasive populations. By exposing both invasion sources to ubiquitous and phylogenetically diverse pathogens (Vibrio spp.), we demonstrate that within a few generations hosts adapted to newly encountered pathogen communities. However, local adaptation only became apparent in selective environments, i.e. at elevated temperatures reflecting patterns of disease outbreaks in natural populations. Resistance against sympatric and allopatric Vibrio spp. strains was dominantly inherited in crosses between both invasion sources, resulting in an overall higher resistance of admixed individuals than pure lines. Therefore, we suggest that a simple genetic resistance mechanism of the host is matched to a common virulence mechanism shared by local Vibrio strains. This combination might have facilitated a fast evolutionary response that can explain another dimension of why invasive species can be so successful in newly invaded ranges.

Identification of a Serum amyloid A gene and the association of SNPs with Vibrio-resistance and growth traits in the clam Meretrix meretrix

Serum amyloid A (SAA), an acute response protein as well as an apolipoprotein, is considered to play crucial roles in both innate immunity and lipid metabolism. In this study, a SAA gene (MmSAA) was identified in the clam Meretrix meretrix. The full length DNA of MmSAA was 1407bp, consisting of three exons and two introns. The distribution of MmSAA in clam tissues was examined with the highest expression in hepatopancreas. In response to the Vibrio parahaemolyticus challenge, MmSAA mRNA showed significantly higher expression at 24 h post-challenge in experimental clams (P < 0.05). Forty-eight single nucleotide polymorphisms (SNPs) in the DNA partial sequence of MmSAA were discovered and examined for their association with Vibrio-resistance and growth traits, respectively. The single SNP association analysis indicated that five single SNPs (g.42, g.72, g.82, g.147 and g.165) were significantly associated with Vibrio-resistance (P < 0.05). Haplotype analysis produced additional support for association with the Chi-square values 6.393 (P = 0.012). Among the five selected SNPs, the effect of a missense mutation (g.82, A → G) was detected by site-directed mutagenesis with fusion expression of protein assay, and the result showed that the recombinant plasmids containing wild-type pET30a-MmSAA had more inhibition effect than the mutant ones on the growth rate of the host bacteria. In addition, four growth traits of the clams in 09G3SPSB population were recorded and the SNP g.176 was found to be significantly associated with the growth traits with the Global score value 0.790 (P = 0.015). Our findings suggested that common genetic variation in MmSAA might contribute to the risk of susceptibility to Vibrio infection and might be associated with the growth traits in the clams M. meretrix, and more works are still needed to validate these SNPs as potential markers for actual selective breeding.

Development of Vibrio spp. infection resistance related SNP markers using multiplex SNaPshot genotyping method in the clam Meretrix meretrix

The clam Meretrix meretrix is a commercially important mollusc species in the coastal areas of South and Southeast Asia. In the present study, large-scale SNPs were genotyped by the Multiplex SNaPshot genotyping method among the stocks of M. meretrix with different Vibrio spp. infection resistance profile. Firstly, the AUTOSNP software was applied to mine SNPs from M. meretrix transcriptome, and 323 SNP loci (including 120 indels) located on 64 contigs were selected based on Uniprot-GO associations. Then, 38 polymorphic SNP loci located on 15 contigs were genotyped successfully in the clam stocks with different resistance to Vibrio parahaemolyticus infection (11-R and 11-S groups). Pearson's Chi-square test was applied to compare the allele and genotype frequency distributions of the SNPs between the different stocks, and seven SNP markers located on three contigs were found to be associated with V. parahaemolyticus infection resistance trait. Haplotype-association analysis showed that six haplotypes had significantly different frequency distributions in 11-S and 11-R (P < 0.05). With selective genotyping between 09-R and 09-C populations, which had different resistance to Vibrio harveyi infection, four out of the seven selected SNPs had significantly different distributions (P < 0.05) and therefore they were considered to be associated with Vibrio spp. infection resistance. Sequence alignments and annotations indicated that the contigs containing the associated SNPs had high similarity to the immune related genes. All these results would be useful for the future marker-assisted selection of M. meretrix strains with high Vibrio spp. infection resistance.

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

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

Expression patterns of an i-type lysozyme in the clam Meretrix meretrix along with larval development

An i-type lysozyme (MmeLys) has been proved to function in immunity of the clam Meretrix meretrix in our previous studies. In this study, the expression patterns of MmeLys mRNA and protein at four chief developmental stages of M. meretrix were analyzed, which was able to provide information about how ontogeny of immunity and, in particular, antibacterial ability occured in the bivalve. The results of real-time PCR and western blot showed that MmeLys expressions were activated in D-veligers and dramatically increased to the highest level in pediveligers. It is proposed that the expression changes at these two stages might be due to the visceral organs changes, which were related to the archenteron formation in D-veligers and the organ-restructuring in pediveligers during metamorphosis. In addition, new methods of whole mount in situ hybridization and whole mount immunofluorescence were applied to identify the MmeLys expression tissues, and these tissues (i.e. hepatopancreas, gill, mantle, mouth, velum and foot) may be involved in the immune function during development of clams. Our study is valuable to a certain extent for exploring the origin of immune functions in clams and provides new methodology for future studies on the immune ontogeny of bivalves.

Three vibrio-resistance related EST-SSR markers revealed by selective genotyping in the clam Meretrix meretrix

The clam Meretrix meretrix is an important commercial bivalve distributed in the coastal areas of South and Southeast Asia. In this study, marker-trait association analyses were performed based on the stock materials of M. meretrix with different vibrio-resistance profile obtained by selective breeding. Forty-eight EST-SSR markers were screened and 27 polymorphic SSRs of them were genotyped in the clam stocks with different resistance to Vibrio parahaemolyticus (11-R and 11-S) and to Vibrio harveyi (09-R and 09-C). Allele frequency distributions of the SSRs among different stocks were compared using Pearson's Chi-square test, and three functional EST-SSR markers (MM959, MM4765 and MM8364) were found to be associated with vibrio-resistance trait. The 140-bp allele of MM959 and 128-bp allele of MM4765 had significantly higher frequencies in resistant groups (11-R and 09-R) than in susceptive/control groups (11-S and 09-C) (P < 0.01), which suggested that the clams carrying these two alleles have stronger resistance against vibrio. Clam individuals of 11-S were divided into three subgroups based on the survival time post-challenge, and the multi-dimensional scaling (MDS) analysis showed that clusters generated by genetic similarity revealed by the three SSR markers were consistent with the three subgroups distinctions. The putative functions of contig959, contig4765 and contig8364 also suggested that the three SSR-involved genes might play important roles in immunity of M. meretrix. All these results supported that EST-SSR markers MM959, MM4765 and MM8364 were associated with vibrio-resistance and would be useful for marker-assisted selection (MAS) in M. meretrix genetic breeding.

Development of molecular resources for an intertidal clam, Sinonovacula constricta, using 454 transcriptome sequencing

Background

The razor clam Sinonovacula constricta is a benthic intertidal bivalve species with important commercial value. Despite its economic importance, knowledge of its transcriptome is scarce. Next generation sequencing technologies offer rapid and efficient tools for generating large numbers of sequences, which can be used to characterize the transcriptome, to develop effective molecular markers and to identify genes associated with growth, a key breeding trait.

Results

Total RNA was isolated from the mantle, gill, liver, siphon, gonad and muscular foot tissues. High-throughput deep sequencing of S. constricta using 454 pyrosequencing technology yielded 859,313 high-quality reads with an average read length of 489 bp. Clustering and assembly of these reads produced 16,323 contigs and 131,346 singletons with average lengths of 1,376 bp and 458 bp, respectively. Based on transcriptome sequencing, 14,615 sequences had significant matches with known genes encoding 147,669 predicted proteins. Subsequently, previously unknown growth-related genes were identified. A total of 13,563 microsatellites (SSRs) and 13,634 high-confidence single nucleotide polymorphism loci (SNPs) were discovered, of which almost half were validated.

Conclusion

De novo sequencing of the razor clam S. constricta transcriptome on the 454 GS FLX platform generated a large number of ESTs. Candidate growth factors and a large number of SSRs and SNPs were identified. These results will impact genetic studies of S. constricta.

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.