Generation and analysis of ESTs from the eastern oyster, Crassostrea virginica Gmelin and identification of microsatellite and SNP markers

Comparative physiology and aquaculture: Toward Omics-enabled improvement of aquatic animal health and sustainable production

Omics-technologies have revolutionized biomedical research over the past two decades, and are now poised to play a transformative role in aquaculture. This article serves as an introduction to a Virtual Special Issue of Comparative Biochemistry and Physiology - Part D: Genomics and Proteomics (CBPD), with the objective to showcase the state-of-the-science for Omics in aquaculture. In this editorial, we describe the role that Omics can play in aquaculture, and provide a synopsis for each of the Special Issue articles that use these technologies to improve aquaculture practices. Current genomic resources available for some aquaculture species are also described. The number of datasets is impressive for species such as Atlantic salmon and rainbow trout, totaling in the thousands (NCBI Gene Expression Omnibus and Sequence Read Archive). We present a conceptual framework that describes how Omics can be leveraged to understand complex responses of aquatic animals in culture for relevant physiological outcomes, such as fecundity, growth, and immunity. Lastly, knowledge gaps and new directions are identified to address current obstacles in aquaculture. Articles in this Special Issue on aquaculture in CBPD highlight the diversity and scope of Omics in aquaculture. As the technology becomes more cost-effective, it is anticipated that genomics, transcriptomics, proteomics, metabolomics and lipidomics will play increasingly important roles in stock diagnostics (e.g. genetics, health, performance). The timing is right, as global concerns are reaching critical levels over food availability/security and water restrictions for humankind.

Population structure and genetic diversity of hatchery stocks as revealed by combined mtDNA fragment sequences in Indian major carp, Catla catla

Catla catla is the second most important Indian major carp due to high growth rate and acceptance to consumers for food value. It is widely cultured in the Indian subcontinent as monoculture or polyculture. In the present study, genetic diversity among hatchery stocks (total 218 samples of catla) collected from different geographical regions of India was examined using mtDNA fragment sequence of Cyt b (306 bp) and D loop (710 bp). High numbers (57) of population specific haplotypes were observed in the present study. The results revealed significant genetic heterogeneity for the sequence data (F_ST = 0.27546, p < .05). Analysis of molecular variance revealed significant genetic differentiation among different catla populations. The information generated in present study could be useful to develop broad genetic base populations of catla.

Construction and evaluation of a high-density SNP array for the Pacific oyster (Crassostrea gigas)

Single nucleotide polymorphisms (SNPs) are widely used in genetics and genomics research. The Pacific oyster (Crassostrea gigas) is an economically and ecologically important marine bivalve, and it possesses one of the highest levels of genomic DNA variation among animal species. Pacific oyster SNPs have been extensively investigated; however, the mechanisms by which these SNPs may be used in a high-throughput, transferable, and economical manner remain to be elucidated. Here, we constructed an oyster 190K SNP array using Affymetrix Axiom genotyping technology. We designed 190,420 SNPs on the chip; these SNPs were selected from 54 million SNPs identified through re-sequencing of 472 Pacific oysters collected in China, Japan, Korea, and Canada. Our genotyping results indicated that 133,984 (70.4%) SNPs were polymorphic and successfully converted on the chip. The SNPs were distributed evenly throughout the oyster genome, located in 3,595 scaffolds with a length of ~509.4 million; the average interval spacing was 4,210 bp. In addition, 111,158 SNPs were distributed in 21,050 coding genes, with an average of 5.3 SNPs per gene. In comparison with genotypes obtained through re-sequencing, ~69% of the converted SNPs had a concordance rate of >0.971; the mean concordance rate was 0.966. Evaluation based on genotypes of full-sib family individuals revealed that the average genotyping accuracy rate was 0.975. Carrying 133 K polymorphic SNPs, our oyster 190K SNP array is the first commercially available high-density SNP chip for mollusks, with the highest throughput. It represents a valuable tool for oyster genome-wide association studies, fine linkage mapping, and population genetics.

Transcriptomic Analysis of Differentially Expressed Genes During Larval Development of Rapana venosa by Digital Gene Expression Profiling

During the life cycle of shellfish, larval development, especially metamorphosis, has a vital influence on the dynamics, distribution, and recruitment of natural populations, as well as seed breeding. Rapana venosa, a carnivorous gastropod, is an important commercial shellfish in China, and is an ecological invader in the United States, Argentina, and France. However, information about the mechanism of its early development is still limited, because research in this area has long suffered from a lack of genomic resources. In this study, 15 digital gene expression (DGE) libraries from five developmental stages of R. venosa were constructed and sequenced on the IIIumina Hi-Sequation 2500 platform. Bioinformaticsanalysis identified numerous differentially and specifically expressed genes, which revealed that genes associated with growth, nervous system, digestive system, immune system, and apoptosis participate in important developmental processes. The functional analysis of differentially expressed genes was further implemented by gene ontology, and Kyoto encyclopedia of genes and genomes enrichment. DGE profiling provided a general picture of the transcriptomic activities during the early development of R. venosa, which may provide interesting hints for further study. Our data represent the first comparative transcriptomic information available for the early development of R. venosa, which is a prerequisite for a better understanding of the physiological traits controlling development.

De novo transcriptome sequencing and analysis of Rapana venosa from six different developmental stages using Hi-seq 2500

The carnivorous whelk Rapana venosa is regarded as a biological invader with strong ecological fitness in the United States, Argentina, France and other countries. R. venosa may seriously damage bivalve resources. Nonetheless, in China, R. venosa is an important commercial species. Larval development, especially metamorphosis, influences the natural population and industrial breeding. However, there are few studies on the early development of R. venosa, and our understanding is further limited by a lack of genomic information. In this study, de novo sequencing was performed to obtain a comprehensive transcriptome profile during early development. A Hi-seq 2500 sequencing run produced 148,737,902 raw reads that were assembled into 1,137,556 unigenes (average length of 619 nucleotides, of which 49,673 could be annotated). The unigenes were assigned to biological processes and functions after annotation in Gene Ontology, eukaryotic Ortholog Groups and Kyoto Encyclopedia of Genes and Genomes. We also identified 93,196 simple sequence repeats among the unigenes. Six unique sequences associated with neuroendocrine function were analyzed by quantitative real-time PCR. Our data represent the first comprehensive transcriptomic resource for R. venosa. Functional annotation of the unigenes involved in various biological processes could stimulate research on the mechanisms of early development in this species. Understanding the mechanism of early development and metamorphosis would benefit antifouling research and aquaculture of R. venosa.

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.

