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

Transcriptional responses of Acropora hyacinthus embryo under the benzo(a)pyrene stress by deep sequencing

Coral embryos are a critical and sensitive period for the early growth and development of coral. Benzo(a)pyrene (BaP) is widely distributed in the ocean and has strong toxicity, but there is little information on the toxic effects to coral embryos exposed to this widespread environmental contaminant. Thus, in this study, we utilized the Illumina Hiseq™ 4000 platform to explore the gene response of Acropora hyacinthus embryos under the BaP stress. A total of 130,042 Unigenes were obtained and analyzed, and approximately 37.67% of those matched with sequences from four different species. In total, 2606 Unigenes were up-regulated, and 3872 Unigenes were down-regulated. After Gene Ontology (GO) annotation, the results show that the "cellular process" and "metabolic process" were leading in the category of biological processes, which the "binding" and "catalytic activity" were the most abundant subcategories in molecular function. Based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, the most differentially expressed genes (DEGs) were enriched, as well as down-regulated in the pathways of oxidative phosphorylation, metabolism of xenobiotics, immune-related genes, apoptosis and human disease genes. At the same time, 388,197 of Single-nucleotide Polymorphisms (SNPs) and 6164 of Simple Sequence Repeats (SSRs) were obtained, which can be served as the richer and more valuable SSRs molecular markers in the future. The results of this study can help to better understand the toxicological mechanism of coral embryo exposed to BaP, and it is also essential for the protection and restoration of coral reef ecosystem in the future.

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.

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.

mRNA-Seq and microarray development for the Grooved Carpet shell clam, Ruditapes decussatus: a functional approach to unravel host-parasite interaction

Background

The Grooved Carpet shell clam Ruditapes decussatus is the autochthonous European clam and the most appreciated from a gastronomic and economic point of view. The production is in decline due to several factors such as Perkinsiosis and habitat invasion and competition by the introduced exotic species, the manila clam Ruditapes philippinarum. After we sequenced R. decussatus transcriptome we have designed an oligo microarray capable of contributing to provide some clues on molecular response of the clam to Perkinsiosis.

Results

A database consisting of 41,119 unique transcripts was constructed, of which 12,479 (30.3%) were annotated by similarity. An oligo-DNA microarray platform was then designed and applied to profile gene expression in R. decussatus heavily infected by Perkinsus olseni. Functional annotation of differentially expressed genes between those two conditionswas performed by gene set enrichment analysis. As expected, microarrays unveil genes related with stress/infectious agents such as hydrolases, proteases and others. The extensive role of innate immune system was also analyzed and effect of parasitosis upon expression of important molecules such as lectins reviewed.

Conclusions

This study represents a first attempt to characterize Ruditapes decussatus transcriptome, an important marine resource for the European aquaculture. The trancriptome sequencing and consequent annotation will increase the available tools and resources for this specie, introducing the possibility of high throughput experiments such as microarrays analysis. In this specific case microarray approach was used to unveil some important aspects of host-parasite interaction between the Carpet shell clam and Perkinsus, two non-model species, highlighting some genes associated with this interaction. Ample information was obtained to identify biological processes significantly enriched among differentially expressed genes in Perkinsus infected versus non-infected gills. An overview on the genes related with the immune system on R. decussatus transcriptome is also reported.

Simple sequence repeats amplification

The technique of SSR amplification is a prerequisite to generate the molecular profiles of various alleles of an individual or genotype. Amplification is the multifold duplication and accumulation of a targeted region which is achieved by polymerase chain reaction. It needs ingredients such as buffer, MgCl2, dNTPs, primers, and DNA polymerase enzyme. The utilization of these essential PCR components in optimal concentrations determines the success of amplification. Thus SSRs, as primers, play an important role in enhancing the amplification and thereby generating the genotype profile. With the advent of technology, fluorophore-labeled primers along with automated capillary electrophoresis system have enhanced the efficiency of detection.

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.

Three EST-SSR markers associated with QTL for the growth of the clam Meretrix meretrix revealed by selective genotyping

The clam Meretrix meretrix is a member of widely cultured, commercially important clams. A marker-trait association analysis was performed using expressed sequence tag (EST) simple sequence repeat (SSR) markers for marker-assisted selection in M. meretrix. Three markers, MM1272, MM2034, and MM7721, were found to be significantly associated with quantitative trait loci (QTLs) controlling shell length (P < 0.0001) in clams of a fast-growing population (JSF) and a control population (JSC). The 144-bp allele of MM1272, the 154-bp allele of MM2034, and the 152- and 165-bp alleles of MM7721 showed a significantly higher frequency in the JSF population (17.65, 36.41, 28.67, and 29.33 %) than in the JSC population (4.65, 8.33, 3.47, and 5.56 %). The three markers showed lower values for the number of alleles and observed heterozygosity as well as a higher proportion of homozygotes in JSF than in JSC population. The three markers have been further confirmed in the high and low tails of another population (09G₃SPSB); similarly, lower values for the number of alleles and observed heterozygosity as well as a higher proportion of homozygotes were found in 09G₃SPSB(H). The putative functions of the three gene fragments containing MM1272, MM2034, and MM7721 also suggested that the three SSR-containing genes might be involved in growth of M. meretrix. All the results suggest that the three EST-SSR markers associated with growth QTLs would be useful for marker-assisted selection in M. meretrix breeding.

Sequencing and characterization of striped venus transcriptome expand resources for clam fishery genetics

Background

The striped venus Chamelea gallina clam fishery is among the oldest and the largest in the Mediterranean Sea, particularly in the inshore waters of northern Adriatic Sea. The high fishing pressure has lead to a strong stock abundance decline, enhanced by several irregular mortality events. The nearly complete lack of molecular characterization limits the available genetic resources for C. gallina. We achieved the first transcriptome of this species with the aim of identifying an informative set of expressed genes, potential markers to assess genetic structure of natural populations and molecular resources for pathogenic contamination detection.

Methodology/Principal Findings

The 454-pyrosequencing of a normalized cDNA library of a pool C. gallina adult individuals yielded 298,494 raw reads. Different steps of reads assembly and filtering produced 36,422 contigs of high quality, one half of which (18,196) were annotated by similarity. A total of 111 microsatellites and 20,377 putative SNPs were identified. A panel of 13 polymorphic transcript-linked microsatellites was developed and their variability assessed in 12 individuals. Remarkably, a scan to search for contamination sequences of infectious origin indicated the presence of several Vibrionales species reported to be among the most frequent clam pathogen's species. Results reported in this study were included in a dedicated database available at http://compgen.bio.unipd.it/chameleabase.

Conclusions/Significance

This study represents the first attempt to sequence and de novo annotate the transcriptome of the clam C. gallina. The availability of this transcriptome opens new perspectives in the study of biochemical and physiological role of gene products and their responses to large and small-scale environmental stress in C. gallina, with high throughput experiments such as custom microarray or targeted re-sequencing. Molecular markers, such as the already optimized EST-linked microsatellites and the discovered SNPs will be useful to estimate effects of demographic processes and to detect minute levels of population structuring.