Comparison of RNA expression profiles based on maize expressed sequence tag frequency analysis and micro-array hybridization

Understanding Ustilago maydis Infection of Multiple Maize Organs

Ustilago maydis is a smut fungus that infects all aerial maize organs, namely, seedling leaves, tassels, and ears. In all organs, tumors are formed by inducing hypertrophy and hyperplasia in actively dividing cells; however, the vast differences in cell types and developmental stages for different parts of the plant requires that U. maydis have both general and organ-specific strategies for infecting maize. In this review, we summarize how the maize–U. maydis interaction can be studied using mutant U. maydis strains to better understand how individual effectors contribute to this interaction, either through general or specific expression in a cell type, tissue, or organ. We also examine how male sterile maize mutants that do not support tumor formation can be used to explore key features of the maize anthers that are required for successful infection. Finally, we discuss key unanswered questions about the maize–U. maydis interaction and how new technologies can potentially be used to answer them.

Large-scale analysis of differential gene expression in coffee genotypes resistant and susceptible to leaf miner-toward the identification of candidate genes for marker assisted-selection

BACKGROUND

A successful development of herbivorous insects into plant tissues depends on coordination of metabolic processes. Plants have evolved complex mechanisms to recognize such attacks, and to trigger a defense response. To understand the transcriptional basis of this response, we compare gene expression profiles of two coffee genotypes, susceptible and resistant to leaf miner (Leucoptera coffella). A total of 22000 EST sequences from the Coffee Genome Database were selected for a microarray analysis. Fluorescence probes were synthesized using mRNA from the infested and non-infested coffee plants. Array hybridization, scanning and data normalization were performed using Nimble Scan® e ArrayStar® platforms. Genes with foldchange values +/-2 were considered differentially expressed. A validation of 18 differentially expressed genes was performed in infected plants using qRT-PCR approach.

RESULTS

The microarray analysis indicated that resistant plants differ in gene expression profile. We identified relevant transcriptional changes in defense strategies before insect attack. Expression changes (>2.00-fold) were found in resistant plants for 2137 genes (1266 up-regulated and 873 down-regulated). Up-regulated genes include those responsible for defense mechanisms, hypersensitive response and genes involved with cellular function and maintenance. Also, our analyses indicated that differential expression profiles between resistant and susceptible genotypes are observed in the absence of leaf-miner, indicating that defense is already build up in resistant plants, as a priming mechanism. Validation of selected genes pointed to four selected genes as suitable candidates for markers in assisted-selection of novel cultivars.

CONCLUSIONS

Our results show evidences that coffee defense responses against leaf-miner attack are balanced with other cellular functions. Also analyses suggest a major metabolic reconfiguration that highlights the complexity of this response.

Developmental genetics of the perianthless flowers and bracts of a paleoherb species, Saururus chinensis

Saururus chinensis is a core member of Saururaceae, a perianthless (lacking petals or sepals) family. Due to its basal phylogenetic position and unusual floral composition, study of this plant family is important for understanding the origin and evolution of perianthless flowers and petaloid bracts among angiosperm species. To isolate genes involved in S. chinensis flower development, subtracted floral cDNA libraries were constructed by using suppression subtractive hybridization (SSH) on transcripts isolated from developing inflorescences and seedling leaves. The subtracted cDNA libraries contained a total of 1,141 ESTs and were used to create cDNA microarrays to analyze transcript profiles of developing inflorescence tissues. Subsequently, qRT-PCR analyses of eight MADS-box transcription factors and in situ hybridizations of two B-class MADS-box transcription factors were performed to verify and extend the cDNA microarray results. Finally, putative phylogenetic relationships within the B-class MADS-box gene family were determined using the discovered S. chinensis B-class genes to compare K-domain sequences with B genes from other basal angiosperms. Two hundred seventy-seven of the 1,141 genes were found to be expressed differentially between S. chinensis inflorescence tissues and seedling leaves, 176 of which were grouped into at least one functional category, including transcription (14.75%), energy (12.59%), metabolism (9.12%), protein-related function (8.99%), and cellular transport (5.76%). qRT-PCR and in situ hybridization of selected MADS-box genes supported our microarray data. Phylogenetic analysis indicated that a total of six B-class MADS-box genes were isolated from S. chinensis. The differential regulation of S. chinensis B-class MADS-box transcription factors likely plays a role during the development of subtending bracts and perianthless flowers. This study contributes to our understanding of inflorescence development in Saururus, and represents an initial step toward understanding the formation of petaloid bracts in this species.

Expressed sequence tags in cultivated peanut (Arachis hypogaea): discovery of genes in seed development and response to Ralstonia solanacearum challenge

Although an important oil crop, peanut has only 162,030 expressed sequence tags (ESTs) publicly available, 86,943 of which are from cultivated plants. More ESTs from cultivated peanuts are needed for isolation of stress-resistant, tissue-specific and developmentally important genes. Here, we generated 63,234 ESTs from our 5 constructed peanut cDNA libraries of Ralstonia solanacearum challenged roots, R. solanacearum challenged leaves, and unchallenged cultured peanut roots, leaves and developing seeds. Among these ESTs, there were 14,547 unique sequences with 7,961 tentative consensus sequences and 6,586 singletons. Putative functions for 47.8 % of the sequences were identified, including transcription factors, tissue-specific genes, genes involved in fatty acid biosynthesis and oil formation regulation, and resistance gene analogue genes. Additionally, differentially expressed genes, including those involved in ethylene and jasmonic acid signal transduction pathways, from both peanut leaves and roots, were identified in R. solanacearum challenged samples. This large expression dataset from different peanut tissues will be a valuable source for marker development and gene expression analysis. It will also be helpful for finding candidate genes for fatty acid synthesis and oil formation regulation as well as for studying mechanisms of interactions between the peanut host and R. solanacearum pathogen.

Large-scale sequencing of normalized full-length cDNA library of soybean seed at different developmental stages and analysis of the gene expression profiles based on ESTs

Although GenBank has now covered over 1,400,000 expressed sequence tags (ESTs) from soybean, most ESTs available to the public have been derived from tissues or environmental conditions rather than developing seeds. It is absolutely necessary for annotating the molecular mechanisms of soybean seed development to analyze completely the gene expression profiles of its immature seed at various stages. Here we have constructed a full-length-enriched cDNA library comprised of a total of 45,408 cDNA clones which cover various stages of soybean seed development. Furthermore, we have sequenced from 5' ends of these clones, 36,656 ESTs were obtained in the present study. These EST sequences could be categorized into 27,982 unigenes, including 22,867 contigs and 5,115 singletons, among which 27,931 could be mapped onto soybean 20 chromosome sequences. Comparative genomic analysis with other plants has revealed that these unigenes include lots of candidate genes specific to dicot, legume and soybean. Approximately 1,789 of these unigenes currently show no homology to known soybean sequences, suggesting that many represent mRNAs specifically expressed in seeds. Novel abundant genes involved in the oil synthesis have been found in this study, may serve as a valuable resource for soybean seed improvement.

Meta Analysis of Gene Expression Data within and Across Species