Functional marker detection and analysis on a comprehensive transcriptome of large yellow croaker by next generation sequencing

Large yellow croaker (Larimichthys crocea) is an important economic fish in China and Eastern Asia. Because of the exhaustive fishing and overdense aquaculture, the wild population and the mariculture of the species are facing serious challenges on germplasm degeneration and susceptibility to infectious disease agents. However, a comprehensive transcriptome from multi-tissues of the species has not been reported and functional molecular markers have not yet been detected and analyzed. In this work, we applied RNA-seq with the Illumina Hiseq2000 platform for a multi-tissue sample of large yellow croaker and assembled the transcriptome into 88,103 transcripts. Of them, 52,782 transcripts have been successfully annotated by nt/nr, InterPro, GO and KEGG database. Comparing with public fish proteins, we have found that 34,576 protein coding transcripts are shared in large yellow croaker with zebrafish, medaka, pufferfish, and stickleback. For functional markers, we have discovered 1,276 polymorphic SSRs and 261, 000 SNPs. The functional impact analysis of SNPs showed that the majority (~75%) of small variants cause synonymous mutations in proteins, followed by variations in 3' UTR region. The functional enrichment analysis illuminated that transcripts involved in DNA bindings, enzyme activities, and signal pathways prominently exhibit less single-nucleotide variants but genes for the constituent of the muscular tissue, the cytoskeleton, and the immunity system contain more frequent SNP mutations, which may reflect the structural and functional selections of the translated proteins. This is the first work for the high-throughput detection and analysis of functional polymorphic SSR and SNP markers in a comprehensive transcriptome of large yellow croaker. Our study provides valuable transcript sequence and functional marker resources for the quantitative trait locus (QTL) identification and molecular selection of the species in the research community.

Transcriptomic analysis of candidate osmoregulatory genes in the eastern oyster Crassostrea virginica

BACKGROUND

The eastern oyster, Crassostrea virginica, is a euryhaline species that can thrive across a wide range of salinities (5-35). As with all estuarine species, individual oysters must be able to regulate their osmotic balance in response to constant temporal variation in salinity. At the population level, recurrent viability selection may be an additional mechanism shaping adaptive osmoregulatory phenotypes at the margins of oyster salinity tolerance. To identify candidate genes for osmoregulation, we sequenced, assembled, and annotated the transcriptome of wild juvenile eastern oysters from 'high' and 'low' salinity regimes. Annotations and candidates were mostly based on the Pacific oyster (Crassostrea gigas) genome sequence so osmoregulatory relevance in C. virginica was explored by testing functional enrichment of genes showing spatially discrete patterns of expression and by quantifying coding sequence divergence.

RESULTS

The assembly of sequence reads and permissive clustering of potentially oversplit alleles resulted in 98,729 reftigs (contigs and singletons). Of these, 50,736 were annotated with 9,307 belonging to a set of candidate osmoregulatory genes identified from the C. gigas genome. A total of 218,777 SNPs (0.0185 SNPs/bp) were identified in annotated reftigs of C. virginica. Amino acid divergence between translations of C. virginica annotated reftigs and C. gigas coding sequence averaged 23.2 % with an average dN/dS ratio of 0.074, suggesting purifying selection on protein sequences. The high and low salinity source oysters each expressed a subset of genes unique to that group, and the functions for these annotated genes were consistent with known molecular mechanisms for osmotic regulation in molluscs.

CONCLUSIONS

Most of the osmoregulatory gene candidates experimentally identified in C. gigas are present in this C. virginica transcriptome. In general these congeners show coding sequence divergence too high to make the C. gigas genome a useful reference for C. virginica bioinformatics. However, strong purifying selection is characteristic of the osmoregulatory candidates so functional annotations are likely to correspond. An initial examination of C. virginica presence/absence expression patterns across the salinity gradient in a single estuary suggests that many of these candidates have expression patterns that co-vary with salinity, consistent with osmoregulatory function in C. virginica.

Subtracted Transcriptome Profile of Tiger Shrimp (Penaeus monodon) That Survived WSSV Challenge

There is increased interest in the development of virus-resistant or improved shrimp stock because production is currently hindered by outbreaks and limited understanding of shrimp defense. Recent advancement now allows for high-throughput molecular studies on shrimp immunity. We used next-generation sequencing (NGS) coupled with suppression subtractive hybridization (SSH) to generate a transcriptome database of genes from tiger shrimp that survived White spot syndrome virus (WSSV) challenge. A total of 9,597 unique sequences were uploaded to NCBI Sequence Read Archive with accession number SRR577080. Sixty-five unique sequences, 6% of the total, were homologous to genes of Penaeus monodon. Genes that were initially related to bacterial infection and environmental stress such as 14-3-3 gene, heat shock protein 90, and calreticulin were also found including a few full-length gene sequences. Initial analysis of the expression of some genes was done. Hemocyanin, ferritin, and fortilin-binding protein exhibited differential expression between survivor and control tiger shrimps. Furthermore, candidate microsatellite markers for brood stock selection were mined and tested. Four trinucleotide and one dinucleotide microsatellites were successfully amplified. The study highlights the advantage of the NGS platform coupled with SSH in terms of gene discovery and marker generation.

SNP discovery from transcriptome of the swimbladder of Takifugu rubripes

Single nucleotide polymorphisms (SNPs) have become the marker of choice for genome-wide association studies in many species. High-throughput sequencing of RNA was developed primarily to analyze global gene expression, while it is an efficient way to discover SNPs from the expressed genes. In this study, we conducted transcriptome sequencing of the swimbladder of Takifugu rubripes using Illumina HiSeq2000 platform to identify gene-associated SNPs in the swimbladder. A total of 30,312,181 unique-mapped-reads were obtained from 44,736,850 raw reads. A total of 62,270 putative SNPs were discovered, which were located in 11,306 expressed genes and 2,246 scaffolds. The average minor allele frequency (MAF) of the SNPs was 0.26. GO and KEGG pathway analysis were conducted to analyze the genes containing SNPs. Validation of selected SNPs revealed that 54% of SNPs (26/48) were true SNPs. The results suggest that RNA-Seq is an efficient and cost-effective approach to discover gene-associated SNPs. In this study, a large number of SNPs were identified and these data will be useful resources for population genetic study, evolution analysis, resource assessment, genetic linkage analysis and genome-wide association studies.

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.

