Identification of host-defense genes and development of microsatellite markers from ESTs of hard clam Meretrix meretrix

Development and validation of genome-wide SSR molecular markers of Tapes dorsatus

BACKGROUND

Tapes dorsatus is an economically important benthic animal in the Beibu Gulf of China. However, the deficiency of microsatellite markers has hindered the study of its genetics. The development of microsatellite markers will provide useful tools for genetic improvement, variety identification, phylogenetic analysis and resource conservation.

METHODS AND RESULTS

Within the genome sequence, 145,008 simple sequence repeats (SSRs) were identified, and 29,691 primer pairs were designed successfully. A total of 100 primer pairs were randomly synthesized for testing, and 93 primers yielded products. Sixty highly polymorphic primers were used to reveal the genetic diversity of 50 T. dorsatus individuals. The average number of alleles (Na) of the population was 10.40; the average number of effective alleles was 6.16, the average expected heterozygosity (He) was 0.82, and the average polymorphic information content was 0.80. The genetic structure of the population was detected, by which the population could be divided into three subpopulations.

CONCLUSION

We identified 145,008 SSRs in the genome of T. dorsatus and designed 29,691 primer pairs in this study. Of 100 synthesized primers, 60 were highly polymorphic and used to reveal the genetic diversity and structure of the population. The SSR markers identified here will provide useful tools and a foundation for genetic diversity, linkage mapping and molecular marker-aided breeding in T. dorsatus.

Proteome Response of Meretrix Bivalves Hepatopancreas Exposed to Paralytic Shellfish Toxins Producing Dinoflagellate Gymnodinium catenatum

Paralytic shellfish toxins (PSTs) contamination of seafood has become a growing global problem. However, the molecular response of bivalves, some of the most popular seafoods, to PSP toxins has seldom been reported and the underlying molecular mechanisms of the interactions between Meretrix meretrix bivalves and PSTs-producing dinoflagellates are scarcely known. This study compared the protein expression profiles between PSP toxin-contaminated and non-PSP toxin contaminated M. meretrix, determined proteome responses and identified potential biomarkers based on feeding experiments. Results showed that the content of total PSP toxins in contaminated bivalves was 40.63 ± 4.08 μg saxitoxin (STX) equivalents per gram, with 95.3% in hepatopancreas, followed by gill (1.82%) and foot (1.79%). According to two-dimensional gel electrophoresis (2-DE), 15 differentially expressed proteins (at least 2-fold difference) between the hepatopancreas of bivalves with and without PSP toxins were detected. Eight of them were successfully identified by MALDI-TOF MS. These were catalase, protein ultraspiracle homolog, G2 and S phase-expression protein, paramyosin, Mn-superoxide dismutase, response regulator receiver domain-containing protein, sarcoplasmic calcium-binding protein and major facilitator superfamily transporters. The differences in the expression levels of the last three proteins involving in cell signaling, structure and membrane transport were 4.2, 5.3 and 4.9-fold, respectively. These proteins could be further developed as potential biomarkers. The other two up-regulated proteins, Mn-superoxide dismutase and catalase, were involved in cell defence mechanisms against oxidative stress, suggesting PSP toxin acts as xenobiotics and poses oxidative stress in bivalves. This study gives insights into the response of bivalves to PSP toxin-producing dinoflagellate at the proteomic level and the potential of using 2-DE to develop specific protein markers in bivalves.

A preliminary attempt to explore the potential functions of a tetraspanin gene (MmTSPAN) in the innate immunity of hard clam Meretrix meretrix: Sequence features and expression profiles

