GEPIS--quantitative gene expression profiling in normal and cancer tissues

MOTIVATION: Expression profiling in diverse tissues is fundamental to understanding gene function as well as therapeutic target identification. The vast collection of expressed sequence tags (ESTs) and the associated tissue source information provides an attractive opportunity for studying gene expression. RESULTS: To facilitate EST-based expression analysis, we developed GEPIS (gene expression profiling in silico), a tool that integrates EST and tissue source information to compute gene expression patterns in a large panel of normal and tumor samples. We found EST-based expression patterns to be consistent with published papers as well as our own experimental results. We also built a GEPIS Regional Atlas that depicts expression characteristics of all genes in a selected genomic region. This program can be adapted for large-scale screening for genes with desirable expression patterns, as illustrated by our large-scale mining for tissue- and tumor-specific genes. AVAILABILITY: The email server version of the GEPIS application is freely available at http://share.gene.com/share/gepis. An interactive version of GEPIS will soon be freely available at http://www.cgl.ucsf.edu/Research/genentech/gepis/. The source code, modules, data and gene lists can be downloaded at http://share.gene.com/share/gepis.     

Census of orthologous genes and self-organizing maps of biologically relevant transcriptional patterns in chickens (Gallus gallus)

The launch of large-scale chicken expressed sequence tags (EST) projects has placed the chicken in the lead for the number of EST sequences in agriculturally important animals. More than 451,000 chicken ESTs derived from over 158 libraries have been deposited in the NCBI dbEST database as of December 2003. But how many genes these ESTs represent and how they are expressed in different chicken tissues/organs remain undetermined. In the present research, we developed a human gene-based strategy for census of chicken orthologous genes and identification of their expression patterns. Among 34,157 human coding genes used in the study, BLAST analysis revealed that 11,066 genes provisionally matched 248,628 chicken ESTs. Based on the average EST abundance of the orthologous genes, the current public repository of chicken ESTs could represent 20,000 provisional genes. Analysis of gene expression in 14 single tissues/organs showed that approximately 15% of genes were expressed exclusively in single tissue/organ whereas the remaining 85% of genes were co-expressed in two or more tissues/organs. A majority (91.15%) of genes expressed in chicken embryos were also expressed at post-hatch stages, indicating that most genes activated in chicken embryos could serve housekeeping functions. Self-organizing maps (SOM) analysis organized 8807 provisional genes in selected chicken tissues into 98 clusters with each cluster being indicative of common regulatory factors and pathways. A total of 969 provisional orthologous genes were identified as preferentially expressed genes (PEGs) in various chicken tissues/organs (LOD>3.0). No doubt, the present study on gene expression patterns will provide insight into dynamics of metabolic pathways and tissue/organ programming and reprogramming in chickens.     

Isolation of developmentally regulated novel genes based on sequence identity and gene expression pattern.

Based on the surmise that a variety of genes might play important roles in embryonic development and tissue differentiation, and that some of them are likely to be expressed in undifferentiated ES cells, we attempted to identify new genes from the ES cell cDNA library. The modified method of expressed sequence tags (ESTs) and the examination of the expression patterns in adult tissues and in vitro differentiated ES cells were utilized in this study. We have isolated and identified several novel cDNA clones with interesting developmental expression pattern. Among the 83 clones randomly chosen, 23 clones (27.7%) have no homology to any sequences in public databases. The rest contain limited or complete sequence homology to the previously reported mammalian genes or ESTs, yet some clones have not been previously identified in the mouse. To examine the expression profile of clones during development and differentiation, sets of slot blots were hybridized with developmental stage specific or tissue specific probes. Out of 40 novel clones tested (21 totally unknown clones and 19 unidentified clones in mouse), most of them were up- or down-regulated as differentiation proceeded, and some clones showed differentiation-stage specific expression profiles. Surprisingly, a majority of genes were also expressed in adult tissues, and some clones even revealed tissue specific expression. These results demonstrate that not only was the strategy we employed in this study quite efficient for screening novel genes, but that the information gained by such studies would also be a useful guide for further analysis of these genes. It also suggests the feasibility of this approach to explore the genomewide network of gene expression during complicated biological processes, such as embryonic development and tissue differentiation.     

ANALYSIS OF EXPRESSED SEQUENCE TAGS FROM THE GREEN ALGA ULVA LINZA (CHLOROPHYTA)1

