Construction of cDNA libraries for highly efficient DNA sequencing from the 3' end of expressed genes

The majority of expressed sequence tag (EST) sequences available today have been derived from the 5' ends of cDNA clones. Obtaining high-quality DNA sequences from the 3' ends of oligo(dT)-primed cDNA on a large scale has been difficult because of slippage of the DNA polymerase enzyme used in direct PCR and cycle sequencing. With the completion of whole genome sequencing for more and more organisms, mRNA 3'-UTR sequences can be particularly useful for clustering large numbers of ESTs for the effective discrimination of individual genes and gene families. We have identified a flaw in the widely used oligo(dT) primers for cDNA synthesis, and here we describe an improved priming approach to effectively synthesize cDNA devoid of homopolymeric nucleotide stretches from mRNA poly(A) tails to enable highly efficient and reliable DNA sequence determination from 3' mRNA ends. Using this method, we produced a rat lung cDNA library and successfully sequenced the 3' ends of 98% of all attempted clones.     

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.     

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.     

Generation of expressed sequence tags from low-CO2 and high-CO2 adapted cells of Chlamydomonas reinhardtii.

To characterize genes whose expression is induced in carbon-stress conditions, 12,969 and 13,450 5'-end expressed sequence tags (ESTs) were generated from cells grown in low-CO2 and high-CO2 conditions of the unicellular green alga, Chlamydomonas reinhardtii. These ESTs were clustered into 4436 and 3566 non-redundant EST groups, respectively. Comparison of their sequences with those of 3433 non-redundant ESTs previously generated from the cells under the standard growth condition indicated that 2665 and 1879 EST groups occurred only in the low-CO2 and high-CO2 populations, respectively. It was also noted that 96.2% and 96.0% of the cDNA species respectively obtained from the low-CO2 and high-CO2 conditions had no similar EST sequence deposited in the public databases. The EST species identified only in the low-CO2 treated cells included genes previously reported to be expressed specifically in low-CO2 acclimatized cells, suggesting that the ESTs generated in this study will be a useful source for analysis of genes related to carbon-stress acclimatization. The sequence information and search results of each clone will appear at the web site: http://www.kazusa.or.jp/en/plant/chlamy/EST/.     

Generation of 7137 non-redundant expressed sequence tags from a legume, Lotus japonicus.

For comprehensive analysis of genes expressed in a model legume, Lotus japonicus, a total of 22,983 5' end expressed sequence tags (ESTs) were accumulated from normalized and size-selected cDNA libraries constructed from young (2 weeks old) plants. The EST sequences were clustered into 7137 non-redundant groups. Similarity search against public non-redundant protein database indicated that 3302 groups showed similarity to genes of known function, 1143 groups to hypothetical genes, and 2692 were novel sequences. Homologues of 5 nodule-specific genes which have been reported in other legume species were contained in the collected ESTs, suggesting that the EST source generated in this study will become a useful tool for identification of genes related to legume-specific biological processes. The sequence data of individual ESTs are available at the web site: http://www.kazusa.or.jp/en/plant/lotus/EST/.     

基于PC/Linux的核酸序列电子延伸系统的构建及其应用

新基因全长cDNA序列的获得常常是分子生物学工作者面临的难题。人类基因组计划及其相关计划的实施导致了大量表达序列标签(EST)的产生。利用一定的生物信息学算法,这些EST序列往往可用来对新基因片段进行延伸。采用Linux操作系统,利用Blast软件和Phrap软件以及EST数据库在微机上构建了EST序列的电子延伸系统,并对来自于人胎肝的11386条EST序列和511条插入片段全长cDNA序列进行了电子延伸,结果显示8373条EST序列和389条插入片段全长cDNA序列得到了程度不等的延伸,部分结果通过RACE实验得到证实。该套系统可高效地、规模化进行EST序列的延伸,可为通过实验获得新基因全长cDNA序列提供重要线索。 Abstract:Normally it is difficult to obtain full-length cDNA sequence of novel genes.More and more expressed sequence tags(ESTs) have been obtained since the start-up of human genome project.Powerful system is badly needed for data mining on these EST sequences.Based on a personal computer coupled with Linux operating system and EST database,the Blast software and Phrap software were used to construct a platform for in silico elongation of ESTs in our lab.The performance was tested using 11386 EST sequences and 511 partial-length cDNA sequences.Results demonstrated that 8373 EST and 389 cDNA sequence were elongated using this system.Thus the platform seems to be a fast way for full-length cDNA sequence cloning of new genes.     

At least 104 nucleotides are transposed from the 5' terminus of the avian sarcoma virus genome to the 5' termini of smaller viral mRNAs.

Cells producing avian sarcoma virus (ASV) contain at least three virus-specific mRNAs, two of which are encoded within the 3' half of the viral genome. Each of these viral RNAs can hybridize with single-stranded DNA(cDNA5') that is complementary to a sequence of 101 nucleotides found at the 5' terminus of the ASV genome, but not within the 3' half of the genome. We proposed previously (Weiss, Varmus and Bishop, 1977) that this nucleotide sequence may be transposed to the 5' termini of viral mRNAs during the genesis of these RNAs. We now substantiate this proposal by reporting the isolation and chemical characterization of the nucleotide sequences complementary to cDNA5' in the genome and mRNAs of the Prague B strain of ASV. We isolated the three identified classes of ASVmRNA (38, 28 and 21S) by molecular hybridization; each class of RNA contained a "capped" oligonucleotide identical to that found at the 5' terminus of the ASV genome. When hybridized with cDNA5', each class of RNA gave rise to RNAase-resistant duplex hybrids that probably encompassed the full extent of cDNA5'. The molar yields of duplex conformed approximately to the number of virus-specific RNA molecules in the initial samples; hence most if not all of the molecules of virus-specific RNA could give rise to the duplexes. The duplexes prepared from the various RNAs all contained the capped oligonucleotide found at the 5' terminus of the viral genome and had identical "fingerprints" when analyzed by two-dimensional fractionation following hydrolysis with RNAase T1. In contrast, RNA representing the 3' half of the ASV genome did not form hybrids with cDNA5'. We conclude that a sequence of more than 100 nucleotides is transposed from the 5' end of the ASV genome to the 5' termini of smaller viral RNAs during the genesis of these RNAs. Transposition of nucleotide sequences during the production of mRNA has now been described for three families of animal viruses and may be a common feature of mRNA biogenesis in eucaryotic cells. The mechanism of transposition, however, and the function of the transposed sequences are not known.     

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.