Tetraspanins belong to the transmembrane 4 superfamily (TM4SF), and play crucial roles in immune responses. In the present study, a novel tetraspanin gene (designated MmTSPAN) was cloned and characterized from the hard clam Meretrix meretrix. The complete cDNA sequence of MmTSPAN contained an open reading frame (ORF) of 816 bp, which encoded a protein of 271 amino acids. MmTSPAN exhibited highly similarity with previously identified tetraspanins from other species. It contained four transmembrane domains (12-35 aa, 69-92 aa, 99-123 aa and 238-261 aa), characteristic CCG motif and four conservative cysteine residues. The mRNA transcripts of MmTSPAN were ubiquitously detectable in all the tested tissues, with the highest expression level in hepatopancreas. Temporal transcriptional levels in the hepatopancreas revealed significant up-regulation of MmTSPAN by Vibrio splendidus stimulation, with a 3.14-fold increase at 6 h compared to the control, and reaching 32.98-fold at 24 h. These results provide useful information for further study of the function of tetraspanin in the innate immune system of M. meretrix, and may offer a new therapeutic target for diseases of M. meretrix.

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

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

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.

De novo assembly, gene annotation, and simple sequence repeat marker development using Illumina paired-end transcriptome sequences in the pearl oyster Pinctada maxima

We analyzed the mantle transcriptome of pearl oyster Pinctada maxima and developed EST-SSR markers using Illumina HiSeq 2000 paired-end sequencing technology. A total of 49,500,748 raw reads were generated. De novo assembly generated 108,704 unigenes with an average length of 407 bp. Sequence similarity search with known proteins or nucleotides revealed that 30,200 (27.78%) and 25,824 (23.76%) consensus sequences were homologous with the sequences in the non-redundant protein and Swiss-Prot databases, respectively, and that 19,701 (18.12%) of these unigenes were possibly involved in approximately 234 known signaling pathways in the Kyoto Encyclopedia of Genes and Genomes database. Ninety one biomineralization-related unigenes were detected. In a cultured stock, 1764 simple sequence repeats were identified and 56 primer pairs were randomly selected and tested. The rate of successful amplification was 68.3%. The developed molecular markers are helpful for further studies on genetic linkage analysis, gene localization, and quantitative trait loci mapping.

Transcriptome analysis of the differences in gene expression between testis and ovary in green mud crab (Scylla paramamosain)

Background

The green mud crab (Scylla paramamosain) is the most prevalent crustacean on the southeast coast of China. The molecular regulatory mechanism of sex determination and gonadal differentiation in this species has received considerable attention in recent years because of the huge differences—both biological and economic—between male and female crabs. In this study, next-generation sequencing technology was used to develop deep-coverage transcriptomic sequencing data for the testis and ovary of S. paramamosain.

Results

A total of 365,116 reads (testis 171,962, ovary 193,154) with an average sequence length of 285 bp were produced from testis and ovary cDNA libraries. After filtering out contaminating reads, the clean reads were assembled, producing a total of 21,791 isotigs and leaving 22,814 reads as singlets. Using the BLASTX program, 3,471 unique sequences (2,275 isotigs and 1,196 singletons) were annotated with known protein sequences from the NCBI non-redundant (Nr) protein sequence database. The Gene Ontology and KEGG (Kyoto Encyclopedia of Genes and Genomes) analyses allowed the 224 unique sequences that were annotated with enzyme code (EC) numbers to be mapped into 174 KEGG pathways. After comparing the ovary and testis libraries, 4,021 gonad-differentially, 10,522 ovary-specifically, and 19,013 testis-specifically expressed genes were identified. Moreover, 33 ovary-specific, 14 testis-specific, and 34 gonad-differential transcripts were confirmed by semi-quantitative PCR and quantitative real-time PCR. In addition, 8,610 putative simple sequence repeats (SSRs) and 23,879 potential single nucleotide polymorphisms (SNPs) were identified.

Conclusion

This is the first large-scale RNA sequencing of S. paramamosain to be reported. We have identified many important functional genes and made a preliminary attempt to construct the regulatory network involved in the gonadal development of crustaceans. The annotated transcriptome data will provide fundamental support for future research into the reproduction biology of S. paramamosain. A large number of candidate SSRs and SNPs were detected, which could be used as genetic markers for population genetics and functional genomics in this species.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-585) contains supplementary material, which is available to authorized users.

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.

