Development of Expressed Sequence Tags from the Pearl Oyster, <Emphasis Type="Italic">Pinctada martensii</Emphasis> 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.     

Generation and analysis of 10 000 ESTs from the half‐smooth tongue sole Cynoglossus semilaevis and identification of microsatellite and SNP markers

Three normalized cDNA libraries were constructed, two of which were constructed from reproductive tissues ovary and testis, and the other one from pooled immune tissues including head kidney, intestine, liver and spleen. A total of 10 542 clones were sequenced generating 10 128 expressed sequence tags (ESTs). Cluster analysis indicated a total of 5808 unique sequences including 1712 contigs and 4096 singletons. A total of 4249 (73%) of the unique ESTs had significant hits to the non‐redundant protein database, 2253 of which were annotated using Gene Ontology (GO) terms. A total of 311 microsatellites (with 246 having sufficient flanking sequences for primer design) and 6294 putative SNPs were identified. These genome resources provide the material basis for future microarray development, marker validation and genetic linkage and QTL analysis.     

Immune gene discovery by expressed sequence tag (EST) analysis of hemocytes in the ridgetail white prawn Exopalaemon carinicauda

The ridgetail white prawn Exopalaemon carinicauda is one of the most important commercial species in eastern China. However, little information of immune genes in E. carinicauda has been reported. To identify distinctive genes associated with immunity, an expressed sequence tag (EST) library was constructed from hemocytes of E. carinicauda. A total of 3411 clones were sequenced, yielding 2853 ESTs and the average sequence length is 436 bp. The cluster and assembly analysis yielded 1053 unique sequences including 329 contigs and 724 singletons. Blast analysis identified 593 (56.3%) of the unique sequences as orthologs of genes from other organisms (E-value < 1e-5). Based on the COG and Gene Ontology (GO), 593 unique sequences were classified. Through comparison with previous studies, 153 genes assembled from 367 ESTs have been identified as possibly involved in defense or immune functions. These genes are categorized into seven categories according to their putative functions in shrimp immune system: antimicrobial peptides, prophenoloxidase activating system, antioxidant defense systems, chaperone proteins, clottable proteins, pattern recognition receptors and other immune-related genes. According to EST abundance, the major immune-related genes were thioredoxin (141, 4.94% of all ESTs) and calmodulin (14, 0.49% of all ESTs). The EST sequences of E. carinicauda hemocytes provide important information of the immune system and lay the groundwork for development of molecular markers related to disease resistance in prawn species.     

EST-based identification of immune-relevant genes from spleen of Indian catfish, Clarias batrachus (Linnaeus, 1758)

A normalized cDNA library from spleen of Indian catfish, Clarias batrachus, was constructed with a redundancy factor of 2.29. A total of 2045 clones from the library were single-pass sequenced, which generated 1937 high quality ESTs with an average read length of approximately 700 bp. Based on sequence similarities, 65 ESTs were found to be associated with immune functions, which were mainly associated with response to stress, response to chemical stimulus, cellular response to stimulus, response to external stimulus, immune response and regulation of response to stimulus. The immune-relevant gene for CD141, thrombomodulin, has been identified in Teleosts for the first time. Six EST-SSRs and three SNPs were found associated with eight immune-relevant genes. These markers associated with important immune genes would be useful for the identification of trait associated alleles for marker-assisted selection. The identification of the putative immune-related genes provides a meaningful framework to understand the Indian catfish immune system and defense mechanisms.     

Development and mapping of EST-derived simple sequence repeat markers for hexaploid wheat.

Expressed sequence tags (ESTs) are a valuable source of molecular markers. To enhance the resolution of an existing linkage map and to identify putative functional polymorphic gene loci in hexaploid wheat (Triticum aestivum L.), over 260,000 ESTs from 5 different grass species were analyzed and 5418 SSR-containing sequences were identified. Using sequence similarity analysis, 156 cross-species superclusters and 138 singletons were used to develop primer pairs, which were then tested on the genomic DNA of barley (Hordeum vulgare), maize (Zea mays), rice (Oryza sativa), and wheat. Three-hundred sixty-eight primer pairs produced PCR amplicons from at least one species and 227 primer pairs amplified DNA from two or more species. EST-SSR sequences containing dinucleotide motifs were significantly more polymorphic (74%) than those containing trinucleotides (56%), and polymorphism was similar for markers in both coding and 5' untranslated (UTR) regions. Out of 112 EST-SSR markers, 90 identified 149 loci that were integrated into a reference wheat genetic map. These loci were distributed on 19 of the 21 wheat chromosomes and were clustered in the distal chromosomal regions. Multiple-loci were detected by 39% of the primer pairs. Of the 90 mapped ESTs, putative functions for 22 were identified using BLASTX queries. In addition, 80 EST-SSR markers (104 loci) were located to chromosomes using nullisomic-tetrasomic lines. The enhanced map from this study provides a basis for comparative mapping using orthologous and PCR-based markers and for identification of expressed genes possibly affecting important traits in wheat.     

