Correction to: The transcriptome landscape of early maize meiosis.

CORRECTION: Following publication of the original article [1], the authors reported that the number of genes overlaying the bar graph in Fig. 3A were incorrectly counted and inserted (i.e. including a title tile, and in inverse order). The corrected numbers are below and match with the listed genes supplied in Additional File: Table S2.

The plant-specific protein FEHLSTART controls male meiotic entry, initializing meiotic synchronization in Arabidopsis.

Meiosis marks the transition from the sporophyte to the gametophyte generation in the life cycle of flowering plants, and creates genetic variations through homologous recombination. In most flowering plants, meiosis is highly synchronized within each anther, which is significant for efficient fertilization. To date, little is known about the molecular mechanisms of entry into meiosis and exit from it, and only a few genes in Arabidopsis have been characterized with a role in regulating meiotic progression. In this study, we report the functional characterization of a plant-specific basic helix-loop-helix (bHLH) protein, FEHLSTART (FST), a defect in which leads to premature meiotic entry and asynchronous meiosis, and results in decreased seed yield. Investigation of the time course of meiosis showed that the onset of leptotene, the first stage of prophase I, frequently occurred earlier in fst-1 than in the wild type. Asynchronous meiosis followed, which could manifest in the disruption of regular spindle structures and symmetric cell divisions in fst-1 mutants during the meiosis I/II transition. In accordance with frequently accelerated meiotic entry, whole-transcriptome analysis of fst-1 anthers undergoing meiosis revealed that 19 circadian rhythm genes were affected and 47 pollen-related genes were prematurely expressed at a higher level. Taken together, we propose that FST is required for normal meiotic entry and the establishment of meiotic synchrony.

The transcriptome landscape of early maize meiosis.

BACKGROUND: A major step in the higher plant life cycle is the decision to leave the mitotic cell cycle and begin the progression through the meiotic cell cycle that leads to the formation of gametes. The molecular mechanisms that regulate this transition and early meiosis remain largely unknown. To gain insight into gene expression features during the initiation of meiotic recombination, we profiled early prophase I meiocytes from maize (Zea mays) using capillary collection to isolate meiocytes, followed by RNA-seq. RESULTS: We detected ~2,000 genes as preferentially expressed during early meiotic prophase, most of them uncharacterized. Functional analysis uncovered the importance of several cellular processes in early meiosis. Processes significantly enriched in isolated meiocytes included proteolysis, protein targeting, chromatin modification and the regulation of redox homeostasis. The most significantly up-regulated processes in meiocytes were processes involved in carbohydrate metabolism. Consistent with this, many mitochondrial genes were up-regulated in meiocytes, including nuclear- and mitochondrial-encoded genes. The data were validated with real-time PCR and in situ hybridization and also used to generate a candidate maize homologue list of known meiotic genes from Arabidopsis. CONCLUSIONS: Taken together, we present a high-resolution analysis of the transcriptome landscape in early meiosis of an important crop plant, providing support for choosing genes for detailed characterization of recombination initiation and regulation of early meiosis. Our data also reveal an important connection between meiotic processes and altered/increased energy production.

Comparative transcriptomics of early meiosis in Arabidopsis and maize.

Though sexually reproductive plants share the same principle and most processes in meiosis, there are distinct features detectable. To address the similarities and differences of early meiosis transcriptomes from the dicot model system Arabidopsis and monocot model system maize, we performed comparative analyses of RNA-seq data of isolated meiocytes, anthers and seedlings from both species separately and via orthologous genes. Overall gene expression showed similarities, such as an increased number of reads mapping to unannotated features, and differences, such as the amount of differentially expressed genes. We detected major similarities and differences in functional annotations of genes up-regulated in meiocytes, which point to conserved features as well as unique features. Transcriptional regulation seems to be quite similar in Arabidopsis and maize, and we could reveal known and novel transcription factors and cis-regulatory elements acting in early meiosis. Taken together, meiosis between Arabidopsis and maize is conserved in many ways, but displays key distinctions that lie in the patterns of gene expression.

Male-specific region of the bovine Y chromosome is gene rich with a high transcriptomic activity in testis development.

