DNA sequence-based "bar codes" for tracking the origins of expressed sequence tags from a maize cDNA library constructed using multiple mRNA sources

To enhance gene discovery, expressed sequence tag (EST) projects often make use of cDNA libraries produced using diverse mixtures of mRNAs. As such, expression data are lost because the origins of the resulting ESTs cannot be determined. Alternatively, multiple libraries can be prepared, each from a more restricted source of mRNAs. Although this approach allows the origins of ESTs to be determined, it requires the production of multiple libraries. A hybrid approach is reported here. A cDNA library was prepared using 21 different pools of maize (Zea mays) mRNAs. DNA sequence "bar codes" were added during first-strand cDNA synthesis to uniquely identify the mRNA source pool from which individual cDNAs were derived. Using a decoding algorithm that included error correction, it was possible to identify the source mRNA pool of more than 97% of the ESTs. The frequency at which a bar code is represented in an EST contig should be proportional to the abundance of the corresponding mRNA in the source pool. Consistent with this, all ESTs derived from several genes (zein and adh1) that are known to be exclusively expressed in kernels or preferentially expressed under anaerobic conditions, respectively, were exclusively tagged with bar codes associated with mRNA pools prepared from kernel and anaerobically treated seedlings, respectively. Hence, by allowing for the retention of expression data, the bar coding of cDNA libraries can enhance the value of EST projects.     

Evaluation of cDNA Libraries from Different Developmental Stages of Schistosoma mansoni for Production of Expressed Sequence Tags (ESTs)

A comparative study of the gene expression profile in differentdevelopmental stages of Schistosoma mansoni has been initiatedbased on the expressed sequence tag(EST) approach. A total of1401 ESTs were generated from seven different cDNA librariesconstructed from four distinct stages of the parasite life cycle.The libraries were first evaluated for their quality for a large-scalecDNA sequencing program. Most of them were shown to have lessthan 20% useless clones and more than 50% new genes. The redundancyof each library was also analyzed, showing that one adult wormcDNA library was composed of a small number of highly frequentgenes. When comparing ESTs from distinct libraries, we coulddetect that most genes were present only in a single library,but others were expressed in more than one developmental stageand may represent housekeeping genes in the parasite. When consideringonly once the genes present in more than one library, a totalof 466 unique genes were obtained, corresponding to 427 newS. mansoni genes. From the total of unique genes, 20.2% wereidentified based on homology with genes from other organisms,8.3% matched S. mansoni characterized genes and 71.5% representunknown genes.     

Analysis of Medicago truncatula nodule expressed sequence tags

Systematic sequencing of expressed sequence tags (ESTs) can give a global picture of the assembly of genes involved in the development and function of organs. Indeterminate nodules representing different stages of the developmental program are especially suited to the study of organogenesis. With the vector lambdaHybriZAP, a cDNA library was constructed from emerging nodules of Medicago truncatula induced by Sinorhizobium meliloti. The 5' ends of 389 cDNA clones were sequenced, then these ESTs were analyzed both by sequence homology search and by studying their expression in roots and nodules. Two hundred fifty-six ESTs exhibited significant similarities to characterized data base entries and 40 of them represented 26 nodulin genes, while 133 had no similarity to sequences with known function. Only 60 out of the 389 cDNA clones corresponded to previously submitted M. truncatula EST sequences. For 117 cDNAs, reverse Northern (RNA) hybridization with root and nodule RNA probes revealed enhanced expression in the nodule, 48 clones are likely to code for novel nodulins, 33 cDNAs are clones of already known nodulin genes, and 36 clones exhibit similarity to other characterized genes. Thus, systematic analysis of the EST sequences and their expression patterns is a powerful way to identify nodule-specific and nodulation-related genes.     

Gene expression profile of human bone marrow stromal cells: high-throughput expressed sequence tag sequencing analysis