De novo assembly and characterization of spotted seal Phoca largha transcriptome using Illumina paired-end sequencing

Spotted seal (Phoca largha) is categorized as a critically endangered species in China. The aim of this study was to investigate spotted seal transcriptome by the approach of Illumina paired-end sequencing technology. We obtained a total of 52,146,394 reads for the mixed tissues of liver and spleen from the spotted seal. The de novo assemblies yielded 354,014 contigs and 178,466 unigenes. In the transcriptome, 193 unigenes were assigned to defense mechanisms. Three unigenes encoded MHC class I and 17 unigenes encoded MHC class II. In addition, bioinformatics analysis revealed a total of 4425 simple sequence repeats (SSRs). Fifty SSRs were randomly selected to validate amplification and determine the degree of polymorphism in the genomic DNA pools. Thirty-five primer pairs successfully amplified the expected DNA fragments and detected significant polymorphism among 28 spotted seal individuals. These results would contribute to the understanding of the genetic makeup of spotted seal transcriptome and provide useful information for functional genomic research in this species.

A consensus microsatellite-based linkage map for the hermaphroditic bay scallop (Argopecten irradians) and its application in size-related QTL analysis

Bay scallop (Argopecten irradians) is one of the most economically important aquaculture species in China. In this study, we constructed a consensus microsatellite-based genetic linkage map with a mapping panel containing two hybrid backcross-like families involving two subspecies of bay scallop, A. i. irradians and A. i. concentricus. One hundred sixty-one microsatellite and one phenotypic (shell color) markers were mapped to 16 linkage groups (LGs), which corresponds to the haploid chromosome number of bay scallop. The sex-specific map was 779.2 cM and 781.6 cM long in female and male, respectively, whereas the sex-averaged map spanned 849.3 cM. The average resolution of integrated map was 5.9 cM/locus and the estimated coverage was 81.3%. The proportion of distorted markers occurred more in the hybrid parents, suggesting that the segregation distortion was possibly resulted from heterospecific interaction between genomes of two subspecies of bay scallop. The overall female-to-male recombination rate was 1.13:1 across all linked markers in common to both parents, and considerable differences in recombination also existed among different parents in both families. Four size-related traits, including shell length (SL), shell height (SH), shell width (SW) and total weight (TW) were measured for quantitative trait loci (QTL) analysis. Three significant and six suggestive QTL were detected on five LGs. Among the three significant QTL, two (qSW-10 and qTW-10, controlling SW and TW, respectively) were mapped on the same region near marker AiAD121 on LG10 and explained 20.5% and 27.7% of the phenotypic variance, while the third (qSH-7, controlling SH) was located on LG7 and accounted for 15.8% of the phenotypic variance. Six suggestive QTL were detected on four different LGs. The linkage map and size-related QTL obtained in this study may facilitate marker-assisted selection (MAS) in bay scallop.

Characterization of the spotted seal Phoca largha transcriptome using Illumina paired-end sequencing and development of SSR markers

Next-generation sequencing provides a powerful new approach for developing functional genomic tools for non-model species. This study aims to investigate the spotted seal (Phoca largha) transcriptome by the approach of Illumina paired-end sequencing technology. We obtained a total of 52,146,394 reads for the mixed tissues of liver and spleen from spotted seal. The de novo assemblies yielded 354,014 contigs and 178,466 unigenes. In addition, bioinformatics analysis revealed a total of 4425 simple sequence repeats (SSRs). Fifty SSRs were randomly selected to validate amplification and determine the degree of polymorphism in the genomic DNA pools. Thirty-five primer pairs successfully amplified expected DNA fragments and detected significant polymorphism among 28 spotted seal individuals. These results contribute to the understanding of the genetic makeup of spotted seal transcriptome and provide useful information for functional genomic research in this species.

Transcriptomic analysis of the clam Meretrix meretrix on different larval stages