The male-specific region of the mammalian Y chromosome (MSY) contains clusters of genes essential for male reproduction. The highly repetitive and degenerative nature of the Y chromosome impedes genomic and transcriptomic characterization. Although the Y chromosome sequence is available for the human, chimpanzee, and macaque, little is known about the annotation and transcriptome of nonprimate MSY. Here, we investigated the transcriptome of the MSY in cattle by direct testis cDNA selection and RNA-seq approaches. The bovine MSY differs radically from the primate Y chromosomes with respect to its structure, gene content, and density. Among the 28 protein-coding genes/families identified on the bovine MSY (12 single- and 16 multicopy genes), 16 are bovid specific. The 1,274 genes identified in this study made the bovine MSY gene density the highest in the genome; in comparison, primate MSYs have only 31-78 genes. Our results, along with the highly transcriptional activities observed from these Y-chromosome genes and 375 additional noncoding RNAs, challenge the widely accepted hypothesis that the MSY is gene poor and transcriptionally inert. The bovine MSY genes are predominantly expressed and are differentially regulated during the testicular development. Synonymous substitution rate analyses of the multicopy MSY genes indicated that two major periods of expansion occurred during the Miocene and Pliocene, contributing to the adaptive radiation of bovids. The massive amplification and vigorous transcription suggest that the MSY serves as a genomic niche regulating male reproduction during bovid expansion.

Analyzing the meiotic transcriptome using isolated meiocytes of Arabidopsis thaliana.

Improved transcriptome sequencing technologies (RNA-seq) have advanced our understanding of the tissue-specific transcriptome landscapes, including those of messenger RNAs, noncoding RNAs and small RNAs. However, transcriptome profiles of plant meiocytes remain challenging due to the lack of efficient methods to enrich meiocytes for the analysis of temporal and spatial gene expression patterns during meiosis. In this chapter, we describe a method to analyze the Arabidopsis meiotic transcriptome using isolated male meiocytes.

Characterization of a set of novel meiotically-active promoters in Arabidopsis.

BACKGROUND: Homologous recombination, together with selection, laid the foundation for traditional plant breeding. The recombination process that takes place during meiotic cell division is crucial for the creation of novel variations of highly desired traits by breeders. Gaining control over this process is important for molecular breeding to achieve more precise, large-scale and quicker plant improvement. As conventional ubiquitous promoters are neither tissue-specific nor efficient in driving gene expression in meiocytes, promoters with high meiotic activities are potential candidates for manipulating the recombination process. So far, only a few meiotically-active promoters have been reported. Recently developed techniques to profile the transcriptome landscape of isolated meiocytes provided the means to discover promoters from genes that are actively expressed in meiosis. RESULTS: In a screen for meiotically-active promoters, we examined ten promoter sequences that are associated with novel meiotic candidate genes. Each promoter was tested by expressing a GFP reporter gene in Arabidopsis. Characterization of regulatory regions revealed that these meiotically-active promoters possessed conserved motifs and motif arrangement. Some of the promoters unite optimal properties which are invaluable for meiosis-directed studies such as delivering specific gene expression in early meiosis I and/or meiosis II. Furthermore, the examination of homologs of the corresponding genes within green plants points to a great potential of applying the information from Arabidopsis to other species, especially crop plants. CONCLUSIONS: We identified ten novel meiotically-active promoters; which, along with their homologs, are prime candidates to specifically drive gene expression during meiosis in plants and can thus provide important tools for meiosis study and crop breeding.

The expansion of the PRAME gene family in Eutheria.

The PRAME gene family belongs to the group of cancer/testis genes whose expression is restricted primarily to the testis and a variety of cancers. The expansion of this gene family as a result of gene duplication has been observed in primates and rodents. We analyzed the PRAME gene family in Eutheria and discovered a novel Y-linked PRAME gene family in bovine, PRAMEY, which underwent amplification after a lineage-specific, autosome-to-Y transposition. Phylogenetic analyses revealed two major evolutionary clades. Clade I containing the amplified PRAMEYs and the unamplified autosomal homologs in cattle and other eutherians is under stronger functional constraints; whereas, Clade II containing the amplified autosomal PRAMEs is under positive selection. Deep-sequencing analysis indicated that eight of the identified 16 PRAMEY loci are active transcriptionally. Compared to the bovine autosomal PRAME that is expressed predominantly in testis, the PRAMEY gene family is expressed exclusively in testis and is up-regulated during testicular maturation. Furthermore, the sense RNA of PRAMEY is expressed specifically whereas the antisense RNA is expressed predominantly in spermatids. This study revealed that the expansion of the PRAME family occurred in both autosomes and sex chromosomes in a lineage-dependent manner. Differential selection forces have shaped the evolution and function of the PRAME family. The positive selection observed on the autosomal PRAMEs (Clade II) may result in their functional diversification in immunity and reproduction. Conversely, selective constraints have operated on the expanded PRAMEYs to preserve their essential function in spermatogenesis.

