Estimating and comparing the rates of gene discovery and expressed sequence tag (EST) frequencies in EST surveys

MOTIVATION: Expressed sequence tag (EST) surveys are an efficient way to characterize large numbers of genes from an organism. The rate of gene discovery in an EST survey depends on the degree of redundancy of the cDNA libraries from which sequences are obtained. However, few statistical methods have been developed to assess and compare redundancies of various libraries from preliminary EST surveys. RESULTS: We consider statistics for the comparison of EST libraries based upon the frequencies with which genes occur in subsamples of reads. These measures are useful in determining which one of several libraries is more likely to yield new genes in future reads and what proportion of additional reads one might want to take from the libraries in order to be likely to obtain new genes. One approach is to compare single sample measures that have been successfully used in species estimation problems, such as coverage of a library, defined as the proportion of the library that is represented in the given sample of reads. Another single library measure is an estimate of the expected number of additional genes that will be found in a new sample of reads. We also propose statistics that jointly use data from all the libraries. Analogous formulas for coverage and the expected numbers of new genes are presented. These measures consider coverage in a single library based upon reads from all libraries and similarly, the expected numbers of new genes that will be discovered by taking reads from all libraries with fixed proportions. Together, the statistics presented provide useful comparative measures for the libraries that can be used to guide sampling from each of the libraries to maximize the rate of gene discovery. Finally, we present tests for whether genes are equally represented or expressed in a set of libraries. Binomial and chi2 tests are presented for gene-by-gene comparisons of expression. Overall tests of the equality of proportional representation are presented and multiple comparisons issues are addressed. These methods can be used to evaluate changes in gene expression reflected in the composition of EST libraries prepared from different tissue types or cells exposed to different environmental conditions. AVAILABILITY: Software will be made available at http://www.mathstat.dal.ca/~tsusko     

蝶兰胚珠的cDNA文库构建及其特异基因表达的研究

以蝶兰(Phalaenopsis“Mt.Kaala”cv SM9108)为材料,分别提取大孢子母细胞时期胚珠和成熟胚珠的PolyA RNA,反转录成cDNA,构建起两个cDNA文库。克隆筛选采用差异杂交法。从上述两个cDNA文库中,各选择一个筛选出的cDNA,对其在植物体不同器官和不同发育时期的胚珠内的表达进行了分析。结果表明该两个cDNA均为胚珠特异,并且分别在胚珠发育的特定时期表达。推测该两个cDNA的表达受胚珠内部的不同因子调控。     

Loss of ovule identity induced by overexpression of the fertilization-related kinase 2 (ScFRK2), a MAPKKK from Solanum chacoense

In order to gain information about protein kinases acting during plant fertilization and embryogenesis, a reverse genetic approach was used to determine the role of protein kinases expressed in reproductive tissues. Two cDNA clones named ScFRK1 and ScFRK2 (Solanum chacoense fertilization-related kinase 1 and 2) were isolated from an expressed sequence tag (EST) library normalized for weakly expressed genes in fertilized ovaries. These showed significant sequence similarities to members of the mitogen-activated protein kinase kinase kinase (MAPKKK) family. RNA gel blot and RNA in situ hybridization analyses confirmed the strong up-regulation of ScFRK2 in ovules after fertilization. In addition, ScFRK2 mRNAs accumulate during early ovule development in the megasporocyte and in the integument of developing ovules. Overexpression of ScFRK2 led to the production of fruits with a severely reduced number of seeds. The seeds that were produced also exhibited developmental retardation. Analysis of ovaries prior to fertilization showed that the seedless phenotype was caused by a homeotic conversion of ovules into carpel-like structures. The present observations are consistent with the role of ScFRK2 in pre- and post-fertilization events. Furthermore, overexpression of ScFRK2 led to changes in the expression of the class D floral homeotic gene ScFBP11, suggesting that the ScFRK2 kinase may interact, directly or indirectly, with the FBP7/11 pathway that directs establishment of ovule identity.     

Progress in maize gene discovery: a project update

The Maize Gene Discovery Project (MGDP) is a 5-year NSF-funded plant genome initiative that began in 1998. The MGDP collaboration involves researchers at six universities from diverse disciplines with the common goal of discovering new maize genes and developing tools for the phenotypic characterization of maize mutants. The project utilizes several approaches: EST sequencing, cDNA microarray production, and the discovery of gene function and genomic sequence through the use of a recombinant Mu1 transposon (RescueMu). Current achievements of the MGDP (NSF 98–72657) include the sequencing of over 120,000 maize ESTs from diverse cDNA libraries, and over 70,000 RescueMu flanking sequences, as well as the cataloguing of mutant seed and cob phenotypes of 23,000 maize ears, 6,200 families of maize seedlings, and 4,000 families of adult maize plants carrying MuDR/Mu and RescueMu insertion alleles. A consolidation of over 24,000 unique sequences from 19 libraries has been made into the first two of the planned set of four "Unigene" microarray slides. In addition, slides for four EST libraries have been produced. These microarray slides, EST clones, library plates of immortalized RescueMu bacterial cultures, and seed are all available online (http://www.zmdb.iastate.edu). The ZmDB website posts periodic assemblies of all maize EST and genomic sequences available from GenBank. ZmDB is also a portal for sequence analysis software designed to aid in gene discovery: MuSeqBox, GeneSeqer, and SplicePredictor . In addition, ZmDB contains links to other plant and genetics websites. Electronic Publication     

Analysis of the transcriptional response to Rice Yellow Mottle Virus infection in<Emphasis Type="Italic"> Oryza sativa indica</Emphasis> and<Emphasis Type="Italic"> japonica</Emphasis> cultivars

Several cDNA libraries were constructed using mRNA isolated from roots, panicles, cell suspensions and leaves of non-stressed Oryza sativa indica (IR64) and japonica (Azucena) plants, from wounded leaves, and from leaves of both cultivars inoculated with Rice Yellow Mottle Virus (RYMV). A total of 5549 cleaned expressed sequence tags (ESTs) were generated from these libraries. They were classified into functional categories on the basis of homology, and analyzed for redundancy within each library. The expression profiles represented by each library revealed great differences between indica and japonica backgrounds. EST frequencies during the early stages of RYMV infection indicated that changes in the expression of genes involved in energy metabolism and photosynthesis are differentially accentuated in susceptible and partially resistant cultivars. Mapping of these ESTs revealed that several co-localize with previously described resistance gene analogs and QTLs (quantitative trait loci).     

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再分别进行聚类分析,而后综合两个聚类分析的结果,由此可发现不同文库间的共同与差异表达基因,可用于特定性状相关的基因功能预测。     

Production of complex nucleic acid libraries using highly parallel in situ oligonucleotide synthesis

Generation of complex libraries of defined nucleic acid sequences can greatly aid the functional analysis of protein and gene function. Previously, such studies relied either on individually synthesized oligonucleotides or on cellular nucleic acids as the starting material. As each method has disadvantages, we have developed a rapid and cost-effective alternative for construction of small-fragment DNA libraries of defined sequences. This approach uses in situ microarray DNA synthesis for generation of complex oligonucleotide populations. These populations can be recovered and either used directly or immortalized by cloning. From a single microarray, a library containing thousands of unique sequences can be generated. As an example of the potential applications of this technology, we have tested the approach for the production of plasmids encoding short hairpin RNAs (shRNAs) targeting numerous human and mouse genes. We achieved high-fidelity clone retrieval with a uniform representation of intended library sequences.