Development of EST SSRs in Finger Millet (Eleusine coracana ssp coracana) and their Transferability to Pearl Millet (Pennisetum glaucum)

EST-SSR markers were developed using sequence information from 1740 expressed sequence tags (ESTs) of finger millet available in the public domain. A set of 31 SSR markers were synthesized based on di, tri, tetra and penta-nucleotide repeat sequences. These were used for PCR analysis of 11 elite germplasm lines of finger millet of Indian and African origin. Out of 31 SSR markers, amplification products were obtained for 17 primer pairs. Of these nine were found polymorphic with two alleles per locus. These 17 SSR primer pairs were also tested for amplification in three varieties of pearl millet (Pennisetum glaucum) and 11 could be transferred to pearl millet. The informative EST SSR markers developed, can be used in finger millet as well as pearl millet genetic improvement projects.     

Simple sequence repeat marker development from bacterial artificial chromosome end sequences and expressed sequence tags of flax (Linum usitatissimum L.)

Flax is an important oilseed crop in North America and is mostly grown as a fibre crop in Europe. As a self-pollinated diploid with a small estimated genome size of ~370 Mb, flax is well suited for fast progress in genomics. In the last few years, important genetic resources have been developed for this crop. Here, we describe the assessment and comparative analyses of 1,506 putative simple sequence repeats (SSRs) of which, 1,164 were derived from BAC-end sequences (BESs) and 342 from expressed sequence tags (ESTs). The SSRs were assessed on a panel of 16 flax accessions with 673 (58 %) and 145 (42 %) primer pairs being polymorphic in the BESs and ESTs, respectively. With 818 novel polymorphic SSR primer pairs reported in this study, the repertoire of available SSRs in flax has more than doubled from the combined total of 508 of all previous reports. Among nucleotide motifs, trinucleotides were the most abundant irrespective of the class, but dinucleotides were the most polymorphic. SSR length was also positively correlated with polymorphism. Two dinucleotide (AT/TA and AG/GA) and two trinucleotide (AAT/ATA/TAA and GAA/AGA/AAG) motifs and their iterations, different from those reported in many other crops, accounted for more than half of all the SSRs and were also more polymorphic (63.4 %) than the rest of the markers (42.7 %). This improved resource promises to be useful in genetic, quantitative trait loci (QTL) and association mapping as well as for anchoring the physical/genetic map with the whole genome shotgun reference sequence of flax.     

Bioinformatic Mining of Type I Microsatellites from Expressed Sequence Tags of Channel Catfish (<Emphasis Type="Italic">Ictalurus punctatus</Emphasis>)

Gene-derived markers are pivotal to the analysis of genome structure, organization, and evolution and necessary for comparative genomics. However, gene-derived markers are relatively difficult to develop. This project utilized the genomic resources of channel catfish expressed sequence tags (ESTs) to identify simple sequence repeats (SSRs), or microsatellites. It took the advantage of ESTs for the establishment of gene identities, and of microsatellites for the acquisition of high polymorphism. When microsatellites are tagged to genes, the microsatellites can then be used as gene markers. A bioinformatic analysis of 43,033 ESTs identified 4855 ESTs containing microsatellites. Cluster analysis indicated that 1312 of these ESTs fell into 569 contigs, and the remaining 3534 ESTs were singletons. A total of 4103 unique microsatellite-containing genes were identified. The dinucleotide CA/TG and GA/TC pairs were the most abundant microsatellites. AT-rich microsatellite types were predominant among trinucleotide and tetranucleotide microsatellites, consistent with our earlier estimation that the catfish genome is highly AT-rich. Our preliminary results indicated that the majority of the identified microsatellites were polymorphic and, therefore, useful for genetic linkage mapping of catfish. Mapping of these gene-derived markers is under way, which will set the foundation for comparative genome analysis in catfish.     