ZNF280BY and ZNF280AY: autosome derived Y-chromosome gene families in Bovidae.

BACKGROUND: Recent progress in exploring the Y-chromosome gene content in humans, mice and cats have suggested that "autosome-to-Y" transposition of the male fertility genes is a recurrent theme during the mammalian Y-chromosome evolution. These transpositions are lineage-dependent. The purpose of this study is to investigate the lineage-specific Y-chromosome genes in bovid. RESULTS: We took a direct testis cDNA selection strategy and discovered two novel gene families, ZNF280BY and ZNF280AY, on the bovine (Bos taurus) Y-chromosome (BTAY), which originated from the transposition of a gene block on the bovine chromosome 17 (BTA17) and subsequently amplified. Approximately 130 active ZNF280BY loci (and ~240 pseudogenes) and ~130 pseudogenized ZNF280AY copies are present over the majority of the male-specific region (MSY). Phylogenetic analysis indicated that both gene families fit with the "birth-and-death" model of evolution. The active ZNF280BY loci share high sequence similarity and comprise three major genomic structures, resulted from insertions/deletions (indels). Assembly of a 1.2 Mb BTAY sequence in the MSY ampliconic region demonstrated that ZNF280BY and ZNF280AY, together with HSFY and TSPY families, constitute the major elements within the repeat units. The ZNF280BY gene family was found to express in different developmental stages of testis with sense RNA detected in all cell types of the seminiferous tubules while the antisense RNA detected only in the spermatids. Deep sequencing of the selected cDNAs revealed that different loci of ZNF280BY were differentially expressed up to 60-fold. Interestingly, different copies of the ZNF280AY pseudogenes were also found to differentially express up to 10-fold. However, expression level of the ZNF280AY pseudogenes was almost 6-fold lower than that of the ZNF280BY genes. ZNF280BY and ZNF280AY gene families are present in bovid, but absent in other mammalian lineages. CONCLUSIONS: ZNF280BY and ZNF280AY are lineage-specific, multi-copy Y-gene families specific to Bovidae, and are derived from the transposition of an autosomal gene block. The temporal and spatial expression patterns of ZNF280BYs in testis suggest a role in spermatogenesis. This study offers insights into the genomic organization of the bovine MSY and gene regulation in spermatogenesis, and provides a model for studying evolution of multi-copy gene families in mammals.

Meiosis-specific gene discovery in plants: RNA-Seq applied to isolated Arabidopsis male meiocytes.

BACKGROUND: Meiosis is a critical process in the reproduction and life cycle of flowering plants in which homologous chromosomes pair, synapse, recombine and segregate. Understanding meiosis will not only advance our knowledge of the mechanisms of genetic recombination, but also has substantial applications in crop improvement. Despite the tremendous progress in the past decade in other model organisms (e.g., Saccharomyces cerevisiae and Drosophila melanogaster), the global identification of meiotic genes in flowering plants has remained a challenge due to the lack of efficient methods to collect pure meiocytes for analyzing the temporal and spatial gene expression patterns during meiosis, and for the sensitive identification and quantitation of novel genes. RESULTS: A high-throughput approach to identify meiosis-specific genes by combining isolated meiocytes, RNA-Seq, bioinformatic and statistical analysis pipelines was developed. By analyzing the studied genes that have a meiosis function, a pipeline for identifying meiosis-specific genes has been defined. More than 1,000 genes that are specifically or preferentially expressed in meiocytes have been identified as candidate meiosis-specific genes. A group of 55 genes that have mitochondrial genome origins and a significant number of transposable element (TE) genes (1,036) were also found to have up-regulated expression levels in meiocytes. CONCLUSION: These findings advance our understanding of meiotic genes, gene expression and regulation, especially the transcript profiles of MGI genes and TE genes, and provide a framework for functional analysis of genes in meiosis.

Reaching the depth of the Chinese hamster ovary cell transcriptome.

The high-throughput DNA sequencing Illumina Solexa GAII platform was employed to characterize the transcriptome of an antibody-producing Chinese hamster ovary (CHO) cell line. More than 55 million sequencing reads were generated and mapped to an existing set of CHO unigenes derived from expressed sequence tags (ESTs), as well as several public sequence databases. A very significant fraction of sequencing reads has not been previously seen. The frequency with which fragments of a unigene were sequenced was taken as an estimate of the abundance level of the corresponding transcripts. A wide dynamic range of transcript abundance levels was observed, spanning six orders of magnitude. However, the distribution of coverage across transcript lengths was found to vary, from relatively uniform to highly variable. This observation suggests that more challenges are yet to be resolved before direct sequencing can be used as a true quantitative measure of transcript level and for differential gene expression analysis. With the depth that high-throughput sequencing methods can reach, one can expect that the entire transcriptome of this industrially important organism will be decoded in the near future.