Pyrosequencing of Haliotis diversicolor transcriptomes: insights into early developmental molluscan gene expression

Background

The abalone Haliotis diversicolor is a good model for study of the settlement and metamorphosis, which are widespread marine ecological phenomena. However, information on the global gene backgrounds and gene expression profiles for the early development of abalones is lacking.

Methodology/Principal Findings

In this study, eight non-normalized and multiplex barcode-labeled transcriptomes were sequenced using a 454 GS system to cover the early developmental stages of the abalone H. diversicolor. The assembly generated 35,415 unigenes, of which 7,566 were assigned GO terms. A global gene expression profile containing 636 scaffolds/contigs was constructed and was proven reliable using qPCR evaluation. It indicated that there may be existing dramatic phase transitions. Bioprocesses were proposed, including the ‘lock system’ in mature eggs, the collagen shells of the trochophore larvae and the development of chambered extracellular matrix (ECM) structures within the earliest postlarvae.

Conclusion

This study globally details the first 454 sequencing data for larval stages of H. diversicolor. A basic analysis of the larval transcriptomes and cluster of the gene expression profile indicates that each stage possesses a batch of specific genes that are indispensable during embryonic development, especially during the two-cell, trochophore and early postlarval stages. These data will provide a fundamental resource for future physiological works on abalones, revealing the mechanisms of settlement and metamorphosis at the molecular level.

Genome-Wide SNP Discovery from Transcriptome of Four Common Carp Strains

Background

Single nucleotide polymorphisms (SNPs) have been used as genetic marker for genome-wide association studies in many species. Gene-associated SNPs could offer sufficient coverage in trait related research and further more could themselves be causative SNPs for traits. Common carp (Cyprinus carpio) is one of the most important aquaculture species in the world accounting for nearly 14% of freshwater aquaculture production. There are various strains of common carp with different economic traits, however, the genetic mechanism underlying the different traits have not been elucidated yet. In this project, we identified a large number of gene-associated SNPs from four strains of common carp using next-generation sequencing.

Results

Transcriptome sequencing of four strains of common carp (mirror carp, purse red carp, Xingguo red carp, Yellow River carp) was performed with Solexa HiSeq2000 platform. De novo assembled transcriptome was used as reference for alignments, and SNP calling was done through BWA and SAMtools. A total of 712,042 Intra-strain SNPs were discovered in four strains, of which 483,276 SNPs for mirror carp, 486,629 SNPs for purse red carp, 478,028 SNPs for Xingguo red carp and 488,281 SNPs for Yellow River carp were discovered, respectively. Besides, 53,893 inter-SNPs were identified. Strain-specific SNPs of four strains were 53,938, 53,866, 48,701, 40,131 in mirror carp, purse red carp, Xingguo red carp and Yellow River carp, respectively. GO and KEGG pathway analysis were done to reveal strain-specific genes affected by strain-specific non-synonymous SNPs. Validation of selected SNPs revealed that 48% percent of SNPs (12 of 25) were tested to be true SNPs.

Conclusions

Transcriptome analysis of common carp using RNA-Seq is a cost-effective way of generating numerous reads for SNP discovery. After validation of identified SNPs, these data will provide a solid base for SNP array designing and genome-wide association studies.

De novo assembly, characterization and functional annotation of pineapple fruit transcriptome through massively parallel sequencing

Background

Pineapple (Ananas comosus var. comosus), is an important tropical non-climacteric fruit with high commercial potential. Understanding the mechanism and processes underlying fruit ripening would enable scientists to enhance the improvement of quality traits such as, flavor, texture, appearance and fruit sweetness. Although, the pineapple is an important fruit, there is insufficient transcriptomic or genomic information that is available in public databases. Application of high throughput transcriptome sequencing to profile the pineapple fruit transcripts is therefore needed.

Methodology/Principal Findings

To facilitate this, we have performed transcriptome sequencing of ripe yellow pineapple fruit flesh using Illumina technology. About 4.7 millions Illumina paired-end reads were generated and assembled using the Velvet de novo assembler. The assembly produced 28,728 unique transcripts with a mean length of approximately 200 bp. Sequence similarity search against non-redundant NCBI database identified a total of 16,932 unique transcripts (58.93%) with significant hits. Out of these, 15,507 unique transcripts were assigned to gene ontology terms. Functional annotation against Kyoto Encyclopedia of Genes and Genomes pathway database identified 13,598 unique transcripts (47.33%) which were mapped to 126 pathways. The assembly revealed many transcripts that were previously unknown.

Conclusions

The unique transcripts derived from this work have rapidly increased of the number of the pineapple fruit mRNA transcripts as it is now available in public databases. This information can be further utilized in gene expression, genomics and other functional genomics studies in pineapple.

Expression of candidate genes related to metabolism, immunity and cellular stress during massive mortality in the American oyster Crassostrea virginica larvae in relation to biochemical and physiological parameters

Quantification of mRNA of genes related to metabolism, immunity and cellular stress was examined in relation to a massive mortality event during the culture of American oyster larvae, Crassostrea virginica which was probably, in regard to previous microbiological analysis, induced by Vibrio infection. To document molecular changes associated with the mortality event, mRNA levels were compared to biochemical and physiological data, previously described in a companion paper. Among the 18 genes studied, comparatively to the antibiotic control, 10 showed a lower relative gene expression when the massive mortality occurred. Six of them are presumed to be related to metabolism, corroborating the metabolic depression associated with the mortality event suggested by biochemical and physiological analyses. Relationships between the regulation of antioxidant enzyme activities, lipid peroxidation, and the mRNA abundance of genes linked to oxidative stress, cytoprotection, and immune response are also discussed. Finally, we observed an increase in the transcript abundance of two genes involved in apoptosis and cell regulation simultaneously with mortality, suggesting that these processes might be linked.

Identification, utilisation and mapping of novel transcriptome-based markers from blackcurrant (Ribes nigrum)