Since the second half of the 1990s, a large number of genome-wide analyses have been described that study gene expression at the transcript level. To this end, two major strategies have been adopted, a first one relying on hybridization techniques such as microarrays, and a second one based on sequencing techniques such as serial analysis of gene expression (SAGE), cDNA-AFLP, and analysis based on expressed sequence tags (ESTs). Despite both types of profiling experiments becoming routine techniques in many research groups, their application remains costly and laborious. As a result, the number of conditions profiled in individual studies is still relatively small and usually varies from only two to few hundreds of samples for the largest experiments. More and more, scientific journals require the deposit of these high throughput experiments in public databases upon publication. Mining the information present in these databases offers molecular biologists the possibility to view their own small-scale analysis in the light of what is already available. However, so far, the richness of the public information remains largely unexploited. Several obstacles such as the correct association between ESTs and microarray probes with the corresponding gene transcript, the incompleteness and inconsistency in the annotation of experimental conditions, and the lack of standardized experimental protocols to generate gene expression data, all impede the successful mining of these data. Here, we review the potential and difficulties of combining publicly available expression data from respectively EST analyses and microarray experiments. With examples from literature, we show how meta-analysis of expression profiling experiments can be used to study expression behavior in a single organism or between organisms, across a wide range of experimental conditions. We also provide an overview of the methods and tools that can aid molecular biologists in exploiting these public data.

Monitoring the expression of maize genes in developing kernels under drought stress using oligo-microarray

Preharvest aflatoxin contamination of grain grown on the US southeastern Coast Plain is provoked and aggravated by abiotic stress. The primary abiotic stress is drought along with high temperatures. The objectives of the present study were to monitor gene expression in developing kernels in response to drought stress and to identify drought-responsive genes for possible use in germplasm assessment. The maize breeding line Tex6 was used, and gene expression profiles were analyzed in developing kernels under drought stress verses well-watered conditions at the stages of 25, 30, 35, 40, 45 d after pollination (DAP) using the 70 mer maize oligo-arrays. A total of 9 573 positive array spots were detected with unique gene IDs, and 7 988 were common in both stressed and well-watered samples. Expression patterns of some genes in several stress response-associated pathways, including abscisic acid, jasmonic acid and phenylalanine ammonia-lyase, were examined, and these specific genes were responsive to drought stress positively. Real-time quantitative polymerase chain reaction validated microarray expression data. The comparison between Tex6 and B73 revealed that there were significant differences in specific gene expression, patterns and levels. Several defense-related genes had been downregulated, even though some defense-related or drought responsive genes were upregulated at the later stages.

A single molecule scaffold for the maize genome

About 85% of the maize genome consists of highly repetitive sequences that are interspersed by low-copy, gene-coding sequences. The maize community has dealt with this genomic complexity by the construction of an integrated genetic and physical map (iMap), but this resource alone was not sufficient for ensuring the quality of the current sequence build. For this purpose, we constructed a genome-wide, high-resolution optical map of the maize inbred line B73 genome containing >91,000 restriction sites (averaging 1 site/∼23 kb) accrued from mapping genomic DNA molecules. Our optical map comprises 66 contigs, averaging 31.88 Mb in size and spanning 91.5% (2,103.93 Mb/∼2,300 Mb) of the maize genome. A new algorithm was created that considered both optical map and unfinished BAC sequence data for placing 60/66 (2,032.42 Mb) optical map contigs onto the maize iMap. The alignment of optical maps against numerous data sources yielded comprehensive results that proved revealing and productive. For example, gaps were uncovered and characterized within the iMap, the FPC (fingerprinted contigs) map, and the chromosome-wide pseudomolecules. Such alignments also suggested amended placements of FPC contigs on the maize genetic map and proactively guided the assembly of chromosome-wide pseudomolecules, especially within complex genomic regions. Lastly, we think that the full integration of B73 optical maps with the maize iMap would greatly facilitate maize sequence finishing efforts that would make it a valuable reference for comparative studies among cereals, or other maize inbred lines and cultivars.

The maize genome contains abundant repeats interspersed by low-copy, gene-coding sequences that make it a challenge to sequence; consequently, current BAC sequence assemblies average 11 contigs per clone. The iMap deals with such complexity by the judicious integration of IBM genetic and B73 physical maps, but the B73 genome structure could differ from the IBM population because of genetic recombination and subsequent rearrangements. Accordingly, we report a genome-wide, high-resolution optical map of maize B73 genome that was constructed from the direct analysis of genomic DNA molecules without using genetic markers. The integration of optical and iMap resources with comparisons to FPC maps enabled a uniquely comprehensive and scalable assessment of a given BAC's sequence assembly, its placement within a FPC contig, and the location of this FPC contig within a chromosome-wide pseudomolecule. As such, the overall utility of the maize optical map for the validation of sequence assemblies has been significant and demonstrates the inherent advantages of single molecule platforms. Construction of the maize optical map represents the first physical map of a eukaryotic genome larger than 400 Mb that was created de novo from individual genomic DNA molecules.

Isolation of expressed sequences from a specific chromosome of Thinopyrum intermedium infected by BYDV

To map important ESTs to specific chromosomes and (or) chromosomal regions is difficult in hexaploid wheat because of its large genome size and serious interference of homoeologous sequences. Large-scale EST sequencing and subsequent chromosome localization are both laborious and time-consuming. The wheat alien addition line TAi-27 contains a pair of chromosomes of Thinopyrum intermedium (Host) Barkworth & D.R. Dewey that carry the resistance gene against barley yellow dwarf virus. In this research, we developed a modified technique based on chromosome microdissection and hybridization-specific amplification to isolate expressed sequences from the alien chromosome of TAi-27 by hybridization between the DNA of the microdissected alien chromosome and cDNA of Th. intermedium infected by barley yellow dwarf virus. Twelve clones were selected, sequenced, and analyzed. Three of them were unknown genes without any hit in the GenBank database and the other nine were highly homologous with ESTs of wheat, barley, and (or) other plants in Gramineae induced by abiotic or biotic stress. The method used in this research to isolate expressed sequences from a specific chromosome has the following advantages: (i) the obtained expressed sequences are larger in size and have 3' end information and (ii) the operation is less complicated. It would be an efficient improved method for genomics and functional genomics research of polyploid plants, especially for EST development and mapping. The obtained expressed sequence data are also informative in understanding the resistance genes on the alien chromosome of TAi-27.

Tissue- and expression level-specific chromatin looping at maize b1 epialleles

This work examines the involvement of chromatin looping in the transcriptional regulation of two epialleles of the maize (Zea mays) b1 gene, B-I and B'. These two epialleles are tissue-specifically regulated and are involved in paramutation. B-I and B' are expressed at high and low levels, respectively. A hepta-repeat approximately 100 kb upstream of the transcription start site (TSS) is required for both paramutation and high b1 expression. Using chromosome conformation capture, we show that the hepta-repeat physically interacts with the TSS region in a tissue- and expression level-specific manner. Multiple repeats are required to stabilize this interaction. High b1 expression is mediated by a multiloop structure; besides the hepta-repeat, other sequence regions physically interact with the TSS as well, and these interactions are epiallele- and expression level-specific. Formaldehyde-assisted isolation of regulatory elements uncovered multiple interacting regions as potentially regulatory.

A joint transcriptomic, proteomic and metabolic analysis of maize endosperm development and starch filling