Management of High-Throughput DNA Sequencing Projects: Alpheus.

High-throughput DNA sequencing has enabled systems biology to begin to address areas in health, agricultural and basic biological research. Concomitant with the opportunities is an absolute necessity to manage significant volumes of high-dimensional and inter-related data and analysis. Alpheus is an analysis pipeline, database and visualization software for use with massively parallel DNA sequencing technologies that feature multi-gigabase throughput characterized by relatively short reads, such as Illumina-Solexa (sequencing-by-synthesis), Roche-454 (pyrosequencing) and Applied Biosystem's SOLiD (sequencing-by-ligation). Alpheus enables alignment to reference sequence(s), detection of variants and enumeration of sequence abundance, including expression levels in transcriptome sequence. Alpheus is able to detect several types of variants, including non-synonymous and synonymous single nucleotide polymorphisms (SNPs), insertions/deletions (indels), premature stop codons, and splice isoforms. Variant detection is aided by the ability to filter variant calls based on consistency, expected allele frequency, sequence quality, coverage, and variant type in order to minimize false positives while maximizing the identification of true positives. Alpheus also enables comparisons of genes with variants between cases and controls or bulk segregant pools. Sequence-based differential expression comparisons can be developed, with data export to SAS JMP Genomics for statistical analysis.

Legume Resources: MtDB and Medicago.Org.