Human bone marrow stromal cells (HBMSC) are pluripotent cells with the potential to differentiate into osteoblasts, chondrocytes, myelosupportive stroma, and marrow adipocytes. We used high-throughput DNA sequencing analysis to generate 4258 single-pass sequencing reactions (known as expressed sequence tags, or ESTs) obtained from the 5' (97) and 3' (4161) ends of human cDNA clones from a HBMSC cDNA library. Our goal was to obtain tag sequences from the maximum number of possible genes and to deposit them in the publicly accessible database for ESTs (dbEST of the National Center for Biotechnology Information). Comparisons of our EST sequencing data with nonredundant human mRNA and protein databases showed that the ESTs represent 1860 gene clusters. The EST sequencing data analysis showed 60 novel genes found only in this cDNA library after BLAST analysis against 3.0 million ESTs in NCBI's dbEST database. The BLAST search also showed the identified ESTs that have close homology to known genes, which suggests that these may be newly recognized members of known gene families. The gene expression profile of this cell type is revealed by analyzing both the frequency with which a message is encountered and the functional categorization of expressed sequences. Comparing an EST sequence with the human genomic sequence database enables assignment of an EST to a specific chromosomal region (a process called digital gene localization) and often enables immediate partial determination of intron/exon boundaries within the genomic structure. It is expected that high-throughput EST sequencing and data mining analysis will greatly promote our understanding of gene expression in these cells and of growth and development of the skeleton.     

Comparison of gene representation between diploid and triploid Paragonimus westermani by expressed sequence tag analyses

Expressed sequence tag (EST) analysis of the diploid and triploid Paragonimus westermani genes was done to have a rapid and informative outlook of the gene-expression profiles of the parasites. Totals of 506 and 505 ESTs were generated from the diploid and triploid P. westermani cDNA libraries. Based on the BLASTx search results of the diploid P. westermani ESTs, 308 (60.9%) matched significantly with formerly identified genes and 198 (39.1%) showed no significant homology in the GenBank database. A similar homology pattern was shown from the triploid EST BLASTx search results with 346 (68.5%) sharing homology with previously identified genes and 159 (31.5%) showing no significant homology. The EST data from both libraries were analyzed and grouped into 9 categories. Comparison of the 2 EST pools revealed high similarities among the categories of the significantly matched genes. Single genes matched repeatedly were also observed in the 2 EST data. Some genes were found that are not yet characterized in P. westermani; these genes were matched by both the diploid and triploid ESTs. Further study of these genes may provide us with more understanding on the parasite's biology and their specific functions in the 2 strains.     

Expressed sequence tags and their application for plant research

Expressed Sequence Tags (ESTs) are short, usually unedited sequences obtained by single-pass sequencing of cDNA clones from any cDNA library. Analyzing and comparing ESTs can provide information on gene expression, function and evolution. Large-scale EST sequencing has become an attractive alternative to plant genome sequencing. Currently, plant EST collections comprise over 3.8 million sequences from about 200 species. They have proved to be a valuable tool for gene discovery and plant metabolism analysis. Several plant-specific EST databases have been created which provide access to sequence data and bioinformatics-based tools for data mining. Searching EST collections allows pre-selection of genes for preparing cDNA arrays, targeted to bring maximum information on specialized processes, like stress response, symbiotic nitrogen fixation etc. Also, ESt-based molecular markers such as SNP, SSR, and indels are fast developing tools for breeders and researchers.     

Using mtDNA Sequences to Estimate SNP Parameters in ESTs

Discovery of single nucleotide polymorphisms (SNPs) requires analysis of redundant sequences such as those available in large public databases. The ability to detect SNPs, especially those of low frequency, is dependent on the depth and scale of the discovery effort. Large numbers of SNPs have been identified by mining large-scale EST surveys and whole genome sequencing projects. These surveys however are subject to ascertainment bias and the inherent errors in large-scale single pass sequencing efforts. For example, the number of steps involved in the construction and sequencing of cDNA libraries make ESTs highly error prone, resulting in an increased frequency of nonvalid SNPs obtained in these surveys. Sequences of mtDNA genes are often incorporated into cDNA libraries as an artifact of the library construction process and are typically either subtracted from cDNA libraries or are considered superfluous when evaluating the information content of EST datasets. Sequences of mtDNA genes provide a unique resource for the analysis of SNP parameters in EST projects. This study uses sequences from four turkey muscle cDNA libraries to demonstrate how mtDNA sequences gleaned from collections of ESTs can be used to estimate SNP parameters and thus help predict the validity of SNPs.     

Using mtDNA sequences to estimate SNP parameters in ESTs