The maize endosperm transcriptome was investigated through cDNA libraries developed at three characteristic stages: (i) lag phase [10 days after pollination (DAP)]; (ii) beginning of storage (14 DAP); and (iii) maximum starch accumulation rate (21 DAP). Expressed sequence tags for 711, 757 and 384 relevant clones, respectively, were obtained and checked manually. The proportion of sequences with no clear function decreased from 35% to 20%, and a large increase in storage protein sequences (i.e. 5% to 38%) was observed from stages (i) to (iii). The remaining major categories included metabolism (11%-13%), transcription-RNA processing-protein synthesis (13%-20%), protein destination (5%-9%), cellular communication (3%-9%) and cell rescue-defence (4%). Good agreement was generally found between category rank in the 10-DAP transcriptome and the recently reported 14-DAP proteome, except that kinases and proteins for RNA processing were not detected in the latter. In the metabolism category, the respiratory pathway transcripts represented the largest proportion (25%-37%), and showed a shift in favour of glycolysis at 21 DAP. At this stage, amino acid metabolism increased to 17%, whereas starch metabolism surprisingly decreased to 7%. A second experiment focused on carbohydrate metabolism by comparing gene expression at three levels (transcripts, proteins and enzyme activities) in relation to substrate or product from 10 to 40 DAP. Here, two distinct patterns were observed: invertases and hexoses were predominant at the beginning, whereas enzyme patterns in the starch pathway, at the three levels, anticipated and paralleled starch accumulation, suggesting that, in most cases, transcriptional control is responsible for the regulation of starch biosynthesis.

Construction of a rice glycosyltransferase phylogenomic database and identification of rice-diverged glycosyltransferases