BACKGROUND

Deep-level second generation sequencing (2GS) technologies are now being applied to non-model species as a viable and favourable alternative to Sanger sequencing. Large-scale SNP discovery was undertaken in blackcurrant (Ribes nigrum L.) using transcriptome-based 2GS 454 sequencing on the parental genotypes of a reference mapping population, to generate large numbers of novel markers for the construction of a high-density linkage map.

RESULTS

Over 700,000 reads were produced, from which a total of 7,000 SNPs were found. A subset of polymorphic SNPs was selected to develop a 384-SNP OPA assay using the Illumina BeadXpress platform. Additionally, the data enabled identification of 3,000 novel EST-SSRs. The selected SNPs and SSRs were validated across diverse Ribes germplasm, including mapping populations and other selected Ribes species.SNP-based maps were developed from two blackcurrant mapping populations, incorporating 48% and 27% of assayed SNPs respectively. A relatively high proportion of visually monomorphic SNPs were investigated further by quantitative trait mapping of theta score outputs from BeadStudio analysis, and this enabled additional SNPs to be placed on the two maps.

CONCLUSIONS

The use of 2GS technology for the development of markers is superior to previously described methods, in both numbers of markers and biological informativeness of those markers. Whilst the numbers of reads and assembled contigs were comparable to similar sized studies of other non-model species, here a high proportion of novel genes were discovered across a wide range of putative function and localisation. The potential utility of markers developed using the 2GS approach in downstream breeding applications is discussed.

Transcriptome analysis in Concholepas concholepas (Gastropoda, Muricidae): mining and characterization of new genomic and molecular markers

The marine gastropod Concholepas concholepas, locally known as the "loco", is the main target species of the benthonic Chilean fisheries. Genetic and genomic tools are necessary to study the genome of this species in order to understand the molecular basis of its development, growth, and other key traits to improve the management strategies and to identify local adaptation to prevent loss of biodiversity. Here, we use pyrosequencing technologies to generate the first transcriptomic database from adult specimens of the loco. After trimming, a total of 140,756 Expressed Sequence Tag sequences were achieved. Clustering and assembly analysis identified 19,219 contigs and 105,435 singleton sequences. BlastN analysis showed a significant identity with Expressed Sequence Tags of different gastropod species available in public databases. Similarly, BlastX results showed that only 895 out of the total 124,654 had significant hits and may represent novel genes for marine gastropods. From this database, simple sequence repeat motifs were also identified and a total of 38 primer pairs were designed and tested to assess their potential as informative markers and to investigate their cross-species amplification in different related gastropod species. This dataset represents the first publicly available 454 data for a marine gastropod endemic to the southeastern Pacific coast, providing a valuable transcriptomic resource for future efforts of gene discovery and development of functional markers in other marine gastropods.

The Littorina sequence database (LSD)--an online resource for genomic data

We present an interactive, searchable expressed sequence tag database for the periwinkle snail Littorina saxatilis, an upcoming model species in evolutionary biology. The database is the result of a hybrid assembly between Sanger and 454 sequences, 1290 and 147,491 sequences respectively. Normalized and non-normalized cDNA was obtained from different ecotypes of L. saxatilis collected in the UK and Sweden. The Littorina sequence database (LSD) contains 26,537 different contigs, of which 2453 showed similarity with annotated proteins in UniProt. Querying the LSD permits the selection of the taxonomic origin of blast hits for each contig, and the search can be restricted to particular taxonomic groups. The database allows access to UniProt annotations, blast output, protein family domains (PFAM) and Gene Ontology. The database will allow users to search for genetic markers and identifying candidate genes or genes for expression analyses. It is open for additional deposition of sequence information for L. saxatilis and other species of the genus Littorina. The LSD is available at http://mbio-serv2.mbioekol.lu.se/Littorina/.

Development of expressed sequence tags from the pearl oyster, Pinctada martensii Dunker

The pearl oyster, Pinctada martensii, is the primary species used for the aquaculture production of marine pearls in China and Japan. Genetic tools and resources are needed to study the genome of this species and to understand the molecular basis of development, growth, host defense, pearl formation, and other important traits. In this study, we developed a set of expressed sequence tags (ESTs) for P. martensii. We constructed cDNA libraries from adult tissues and sequenced 7,128 ESTs. Clustering analysis identified 788 contigs (covering 5,769 ESTs) and 1,351 singletons, yielding a total of 2,139 unique genes. Of these unique genes, only 935 had significant (E-value ≤ 0.005) hits in GenBank, and the remaining 1,204 (56.3%) were novel. Most of the known genes are related to cellular structure, protein binding, and metabolic processes. Putative host-defense genes (86) were identified including C-type lectin, ferritin, polyubiquitin, proteases, protease inhibitors, scavenger receptors, heat shock proteins, and RAS oncogenes. The EST sequences developed in this study provide a valuable resource for future efforts on gene identification, marker development, and studies on molecular mechanism of host defense in pearl oysters.

Normalization of full-length-enriched cDNA

A well-recognized obstacle to efficient high-throughput analysis of cDNA libraries is the differential abundance of various transcripts in any particular cell type. Decreasing the prevalence of clones representing abundant transcripts before sequencing, using cDNA normalization, may significantly increase the efficacy of random sequencing and is essential for rare gene discovery. Duplex-specific nuclease (DSN) normalization allows the generation of normalized full-length-enriched cDNA libraries to permit a high gene discovery rate. The method is based on the unique properties of DSN from the Kamchatka crab and involves denaturation-reassociation of cDNA, degradation of the ds-fraction formed by abundant transcripts by DSN, and PCR amplification of the remaining ss-DNA fraction. The method has been evaluated in various plant and animal models.

Transcriptome and proteome analysis of Pinctada margaritifera calcifying mantle and shell: focus on biomineralization

BACKGROUND

The shell of the pearl-producing bivalve Pinctada margaritifera is composed of an organic cell-free matrix that plays a key role in the dynamic process of biologically-controlled biomineralization. In order to increase genomic resources and identify shell matrix proteins implicated in biomineralization in P. margaritifera, high-throughput Expressed Sequence Tag (EST) pyrosequencing was undertaken on the calcifying mantle, combined with a proteomic analysis of the shell.

RESULTS