The clam Meretrix meretrix (Mollusca: Bivalvia) is an important commercial species in China. The deficiency of genomic and transcriptomic data is becoming the bottleneck of further researches on its complex and unique developmental processes. To improve this situation, in this study, a large-scale RNA sequencing was conducted on M. meretrix on larval stages. In particular, mRNAs of trochophore, D-veliger, pediveliger, and postlarva were purified, reverse transcribed, and sequenced through 454 sequencing technology. A total of 704,671 reads were obtained and assembled into 124,737 unique sequences (35,205 contigs and 89,532 singletons). Further analysis showed that 118,075 (94.66%) of these sequences were low-expression-level transcripts. Fifteen thousand two hundred fifteen (12.20%) of the unique sequences were annotated by searching against Uniprot Protein Knowledgebase, while the others (109,522, 87.80%) were left as novel sequences. Gene ontology analysis of the annotated sequences showed that most of them were assigned to certain gene ontology terms. By analyzing the depth of each unique sequence, a preliminary quantification analysis was conducted. An amount of sequences that showed a dramatic transcript discrepancy among the four larval stages were screened, which were related to development, growth, shell formation, and immune responses etc. As the first attempt on large-scale RNA sequencing of marine bivalve larvae, this work would enrich the knowledge of larval development of marine bivalves and provide fundamental support for further researches.

Mining of EST-SSR markers in clam Meretrix meretrix larvae from 454 shotgun transcriptome

A total of 2,970 EST-SSRs (2.38%) were identified by transcriptome sequencing of clam Meretrix meretrix (751,970 reads, ~310.82 Mbp), using 454 Genome Sequencer FLX next-generation sequencing platform. Dinucleotide SSR was the dominant repeat type (40.2%), followed by trinucleotide (37.8%), tetranuleotide (12.0%) and pentanucleotide (2.0%) SSR. The dominant repeat motif was AT (71.3%) in the dinucleotide SSR type and AAC (45.6%) in the trinucleotide SSR type. Nearly 79% of all microsatellites had flanking sequences suitable for PCR primer design. Half of PAL were found to be polymorphic in a subset of 40 primer pairs randomly selected. Specifically, the density of dinucleotide, trinucleotide and tetranucleotide repeats showed significant variation among four development stages (trochophore, D-veliger, pediveliger and postlarva). The results suggested that dinucleotide, trinucleotide and tetranucleotide SSRs may play an important role in contributing to the different expression profiles in larval stages.

Major histocompatibility complex class IIA and IIB genes of the spotted halibut Verasper variegatus: genomic structure, molecular polymorphism, and expression analysis

The major histocompatibility complex (MHC) is a large genomic region characterized by extremely high polymorphism and its association with resistance/susceptibility to disease in vertebrates. In this study, the full lengths of MHC IIA and IIB cDNA were obtained from spotted halibut (Verasper variegates) by expressed sequence tag (EST) and rapid amplification of cDNA ends (RACE) approaches. The genomic structure, molecular polymorphism, and expression patterns were examined to study MHC II gene functions in fish. As in other teleosts, the genomic structure of the spotted halibut MHC IIA contained 4 exons and 3 introns. The deduced amino acid sequence of the class IIA molecule shared 28-79% similarity with those of teleosts and mammals. Nine class IIA alleles were identified from five individuals. Three alleles originating from a single individual suggested the existence of at least two class IIA loci in the genome. Six exons and 5 introns were identified from spotted halibut MHC IIB, and the deduced amino acid sequence shared 33-79% similarity with those of teleosts and mammals. Twelve alleles were identified, among which five were observed in a single individual, which suggested at least three class IIB loci. Quantitative real-time PCR analysis revealed the presence of class IIA and IIB transcripts in nine normal tissues with high expression level in kidney and gill. Furthermore, MHC IIA and IIB are probably two candidates of immune molecules involved in the acute-phase response in spotted halibut, because their transcriptional levels were significantly up-regulated in blood and liver after bacterial challenge.

Molecular characterization and expression analysis of an inhibitor of NF-κB (IκB) from Asiatic hard clam Meretrix meretrix