Potential Indicators of Stress Response Identified by Expressed Sequence Tag Analysis of Hemocytes and Embryos from the American Oyster, Crassostrea virginica

A pilot program was initiated to identify genes from the American oyster, Crassostrea virginica, that are potentially involved in the stress response for use as bioindicators of exposure to environmental pollutants and to toxic and infectious agents. A PCR-based method was used to construct cDNA libraries from pooled embryos and the hemocytes of a single individual. A total of 998 randomly selected clones (expressed sequence tags, ESTs) were sequenced. Approximately 40% of the ESTs are novel sequences. Several potential biomarkers identified include an antimicrobial peptide, recognition molecules (lectin receptors), proteinases and proteinase inhibitors, and a novel metallothionein. Diversity analysis shows that 363 and 286 unique genes were identified from the hemocyte and embryo libraries, respectively, indicating that full-scale EST collection is a valuable approach for the discovery of new genes of potential significance in the molluscan stress response. Received May 11, 2001; accepted July 16, 2001     

Variation in the <Emphasis Type="Italic">COI</Emphasis> gene of the freshwater pearl mussel <Emphasis Type="Italic">Margaritifera margaritifera</Emphasis> from River Vuokkijoki

The freshwater pearl mussel Margaritifera margaritifera L. is one of the most endangered freshwater mussels in the world. Effective conservation of threatened species requires not only ecological, but also genetic information from the target species and populations. Since low genetic diversity can reduce the ability of a species to adapt to environmental changes, maintaining genetic diversity has been identified as one of the key elements in successful conservation programs. We examined genetic variation of the freshwater pearl mussel from the River Vuokkijoki, Karelia, Russia. We sequenced a fragment of the cytochrome c oxidase subunit I gene (COI) from 22 individuals and compared the data to 32 previously published COI sequences available in GenBank. We identified 10 different COI haplotypes in the sequenced samples, three of which had not been previously reported. Our results show that the River Vuokkijoki has high genetic diversity and suggest that the colonization of this northern freshwater pearl mussel population might have occurred from multiple and even distant refugia. Therefore, the freshwater pearl mussel population of the River Vuokkijoki is valuable for the conservation of the whole species.     

Mining and charactering microsatellites in Cucumis melo expressed sequence tags from sequence database

Simple sequence repeats (SSRs) derived from expressed sequence tags (ESTs) are valuable markers because they represent transcribed regions and often have putative functions. We mined and characterized microsatellites in melon ESTs. Three hundred and eighty‐three SSR loci were identified in 309 of 3188 unigenes assembled by 5747 EST and mRNA sequences in GenBank with occurring frequency of 1/4.7 kb. Twenty‐two polymorphic EST‐SSR markers were developed with the mean allele number of 2.9 per locus and mean expected heterozygosity of 0.442. Amplification products were also detected by 15 pairs of primer in Cucumis sativus. Those informative EST‐SSR markers can be used in melon genetic improvement projects.     

Development of Expressed Sequence Tags from the Bay Scallop, Argopecten irradians irradians

The bay scallop, Argopecten irradians irradians, introduced from North America, has become one of the most important aquaculture species in China. Inan effort to identify scallop genes involved in host defense, a high-quality cDNA library was constructed from whole body tissues of the bay scallop. A total of 5828 successful sequencing reactions yielded 4995 expressed sequence tags (ESTs) longer than 100 bp. Cluster and assembly analyses of the ESTs identified 637 contigs (consisting of 2853 sequences) and 2142 singletons, totaling 2779 unique sequences. Basic Local Alignment Search Tool (BLAST) analysis showed that the majority (73%) of the unique sequences had no significant homology (E-value ≤ 0.005) to sequences in GenBank. Among the 748 sequences with significant GenBank matches, 160 (21.4%) were for genes related to metabolism, 131 (17.5%) for cell/organism defense, 124 (16.6%) for gene/protein expression, 83 (11.1%) for cell structure/motility, 70 (9.4%) for cell signaling/communication, 17 (2.3%) for cell division, and 163 (21.8%) matched to genes of unknown functions. The list of host-defense genes included many genes with known and important roles in innate defense such as lectins, defensins, proteases, protease inhibitors, heat shock proteins, antioxidants, and Toll-like receptors. The study provides a significant number of ESTs for gene discovery and candidate genes for studying host defense in scallops and other molluscs.     