Glycosyltransferases (GTs; EC 2.4.x.y) constitute a large group of enzymes that form glycosidic bonds through transfer of sugars from activated donor molecules to acceptor molecules. GTs are critical to the biosynthesis of plant cell walls, among other diverse functions. Based on the Carbohydrate-Active enZymes (CAZy) database and sequence similarity searches, we have identified 609 potential GT genes (loci) corresponding to 769 transcripts (gene models) in rice (Oryza sativa), the reference monocotyledonous species. Using domain composition and sequence similarity, these rice GTs were classified into 40 CAZy families plus an additional unknown class. We found that two Pfam domains of unknown function, PF04577 and PF04646, are associated with GT families GT61 and GT31, respectively. To facilitate functional analysis of this important and large gene family, we created a phylogenomic Rice GT Database (http://ricephylogenomics.ucdavis.edu/cellwalls/gt/). Through the database, several classes of functional genomic data, including mutant lines and gene expression data, can be displayed for each rice GT in the context of a phylogenetic tree, allowing for comparative analysis both within and between GT families. Comprehensive digital expression analysis of public gene expression data revealed that most ( approximately 80%) rice GTs are expressed. Based on analysis with Inparanoid, we identified 282 'rice-diverged' GTs that lack orthologs in sequenced dicots (Arabidopsis thaliana, Populus tricocarpa, Medicago truncatula, and Ricinus communis). Combining these analyses, we identified 33 rice-diverged GT genes (45 gene models) that are highly expressed in above-ground, vegetative tissues. From the literature and this analysis, 21 of these loci are excellent targets for functional examination toward understanding and manipulating grass cell wall qualities. Study of the remainder may reveal aspects of hormone and protein metabolism that are critical for rice biology. This list of 33 genes and the Rice GT Database will facilitate the study of GTs and cell wall synthesis in rice and other plants.

Gene expression profiling during asexual development of the late blight pathogen Phytophthora infestans reveals a highly dynamic transcriptome

Much of the pathogenic success of Phytophthora infestans, the potato and tomato late blight agent, relies on its ability to generate from mycelia large amounts of sporangia, which release zoospores that encyst and form infection structures. To better understand these stages, Affymetrix GeneChips based on 15,650 unigenes were designed and used to profile the life cycle. Approximately half of P. infestans genes were found to exhibit significant differential expression between developmental transitions, with approximately (1)/(10) being stage-specific and most changes occurring during zoosporogenesis. Quantitative reverse-transcription polymerase chain reaction assays confirmed the robustness of the array results and showed that similar patterns of differential expression were obtained regardless of whether hyphae were from laboratory media or infected tomato. Differentially expressed genes encode potential cellular regulators, especially protein kinases; metabolic enzymes such as those involved in glycolysis, gluconeogenesis, or the biosynthesis of amino acids or lipids; regulators of DNA synthesis; structural proteins, including predicted flagellar proteins; and pathogenicity factors, including cell-wall-degrading enzymes, RXLR effector proteins, and enzymes protecting against plant defense responses. Curiously, some stage-specific transcripts do not appear to encode functional proteins. These findings reveal many new aspects of oomycete biology, as well as potential targets for crop protection chemicals.

Rapid EST isolation from chromosome 1R of rye

BACKGROUND

To obtain important expressed sequence tags (ESTs) located on specific chromosomes is currently difficult. Construction of single-chromosome EST library could be an efficient strategy to isolate important ESTs located on specific chromosomes. In this research we developed a method to rapidly isolate ESTs from chromosome 1R of rye by combining the techniques of chromosome microdissection with hybrid specific amplification (HSA).

RESULTS

Chromosome 1R was isolated by a glass needle and digested with proteinase K (PK). The DNA of chromosome 1R was amplified by two rounds of PCR using a degenerated oligonucleotide 6-MW sequence with a Sau3AI digestion site as the primer. The PCR product was digested with Sau3AI and linked with adaptor HSA1, then hybridized with the Sau3AI digested cDNA with adaptor HSA2 of rye leaves with and without salicylic acid (SA) treatment, respectively. The hybridized DNA fragments were recovered by the HSA method and cloned into pMD18-T vector. The cloned inserts were released by PCR using the partial sequences in HSA1 and HSA2 as the primers and then sequenced. Of the 94 ESTs obtained and analyzed, 6 were known sequences located on rye chromosome 1R or on homologous group 1 chromosomes of wheat; all of them were highly homologous with ESTs of wheat, barley and/or other plants in Gramineae, some of which were induced by abiotic or biotic stresses. Isolated in this research were 22 ESTs with unknown functions, probably representing some new genes on rye chromosome 1R.

CONCLUSION

We developed a new method to rapidly clone chromosome-specific ESTs from chromosome 1R of rye. The information reported here should be useful for cloning and investigating the new genes found on chromosome 1R.

Isolation and characterization of expressed sequence tags (ESTs) from subtracted cDNA libraries of Pennisetum glaucum seedlings

Pearl millet (Pennisetum glaucum), used as forage and grain crop is a stress tolerant species. Here we identify differentially regulated transcripts in response to abiotic (salinity, drought and cold) stresses from subtracted cDNA libraries by single-pass sequencing of cDNA clones. A total of 2,494 EST sequences were clustered and assembled into a collection of 1,850 unique sequences with 224 contigs and 1,626 singleton sequences. By sequence comparisons the putative functions of many ESTs could be assigned. Genes with stress related functions include those involved in cellular defense against abiotic stresses and transcripts for proteins involved in stress response signaling and transcription in addition to ESTs encoding unknown functions. These provide new candidate genes for investigation to elucidate their role in abiotic stress. The relative mRNA abundance of 38 selected genes, quantified using real time quantitative RT-PCR, demonstrated the existence of a complex gene regulatory network that differentially modulates gene expression in a kinetics-specific manner in response to different abiotic stresses. Notably, housekeeping and non-target genes were effectively reduced in these subtracted cDNA libraries constructed. These EST sequences are a rich source of stress-related genes and reveal a major part of the stress-response transcriptome that will provide the foundation for further studies into understanding Pennisetum's adaptability to harsh environmental conditions.

Curated genome annotation of Oryza sativa ssp. japonica and comparative genome analysis with Arabidopsis thaliana

We present here the annotation of the complete genome of rice Oryza sativa L. ssp. japonica cultivar Nipponbare. All functional annotations for proteins and non-protein-coding RNA (npRNA) candidates were manually curated. Functions were identified or inferred in 19,969 (70%) of the proteins, and 131 possible npRNAs (including 58 antisense transcripts) were found. Almost 5000 annotated protein-coding genes were found to be disrupted in insertional mutant lines, which will accelerate future experimental validation of the annotations. The rice loci were determined by using cDNA sequences obtained from rice and other representative cereals. Our conservative estimate based on these loci and an extrapolation suggested that the gene number of rice is approximately 32,000, which is smaller than previous estimates. We conducted comparative analyses between rice and Arabidopsis thaliana and found that both genomes possessed several lineage-specific genes, which might account for the observed differences between these species, while they had similar sets of predicted functional domains among the protein sequences. A system to control translational efficiency seems to be conserved across large evolutionary distances. Moreover, the evolutionary process of protein-coding genes was examined. Our results suggest that natural selection may have played a role for duplicated genes in both species, so that duplication was suppressed or favored in a manner that depended on the function of a gene.

Identification of stress-responsive genes in an indica rice (Oryza sativa L.) using ESTs generated from drought-stressed seedlings

The impacts of drought on plant growth and development limit cereal crop production worldwide. Rice (Oryza sativa) productivity and production is severely affected due to recurrent droughts in almost all agroecological zones. With the advent of molecular and genomic technologies, emphasis is now placed on understanding the mechanisms of genetic control of the drought-stress response. In order to identify genes associated with water-stress response in rice, ESTs generated from a normalized cDNA library, constructed from drought-stressed leaf tissue of an indica cultivar, Nagina 22 were used. Analysis of 7794 cDNA sequences led to the identification of 5815 rice ESTs. Of these, 334 exhibited no significant sequence homology with any rice ESTs or full-length cDNAs in public databases, indicating that these transcripts are enriched during drought stress. Analysis of these 5815 ESTs led to the identification of 1677 unique sequences. To characterize this drought transcriptome further and to identify candidate genes associated with the drought-stress response, the rice data were compared with those for abiotic stress-induced sequences obtained from expression profiling studies in Arabidopsis, barley, maize, and rice. This comparative analysis identified 589 putative stress-responsive genes (SRGs) that are shared by these diverse plant species. Further, the identified leaf SRGs were compared to expression profiles for a drought-stressed rice panicle library to identify common sequences. Significantly, 125 genes were found to be expressed under drought stress in both tissues. The functional classification of these 125 genes showed that a majority of them are associated with cellular metabolism, signal transduction, and transcriptional regulation.

DNA Microarray Analyses in Higher Plants

DNA microarrays were originally devised and described as a convenient technology for the global analysis of plant gene expression. Over the past decade, their use has expanded enormously to cover all kingdoms of living organisms. At the same time, the scope of applications of microarrays has increased beyond expression analyses, with plant genomics playing a leadership role in the on-going development of this technology. As the field has matured, the rate-limiting step has moved from that of the technical process of data generation to that of data analysis. We currently face major problems in dealing with the accumulating datasets, not simply with respect to how to archive, access, and process the huge amounts of data that have been and are being produced, but also in determining the relative quality of the different datasets. A major recognized concern is the appropriate use of statistical design in microarray experiments, without which the datasets are rendered useless. A vigorous area of current research involves the development of novel statistical tools specifically for microarray experiments. This article describes, in a necessarily selective manner, the types of platforms currently employed in microarray research and provides an overview of recent activities using these platforms in plant biology.

A microarray analysis of gene expression in the free-living stages of the parasitic nematode Strongyloides ratti

BACKGROUND

The nematode Strongyloides ratti has two adult phases in its lifecycle: one obligate, female and parasitic and one facultative, dioecious and free-living. The molecular control of the development of this free-living generation remains to be elucidated.

RESULTS

We have constructed an S. ratti cDNA microarray and used it to interrogate changes in gene expression during the free-living phase of the S. ratti life-cycle. We have found very extensive differences in gene expression between first-stage larvae (L1) passed in faeces and infective L3s preparing to infect hosts. In L1 stages there was comparatively greater expression of genes involved in growth. We have also compared gene expression in L2 stages destined to develop directly into infective L3s with those destined to develop indirectly into free-living adults. This revealed relatively small differences in gene expression. We find little evidence for the conservation of transcription profiles between S. ratti and S. stercoralis or C. elegans.

CONCLUSION

This is the first multi-gene study of gene expression in S. ratti. This has shown that robust data can be generated, with consistent measures of expression within computationally determined clusters and contigs. We find inconsistencies between EST representation data and microarray hybridization data in the identification of genes with stage-specific expression and highly expressed genes. Many of the genes whose expression is significantly different between L1 and iL3s stages are unknown beyond alignments to predicted genes. This highlights the forthcoming challenge in actually determining the role of these genes in the life of S. ratti.

Organization and variability of the maize genome

With a size approximating that of the human genome, the maize genome is about to become the largest plant genome yet sequenced. Contributing to that size are a whole-genome duplication event and a retrotransposition explosion that produced a large amount of repetitive DNA. This DNA is greatly under-represented in cDNA collections, so analysis of the maize transcriptome has been an expedient way of assessing the gene content of maize. Over 2 million maize cDNA sequences are now available, making maize the third most widely studied organism, behind mouse and man. To date, the sequencing of large-sized DNA clones has been largely driven by the genetic interests of different investigators. The recent construction of a physical map that is anchored to the genetic map will aid immensely in the maize genome-sequencing effort. However, studies showing that the repetitive DNA component is highly polymorphic among maize inbred lines point to the need to sample vertically a few specific regions of the genome to evaluate the extent and importance of this variability.

Changes in gene expression in maize kernel in response to water and salt stress

Increasing pressure on limited water resources for agriculture, together with the global temperature increase, highlight the importance of breeding for drought-tolerant cultivars. A better understanding of the molecular nature of drought stress can be expected through the use of genomics approaches. Here, a macroarray of approximately 2500 maize cDNAs was used for determining transcript changes during water- and salt-stress treatments of developing kernels at 15 days after pollination. Normalization of relative transcript abundances was carried out using a human nebulin control sequence. The proportions of transcripts that changed significantly in abundance upon treatment (>2-fold compared to the control) were determined; 1.5% of the sequences examined were up-regulated by high salinity and 1% by water stress. Both stresses induced 0.8% of the sequences. These include genes involved in various stress responses: abiotic, wounding and pathogen attack (abscisic acid response binding factor, glycine and proline-rich proteins, pathogenesis-related proteins, etc.). The proportion of down-regulated genes was higher than that for up-regulated genes for water stress (3.2%) and lower for salt stress (0.7%), although only eight genes, predominantly involved in energy generation, were down-regulated in both stress conditions. Co-expression of genes of unknown function under defined conditions may help in elucidating their roles in coordinating stress responses.

Gene expression profiling in rice young panicle and vegetative organs and identification of panicle-specific genes through known gene functions

In rice, at the stage from pistil and stamen primordia formation to microsporocyte meiosis, the young panicle organs (YPO) make a great contribution to grain productivity. This period corresponds to the onset of meiosis and marks the transition from vegetative to reproductive stages. By comparing gene expression profiling of YPO with that of rice aerial vegetative organs (AVO), it is possible to gain further molecular insight into this period that is developmentally and functionally important. In this report, a total of 92,582 high-quality ESTs from 5'-end sequencing, including 44,247 from YPO and 48,335 from AVO, were obtained and classified. There were 12,884 (29.12%) ESTs from YPO and 16,304 (33.73%) ESTs from AVO matched to known genes, which generated 1,667 and 2,172 known genes, respectively, after integration of these ESTs. From the functions of known homologous genes, we identified some tissue- and developmental-stage-specified genes in YPO. The expression of these genes clearly reflected the unique functional characteristics of YPO. Furthermore, we estimated that there are about 10,000 mRNAs specifically expressed in rice YPO.

Sorghum expressed sequence tags identify signature genes for drought, pathogenesis, and skotomorphogenesis from a milestone set of 16,801 unique transcripts

Improved knowledge of the sorghum transcriptome will enhance basic understanding of how plants respond to stresses and serve as a source of genes of value to agriculture. Toward this goal, Sorghum bicolor L. Moench cDNA libraries were prepared from light- and dark-grown seedlings, drought-stressed plants, Colletotrichum-infected seedlings and plants, ovaries, embryos, and immature panicles. Other libraries were prepared with meristems from Sorghum propinquum (Kunth) Hitchc. that had been photoperiodically induced to flower, and with rhizomes from S. propinquum and johnsongrass (Sorghum halepense L. Pers.). A total of 117,682 expressed sequence tags (ESTs) were obtained representing both 3' and 5' sequences from about half that number of cDNA clones. A total of 16,801 unique transcripts, representing tentative UniScripts (TUs), were identified from 55,783 3' ESTs. Of these TUs, 9,032 are represented by two or more ESTs. Collectively, these libraries were predicted to contain a total of approximately 31,000 TUs. Individual libraries, however, were predicted to contain no more than about 6,000 to 9,000, with the exception of light-grown seedlings, which yielded an estimate of close to 13,000. In addition, each library exhibits about the same level of complexity with respect to both the number of TUs preferentially expressed in that library and the frequency with which two or more ESTs is found in only that library. These results indicate that the sorghum genome is expressed in highly selective fashion in the individual organs and in response to the environmental conditions surveyed here. Close to 2,000 differentially expressed TUs were identified among the cDNA libraries examined, of which 775 were differentially expressed at a confidence level of 98%. From these 775 TUs, signature genes were identified defining drought, Colletotrichum infection, skotomorphogenesis (etiolation), ovary, immature panicle, and embryo.

Comparative analyses of six solanaceous transcriptomes reveal a high degree of sequence conservation and species-specific transcripts

Background

The Solanaceae is a family of closely related species with diverse phenotypes that have been exploited for agronomic purposes. Previous studies involving a small number of genes suggested sequence conservation across the Solanaceae. The availability of large collections of Expressed Sequence Tags (ESTs) for the Solanaceae now provides the opportunity to assess sequence conservation and divergence on a genomic scale.

Results

All available ESTs and Expressed Transcripts (ETs), 449,224 sequences for six Solanaceae species (potato, tomato, pepper, petunia, tobacco and Nicotiana benthamiana), were clustered and assembled into gene indices. Examination of gene ontologies revealed that the transcripts within the gene indices encode a similar suite of biological processes. Although the ESTs and ETs were derived from a variety of tissues, 55–81% of the sequences had significant similarity at the nucleotide level with sequences among the six species. Putative orthologs could be identified for 28–58% of the sequences. This high degree of sequence conservation was supported by expression profiling using heterologous hybridizations to potato cDNA arrays that showed similar expression patterns in mature leaves for all six solanaceous species. 16–19% of the transcripts within the six Solanaceae gene indices did not have matches among Solanaceae, Arabidopsis, rice or 21 other plant gene indices.

Conclusion

Results from this genome scale analysis confirmed a high level of sequence conservation at the nucleotide level of the coding sequence among Solanaceae. Additionally, the results indicated that part of the Solanaceae transcriptome is likely to be unique for each species.

Reduced representation sequencing: a success in maize and a promise for other plant genomes

Plant, and particularly cereal genomes, are challenging to sequence due to their large size and high repetitive DNA content. Gene-enrichment strategies are alternative or complementary approaches to complete genome sequencing that yield, rapidly and inexpensively, useful sequence data from large and complex genomes. The maize genome is large (2.7 Gbp) and contains large amounts of conserved repetitive elements. Furthermore, the high allelic diversity found between maize inbred lines may necessitate sequencing several inbred lines in order to recover the maize "gene pool". Two gene-enrichment approaches, methylation filtration (MF) and high C(o)t (HC) sequencing have been tested in maize and their ability to sample the gene space has been examined. Combined with other genomic sequencing strategies, gene-enriched genomic sequencing is a practical way to examine the maize gene pool, to order and orient the genic sequences on the genome, and to enable investigation of gene content of other complex plant genomes.

Endosperm-preferred expression of maize genes as revealed by transcriptome-wide analysis of expressed sequence tags

The transcriptome-wide endosperm-preferred expression of maize genes was addressed by analyzing a large database of expressed sequence tags (ESTs). We generated 30,531 high quality sequence-reads from the 5'-ends of cDNA libraries from maize endosperm harvested at 10, 15, and 20 days after pollination. A further 196,900 maize sequence-reads retrieved from public databases were added to this endosperm collection to generate MAIZEST, a database with tools for data storage and analysis. MAIZEST contains 227,431 ESTs, one third of which represents developing endosperm and the remaining two-thirds represent transcripts from 49 cDNA libraries constructed from different organs and tissues. Assembling the MAIZEST ESTs generated 29,206 putative transcripts, of which a set of 4032 assembled sequences was composed exclusively of sequences derived from endosperm cDNA libraries. After sequence analysis using overlapping parameters, a sub-set of 2403 assembled sequences was functionally annotated and revealed a wide variety of putative new genes involved in endosperm development and metabolism.

Profile and analysis of gene expression changes during early development in germinating spores of Ceratopteris richardii

Analysis of an expressed sequence tag library with more than 5,000 sequences from spores of the fern Ceratopteris richardii reveals that more than 3,900 of them represent distinct genes, and almost 70% of these have significant similarity to Arabidopsis (Arabidopsis thaliana) genes. Eight genes are common between three very different dormant plant systems, Ceratopteris spores, Arabidopsis seeds, and Arabidopsis pollen. We evaluated the pattern of mRNA abundance over the first 48 h of spore development using a microarray of cDNAs representing 3,207 distinct genes of C. richardii and determined the relative levels of RNA abundance for 3,143 of these genes using a Bayesian method of statistical analysis. More than 900 of them (29%) show a significant change between any of the five time points analyzed, and these have been annotated based on their sequence similarity with the Arabidopsis proteome. Novel data arising from these analyses identify genes likely to be critical for the germination and subsequent early development of diverse cells and tissues emerging from dormancy.

Identification by suppression subtractive hybridization of genes that are differentially expressed between near-isogenic maize lines in association with sugarcane mosaic virus resistance

The molecular mechanisms underlying the development and progression of sugarcane mosaic virus (SCMV) infection in maize are poorly understood. A study of differential expression was conducted to identify genes involved in resistance to SCMV. In this study, we combined suppression subtractive hybridization and macroarray hybridization to identify genes that are differently expressed in the near isogenic lines F7+ (SCMV resistant) and F7 (susceptible). Altogether, 302 differentially expressed genes were identified in four comparisons based on constitutively expressed and inducible genes, and on compatible and incompatible plant-virus interactions. Apart from genes related to metabolism, most of the functionally classified genes identified belonged to three pathogenesis-related categories: cell rescue, defense, cell death and ageing, signal transduction, and transcription. The latter three groups accounted for 56-66% of the genes classified. Some 19% (60 of 302) of the identified genes had previously been assigned to 29 bins distributed over all ten maize chromosomes. Among the mapped genes, 31% (18 of 58) were located within the Scmv2 and Scmv1 regions on chromosomes 3 and 6, respectively, which have been associated with resistance to SCMV. Promising candidate genes for Scmv1 have been identified, such as AA661457 (receptor-like kinase Xa21-binding protein 3). The implications of the genomic distribution of differentially expressed genes identified by this isogenic comparison are discussed in the context of breeding for resistance.

Comparative analysis of expressed sequence tags from different organs of Vitis vinifera L

Expressed sequence tags (ESTs) are providing a valuable approach to sampling organism-expressed genomes, especially when studying large genomes such as those of many plants. We report on the comparison of 8,647 ESTs generated from six different grape (Vitis vinifera L.) organs: berry, root, leaf, bud, shoot and inflorescence. Clustering and assembly of these ESTs resulted in 4,203 unique sequences and revealed that at this level of EST sampling, each organ shares a low percentage of transcripts with the others. To define organ relationships based on EST counts, we calculated a distance matrix of pairwise correlation coefficients between the libraries which indicated bud, inflorescence and shoot as a group distinct from the other organs considered in this study. A putative function was identified for about 85% of the unique sequences. By assigning them to specific functional classes, we were able to highlight strong differences between organs in the metabolism, protein biosynthesis and photosynthesis categories. This grape EST collection has also proven to be a valuable source for the development of 'functional' simple sequence repeats (SSRs) markers: a total of 405 SSRs have been identified. EST sequences and annotation results have been organised in the IASMA-grape database, freely available at the address http://genomics.iasma.it.

Identification of genes with a biocontrol function in Trichoderma harzianum mycelium using the expressed sequence tag approach

Species of Trichoderma are commercially applied as biological control agents against plant fungal pathogens based on different mechanisms, such as the production of antifungal metabolites, competition for space and nutrients and mycoparasitism. But the integrated biocontrol mechanism of Trichoderma harzianum is not well explored at the genetic level. This study aimed to initiate the preliminary development of an expressed sequence tag (EST) database for T. harzianum and thereby gain potentially useful information on Trichoderma gene sequences in order to elucidate the integrated biocontrol mechanism. Partial sequencing of anonymous cDNA clones is a widely used technique for gene identification. A directional cDNA library has been constructed from mycelium of T. harzianum and 3298 clones have been randomly selected, subjected to single-pass sequencing from the 5' end of the vector, and identified by sequence similarity searches against gene sequences in international databases. Of the 3298 mycelium clones, 2174 exhibit similarity to known genes and 451 to known ESTs, while 673 represent novel gene sequences. Analysis of the identified clones indicated sequence similarity to a broad diversity of genes encoding proteins such as enzymes, structural proteins, and regulatory factors. A significant proportion of genes identified in the mycelium were involved in processes related to mycoparasitism and fungicidal metabolites, as would be expected in biocontrol fungus. These results present the successful application of EST analysis in T. harzianum and provide a preliminary indication of gene expression in mycelium.

A chromosome bin map of 16,000 expressed sequence tag loci and distribution of genes among the three genomes of polyploid wheat

Because of the huge size of the common wheat (Triticum aestivum L., 2n = 6x = 42, AABBDD) genome of 17,300 Mb, sequencing and mapping of the expressed portion is a logical first step for gene discovery. Here we report mapping of 7104 expressed sequence tag (EST) unigenes by Southern hybridization into a chromosome bin map using a set of wheat aneuploids and deletion stocks. Each EST detected a mean of 4.8 restriction fragments and 2.8 loci. More loci were mapped in the B genome (5774) than in the A (5173) or D (5146) genomes. The EST density was significantly higher for the D genome than for the A or B. In general, EST density increased relative to the physical distance from the centromere. The majority of EST-dense regions are in the distal parts of chromosomes. Most of the agronomically important genes are located in EST-dense regions. The chromosome bin map of ESTs is a unique resource for SNP analysis, comparative mapping, structural and functional analysis, and polyploid evolution, as well as providing a framework for constructing a sequence-ready, BAC-contig map of the wheat genome.

Comparative analysis of early embryonic sunflower cDNA libraries

To gain information concerning cell functions and activities during sunflower embryogenesis, an expressed sequence tag (EST) approach was used to analyse gene expression in the early stages of sunflower embryos development. Confocal microscopy observations of whole-mounted embryos allowed us to identify precisely the major steps of the zygotic embryonic development. A time-course analysis was then employed to collect the embryonic material. Three cDNA libraries were constructed from microdissected embryos, and three other cDNA libraries were created using a classical day after pollination schedule. A total of 7106 ESTs were produced and assembled. The total number of putative different genes represents about 43.1 (3064 tentative contigs and singlets) of the analysed sequences. The unigenes that showed similarity to proteins with known or predicted functions (50.3) were classified into 15 different functional categories. The functional profiles were found to be quite similar for all studied embryo stages but statistical analysis revealed that successive and coordinate sets of genes are expressed at each embryonic stage. The analysis allowed us to identify abundant and differentially expressed genes at the early stages of embryos development as well as some putatively interesting genes, showing strong similarities with genes playing key roles in plant and animal embryogenesis. The data presented in this study not only provide a first global overview of the genes expression profile during sunflower embryogenesis but also represent an original and valuable tool for developmental genomics studies on exalbuminous dicots.

A unigene catalogue of 5700 expressed genes in cassava

Two economically important characters, starch content and cassava bacterial blight resistance, were targeted to generate a large collection of cassava ESTs. Two libraries were constructed from cassava root tissues of varieties with high and low starch contents. Other libraries were constructed from plant tissues challenged by the pathogen Xanthomonas axonopodis pv.manihotis. We report here the single pass sequencing of 11,954 cDNA clones from the 5' ends, including 111 from the 3' ends. Cluster analysis permitted the identification of a unigene set of 5,700 sequences. Sequence analyses permitted the assignment of a putative functional category for 37% of sequences whereas approximately 16% sequences did not show any significant similarity with other proteins present in the database and therefore can be considered as cassava specific genes. A group of genes belonging to a large multigene family was identified. We characterize a set of genes detected only in infected libraries putatively involved in the defense response to pathogen infection. By comparing two libraries obtained from cultivars contrasting in their starch content a group of genes associated to starch biosynthesis and differentially expressed was identified. This is the first large cassava EST resource developed today and publicly available thus making a significant contribution to genomic knowledge of cassava.

Large-scale analysis of the barley transcriptome based on expressed sequence tags

To provide resources for barley genomics, 110,981 expressed sequence tags (ESTs) were generated from 22 cDNA libraries representing tissues at various developmental stages. This EST collection corresponds to approximately one-third of the 380,000 publicly available barley ESTs. Clustering and assembly resulted in 14,151 tentative consensi (TCs) and 11 073 singletons, altogether representing 25 224 putatively unique sequences. Of these, 17.5% showed no significant similarity to other barley ESTs present in dbEST. More than 41% of all barley genes are supposed to belong to multigene families and approximately 4% of the barley genes undergo alternative splicing. Based on the functional annotation of the set of unique sequences, the functional category 'Energy' was further analysed to reveal tissue- and stage-specific differences in gene expression. Hierarchical clustering of 362 differentially expressed TCs resulted in the identification of seven major clusters. The clusters reflect biochemical pathways predominantly activated in specific tissues and at various developmental stages. During seed germination glycolysis could be identified as the most predominant biochemical pathway. Germination-specific glycolysis is characterized by the coordinated expression of phosphoenolpyruvate carboxylase and phosphoenolpyruvate carboxykinase, whose antagonistic actions possibly regulate the flux of amino acids into protein biosynthesis and gluconeogenesis respectively. The expression of defence-related and antioxidant genes during germination might be controlled by the ethylene-signalling pathway as concluded from the coordinated expression of those genes and the transcription factors (TF) EIN3 and EREBPG. Moreover, because of their predominant expression in germinating seeds, TF of the AP2 and MYB type are presumably major regulators of germination.

A chromosome bin map of 2148 expressed sequence tag loci of wheat homoeologous group 7

The objectives of this study were to develop a high-density chromosome bin map of homoeologous group 7 in hexaploid wheat (Triticum aestivum L.), to identify gene distribution in these chromosomes, and to perform comparative studies of wheat with rice and barley. We mapped 2148 loci from 919 EST clones onto group 7 chromosomes of wheat. In the majority of cases the numbers of loci were significantly lower in the centromeric regions and tended to increase in the distal regions. The level of duplicated loci in this group was 24% with most of these loci being localized toward the distal regions. One hundred nineteen EST probes that hybridized to three fragments and mapped to the three group 7 chromosomes were designated landmark probes and were used to construct a consensus homoeologous group 7 map. An additional 49 probes that mapped to 7AS, 7DS, and the ancestral translocated segment involving 7BS also were designated landmarks. Landmark probe orders and comparative maps of wheat, rice, and barley were produced on the basis of corresponding rice BAC/PAC and genetic markers that mapped on chromosomes 6 and 8 of rice. Identification of landmark ESTs and development of consensus maps may provide a framework of conserved coding regions predating the evolution of wheat genomes.

Chromosome bin map of expressed sequence tags in homoeologous group 1 of hexaploid wheat and homoeology with rice and Arabidopsis

A total of 944 expressed sequence tags (ESTs) generated 2212 EST loci mapped to homoeologous group 1 chromosomes in hexaploid wheat (Triticum aestivum L.). EST deletion maps and the consensus map of group 1 chromosomes were constructed to show EST distribution. EST loci were unevenly distributed among chromosomes 1A, 1B, and 1D with 660, 826, and 726, respectively. The number of EST loci was greater on the long arms than on the short arms for all three chromosomes. The distribution of ESTs along chromosome arms was nonrandom with EST clusters occurring in the distal regions of short arms and middle regions of long arms. Duplications of group 1 ESTs in other homoeologous groups occurred at a rate of 35.5%. Seventy-five percent of wheat chromosome 1 ESTs had significant matches with rice sequences (E < or = e(-10)), where large regions of conservation occurred between wheat consensus chromosome 1 and rice chromosome 5 and between the proximal portion of the long arm of wheat consensus chromosome 1 and rice chromosome 10. Only 9.5% of group 1 ESTs showed significant matches to Arabidopsis genome sequences. The results presented are useful for gene mapping and evolutionary and comparative genomics of grasses.

Utility of different gene enrichment approaches toward identifying and sequencing the maize gene space

Maize (Zea mays) possesses a large, highly repetitive genome, and subsequently a number of reduced-representation sequencing approaches have been used to try and enrich for gene space while eluding difficulties associated with repetitive DNA. This article documents the ability of publicly available maize expressed sequence tag and Genome Survey Sequences (GSSs; many of which were isolated through the use of reduced representation techniques) to recognize and provide coverage of 78 maize full-length cDNAs (FLCs). All 78 FLCs in the dataset were identified by at least three GSSs, indicating that the majority of maize genes have been identified by at least one currently available GSS. Both methyl-filtration and high-Cot enrichment methods provided a 7- to 8-fold increase in gene discovery rates as compared to random sequencing. The available maize GSSs aligned to 75% of the FLC nucleotides used to perform searches, while the expressed sequence tag sequences aligned to 73% of the nucleotides. Our data suggest that at least approximately 95% of maize genes have been tagged by at least one GSS. While the GSSs are very effective for gene identification, relatively few (18%) of the FLCs are completely represented by GSSs. Analysis of the overlap of coverage and bias due to position within a gene suggest that RescueMu, methyl-filtration, and high-Cot methods are at least partially nonredundant.

Cytoskeletal proteins are coordinately increased in maize genotypes with high levels of eEF1A

The opaque2 (o2) mutation increases the Lys content of maize (Zea mays) endosperm by reducing the synthesis of zein storage proteins and increasing the accumulation of other types of cellular proteins. Elongation factor 1A (eEF1A) is one of these proteins, and its concentration is highly correlated with the amount of other Lys-containing proteins in the endosperm. We investigated the basis for this relationship by comparing patterns of protein accumulation and gene expression between a high (Oh51Ao2) and a low (Oh545o2) eEF1A inbred, as well as between high and low eEF1A recombinant inbred lines obtained from their cross. The content of alpha-zein and several cytoskeletal proteins was measured in high and low eEF1A inbred lines, and the levels of these proteins were found to correlate with that of eEF1A. To extend this analysis, we used an endosperm expressed sequence tag microarray to examine steady-state levels of RNA transcripts in developing endosperm of these genotypes. We identified about 120 genes coordinately regulated in association with eEF1A content. These genes encode proteins involved in several biological structures and processes, including the actin cytoskeleton, the endoplasmic reticulum, and the protein synthesis apparatus. Thus, higher levels of eEF1A in o2 mutants may be related to a more extensive cytoskeletal network surrounding the rough endoplasmic reticulum and increased synthesis of cytoskeleton-associated proteins, all of which contribute significantly to the Lys content of the endosperm.

Anchoring 9,371 maize expressed sequence tagged unigenes to the bacterial artificial chromosome contig map by two-dimensional overgo hybridization

Our goal is to construct a robust physical map for maize (Zea mays) comprehensively integrated with the genetic map. We have used a two-dimensional 24 x 24 overgo pooling strategy to anchor maize expressed sequence tagged (EST) unigenes to 165,888 bacterial artificial chromosomes (BACs) on high-density filters. A set of 70,716 public maize ESTs seeded derivation of 10,723 EST unigene assemblies. From these assemblies, 10,642 overgo sequences of 40 bp were applied as hybridization probes. BAC addresses were obtained for 9,371 overgo probes, representing an 88% success rate. More than 96% of the successful overgo probes identified two or more BACs, while 5% identified more than 50 BACs. The majority of BACs identified (79%) were hybridized with one or two overgos. A small number of BACs hybridized with eight or more overgos, suggesting that these BACs must be gene rich. Approximately 5,670 overgos identified BACs assembled within one contig, indicating that these probes are highly locus specific. A total of 1,795 megabases (Mb; 87%) of the total 2,050 Mb in BAC contigs were associated with one or more overgos, which are serving as sequence-tagged sites for single nucleotide polymorphism development. Overgo density ranged from less than one overgo per megabase to greater than 20 overgos per megabase. The majority of contigs (52%) hit by overgos contained three to nine overgos per megabase. Analysis of approximately 1,022 Mb of genetically anchored BAC contigs indicates that 9,003 of the total 13,900 overgo-contig sites are genetically anchored. Our results indicate overgos are a powerful approach for generating gene-specific hybridization probes that are facilitating the assembly of an integrated genetic and physical map for maize.

Characterization of soybean genomic features by analysis of its expressed sequence tags

We analyzed 314,254 soybean expressed sequence tags (ESTs), including 29,540 from our laboratory and 284,714 from GenBank. These ESTs were assembled into 56,147 unigenes. About 76.92% of the unigenes were homologous to genes from Arabidopsis thaliana ( Arabidopsis). The putative products of these unigenes were annotated according to their homology with the categorized proteins of Arabidopsis. Genes corresponding to cell growth and/or maintenance, enzymes and cell communication belonged to the slow-evolving class, whereas genes related to transcription regulation, cell, binding and death appeared to be fast-evolving. Soybean unigenes with no match to genes within the Arabidopsis genome were identified as soybean-specific genes. These genes were mainly involved in nodule development and the synthesis of seed storage proteins. In addition, we also identified 61 genes regulated by salicylic acid, 1,322 transcription factor genes and 326 disease resistance-like genes from soybean unigenes. SSR analysis showed that the soybean genome was more complex than the Arabidopsis and the Medicago truncatula genomes. GC content in soybean unigene sequences is similar to that in Arabidopsis and M. truncatula. Furthermore, the combined analysis of the EST database and the BAC-contig sequences revealed that the total gene number in the soybean genome is about 63,501.

Identification of differentially expressed transcripts from maturing stem of sugarcane by in silico analysis of stem expressed sequence tags and gene expression profiling

Sugarcane accumulates high concentrations of sucrose in the mature stem and a number of physiological processes on-going in maturing stem tissue both directly and indirectly allow this process. To identify transcripts that are associated with stem maturation, we compared patterns of gene expression in maturing and immature stem tissue by expression profiling and bioinformatic analysis of sets of stem ESTs. This study complements a previous study of gene expression associated directly with sugar metabolism in sugarcane. A survey of sequences derived from stem tissue identified an abundance of several classes of sequence that are associated with fibre biosynthesis in the maturing stem. A combination of EST analyses and microarray hybridization revealed that genes encoding homologues of the dirigent protein, a protein that assists in the stereospecificity of lignin assembly, were the most abundant and most strongly differentially expressed transcripts in maturing stem tissue. There was also evidence of coordinated expression of other categories of fibre biosynthesis and putative defence- and stress-related transcripts in the maturing stem. This study has demonstrated the utility of genomic approaches using large-scale EST acquisition and microarray hybridization techniques to highlight the very significant transcriptional investment the maturing stem of sugarcane has placed in fibre biosynthesis and stress tolerance, in addition to its already well-documented role in sugar accumulation.

PlantGDB, plant genome database and analysis tools

PlantGDB (http://www.plantgdb.org/) is a database of molecular sequence data for all plant species with significant sequencing efforts. The database organizes EST sequences into contigs that represent tentative unique genes. Contigs are annotated and, whenever possible, linked to their respective genomic DNA. Genome sequence fragments are assembled similarly. The goal of the PlantGDB web site is to establish the basis for identifying sets of genes common to all plants or specific to particular species by integrating a number of bioinformatics tools that facilitate gene prediction and cross- species comparisons. For species with large-scale genome sequencing efforts, PlantGDB provides genome browsing capabilities that integrate all available EST and cDNA evidence for current gene models (for Arabidopsis thaliana, see the AtGDB site at http://www.plantgdb.org/AtGDB/).

Comparison of gene representation in midguts from two phytophagous insects, Bombyx mori and Ips pini, using expressed sequence tags

Midgut proteins may provide new molecular targets for insect control. This could be particularly important for some pests, such as pine bark beetles, which are difficult to control by conventional methods. Expressed sequence tags (ESTs) provide information about the activity of a particular tissue, and, in the case of pest insects, may quickly identify potential targets. We present here an EST project representing 574 tentative unique genes (TUGs) expressed in the midgut of the male pine engraver beetle, Ips pini. This tissue uses the mevalonate pathway to produce the monoterpenoid pheromone component, ipsdienol, de novo in response to juvenile hormone (JH) III. Comparison of our ESTs with those previously isolated from larval silkmoth (Bombyx mori) midguts revealed interesting similarities and differences in gene representation that correlate with the conserved and divergent functions of these two tissues. For example, seven mevalonate pathway genes were represented in the I. pini ESTs, while none were found from B. mori. This type of comparison may assist the identification of species-specific targets for future pest control strategies.

Conserved subgroups and developmental regulation in the monocot rop gene family

Rop small GTPases are plant-specific signaling proteins with roles in pollen and vegetative cell growth, abscisic acid signal transduction, stress responses, and pathogen resistance. We have characterized the rop family in the monocots maize (Zea mays) and rice (Oryza sativa). The maize genome contains at least nine expressed rops, and the fully sequenced rice genome has seven. Based on phylogenetic analyses of all available Rops, the family can be subdivided into four groups that predate the divergence of monocots and dicots; at least three have been maintained in both lineages. However, the Rop family has evolved differently in the two lineages, with each exhibiting apparent expansion in different groups. These analyses, together with genetic mapping and identification of conserved non-coding sequences, predict orthology for specific rice and maize rops. We also identified consensus protein sequence elements specific to each Rop group. A survey of ROP-mRNA expression in maize, based on multiplex reverse transcriptase-polymerase chain reaction and a massively parallel signature sequencing database, showed significant spatial and temporal overlap of the nine transcripts, with high levels of all nine in tissues in which cells are actively dividing and expanding. However, only a subset of rops was highly expressed in mature leaves and pollen. Intriguingly, the grouping of maize rops based on hierarchical clustering of expression profiles was remarkably similar to that obtained by phylogenetic analysis. We hypothesize that the Rop groups represent classes with distinct functions, which are specified by the unique protein sequence elements in each group and by their distinct expression patterns.

Maize ROP2 GTPase Provides a Competitive Advantage to the Male Gametophyte

Rop GTPases have been implicated in the regulation of plant signal transduction and cell morphogenesis. To explore ROP2 function in maize, we isolated five Mutator transposon insertions (rop2::Mu alleles). Transmission frequency through the male gametophyte, but not the female, was lower than expected in three of the rop2::Mu mutants. These three alleles formed an allelic series on the basis of the relative transmission rate of each when crossed as trans-heterozygotes. A dramatic reduction in the level of ROP2-mRNA in pollen was associated with the three alleles causing a transmission defect, whereas a rop2::Mu allele that did not result in a defect had wild-type transcript levels, thus confirming that mutation of rop2 causes the mutant phenotype. These data strongly support a role for rop2 in male gametophyte function, perhaps surprisingly, given the expression in pollen of the nearly identical duplicate gene rop9. However, the transmission defect was apparent only when a rop2::Mu heterozygote was used as the pollen donor or when a mixture of wild-type and homozygous mutant pollen was used. Thus, mutant pollen is at a competitive disadvantage compared to wild-type pollen, although mutant pollen grains lacked an obvious cellular defect. Our data demonstrate the importance in vivo of a specific Rop, rop2, in the male gametophyte.

A snapshot of the low temperature stress transcriptome of developing rice seedlings (Oryza sativa L.) via ESTs from subtracted cDNA library

Rice (Oryza sativa L.) is sensitive to chilling particularly during early seedling development. Given the biochemical complexity of tolerance mechanisms, genetic potential for this trait depends on highly coordinated expression of many genes. We used a simple cDNA subtraction strategy to develop Expressed Sequence Tags (ESTs) that represent an important subset of cold stress-upregulated genes. The 3,084 subtracted cDNA clones represent a total of 1,967 unigenes from 1,354 singletons and 613 contigs. As expected in the developing seedlings, genes involved in basic cellular processes, i.e., metabolism, growth and development, protein synthesis, folding and destination, cellular transport, cell division and DNA replication were widely represented. Genes with stress-related and regulatory functions comprised 23.17% of the total ESTs. These categories included proteins with known function in cellular defenses against abiotic (drought, cold and salinity) and biotic (pathogen) stresses, and proteins involved in developmental and stress response signalling and transcription. Based on the types of genes represented, tolerance mechanisms rely on precise integration of developmental processes with stress-related responses. A large fraction of the ESTs (38.7%) represents unknown proteins. This EST library is a rich source of cold stress-related genes, and supplements for other publicly available libraries for comprehensive analysis of the stress-response transcriptome.