Discovery of single nucleotide polymorphisms (SNPs) requires analysis of redundant sequences such as those available in large public databases. The ability to detect SNPs, especially those of low frequency, is dependent on the depth and scale of the discovery effort. Large numbers of SNPs have been identified by mining large-scale EST surveys and whole genome sequencing projects. These surveys however are subject to ascertainment bias and the inherent errors in large-scale single pass sequencing efforts. For example, the number of steps involved in the construction and sequencing of cDNA libraries make ESTs highly error prone, resulting in an increased frequency of nonvalid SNPs obtained in these surveys. Sequences of mtDNA genes are often incorporated into cDNA libraries as an artifact of the library construction process and are typically either subtracted from cDNA libraries or are considered superfluous when evaluating the information content of EST datasets. Sequences of mtDNA genes provide a unique resource for the analysis of SNP parameters in EST projects. This study uses sequences from four turkey muscle cDNA libraries to demonstrate how mtDNA sequences gleaned from collections of ESTs can be used to estimate SNP parameters and thus help predict the validity of SNPs.     

Correlated clustering and virtual display of gene expression patterns in the wheat life cycle by large-scale statistical analyses of expressed sequence tags

Compared to rice, wheat exhibits characteristic growth habits and contains complex genome constituents. To assess global changes in gene expression patterns in the wheat life cycle, we conducted large-scale analysis of expressed sequence tags (ESTs) in common wheat. Ten wheat tissues were used to construct cDNA libraries: crown and root from 14-day-old seedlings; spikelet from early and late flowering stages; spike at the booting stage, heading date and flowering date; pistil at the heading date; and seeds at 10 and 30 days post-anthesis. Several thousand colonies were randomly selected from each of these 10 cDNA libraries and sequenced from both 5' and 3' ends. Consequently, a total of 116 232 sequences were accumulated and classified into 25 971 contigs based on sequence homology. By computing abundantly expressed ESTs, correlated expression patterns of genes across the tissues were identified. Furthermore, relationships of gene expression profiles among the 10 wheat tissues were inferred from global gene expression patterns. Genes with similar functions were grouped with one another by clustering gene expression profiles. This technique might enable estimation of the functions of anonymous genes. Multidimensional analysis of EST data that is analogous to the microarray experiments may offer new approaches to functional genomics of plants.     

Gene identification and expression analysis of 86,136 Expressed Sequence Tags (EST) from the rice genome

Expressed Sequence Tag (EST) analysis has pioneered genome-wide gene discovery and expression profiling. In order to establish a gene expression index in the rice cultivar indica, we sequenced and analyzed 86,136 ESTs from nine rice cDNA libraries from the super hybrid cultivar LYP9 and its parental cultivars. We assembled these ESTs into 13,232 contigs and leave 8,976 singletons. Overall, 7,497 sequences were found similar to the existing sequences in GenBank and 14,711 are novel. These sequences are classified by molecular function, biological process and pathways according to the Gene Ontology. We compared our sequenced ESTs with the publicly available 95,000 ESTs from japonica, and found little sequence variation, despite the large difference between genome sequences. We then assembled the combined 173,000 rice ESTs for further analysis. Using the pooled ESTs, we compared gene expression in metabolism pathway between rice and Arabidopsis according to KEGG. We further profiled gene expression pattern     

Tissue expression map of a large number of expressed sequence tags and its application to in silico screening of stress response genes in common wheat

In order to assess global changes in gene expression patterns in stress-induced tissues, we conducted large-scale analysis of expressed sequence tags (ESTs) in common wheat. Twenty-one cDNA libraries derived from stress-induced tissues, such as callus, as well as liquid cultures and abiotic stress conditions (temperature treatment, desiccation, photoperiod, moisture and ABA) were constructed. Several thousand colonies were randomly selected from each of these 21 cDNA libraries and sequenced from both the 5′ and 3′ ends. By computing abundantly expressed ESTs, correlated expression patterns of genes across the tissues were monitored. Furthermore, the relationships between gene expression profiles among the stress-induced tissues were inferred from the gene expression patterns. Multi-dimensional analysis of EST data is analogous to microarray experiments. As an example, genes specifically induced and/or suppressed by cold acclimation and heat-shock treatments were selected in silico. Four hundred and ninety genes showing fivefold induction or 218 genes for suppression in comparison to the control expression level were selected. These selected genes were annotated with the BLAST search. Furthermore, gene ontology was conducted for these genes with the InterPro search. Because genes regulated in response to temperature treatment were successfully selected, this method can be applied to other stress-treated tissues. Then, the method was applied to screen genes in response to abiotic stresses such as drought and ABA treatments. In silico selection of screened genes from virtual display should provide a powerful tool for functional plant genomics.Electronic Supplementary Material Supplementary material is available to authorised users in the online version of this article at .     