We report the functional analysis of 276 738 sequences, leading to the constitution of an unprecedented catalog of 82 P. margaritifera biomineralization-related mantle protein sequences. Components of the current "chitin-silk fibroin gel-acidic macromolecule" model of biomineralization processes were found, in particular a homolog of a biomineralization protein (Pif-177) recently discovered in P. fucata. Among these sequences, we could show the localization of two other biomineralization protein transcripts, pmarg-aspein and pmarg-pearlin, in two distinct areas of the outer mantle epithelium, suggesting their implication in calcite and aragonite formation. Finally, by combining the EST approach with a proteomic mass spectrometry analysis of proteins isolated from the P. margaritifera shell organic matrix, we demonstrated the presence of 30 sequences containing almost all of the shell proteins that have been previously described from shell matrix protein analyses of the Pinctada genus. The integration of these two methods allowed the global composition of biomineralizing tissue and calcified structures to be examined in tandem for the first time.

CONCLUSIONS

This EST study made on the calcifying tissue of P. margaritifera is the first description of pyrosequencing on a pearl-producing bivalve species. Our results provide direct evidence that our EST data set covers most of the diversity of the matrix protein of P. margaritifera shell, but also that the mantle transcripts encode proteins present in P. margaritifera shell, hence demonstrating their implication in shell formation. Combining transcriptomic and proteomic approaches is therefore a powerful way to identify proteins involved in biomineralization. Data generated in this study supply the most comprehensive list of biomineralization-related sequences presently available among protostomian species, and represent a major breakthrough in the field of molluskan biomineralization.

Identification of immune genes of the Agamaki clam (Sinonovacula constricta) by sequencing and bioinformatic analysis of ESTs

The Agamaki clam (Sinonovacula constricta) is an economically important shellfish in Asia. However, genomic research on this species is still in its infancy, and genomic resources are largely unavailable. The objective of this study was to generate expressed sequence tags (ESTs) from a normalized liver complementary DNA library and to identify genes that function in immune defense. A total of 5,296 ESTs were sequenced, from which 540 contigs and 3,473 singletons were identified. BLAST homology analysis indicated that only 20.7% of these ESTs were homologues of known genes while the remaining 79.3% appeared to be novel sequences. Based on sequence similarities, 43 putative immune genes were identified that encode proteases and protease inhibitors, adhesive proteins, stress proteins, lysosomal enzymes, and signal transduction regulators. Our study thus provides both a large collection of novel transcripts and a detailed annotation of immune genes for an important bivalve species.

The first set of EST resource for gene discovery and marker development in pigeonpea (Cajanus cajan L.)

BACKGROUND

Pigeonpea (Cajanus cajan (L.) Millsp) is one of the major grain legume crops of the tropics and subtropics, but biotic stresses [Fusarium wilt (FW), sterility mosaic disease (SMD), etc.] are serious challenges for sustainable crop production. Modern genomic tools such as molecular markers and candidate genes associated with resistance to these stresses offer the possibility of facilitating pigeonpea breeding for improving biotic stress resistance. Availability of limited genomic resources, however, is a serious bottleneck to undertake molecular breeding in pigeonpea to develop superior genotypes with enhanced resistance to above mentioned biotic stresses. With an objective of enhancing genomic resources in pigeonpea, this study reports generation and analysis of comprehensive resource of FW- and SMD- responsive expressed sequence tags (ESTs).

RESULTS