嗜水气单胞菌诱导的池蝶蚌血细胞cDNA 文库的构建和亲环蛋白基因序列的初步分析

实验利用灭活的嗜水气单胞菌(Aeromonas hydrophila)诱导处于四龄池蝶蚌(Hyriopsis schlegelii)14h,将诱导后的池蝶蚌血细胞的总RNA 进行逆转录, 用LD-PCR 法合成双链cDNA, 从而首次成功构建池蝶蚌血细胞的全长cDNA 文库。原始文库的滴度为4106 cfu/cm3, 重组率为90%, 扩增后文库的滴度为3.55109pfu/mL。目前文库已随机测序672 个样品, 将所得双向序列进行拼接, 去除载体, 并多序列比对去除重复序列后, 发现436 条为已知功能序列, 其余为未知功能序列。序列中最小长度270 bp, 最大长度为2153 bp, 平均大小608.6 bp, 表明插入片段大小理想。从文库中筛选获得免疫相关基因池蝶蚌亲环蛋白A(HsCyp A)全长基因并进行序列分析。结果显示, HsCyp A 全长1229 bp, 序列包括52 bp 的5非编码区、495 bp 的开放阅读框、682 bp 的3非编码区和29 bp 的poly(A)尾, 没有明显的加尾信号。对Cyp A 氨基酸序列二级结构进行了较详细的分析并进行了三维建模, 同时构建了其系统进化树, 分析表明亲环蛋白家族是一个在进化上非常保守的蛋白家族。综合分析, Cyp A 在水生动物中不仅仅只是一种组成型蛋白, 而是可能在病原感染防御中发挥重要作用。       

EST derived PCR-based markers for functional gene homologues in cotton.

We investigated the utility of the Gossypium arboreum EST sequences in the GenBank database for developing PCR-based markers targeting known-function genes in cultivated tetraploid cottons, G. hirsutum and G. barbadense. Four hundred sixty-five randomly selected ESTs from this library were subjected to BLASTn search against all GenBank databases, of which putative function was assigned to 93 ESTs based on high nucleotide homology to previously studied genes. PCR primers were synthesized for 89 of the known-function ESTs. A total of 57 primer pairs amplified G. arboreum genomic DNA, but only 39 amplified in G. hirsutum and G. barbadense, suggesting that sequence divergence may be a factor causing non-amplification for some sites. DNA sequence analysis showed that most primer pairs were targeting the expected homologous loci. While the amplified products that were of larger size than the corresponding EST sequences contain introns, the primer pairs with a smaller amplicon than predicted from the flanking EST sequences did not amplify the expected orthologous gene sequences. Among the 39 primer pairs that amplified tetraploid cotton DNA, 3 detected amplicon size polymorphisms and 10 detected polymorphisms after digestion with one of six restriction enzymes. Ten of the polymorphic loci were subsequently mapped to an anchor RFLP map. Digestion of PCR-amplified sequences offers one means by which cotton genes can be mapped to their chromosomal locations more quickly and economically than by RFLP analysis.     

Development and mapping of functional expressed sequence tag-derived simple sequence repeat markers in a rubber tree RRIM600 × PB217 population

Sets of polymorphic expressed sequence tag–simple sequence repeat (EST-SSR) markers from the rubber tree (Hevea brasiliensis) have been published by many researchers, but none has been specifically developed to study latex and wood yield traits. In this study, a total 10,321 rubber tree EST sequences, generated from suppression subtractive hybridization-cDNA libraries of bark and latex of high- and low-yielding clones, were used as sources for SSR searching. A total of 432 EST-SSR loci were identified and it was possible to design primer pairs for a subset of 298 EST-SSRs. The highest proportion of EST-SSRs was represented by dinucleotide repeats (46.6 %), followed by trinucleotide repeats (44.3 %). Based on BLASTX analysis, 234 ESTs (80 %) showed similarity to genes in NCBI databases and could be divided into 120 putative proteins with known function and 114 unknown proteins. To enhance the resolution of an existing linkage map from previous work on a rubber tree RRIM600 × PB217 population, 69 EST-SSR markers from the above set were tested to be integrated into the reference genetic map. The enriched map of 18 linkage groups spanned 2054.2 cM in length, showed an average genetic distance of 4.3 cM between adjacent markers, and included 63 new EST-SSR markers. The enhanced map from this study provides a basis for comparative mapping using PCR-based markers and identification of expressed genes possibly affecting important traits of interest.