The TIGR Rice Genome Annotation Resource: improvements and new features.

In The Institute for Genomic Research Rice Genome Annotation project (http://rice.tigr.org), we have continued to update the rice genome sequence with new data and improve the quality of the annotation. In our current release of annotation (Release 4.0; January 12, 2006), we have identified 42,653 non-transposable element-related genes encoding 49,472 gene models as a result of the detection of alternative splicing. We have refined our identification methods for transposable element-related genes resulting in 13,237 genes that are related to transposable elements. Through incorporation of multiple transcript and proteomic expression data sets, we have been able to annotate 24 799 genes (31,739 gene models), representing approximately 50% of the total gene models, as expressed in the rice genome. All structural and functional annotation is viewable through our Rice Genome Browser which currently supports 59 tracks. Enhanced data access is available through web interfaces, FTP downloads and a Data Extractor tool developed in order to support discrete dataset downloads.

Physiogenomic resources for rat models of heart, lung and blood disorders.

Cardiovascular disorders are influenced by genetic and environmental factors. The TIGR rodent expression web-based resource (TREX) contains over 2,200 microarray hybridizations, involving over 800 animals from 18 different rat strains. These strains comprise genetically diverse parental animals and a panel of chromosomal substitution strains derived by introgressing individual chromosomes from normotensive Brown Norway (BN/NHsdMcwi) rats into the background of Dahl salt sensitive (SS/JrHsdMcwi) rats. The profiles document gene-expression changes in both genders, four tissues (heart, lung, liver, kidney) and two environmental conditions (normoxia, hypoxia). This translates into almost 400 high-quality direct comparisons (not including replicates) and over 100,000 pairwise comparisons. As each individual chromosomal substitution strain represents on average less than a 5% change from the parental genome, consomic strains provide a useful mechanism to dissect complex traits and identify causative genes. We performed a variety of data-mining manipulations on the profiles and used complementary physiological data from the PhysGen resource to demonstrate how TREX can be used by the cardiovascular community for hypothesis generation.

Standardizing global gene expression analysis between laboratories and across platforms.

To facilitate collaborative research efforts between multi-investigator teams using DNA microarrays, we identified sources of error and data variability between laboratories and across microarray platforms, and methods to accommodate this variability. RNA expression data were generated in seven laboratories, which compared two standard RNA samples using 12 microarray platforms. At least two standard microarray types (one spotted, one commercial) were used by all laboratories. Reproducibility for most platforms within any laboratory was typically good, but reproducibility between platforms and across laboratories was generally poor. Reproducibility between laboratories increased markedly when standardized protocols were implemented for RNA labeling, hybridization, microarray processing, data acquisition and data normalization. Reproducibility was highest when analysis was based on biological themes defined by enriched Gene Ontology (GO) categories. These findings indicate that microarray results can be comparable across multiple laboratories, especially when a common platform and set of procedures are used.

Iterative microarray and RNA interference-based interrogation of the SRC-induced invasive phenotype.

Src kinase has long been recognized as a factor in the progression of colorectal cancer and seems to play a specific role in the development of the metastatic phenotype. In spite of numerous studies conducted to elucidate the exact role of Src in cancer progression, downstream targets of Src remain poorly understood. Gene expression profiling has permitted the identification of large sets of genes that may be functionally interrelated but it is often unclear as to which molecular pathways they belong. Here we have developed an iterative approach to experimentally reconstruct a network of gene activity regulated by Src and contributing to the invasive phenotype. Our strategy uses a combination of phenotypic anchoring of gene expression profiles and loss-of-function screening by way of RNA-mediated interference. Using a panel of human colon cancer cell lines exhibiting differential Src-specific activity and invasivity, we identify the first two levels of gene transcription responsible for the invasive phenotype, where first-tier genes are controlled by Src activity and the second-tier genes are under the influence of the first tier. Specifically, perturbation of first-tier gene activity by either pharmacologic inhibition of Src or RNA-mediated interference-directed knockdown leads to a loss of invasivity and decline of second-tier gene activity. The targeting of first-tier genes may be bypassed altogether because knockdown of second-tier genes led to a similar loss of invasive potential. In this manner, numerous members of a "transcriptional cascade" pathway for metastatic activity have been identified and functionally validated.

Autocrine activation of an osteopontin-CD44-Rac pathway enhances invasion and transformation by H-RasV12.

Activated forms of Ras family members are prevalent in many cancers where Ras mutants transduce signals essential for transformation, angiogenesis, invasion and metastasis. As a cancer progression model, we used NIH3T3 cells to explore the mechanism of Ras-induced tumorigenesis. Ras family mutants H-RasV12 and Rit79L strongly induced foci formation, while Rho family mutants RhoA-QL, Rac1-QL and Cdc42-QL were less effective. A comparison of downstream transcriptional targets of Ras and Rho family members using a 26 383 element cDNA microarray revealed that the osteopontin (OPN) gene exhibited the best correlation between magnitude of gene expression change and level of foci formation (r=0.96, P<0.001). In association with H-RasV12- and Rit79L-mediated transformation, foci secreted OPN protein and upregulated the OPN receptor CD44, suggesting the novel initiation of an aberrant OPN-CD44-Rac autocrine pathway. In support of this were the following observations. First, RGD-deficient OPN protein-binding activity was present in H-RasV12-transformed cells but not in control cells, and binding activity was inhibited by the CD44 blocking antibody. Second, foci formation, cell invasion and Rac activity were induced by H-RasV12 and inhibited by the CD44 blocking antibody. Third, foci formation by H-RasV12 was substantially reduced by a short interfering RNA (siRNA) specifically targeting OPN expression for knockdown. Fourth, H-RasV12-mediated transformation was not blocked by the GRGDS peptide, suggesting that OPN effects were not mediated by the integrins. Lastly, OPN knockdown affected the downstream expression of 160 '2nd tier' genes, and at least a subset of these genes appears to be involved in transformation. Indeed, four genes were selected for knockdown, each resulting in a disruption of foci formation and/or invasion. These results underscore the role of aberrant autocrine signaling and transcriptional networking during tumorigenesis.

The effects of temperature reduction on gene expression and oxidative stress in skeletal muscle from adult zebrafish.

Longevity is inversely proportional to ambient temperature in ectothermic organisms such as fish. However, the mechanism by which reducing temperature over a physiological range increases life span is not known and available data are derived primarily from invertebrates. With a rodent-like longevity and abundant biological resources, the zebrafish is an ideal vertebrate ectothermic model in which to investigate this phenomenon. As an initial approach, the effects of a year-long 10 degrees C reduction in water temperature on global gene expression in tail skeletal muscle from adult zebrafish were determined using an oligonucleotide microarray representing 15,512 genes. Expression levels for approximately 600 genes were up-regulated by 1.7-fold or greater by the reduction in temperature, while a similar number of transcripts were down regulated by more than 1.7-fold. Using gene ontology (GO) classifications for molecular function, two functional groups, "oxygen and reactive oxygen species metabolism" and "response to oxidative stress," were found to be overrepresented among up-regulated genes. Transcripts levels for the genes in these two categories were increased by temperature reduction (TR). However, temperature reduction did not suppress lipid peroxidation potential, protein carbonyl content, or 8-oxoguanine level. Additional studies will be required to further delineate the role of altered gene expression and oxidative stress on the longevity-promoting effects of temperature reduction.

Microarray gene expression profiling during the segmentation phase of zebrafish development.

We analyzed 15,512 unique transcripts from wild-type Danio rerio using a long oligonucleotide microarray containing >16,000 65-mers probes. Total RNA was isolated from staged embryos at 2 h intervals over a 24-h period. On average, at any given time point, 27% of the probe set detected corresponding transcripts in embryonic RNA. There were two predominant patterns in the nearly 4000 genes that changed expression in at least one time point during the first 24 hpf. At 12 hpf, we detected 420 up-regulated and 386 down-regulated genes. By 24 hpf, the number of up- and down-regulated genes had increased to 954 and 766, respectively. While the majority of these genes maintained their new level of expression for the duration of the time course, we identified five genes with phasic regulation over the 24-h time course. Two of these genes, germ cell nuclear factor and mesogenin, have been identified as being expressed during gastrulation (5 1/4 to 10 h postfertilization) and subsequently repressed. A cluster containing 36 distinct ribosomal proteins was up-regulated at 12 h, indicating a capability for de novo protein synthesis during and after this stage. Twenty-three muscle-specific genes were up-regulated late during the initial 24 hpf, corresponding to the development and differentiation of the somites.

Identification of H-Ras, RhoA, Rac1 and Cdc42 responsive genes.

The superfamily of small GTP-binding proteins has expanded dramatically in recent years. The Ras family has long been associated with signaling pathways contributing to normal and aberrant cell growth, while Rho-related protein function is to integrate extracellular signals with specific targets regulating cell morphology, cell aggregation, tissue polarity, cell motility and cytokinesis. Recent findings suggest that certain Rho proteins, including RhoA, Rac1 and Cdc42, can also play a role in signal transduction to the nucleus and cell growth control. However, the nature of the genes regulated by Ras and Rho GTPases, as well as their contribution to their numerous biological effects is still largely unknown. To approach these questions, we investigated the global gene expression pattern induced by activated forms of H-Ras, RhoA, Rac1 and Cdc42 using cDNA microarrays comprising 19 117 unique elements. Using this approach, we identified 1184 genes that were up- or downregulated by at least twofold. Hierarchical cluster analysis revealed the existence of patterns of gene regulation both unique and common to H-Ras V12, RhoA QL, Rac1 QL and Cdc42 QL activation. For example, H-Ras V12 upregulated osteopontin and Akt 1, and H-Ras and RhoA stimulated cyclin G1, cyclin-dependent kinase 8, cyclin A2 and HMGI-C, while Rac1 QL and Cdc42 QL upregulated extracellular matrix and cell adhesion proteins such as alpha-actinin 4, procollagen type I and V and neuropilin. Furthermore, H-Ras V12 downregulated by >eightfold 52 genes compared to only three genes by RhoA QL, Rac1 QL and Cdc42 QL. These results provide key information to begin unraveling the complexity of the molecular mechanisms underlying the transforming potential of Ras and Rho proteins, as well as the numerous morphological and cell cycle effects induced by these small GTPases.64

Assessing unmodified 70-mer oligonucleotide probe performance on glass-slide microarrays.

BACKGROUND: Long oligonucleotide microarrays are potentially more cost- and management-efficient than cDNA microarrays, but there is little information on the relative performance of these two probe types. The feasibility of using unmodified oligonucleotides to accurately measure changes in gene expression is also unclear. RESULTS: Unmodified sense and antisense 70-mer oligonucleotides representing 75 known rat genes and 10 Arabidopsis control genes were synthesized, printed and UV cross-linked onto glass slides. Printed alongside were PCR-amplified cDNA clones corresponding to the same genes, enabling us to compare the two probe types simultaneously. Our study was designed to evaluate the mRNA profiles of heart and brain, along with Arabidopsis cRNA spiked into the labeling reaction at different relative copy number. Hybridization signal intensity did not correlate with probe type but depended on the extent of UV irradiation. To determine the effect of oligonucleotide concentration on hybridization signal, 70-mers were serially diluted. No significant change in gene-expression ratio or loss in hybridization signal was detected, even at the lowest concentration tested (6.25 microm). In many instances, signal intensity actually increased with decreasing concentration. The correlation coefficient between oligonucleotide and cDNA probes for identifying differentially expressed genes was 0.80, with an average coefficient of variation of 13.4%. Approximately 8% of the genes showed discordant results with the two probe types, and in each case the cDNA results were more accurate, as determined by real-time PCR. CONCLUSIONS: Microarrays of UV cross-linked unmodified oligonucleotides provided sensitive and specific measurements for most of the genes studied.

Identification of Src transformation fingerprint in human colon cancer.

We used a classical rodent model of transformation to understand the transcriptional processes, and hence the molecular and cellular events a given cell undergoes when progressing from a normal to a transformed phenotype. Src activation is evident in 80% of human colon cancer, yet the myriad of cellular processes effected at the level of gene expression has yet to be fully documented. We identified a Src 'transformation fingerprint' within the gene expression profiles of Src-transformed rat 3Y1 fibroblasts demonstrating a progression in transformation characteristics. To evaluate the role of this gene set in human cancer development and progression, we extracted the orthologous genes present on the Affymetrix Hu95A GeneChip (12k named genes) and compared expression profiles between the Src-induced rodent cell line model of transformation and staged colon tumors where Src is known to be activated. A similar gene expression pattern between the cell line model and staged colon tumors for components of the cell cycle, cytoskeletal associated proteins, transcription factors and lysosomal proteins suggests the need for co-regulation of several cellular processes in the progression of cancer. Genes not previously implicated in tumorigenesis were detected, as well as a set of 14 novel, highly conserved genes with here-to-fore unknown function. These studies define a set of transformation associated genes whose up-regulation has implications for understanding Src mediated transformation and strengthens the role of Src in the development and progression of human colon cancer. Supportive Supplemental Data can be viewed at http://pga.tigr.org/PGApubs.shtml.