一个新的睾丸特异表达基因的克隆及功能的初步分析

运用“数据库消减杂交”(digital differential display)方法来筛选人类睾丸特异表达新基因,获得了有差异显示的代表新基因的克隆重叠群。挑选其中一个克隆重叠群HS．326528进行多组织RT—PCR,初步获得该重叠群在人睾丸中有高表达。从该重叠群的IMAGE出发,采用生物信息学的方法快速克隆了一个人类新基因的全长cDNA序列,其全长1044bp,开放阅读框为214～529bp,定位于15q26．2,编码由105个氨基酸组成、分子量为11．7kD、等电点为10．09的一个碱性蛋白,该蛋白与已知蛋白无明显的同源性,克隆实验证明该基因的阅读框完全正确,RT—PCR和Northern blot显示该基因在人类睾丸中特异表达,实时PCR结果表明：该基因在成人睾丸中高表达,在精子中有中度表达,在胚胎睾丸中低表达,推测该基因与精子的生成有关,命名为SRG8(homo sapiens spermatogenesis—related gene 8)(GenBank登录号：AY489187),该基因编码的蛋白定位于细胞核。流式结果分析表明,SRG8基因能够促使HeLa细胞由S期向G2期的转变,从而加速细胞的分裂。这些结果表明SRG8基因可能在睾丸的发育及精子的形成过程中起重要的作用。     

Use of a large-scale Triticeae expressed sequence tag resource to reveal gene expression profiles in hexaploid wheat (Triticum aestivum L.).

The US Wheat Genome Project, funded by the National Science Foundation, developed the first large public Triticeae expressed sequence tag (EST) resource. Altogether, 116,272 ESTs were produced, comprising 100,674 5' ESTs and 15 598 3' ESTs. These ESTs were derived from 42 cDNA libraries, which were created from hexaploid bread wheat (Triticum aestivum L.) and its close relatives, including diploid wheat (T. monococcum L. and Aegilops speltoides L.), tetraploid wheat (T. turgidum L.), and rye (Secale cereale L.), using tissues collected from various stages of plant growth and development and under diverse regimes of abiotic and biotic stress treatments. ESTs were assembled into 18,876 contigs and 23,034 singletons, or 41,910 wheat unigenes. Over 90% of the contigs contained fewer than 10 EST members, implying that the ESTs represented a diverse selection of genes and that genes expressed at low and moderate to high levels were well sampled. Statistical methods were used to study the correlation of gene expression patterns, based on the ESTs clustered in the 1536 contigs that contained at least 10 5' EST members and thus representing the most abundant genes expressed in wheat. Analysis further identified genes in wheat that were significantly upregulated (p < 0.05) in tissues under various abiotic stresses when compared with control tissues. Though the function annotation cannot be assigned for many of these genes, it is likely that they play a role associated with the stress response. This study predicted the possible functionality for 4% of total wheat unigenes, which leaves the remaining 96% with their functional roles and expression patterns largely unknown. Nonetheless, the EST data generated in this project provide a diverse and rich source for gene discovery in wheat.     

cDNA(EST)文库的高通量生物信息学分析体系的构建与应用

应用生物信息学方法,构建了一套针对cDNA或EST文库的高通量、自动化分析体系,CLASP(cDNA Library Analysis SystemPrimary)。CLASP基于Linux操作系统,主要由Perl程序构成。它以cDNA文库(ESTs)序列为分析对象,具有自动查找序列同源基因并进行染色体定位(包括细胞遗传学定位和SIS定位)、EST自动延伸等功能;并对不同来源序列进行聚类分析。应用该体系对3对肺癌相关抑制性消减杂交(SSH)cDNA文库进行了分析。结果在所有3对文库的2083条EST中有1492条找到了同源基因,其中1365条得到染色体定位。对所余591条未知基因的EST进行了电子延伸,其中有214条EST得到不同程度的延伸。对上述cDNA文库中已知基因的EST以及电子延伸后的EST再分别进行聚类分析,而后综合两个聚类分析的结果,由此可发现不同文库间的共同与差异表达基因,可用于特定性状相关的基因功能预测。