To identify the genes and gene functions that underlie key aspects of legume biology, researchers have selected the cool season legume Medicago truncatula as a model system for legume research. The mission of the M. truncatula Consortium is to promote unrestricted sharing of data and information that are provided by Medicago research groups worldwide. Through integration of a variety of data and tools, the medicago.org site intends to facilitate progress in the fields of structural, comparative, and functional genomics. To this goal, and as a consortium partner, the Center for Computational Genomics and Bioinformatics (CCGB) at the University of Minnesota has developed MtDB2.0, the M. truncatula database version 2.0. The MtDB2.0 database is the first step toward the global integration of M. truncatula genomic, genetic, and biological information. MtDB2.0 is a relational database that integrates M. truncatula transcriptome data and provides a wide range of user-defined data mining options. The database is interrogated through a series of interfaces, with 58 options grouped into two filters. Sequence identifiers from all public M. truncatula sites [e.g., IDs from GenBank, CCGB, The Institute for Genomic Research (TIGR), National Center for Genome Resources (NCGR), and I'Institut National de la Recherche Agronomique (INRA)] are fully cross-referenced to facilitate comparisons between different sites, and hypertext links to the appropriate database records are provided for all queries' results. MtDB's goal is to provide researchers with the means to quickly and independently identify sequences that match specific research interests based on user-defined criteria. MtDB2.0 offers unrestricted access to advanced and powerful querying tools unmatched by any other public databases. Structurad Query Language (SQL)-encoded queries with a Java-based Web user interface, incorporate different filtering that allow sophisticated data mining of the expressed sequence tag sequencing project results, including the CCGB M. truncatula Unigene set generated with the Phrap assembler. The underlying database and query software have been designed for ease of updates and portability to other model organisms. Public access to the database is at http://www.medicago.org/MtDB.

ForestTreeDB: a database dedicated to the mining of tree transcriptomes.

ForestTreeDB is intended as a resource that centralizes large-scale expressed sequence tag (EST) sequencing results from several tree species (http://foresttree.org/ftdb). It currently encompasses 344,878 quality sequences from 68 libraries, from diverse organs of conifer and hybrid poplar trees. It utilizes the Nimbus data model to provide a hosting system for multiple projects, and uses object-relational mapping APIs in Java and Perl for data accesses within an Oracle database designed to be scalable, maintainable and extendable. Transcriptome builds or unigene sets occupy the focal point of the system. Several of the five current species-specific unigenes were used to design microarrays and SNP resources. The ForestTreeDB web application provides the means for multiple combination database queries. It presents the user with a list of discrete queries to retrieve and download large EST datasets or sequences from precompiled unigene assemblies. Functional annotation assignment is not trivial in conifers which are distantly related to angiosperm model plants. Optimal annotations are achieved through database queries that integrate results from several procedures based open-source tools. ForestTreeDB aims to facilitate sequence mining of coherent annotations in multiple species to support comparative genomic approaches. We plan to continuously enrich ForestTreeDB with other resources through collaborations with other genomic projects.

Legume genome evolution viewed through the Medicago truncatula and Lotus japonicus genomes.

Genome sequencing of the model legumes, Medicago truncatula and Lotus japonicus, provides an opportunity for large-scale sequence-based comparison of two genomes in the same plant family. Here we report synteny comparisons between these species, including details about chromosome relationships, large-scale synteny blocks, microsynteny within blocks, and genome regions lacking clear correspondence. The Lotus and Medicago genomes share a minimum of 10 large-scale synteny blocks, each with substantial collinearity and frequently extending the length of whole chromosome arms. The proportion of genes syntenic and collinear within each synteny block is relatively homogeneous. Medicago-Lotus comparisons also indicate similar and largely homogeneous gene densities, although gene-containing regions in Mt occupy 20-30% more space than Lj counterparts, primarily because of larger numbers of Mt retrotransposons. Because the interpretation of genome comparisons is complicated by large-scale genome duplications, we describe synteny, synonymous substitutions and phylogenetic analyses to identify and date a probable whole-genome duplication event. There is no direct evidence for any recent large-scale genome duplication in either Medicago or Lotus but instead a duplication predating speciation. Phylogenetic comparisons place this duplication within the Rosid I clade, clearly after the split between legumes and Salicaceae (poplar).

Generation, annotation, analysis and database integration of 16,500 white spruce EST clusters.

BACKGROUND: The sequencing and analysis of ESTs is for now the only practical approach for large-scale gene discovery and annotation in conifers because their very large genomes are unlikely to be sequenced in the near future. Our objective was to produce extensive collections of ESTs and cDNA clones to support manufacture of cDNA microarrays and gene discovery in white spruce (Picea glauca [Moench] Voss). RESULTS: We produced 16 cDNA libraries from different tissues and a variety of treatments, and partially sequenced 50,000 cDNA clones. High quality 3' and 5' reads were assembled into 16,578 consensus sequences, 45% of which represented full length inserts. Consensus sequences derived from 5' and 3' reads of the same cDNA clone were linked to define 14,471 transcripts. A large proportion (84%) of the spruce sequences matched a pine sequence, but only 68% of the spruce transcripts had homologs in Arabidopsis or rice. Nearly all the sequences that matched the Populus trichocarpa genome (the only sequenced tree genome) also matched rice or Arabidopsis genomes. We used several sequence similarity search approaches for assignment of putative functions, including blast searches against general and specialized databases (transcription factors, cell wall related proteins), Gene Ontology term assignation and Hidden Markov Model searches against PFAM protein families and domains. In total, 70% of the spruce transcripts displayed matches to proteins of known or unknown function in the Uniref100 database (blastx e-value < 1e-10). We identified multigenic families that appeared larger in spruce than in the Arabidopsis or rice genomes. Detailed analysis of translationally controlled tumour proteins and S-adenosylmethionine synthetase families confirmed a twofold size difference. Sequences and annotations were organized in a dedicated database, SpruceDB. Several search tools were developed to mine the data either based on their occurrence in the cDNA libraries or on functional annotations. CONCLUSION: This report illustrates specific approaches for large-scale gene discovery and annotation in an organism that is very distantly related to any of the fully sequenced genomes. The ArboreaSet sequences and cDNA clones represent a valuable resource for investigations ranging from plant comparative genomics to applied conifer genetics.

Databases and information integration for the Medicago truncatula genome and transcriptome.

An international consortium is sequencing the euchromatic genespace of Medicago truncatula. Extensive bioinformatic and database resources support the marker-anchored bacterial artificial chromosome (BAC) sequencing strategy. Existing physical and genetic maps and deep BAC-end sequencing help to guide the sequencing effort, while EST databases provide essential resources for genome annotation as well as transcriptome characterization and microarray design. Finished BAC sequences are joined into overlapping sequence assemblies and undergo an automated annotation process that integrates ab initio predictions with EST, protein, and other recognizable features. Because of the sequencing project's international and collaborative nature, data production, storage, and visualization tools are broadly distributed. This paper describes databases and Web resources for the project, which provide support for physical and genetic maps, genome sequence assembly, gene prediction, and integration of EST data. A central project Web site at medicago.org/genome provides access to genome viewers and other resources project-wide, including an Ensembl implementation at medicago.org, physical map and marker resources at mtgenome.ucdavis.edu, and genome viewers at the University of Oklahoma (www.genome.ou.edu), the Institute for Genomic Research (www.tigr.org), and Munich Information for Protein Sequences Center (mips.gsf.de).

Microarrays for global expression constructed with a low redundancy set of 27,500 sequenced cDNAs representing an array of developmental stages and physiological conditions of the soybean plant.

BACKGROUND: Microarrays are an important tool with which to examine coordinated gene expression. Soybean (Glycine max) is one of the most economically valuable crop species in the world food supply. In order to accelerate both gene discovery as well as hypothesis-driven research in soybean, global expression resources needed to be developed. The applications of microarray for determining patterns of expression in different tissues or during conditional treatments by dual labeling of the mRNAs are unlimited. In addition, discovery of the molecular basis of traits through examination of naturally occurring variation in hundreds of mutant lines could be enhanced by the construction and use of soybean cDNA microarrays. RESULTS: We report the construction and analysis of a low redundancy 'unigene' set of 27,513 clones that represent a variety of soybean cDNA libraries made from a wide array of source tissue and organ systems, developmental stages, and stress or pathogen-challenged plants. The set was assembled from the 5' sequence data of the cDNA clones using cluster analysis programs. The selected clones were then physically reracked and sequenced at the 3' end. In order to increase gene discovery from immature cotyledon libraries that contain abundant mRNAs representing storage protein gene families, we utilized a high density filter normalization approach to preferentially select more weakly expressed cDNAs. All 27,513 cDNA inserts were amplified by polymerase chain reaction. The amplified products, along with some repetitively spotted control or 'choice' clones, were used to produce three 9,728-element microarrays that have been used to examine tissue specific gene expression and global expression in mutant isolines. CONCLUSIONS: Global expression studies will be greatly aided by the availability of the sequence-validated and low redundancy cDNA sets described in this report. These cDNAs and ESTs represent a wide array of developmental stages and physiological conditions of the soybean plant. We also demonstrate that the quality of the data from the soybean cDNA microarrays is sufficiently reliable to examine isogenic lines that differ with respect to a mutant phenotype and thereby to define a small list of candidate genes potentially encoding or modulated by the mutant phenotype.

wCLUTO: a Web-enabled clustering toolkit.

As structural and functional genomics efforts provide the biological community with ever-broadening sets of interrelated data, the need to explore such complex information for subtle relationships expands. We present wCLUTO, a Web-enabled version of the stand-alone application CLUTO, designed to apply clustering methods to genomic information. Its first application is focused on the clustering transcriptome data from microarrays. Data can be uploaded by the user into the clustering tool, a choice of several clustering methods can be made and configured, and data are presented to the user in a variety of visual formats, including a three-dimensional "mountain" view of the clusters. Parameters can be explored to rapidly examine a variety of clustering results, and the resulting clusters can be downloaded either for manipulation by other programs or to be saved in a format for publication.

Novel gene expression profiles define the metabolic and physiological processes characteristic of wood and its extractive formation in a hardwood tree species, Robinia pseudoacacia.

Wood is of critical importance to humans as a primary feedstock for biofuel, fiber, solid wood products, and various natural compounds including pharmaceuticals. The trunk wood of most tree species has two distinctly different regions: sapwood and heartwood. In addition to the major constituents, wood contains extraneous chemicals that can be removed by extraction with various solvents. The composition and the content of the extractives vary depending on such factors as, species, growth conditions, and time of year when the tree is cut. Despite the great commercial and keen scientific interest, little is known about the tree-specific biology of the formation of heartwood and its extractives. In order to gain insight on the molecular regulations of heartwood and its extractive formation, we carried out global examination of gene expression profiles across the trunk wood of black locust (Robinia pseudoacacia L.) trees. Of the 2,915 expressed sequenced tags (ESTs) that were generated and analyzed in the current study, 55.3% showed no match to known sequences. Cluster analysis of the ESTs identified a total of 2278 unigene sets, which were used to construct cDNA microarrays. Microarray hybridization analyses were then performed to survey the changes in gene expression profiles of trunk wood. The gene expression profiles of wood formation differ according to the region of trunk wood sampled, with highly expressed genes defining the metabolic and physiological processes characteristic of each region. For example, the gene encoding sugar transport had the highest expression in the sapwood, while the structural genes for flavonoid biosynthesis were up-regulated in the sapwood-heartwood transition zone. This analysis also established the expression patterns of 341 previously unknown genes.