A total of 16 cDNA libraries were constructed from four pigeonpea genotypes that are resistant and susceptible to FW ('ICPL 20102' and 'ICP 2376') and SMD ('ICP 7035' and 'TTB 7') and a total of 9,888 (9,468 high quality) ESTs were generated and deposited in dbEST of GenBank under accession numbers GR463974 to GR473857 and GR958228 to GR958231. Clustering and assembly analyses of these ESTs resulted into 4,557 unique sequences (unigenes) including 697 contigs and 3,860 singletons. BLASTN analysis of 4,557 unigenes showed a significant identity with ESTs of different legumes (23.2-60.3%), rice (28.3%), Arabidopsis (33.7%) and poplar (35.4%). As expected, pigeonpea ESTs are more closely related to soybean (60.3%) and cowpea ESTs (43.6%) than other plant ESTs. Similarly, BLASTX similarity results showed that only 1,603 (35.1%) out of 4,557 total unigenes correspond to known proteins in the UniProt database (or= 5 sequences detected 102 single nucleotide polymorphisms (SNPs) in 37 contigs. As an example, a set of 10 contigs were used for confirming in silico predicted SNPs in a set of four genotypes using wet lab experiments. Occurrence of SNPs were confirmed for all the 6 contigs for which scorable and sequenceable amplicons were generated. PCR amplicons were not obtained in case of 4 contigs. Recognition sites for restriction enzymes were identified for 102 SNPs in 37 contigs that indicates possibility of assaying SNPs in 37 genes using cleaved amplified polymorphic sequences (CAPS) assay.

CONCLUSION

The pigeonpea EST dataset generated here provides a transcriptomic resource for gene discovery and development of functional markers associated with biotic stress resistance. Sequence analyses of this dataset have showed conservation of a considerable number of pigeonpea transcripts across legume and model plant species analysed as well as some putative pigeonpea specific genes. Validation of identified biotic stress responsive genes should provide candidate genes for allele mining as well as candidate markers for molecular breeding.

Development of gene-targeted SNP markers for genomic mapping in Pacific abalone Haliotis discus hannai Ino

Useful and novel DNA markers are needed for aquaculture genetics and breeding. In this study, we report the discovery and development of gene-targeted single nucleotide polymorphisms (SNPs) for genomic mapping in the Pacific abalone Haliotis discus hannai Ino. Single EST or EST-contigs from 66 genes that had positive BLASTx matches (E-value ≤ 1e-8) were used for polymerase chain reaction (PCR) amplification. PCR products from the two parents of one mapping family were directly sequenced, and 83 SNP loci were found from 17 genes. Allele-specific PCR (AS-PCR) was developed and optimized for genotyping of 11 SNP loci in 120 progeny of the mapping family. Nine of the loci conformed to the expected Mendelian ratio of 1:1 based on the χ(2) test (P > 0.05) and could potentially be used for linkage map construction. Our data also indicate that the sequencing of two parents may be a practical strategy for the discovery of informative SNPs for linkage mapping in a particular mapping population.

Identification of CYP genes in Mytilus (mussel) and Crassostrea (oyster) species: first approach to the full complement of cytochrome P450 genes in bivalves

Understanding the fate and effects of organic chemicals in animals requires knowledge of cytochrome P450 (CYP) genes, which thus far are poorly known in bivalve mollusks. We searched for CYP sequences in EST databases for Mytilus and Crassostrea species, lophotrochozoan representatives of the protostomes. From ESTs averaging ca. 924bp, we identified 58 CYP genes in Mytilus californianus and 39 CYP genes in Crassostrea gigas. The sequences fell in all known animal CYP clans, and collectively they clustered in phylogenetic analysis with vertebrate CYP families 1, 2, 3, 4, 17, 20, 26 and 27. As in deuterostomes, a majority of the sequences fell in Clan 2. The CYP sequences found thus far in bivalves suggest a diversity consistent with that found in many other animal species. The present description of mollusk genes provides the overall framework for classification of any additional bivalve sequences. The sequences identified also will be useful in obtaining full-length sequences and in designing primers for analysis of expression of mussel and oyster CYP genes, or for recombinant protein expression to identify potential substrates for the bivalve CYP proteins, and understand their roles in xenobiotic detoxification and physiology of bivalves.

Generation and analysis of a 29,745 unique Expressed Sequence Tags from the Pacific oyster (Crassostrea gigas) assembled into a publicly accessible database: the GigasDatabase

BACKGROUND

Although bivalves are among the most-studied marine organisms because of their ecological role and economic importance, very little information is available on the genome sequences of oyster species. This report documents three large-scale cDNA sequencing projects for the Pacific oyster Crassostrea gigas initiated to provide a large number of expressed sequence tags that were subsequently compiled in a publicly accessible database. This resource allowed for the identification of a large number of transcripts and provides valuable information for ongoing investigations of tissue-specific and stimulus-dependant gene expression patterns. These data are crucial for constructing comprehensive DNA microarrays, identifying single nucleotide polymorphisms and microsatellites in coding regions, and for identifying genes when the entire genome sequence of C. gigas becomes available.

DESCRIPTION

In the present paper, we report the production of 40,845 high-quality ESTs that identify 29,745 unique transcribed sequences consisting of 7,940 contigs and 21,805 singletons. All of these new sequences, together with existing public sequence data, have been compiled into a publicly-available Website http://public-contigbrowser.sigenae.org:9090/Crassostrea_gigas/index.html. Approximately 43% of the unique ESTs had significant matches against the SwissProt database and 27% were annotated using Gene Ontology terms. In addition, we identified a total of 208 in silico microsatellites from the ESTs, with 173 having sufficient flanking sequence for primer design. We also identified a total of 7,530 putative in silico, single-nucleotide polymorphisms using existing and newly-generated EST resources for the Pacific oyster.

CONCLUSION

A publicly-available database has been populated with 29,745 unique sequences for the Pacific oyster Crassostrea gigas. The database provides many tools to search cleaned and assembled ESTs. The user may input and submit several filters, such as protein or nucleotide hits, to select and download relevant elements. This database constitutes one of the most developed genomic resources accessible among Lophotrochozoans, an orphan clade of bilateral animals. These data will accelerate the development of both genomics and genetics in a commercially-important species with the highest annual, commercial production of any aquatic organism.

Identification of genes involved in immune response, microsatellite, and SNP markers from expressed sequence tags generated from hemocytes of freshwater pearl mussel (Hyriopsis cumingii)

Triangle sail mussel (Hyriopsis cumingii) is the most important mussel species commercially exploited for freshwater pearl production in China. However, its genome research is still at the infantry. Genomic resources for this species are largely not available. The objectives of this study was to generate expressed sequence tags from a hemocyte cDNA library, to identify genes involved in defense mechanisms, and to identify polymorphic markers from the expressed sequence tag (EST) resources for genetic analysis. A total of 5,290 ESTs were sequenced, obtaining 481 contigs and 1,165 singletons. BLAST similarity analysis indicated almost half (46.5%) of these ESTs were homologs of known genes while 53.5% were transcripts of unknown identities. Based on sequence similarities, 50 genes were identified as putative genes involved in immune and defense functions such as hemocyte immune process, stress proteins, adhesive proteins, proteases and protease regulators, antimicrobial peptides, lysosomal enzymes, cell apoptosis, and cell cycle proteins. A total of 201 microsatellites were identified from these ESTs, with 31 having sufficient flanking sequences for primer design. Polymerase chain reaction amplification was successful for 18 primer pairs and 14 of them were polymorphic. A total of 987 putative single nucleotide polymorphisms were identified including 204 transitions, 611 transversions, and 172 indels; 12 of them were involved in nine genes of defense mechanisms. These resources provide the material basis for future marker validation and genetic linkage and quantitative trait loci analysis in the freshwater pearl mussel.

Expressed sequence tags from the zhikong scallop (Chlamys farreri): discovery and annotation of host-defense genes

A high-quality cDNA library was constructed from whole body tissues of the zhikong scallop, Chlamys farreri, challenged by Listonella anguillarum. A total of 5720 clones were sequenced, yielding 5123 expressed sequence tags (ESTs). Among the 3326 unique genes identified, 2289 (69%) genes had no significant (E-value < 1e-5) matches to known sequences in public databases and 194 (6%) matched proteins of unknown functions. The remaining 843 (25%) genes that exhibited homology with genes of known functions, showed broad involvement in metabolic processes (31%), cell structure and motility (20%), gene and protein expression (12%), cell signaling and cell communication (8%), cell division (4%), and notably, 25% of those genes were related to immune function. They included stress response genes, complement-like genes, proteinase and proteinase inhibitors, immune recognition receptors and immune effectors. The EST collection obtained in this study provides a useful resource for gene discovery and especially for the identification of host-defense genes and systems in scallops and other molluscs.

Functional annotation and analysis of expressed sequence tags from the hepatopancreas of mitten crab (Eriocheir sinensis)

Due to its popularity as a traditional food, intensive harvesting of the mitten crab (Eriocheir sinensis) is common and has lead to an increase in disease incidence, resulting in catastrophic losses to crab aquaculture. The hepatopancreas of E. sinensis is not only an important digestive organ but also an indispensable immune organ. We constructed a nonnormalized cDNA library from the hepatopancreas of E. sinensis and acquired 3,297 high-quality expressed sequence tags representing 1,178 unigenes. More than half of these unigenes were novel genes for this species; the remaining had homologs in public databases, which is of great importance for future functional research. We also investigated the association of these genes with immune processes for insight into one of the main functions of the hepatopancreas besides metabolism. Despite the relatively low sampling scalar of our cDNA library, we were able to demonstrate several important properties of the hepatopancreatic transcriptome and identified numerous genes that were closely associated with immune responses. These results might serve as the basis for an in-depth genomics study of E. sinensis, including transcriptome analysis, physical mapping, and whole genome sequencing.

Transcriptomic analysis of the entomopathogenic nematode Heterorhabditis bacteriophora TTO1

BACKGROUND

The entomopathogenic nematode Heterorhabditis bacteriophora and its symbiotic bacterium, Photorhabdus luminescens, are important biological control agents of insect pests. This nematode-bacterium-insect association represents an emerging tripartite model for research on mutualistic and parasitic symbioses. Elucidation of mechanisms underlying these biological processes may serve as a foundation for improving the biological control potential of the nematode-bacterium complex. This large-scale expressed sequence tag (EST) analysis effort enables gene discovery and development of microsatellite markers. These ESTs will also aid in the annotation of the upcoming complete genome sequence of H. bacteriophora.

RESULTS

A total of 31,485 high quality ESTs were generated from cDNA libraries of the adult H. bacteriophora TTO1 strain. Cluster analysis revealed the presence of 3,051 contigs and 7,835 singletons, representing 10,886 distinct EST sequences. About 72% of the distinct EST sequences had significant matches (E value < 1e-5) to proteins in GenBank's non-redundant (nr) and Wormpep190 databases. We have identified 12 ESTs corresponding to 8 genes potentially involved in RNA interference, 22 ESTs corresponding to 14 genes potentially involved in dauer-related processes, and 51 ESTs corresponding to 27 genes potentially involved in defense and stress responses. Comparison to ESTs and proteins of free-living nematodes led to the identification of 554 parasitic nematode-specific ESTs in H. bacteriophora, among which are those encoding F-box-like/WD-repeat protein theromacin, Bax inhibitor-1-like protein, and PAZ domain containing protein. Gene Ontology terms were assigned to 6,685 of the 10,886 ESTs. A total of 168 microsatellite loci were identified with primers designable for 141 loci.

CONCLUSION

A total of 10,886 distinct EST sequences were identified from adult H. bacteriophora cDNA libraries. BLAST searches revealed ESTs potentially involved in parasitism, RNA interference, defense responses, stress responses, and dauer-related processes. The putative microsatellite markers identified in H. bacteriophora ESTs will enable genetic mapping and population genetic studies. These genomic resources provide the material base necessary for genome annotation, microarray development, and in-depth gene functional analysis.

A simple and efficient method for isolating polymorphic microsatellites from cDNA

Background

Microsatellites in cDNA are useful as molecular markers because they represent transcribed genes and can be used as anchor markers for linkage and comparative mapping, as well as for studying genome evolution. Microsatellites in cDNA can be detected in existing ESTs by data mining. However, in most fish species, no ESTs are available or the number of ESTs is limited, although fishes represent half of the vertebrates on the earth. We developed a simple and efficient method for isolation of microsatellites from cDNA in fish.

Results

The method included normalization of 150 ng cDNA using 0.5 U duplex-specific nuclease (DSN) at 65°C for 30 min, enrichment of microsatellites using biotinylated oligonucleotides and magnetic field, and directional cloning of cDNA into a vector. We tested this method to enrich CA- and GA-microsatellites from cDNA of Asian seabass, and demonstrated that enrichment of microsatellites from normalized cDNA could increased the efficiency of microsatellite isolation over 30 times as compared to direct sequencing of clones from cDNA libraries. One hundred and thirty-nine (36.2%) out of 384 clones from normalized cDNA contained microsatellites. Unique microsatellite sequences accounted for 23.6% (91/384) of sequenced clones. Sixty microsatellites isolated from cDNA were characterized, and 41 were polymorphic. The average allele number of the 41 microsatellites was 4.85 ± 0.54, while the expected heterozygosity was 0.56 ± 0.03. All the isolated microsatellites inherited in a Mendelian pattern.

Conclusion

Normalization of cDNA substantially increased the efficiency of enrichment of microsatellites from cDNA. The described method for isolation of microsatellites from cDNA has the potential to be applied to a wide range of fish species. The microsatellites isolated from cDNA could be useful for linkage and comparative mapping, as well as for studying genome evolution.

MytiBase: a knowledgebase of mussel (M. galloprovincialis) transcribed sequences

BACKGROUND

Although bivalves are among the most studied marine organisms due to their ecological role, economic importance and use in pollution biomonitoring, very little information is available on the genome sequences of mussels. This study reports the functional analysis of a large-scale Expressed Sequence Tag (EST) sequencing from different tissues of Mytilus galloprovincialis (the Mediterranean mussel) challenged with toxic pollutants, temperature and potentially pathogenic bacteria.

RESULTS

We have constructed and sequenced seventeen cDNA libraries from different Mediterranean mussel tissues: gills, digestive gland, foot, anterior and posterior adductor muscle, mantle and haemocytes. A total of 24,939 clones were sequenced from these libraries generating 18,788 high-quality ESTs which were assembled into 2,446 overlapping clusters and 4,666 singletons resulting in a total of 7,112 non-redundant sequences. In particular, a high-quality normalized cDNA library (Nor01) was constructed as determined by the high rate of gene discovery (65.6%). Bioinformatic screening of the non-redundant M. galloprovincialis sequences identified 159 microsatellite-containing ESTs. Clusters, consensuses, related similarities and gene ontology searches have been organized in a dedicated, searchable database http://mussel.cribi.unipd.it.

CONCLUSION

We defined the first species-specific catalogue of M. galloprovincialis ESTs including 7,112 unique transcribed sequences. Putative microsatellite markers were identified. This annotated catalogue represents a valuable platform for expression studies, marker validation and genetic linkage analysis for investigations in the biology of Mediterranean mussels.

Analysis of genes isolated from plated hemocytes of the Pacific oyster, Crassostreas gigas

A complementary deoxyribonucleic acid library was constructed from hemocytes of Crassostrea gigas that had been plated on poly-lysine plates for 24 h. From this library, 2,198 expressed sequence tags (ESTs) of greater than or equal to 100 bp were generated and analyzed. A large number of genes that potentially could be involved in the physiology of the oyster hemocyte were uncovered. They included proteins involved in cytoskeleton rearrangement, proteases and antiproteases, regulators of transcription and translation, cell death regulators, receptors and their associated protein factors, lectins, signal transduction proteins, and enzymes involved in eicosanoid and steroid synthesis and xenobiotic metabolism. Based on their relationship with innate immunity, the expression of selected genes was analyzed by quantitative polymerase chain reaction in gills from bacterial-challenged oysters. Several genes observed in the library were significantly upregulated by bacterial challenge including interleukin 17, astacin, cystatin B, the EP4 receptor for prostaglandin E, the ectodysplasin receptor, c-jun, and the p100 subunit of nuclear factor-kB. Using a similar approach, we have been analyzing the genes expressed in trout macrophages. While there are significant differences between the types of genes present in vertebrate macrophages compared with oyster hemocytes, there are some striking similarities including proteins involved in cytoskeletal rearrangement, proteases and antiproteases, and genes involved in certain signal transduction pathways underlying immune processes such as phagocytosis. Finally, C. virginica homologs of some of the C. gigas genes uncovered in the ESTs were obtained by aligning the ESTs reported here, against the assembled C. virginica ESTs at the National Center for Biotechnology Information.

Expressed sequence tags (ESTs) from immune tissues of turbot (Scophthalmus maximus) challenged with pathogens

Background

The turbot (Scophthalmus maximus; Scophthalmidae; Pleuronectiformes) is a flatfish species of great relevance for marine aquaculture in Europe. In contrast to other cultured flatfish, very few genomic resources are available in this species. Aeromonas salmonicida and Philasterides dicentrarchi are two pathogens that affect turbot culture causing serious economic losses to the turbot industry. Little is known about the molecular mechanisms for disease resistance and host-pathogen interactions in this species. In this work, thousands of ESTs for functional genomic studies and potential markers linked to ESTs for mapping (microsatellites and single nucleotide polymorphisms (SNPs)) are provided. This information enabled us to obtain a preliminary view of regulated genes in response to these pathogens and it constitutes the basis for subsequent and more accurate microarray analysis.

Results

A total of 12584 cDNAs partially sequenced from three different cDNA libraries of turbot (Scophthalmus maximus) infected with Aeromonas salmonicida, Philasterides dicentrarchi and from healthy fish were analyzed. Three immune-relevant tissues (liver, spleen and head kidney) were sampled at several time points in the infection process for library construction. The sequences were processed into 9256 high-quality sequences, which constituted the source for the turbot EST database. Clustering and assembly of these sequences, revealed 3482 different putative transcripts, 1073 contigs and 2409 singletons. BLAST searches with public databases detected significant similarity (e-value ≤ 1e-5) in 1766 (50.7%) sequences and 816 of them (23.4%) could be functionally annotated. Two hundred three of these genes (24.9%), encoding for defence/immune-related proteins, were mostly identified for the first time in turbot. Some ESTs showed significant differences in the number of transcripts when comparing the three libraries, suggesting regulation in response to these pathogens. A total of 191 microsatellites, with 104 having sufficient flanking sequences for primer design, and 1158 putative SNPs were identified from these EST resources in turbot.

Conclusion

A collection of 9256 high-quality ESTs was generated representing 3482 unique turbot sequences. A large proportion of defence/immune-related genes were identified, many of them regulated in response to specific pathogens. Putative microsatellites and SNPs were identified. These genome resources constitute the basis to develop a microarray for functional genomics studies and marker validation for genetic linkage and QTL analysis in turbot.

Normalization of full-length enriched cDNA

Analysis of rare messages in cDNA libraries is extremely difficult due to the substantial variations in the abundance of different transcripts in cells and tissues. Therefore, for rare transcript searches and analyses, the generation of equalized (normalized) cDNA is essential. Several cDNA normalization methods have been developed since 1990. A number of these methods have been optimized for the normalization of full-length enriched cDNA, and used in various applications, including transcriptome analysis and functional screening of cDNA libraries. One such procedure (named DSN-normalization) is based on the unique properties of duplex-specific nuclease (DSN) from kamchatka crab and allows the generation of normalized cDNA libraries with a high gene discovery rate.

Generation and analysis of expressed sequence tags from the ciliate protozoan parasite Ichthyophthirius multifiliis

BACKGROUND

The ciliate protozoan Ichthyophthirius multifiliis (Ich) is an important parasite of freshwater fish that causes 'white spot disease' leading to significant losses. A genomic resource for large-scale studies of this parasite has been lacking. To study gene expression involved in Ich pathogenesis and virulence, our goal was to generate expressed sequence tags (ESTs) for the development of a powerful microarray platform for the analysis of global gene expression in this species. Here, we initiated a project to sequence and analyze over 10,000 ESTs.

RESULTS

We sequenced 10,368 EST clones using a normalized cDNA library made from pooled samples of the trophont, tomont, and theront life-cycle stages, and generated 9,769 sequences (94.2% success rate). Post-sequencing processing led to 8,432 high quality sequences. Clustering analysis of these ESTs allowed identification of 4,706 unique sequences containing 976 contigs and 3,730 singletons. These unique sequences represent over two million base pairs (~10% of Plasmodium falciparum genome, a phylogenetically related protozoan). BLASTX searches produced 2,518 significant (E-value < 10-5) hits and further Gene Ontology (GO) analysis annotated 1,008 of these genes. The ESTs were analyzed comparatively against the genomes of the related protozoa Tetrahymena thermophila and P. falciparum, allowing putative identification of additional genes. All the EST sequences were deposited by dbEST in GenBank (GenBank: EG957858-EG966289). Gene discovery and annotations are presented and discussed.

CONCLUSION

This set of ESTs represents a significant proportion of the Ich transcriptome, and provides a material basis for the development of microarrays useful for gene expression studies concerning Ich development, pathogenesis, and virulence.