There is a general lack of genomic information available for chlorophyte seaweed genera such as Ulva, and in particular there is no information concerning the genes that contribute to adhesion and cell wall biosynthesis for this organism. Partial sequencing of cDNA libraries to generate expressed sequence tags (ESTs) is an effective means of gene discovery and characterization of expression patterns. In this study, a cDNA library was created from sporulating tissue of Ulva linza L. Initially, 650 ESTs were randomly selected from a cDNA library and sequenced from their 5′ ends to obtain an indication of the level of redundancy of the library (21%). The library was normalized to enrich for rarer sequences, and a further 1920 ESTs were sequenced. These sequences were subjected to contig assembly that resulted in a unigene set of approximately 1104 ESTs. Forty‐eight percent of these sequences exhibited significant similarity to sequences in the databases. Phylogenetic comparisons are made between selected sequences with similarity in the databases to proteins involved in aspects of extracellular matrix/cell wall assembly and adhesion.     

A 1.5-Mb physical map of the hidrotic ectodermal dysplasia (Clouston syndrome) gene region on human chromosome 13q11

The HED (hidrotic ectodermal dysplasia) or Clouston syndrome gene (named ED2) has been mapped to the pericentromeric region of chromosome 13 (13q11) to a 2.4-cM interval flanked by markers D13S1828 and D13S1830. We have developed a BAC/PAC-based contig map of this region. This contig, comprising 23 clones and spanning 1.5 Mb, was established by mapping of 27 BAC/PAC end-derived STSs, 11 known polymorphic markers, 2 previously mapped genes, and 14 ESTs. The genomic clone overlaps were confirmed by restriction fragment fingerprint analysis. This contig provides the basis for genomic sequencing and gene identification in the ED2 critical region. Of the 14 ESTs mapped to the contig, 6 show homology to human genes and 8 appear to be novel. Expression patterns of the genes/ESTs were tested by Northern blot and RT-PCR. Full characterization of some of these genes, as well as the novel ESTs, will be useful in assessing their involvement in the HED/Clouston syndrome.     

Using ESTs to improve the accuracy of de novo gene prediction

Background  

ESTs are a tremendous resource for determining the exon-intron structures of genes, but even extensive EST sequencing tends to leave many exons and genes untouched. Gene prediction systems based exclusively on EST alignments miss these exons and genes, leading to poor sensitivity. De novo gene prediction systems, which ignore ESTs in favor of genomic sequence, can predict such "untouched" exons, but they are less accurate when predicting exons to which ESTs align. TWINSCAN is the most accurate de novo gene finder available for nematodes and N-SCAN is the most accurate for mammals, as measured by exact CDS gene prediction and exact exon prediction.     

cDNA-AFLP analysis on bolting or flowering of flowering Chinese cabbage and molecular characteristics of BrcuDFR-like/BrcuAXS gene

The molecular basis of flower bud differentiation in flowering Chinese cabbage (Brassica rapa L. ssp. Chinensis var. utilis Tsen et Lee) was studied in this work. Samples were taken from two varieties, the early-blooming "Youqin 49" and the late-blooming "Youqingtiancaixin 80", at five different developmental stages and studied via cDNA-AFLP. Nineteen expression sequence tags (ESTs) associated with bolting or flowering were isolated and cloned. Blast results indicated that 15 ESTs were involved in the synthesis of anthocayanins, photosynthesis, signal transduction, and phytochrome production. Two ESTs had high similarity to hypothetical proteins with unknown function. Two other ESTs shared no similarity to any sequence in the NCBI database and potentially may be newly identified genes. The deduced amino acid sequences of EST amplified by primer A6T4 or A8T4 had high similarity to both dihydroflavonol reductase (DFR) and UDP-D: -apiose/UDP-D: -xylose synthase (AXS), thus was named BrcuDFR-like/BrcuAXS. Using the cDNA sequence, a putative BrcuDFR-like/BrcuAXS gene was cloned and characterized from flowering Chinese cabbage via rapid amplification of cDNA ends (RACE). The full-length cDNA has 1332?bp with an open frame of 919?bp which codes for a polypeptide of 313 amino acids. The corresponding genome sequence is 2,046?bp. Comparison of cDNA and its corresponding genomic sequence indicates that BrcuDFR-like/BrcuAXS contains 9 exons and 8 introns. The temporal expression patterns indicated the gene is more likely to encode the DFR protein, which catalyzes the synthesis of anthocayanins, than UDP-D: -apiose/UDP-D: -xylose synthase (AXS), which catalyzes the conversion of UDP-D: -glucuronate to a mixture of UDP-D: -apiose and UDP-D: -xylose. Further work is needed to determine what role BrcuDFR-like/BrcuAXS plays during floral organ development.     

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

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