Inhibitor of NF-κB (IκB) is an important member of Rel/NF-κB signaling pathway, which is an important mediator of immune responses in innate immune system. In this study, the IκB cDNA of hard clam Meretrix meretrix (designated as Mm-IκB) was cloned and characterized. The full-length cDNA of Mm-IκB was of 2098 bp, containing a 5' untranslated region (UTR) of 123 bp, a 3' UTR of 810 bp with a poly (A) tail, and an open reading frame (ORF) of 1164 bp encoding a polypeptide of 387 amino acids. The high similarity of Mm-IκB with other IκBs from invertebrates indicated that Mm-IκB should be a member of IκB family. Similar to most IκBs, Mm-IκB possessed all conserved features critical for the fundamental structure and function of IκBs, such as five ankyrin repeats and a conserved degradation motif (DS(44)RYSS(48)). Two PEST domains and a phosphorylation site motif (S(367)EEE(370)) at the C-terminus of Mm-IκB were identified. By quantitative real-time RT-PCR analysis, mRNA level of Mm-IκB was found to be most abundantly expressed in the tissues of mantle, gill and hepatopancreas, weakly expressed in muscle, foot and haemocyte. The Mm-IκB gene expression was significantly up-regulated at 24 h in haemocyte and at 12 h in gill after Vibrio anguillarum challenge, respectively. The results suggested the involvement of Mm-IκB in response against bacterial infection and further highlighted its functional importance in the immune system of M. meretrix.

Identification and expression of a putative LPS-induced TNF-α factor from Asiatic hard clam Meretrix meretrix

LPS-induced TNF-α (LITAF) is a novel transcriptional factor that mediates the expression of inflammatory cytokines in LPS-induced processes. In the present study, the full-length cDNA encoding LITAF (designated as Mm-LITAF) was identified from Asiatic hard clam, Meretrix meretrix, by expressed sequence tag and rapid amplification of cDNA ends (RACE) approaches. The full-length cDNA of Mm-LITAF was 1653 bp, consisting of a 5' untranslated region (UTR) of 91 bp, a 3'UTR of 1166 bp with one cytokine RNA instability motif (ATTTA) and one polyadenylation signal (AATAAA), and an open reading frame (ORF) of 396 bp encoding a polypeptide of 131 amino acids with a theoretical isoelectric point of 7.49, and predicted molecular weight of 14.47 kDa. The deduced amino acid of Mm-LITAF shared 29-63% similarity with the LITAFs from other species, indicating that Mm-LITAF should be a member of the LITAF family. Two highly conserved CXXC motifs forming a compact Zn(2+)-binding structure were also identified in Mm-LITAF. A quantitative reverse transcriptase real-time PCR (qRT-PCR) assay was developed to assess the expression of Mm-LITAF mRNA in different tissues, and the temporal expression of Mm-LITAF in clams challenged with Vibrio anguillarum. The mRNA transcript of Mm-LITAF could be detected in all the examined tissues with the highest expression level in the gill. Mm-LITAF expression was up-regulated significantly at 16 h in the gill and at 8 h in haemocytes after bacterial challenge, respectively. These results suggest that the Mm-LITAF is a constitutive and inducible acute-phase protein that perhaps involved in the innate immune response of hard clam.

Sixteen polymorphic simple sequence repeat markers from expressed sequence tags of the Chinese Mitten Crab Eriocheir sinensis

The Chinese mitten crab (Eriocheir sinensis) is an economically important aquaculture species in China. In this study, we developed and evaluated simple sequence repeat markers from expressed sequence tags of E. sinensis. Among the 40 wild E. sinensis individuals tested, 16 loci were polymorphic. The number of alleles per locus ranged from two to ten. The observed heterozygosity ranged from 0.0667 to 0.9667, whereas the expected heterozygosity ranged from 0.0661 to 0.9051. These markers have the potential for use in genetic studies of population structure and intraspecific variation in E. sinensis.