Systemic analysis of microarray data sets towards identifying genes regulating root development

Nitrogen and water are essential for plant growth and development. In this study, we designed experiments to produce gene expression data of poplar roots under nitrogen starvation and water deprivation conditions. We found low concentration of nitrogen led first to increased root elongation followed by lateral root proliferation and eventually increased root biomass. To identify genes regulating root growth and development under nitrogen starvation and water deprivation, we designed a series of data analysis procedures, through which, we have successfully identified biologically important genes. Differentially Expressed Genes (DEGs) analysis identified the genes that are differentially expressed under nitrogen starvation or drought. Protein domain enrichment analysis identified enriched themes (in same domains) that are highly interactive during the treatment. Gene Ontology (GO) enrichment analysis allowed us to identify biological process changed during nitrogen starvation. Based on the above analyses, we examined the local Gene Regulatory Network (GRN) and identified a number of transcription factors. After testing, one of them is a high hierarchically ranked transcription factor that affects root growth under nitrogen starvation. It is very tedious and time-consuming to analyze gene expression data. To avoid doing analysis manually, we attempt to automate a computational pipeline that now can be used for identification of DEGs and protein domain analysis in a single run. It is implemented in scripts of Perl and R.

Advanced glycation end products regulate extracellular matrix protein and protease expression by human glomerular mesangial cells

Advanced glycation end products (AGEs) may play a role in the pathogenesis of diabetic nephropathy, by modulating extracellular matrix turnover. AGEs are known to activate specific membrane receptors, including the receptor for AGE (RAGE). In the present study, we analyzed the various receptors for AGEs expressed by human mesangial cells and we studied the effects of glycated albumin and of carboxymethyl lysine on matrix protein and remodelling enzyme synthesis. Membrane RAGE expression was confirmed by FACS analysis. Microarray methods, RT-PCR, and Northern blot analysis were used to detect and confirm specific gene induction. Zymographic analysis and ELISA were used to measure the induction of tPA and PAI-1. We show herein that cultured human mesangial cells express AGE receptor type 1, type 2 and type 3 and RAGE. AGEs (200 microg/ml) induced at least a 2-fold increase in mRNA for 10 genes involved in ECM remodelling, including tPA, PAI-1 and TIMP-3. The increase in tPA synthesis was confirmed by fibrin zymography. The stimulation of PAI-1 synthesis was confirmed by ELISA. AGEs increased PAI-1 mRNA through a signalling pathway involving reactive oxygen species, the MAP kinases ERK-1/ERK-2 and the nuclear transcription factor NF-kappaB, but not AP-1. Carboxymethyl lysine (CML, 5 microM), which is a RAGE ligand, also stimulated PAI-1 synthesis by mesangial cells. In addition, a blocking anti-RAGE antibody partially inhibited the AGE-stimulated gene expression and decreased the PAI-1 accumulation induced by AGEs and by CML. Inhibition of AGE receptors or neutralization of the protease inhibitors TIMP-3 and PAI-1 could represent an important new therapeutic strategy for diabetic nephropathy.

Expression Analysis of the Theileria parva Subtelomere-Encoded Variable Secreted Protein Gene Family

Background The intracellular protozoan parasite Theileria parva transforms bovine lymphocytes inducing uncontrolled proliferation. Proteins released from the parasite are assumed to contribute to phenotypic changes of the host cell and parasite persistence. With 85 members, genes encoding subtelomeric variable secreted proteins (SVSPs) form the largest gene family in T. parva. The majority of SVSPs contain predicted signal peptides, suggesting secretion into the host cell cytoplasm. Methodology/Principal Findings We analysed SVSP expression in T. parva-transformed cell lines established in vitro by infection of T or B lymphocytes with cloned T. parva parasites. Microarray and quantitative real-time PCR analysis revealed mRNA expression for a wide range of SVSP genes. The pattern of mRNA expression was largely defined by the parasite genotype and not by host background or cell type, and found to be relatively stable in vitro over a period of two months. Interestingly, immunofluorescence analysis carried out on cell lines established from a cloned parasite showed that expression of a single SVSP encoded by TP03_0882 is limited to only a small percentage of parasites. Epitope-tagged TP03_0882 expressed in mammalian cells was found to translocate into the nucleus, a process that could be attributed to two different nuclear localisation signals. Conclusions Our analysis reveals a complex pattern of Theileria SVSP mRNA expression, which depends on the parasite genotype. Whereas in cell lines established from a cloned parasite transcripts can be found corresponding to a wide range of SVSP genes, only a minority of parasites appear to express a particular SVSP protein. The fact that a number of SVSPs contain functional nuclear localisation signals suggests that proteins released from the parasite could contribute to phenotypic changes of the host cell. This initial characterisation will facilitate future studies on the regulation of SVSP gene expression and the potential biological role of these enigmatic proteins.

Loss of Raf-1 kinase inhibitor protein (RKIP) is strongly associated with high-grade tumor budding and correlates with an aggressive phenotype in pancreatic ductal adenocarcinoma (PDAC)

BACKGROUND Raf-1 kinase inhibitor protein (RKIP) has emerged as a significant metastatic suppressor in a variety of human cancers and is known to inhibit Ras/Raf/MEK/ERK signaling. By suppressing the activation of the NFkB/SNAIL circuit, RKIP can regulate the induction of epithelial-mesenchymal transition (EMT). The aim of this study was to evaluate RKIP expression and to determine its association with clinicopathological features, including EMT in form of tumor budding in pancreatic ductal adenocarcinoma (PDAC). METHODS Staining for RKIP was performed on a multipunch Tissue Microarray (TMA) of 114 well-characterized PDACs with clinico-pathological, follow-up and adjuvant therapy information. RKIP-expression was assessed separately in the main tumor body and in the tumor buds. Another 3 TMAs containing normal pancreatic tissue, precursor lesions (Pancreatic Intraepithelial Neoplasia, PanINs) and matched lymph node metastases were stained in parallel. Cut-off values were calculated by receiver operating characteristic (ROC) curve analysis. RESULTS We found a significant progressive loss of RKIP expression between normal pancreatic ductal epithelia (average: 74%), precursor lesions (PanINs; average: 37%), PDAC (average 20%) and lymph node metastases (average 8%, p<0.0001). RKIP expression was significantly lower in tumor buds (average: 6%) compared to the main tumor body (average 20%; p<0.005). RKIP loss in the tumor body was marginally associated with advanced T-stage (p=0.0599) as well as high-grade peritumoral (p=0.0048) and intratumoral budding (p=0.0373). RKIP loss in the buds showed a clear association with advanced T stage (p=0.0089). CONCLUSIONS The progressive loss of RKIP seems to play a major role in the neoplastic transformation of pancreas, correlates with aggressive features in PDAC and is associated with the presence of EMT in form of tumor budding.

Loss of tumor suppressor mir-203 mediates overexpression of LIM and SH3 Protein 1 (LASP1) in high-risk prostate cancer thereby increasing cell proliferation and migration

Several studies have linked overexpression of the LIM and SH3 domain protein 1 (LASP1) to progression of breast, colon, liver, and bladder cancer. However, its expression pattern and role in human prostate cancer (PCa) remained largely undefined. Analysis of published microarray data revealed a significant overexpression of LASP1 in PCa metastases compared to parental primary tumors and normal prostate epithelial cells. Subsequent gene-set enrichment analysis comparing LASP1-high and -low PCa identified an association of LASP1 with genes involved in locomotory behavior and chemokine signaling. These bioinformatic predictions were confirmed in vitro as the inducible short hairpin RNA-mediated LASP1 knockdown impaired migration and proliferation in LNCaP prostate cancer cells. By immunohistochemical staining and semi-quantitative image analysis of whole tissue sections we found an enhanced expression of LASP1 in primary PCa and lymph node metastases over benign prostatic hyperplasia. Strong cytosolic and nuclear LASP1 immunoreactivity correlated with PSA progression. Conversely, qRT-PCR analyses for mir-203, which is a known translational suppressor of LASP1 in matched RNA samples revealed an inverse correlation of LASP1 protein and mir-203 expression. Collectively, our results suggest that loss of mir-203 expression and thus uncontrolled LASP1 overexpression might drive progression of PCa.

Common target genes of palatal and gingival fibroblasts for EMD: the microarray approach.

BACKGROUND AND OBJECTIVE Connective tissue grafts are frequently applied, together with Emdogain(®) , for root coverage. However, it is unknown whether fibroblasts from the gingiva and from the palate respond similarly to Emdogain. The aim of this study was therefore to evaluate the effect of Emdogain(®) on fibroblasts from palatal and gingival connective tissue using a genome-wide microarray approach. MATERIAL AND METHODS Human palatal and gingival fibroblasts were exposed to Emdogain(®) and RNA was subjected to microarray analysis followed by gene ontology screening with Database for Annotation, Visualization and Integrated Discovery functional annotation clustering, Kyoto Encyclopedia of Genes and Genomes pathway analysis and the Search Tool for the Retrieval of Interacting Genes/Proteins functional protein association network. Microarray results were confirmed by quantitative RT-PCR analysis. RESULTS The transcription levels of 106 genes were up-/down-regulated by at least five-fold in both gingival and palatal fibroblasts upon exposure to Emdogain(®) . Gene ontology screening assigned the respective genes into 118 biological processes, six cellular components, eight molecular functions and five pathways. Among the striking patterns observed were the changing expression of ligands targeting the transforming growth factor-beta and gp130 receptor family as well as the transition of mesenchymal epithelial cells. Moreover, Emdogain(®) caused changes in expression of receptors for chemokines, lipids and hormones, and for transcription factors such as SMAD3, peroxisome proliferator-activated receptor gamma and those of the ETS family. CONCLUSION The present data suggest that Emdogain(®) causes substantial alterations in gene expression, with similar patterns observed in palatal and gingival fibroblasts.

Salivary Pellicle vs Whole Saliva: The Response of Oral Fibroblasts Based on a Genome-Wide Microarray.

PURPOSE Whole saliva comprises components of the salivary pellicle that spontaneously forms on surfaces of implants and teeth. However, there are no studies that functionally link the salivary pellicle with a possible change in gene expression. MATERIALS AND METHODS This study examined the genetic response of oral fibroblasts exposed to the salivary pellicle and whole saliva. Oral fibroblasts were seeded onto a salivary pellicle and the respective untreated surface. Oral fibroblasts were also exposed to freshly harvested sterile-filtered whole saliva. A genome-wide microarray of oral fibroblasts was performed, followed by gene ontology screening with DAVID functional annotation clustering, KEGG pathway analysis, and the STRING functional protein association network. RESULTS Exposure of oral fibroblasts to saliva caused 61 genes to be differentially expressed (P < .05). Gene ontology screening assigned the respective genes into 262 biologic processes, 3 cellular components, 13 molecular functions, and 7 pathways. Most remarkable was the enrichment in the inflammatory response. None of the genes regulated by whole saliva was significantly changed when cells were placed onto a salivary pellicle. CONCLUSION The salivary pellicle per se does not provoke a significant inflammatory response of oral fibroblasts in vitro, whereas sterile-filtered whole saliva does produce a strong inflammatory response.

Retinoic acid receptor beta and angiopoietin-like protein 1 are involved in the regulation of human androgen biosynthesis

Androgens are essential for sexual development and reproduction. However, androgen regulation in health and disease is poorly understood. We showed that human adrenocortical H295R cells grown under starvation conditions acquire a hyperandrogenic steroid profile with changes in steroid metabolizing enzymes HSD3B2 and CYP17A1 essential for androgen production. Here we studied the regulatory mechanisms underlying androgen production in starved H295R cells. Microarray expression profiling of normal versus starved H295R cells revealed fourteen differentially expressed genes; HSD3B2, HSD3B1, CYP21A2, RARB, ASS1, CFI, ASCL1 and ENC1 play a role in steroid and energy metabolism and ANGPTL1, PLK2, DUSP6, DUSP10 and FREM2 are involved in signal transduction. We discovered two new gene networks around RARB and ANGPTL1, and show how they regulate androgen biosynthesis. Transcription factor RARB stimulated the promoters of genes involved in androgen production (StAR, CYP17A1 and HSD3B2) and enhanced androstenedione production. For HSD3B2 regulation RARB worked in cooperation with Nur77. Secretory protein ANGPTL1 modulated CYP17A1 and DUSP6 expression by inducing ERK1/2 phosphorylation. By contrast, our studies revealed no evidence for hormones or cell cycle involvement in regulating androgen biosynthesis. In summary, these studies establish a firm role for RARB and ANGPTL1 in the regulation of androgen production in H295R cells.

Loss of activator protein-2alpha results in overexpression of the thrombin receptor and correlates with the malignant phenotype of human melanoma

Increasing evidence demonstrates that the thrombin receptor (protease activated receptor-1, PAR-1) plays a major role in tumor invasion and contributes to the metastatic phenotype of human melanoma. We demonstrate that the metastatic potential of human melanoma cells correlates with overexpression of PAR-1. The promoter of the PAR-1 gene contains multiple putative AP-2 and Sp1 consensus elements. We provide evidence that an inverse correlation exists between the expression of AP-2 and the expression of PAR-1 in human melanoma cells. Re-expression of AP-2 in WM266-4 melanoma cells (AP-2 negative) resulted in decreased mRNA and protein expression of PAR-1 and significantly reduced the tumor potential in nude mice. ChIP analysis of the PAR-1 promoter regions bp −365 to −329 (complex 1) and bp −206 to −180 (complex 2) demonstrates that in metastatic cells Sp1 is predominantly binding to the PAR-1 promoter, while in nonmetastatic cells AP-2 is bound. In vitro analysis of complex 1 demonstrates that AP-2 and Sp1 bind to this region in a mutually exclusive manner. Transfection experiments with full-length and progressive deletions of the PAR-1 promoter luciferase constructs demonstrated that metastatic cells had increased promoter activity compared to low and nonmetastatic melanoma cells. Our data shows that exogenous AP-2 expression decreased promoter activity, while transient expression of Sp1 further activated expression of the reporter gene. Mutational analysis of complex 1 within PAR-1 luciferase constructs further demonstrates that the regulation of PAR-1 is mediated through interactions with AP-2 and Sp1. Moreover, loss of AP-2 in metastatic cells alters the AP-2 to Sp1 ratio and DNA-binding activity resulting in overexpression of PAR-1. In addition, we evaluated the expression of AP-2 and PAR-1 utilizing a tissue microarray of 93 melanocytic lesions spanning from benign nevi to melanoma metastasis. We report loss of AP-2 expression in malignant tumors compared to benign tissue while PAR-1 was expressed more often in metastatic melanoma cells than in benign melanocytes. We propose that loss of AP-2 results in increased expression of PAR-1, which in turn results in upregulation of gene products that contribute to the metastatic phenotype of melanoma. ^

A mechanism of insulin-like growth factor binding protein 2-induced cell mobility in glioma

Insulin-like growth factor binding protein 2 (IGFBP2) is a protein known to be overexpressed in a majority of glioblastoma multiforme (GBM) tumors. While it is known the IGFBP2 is involved in promoting GBM tumor cell invasion, no mechanism exists for how the protein is involved in signal transduction pathways leading to enhanced cell invasion. ^ We follow up on preliminary microarray data on IGFBP2-overexpressing GBM cells and protein sequence analysis of IGFBP2 in generating the hypothesis that IGFBP2 interacts with integnn α5 in regulating cell mobility. Microarray data showing upregulation of integrin α5 by IGFBP2 is validated and evidence of protein-protein interaction between IGFBP2 and integrin α5 is found. The exact binding domain on IGFBP2 responsible for its interaction with integrin α5 is also determined, confirming our initial findings and reaffirming that the IGFBP2/integrin α5 interaction is specific. Disruption of this interaction resulted in attenuation of IGFBP2-enhanced cell mobility. Further, we found that cell mobility is only enhanced when IGFBP2 and integrin α5 are both overexpressed and able to interact with each other. ^ We also determined fibronectin to be a critical player in the activation of the IGFBP2/integrin α5 pathway. The activation of this pathway appears to be progressive and initiates once GBM cells have sufficiently established anchorage. ^

Role of activation of the protein tyrosine kinase, Src, in colorectal cancer chemoresistance

Advances in therapy for colorectal cancer have been hampered by development of resistance to chemotherapy. The Src family of protein tyrosine kinases has been associated with colorectal cancer development and progression. Activation of the prototypic member of the family, Src, occurs in advanced colorectal cancer and is associated with a worse outcome. This work tests the hypotheses that Src activation contributes to chemoresistance in some colon tumors and that this resistance can be overcome by use of Src inhibitors. The aims of the proposal were to (1) determine if constitutive Src activation is sufficient to induce oxaliplatin resistance; (2) evaluate the role of reactive oxygen species (ROS) in the activation of Src after oxaliplatin treatment; (3) determine the frequency of Src activation in liver metastases after oxaliplatin treatment; and (4) evaluate the safety, preliminary efficacy, and pharmacodynamics of the combination of dasatinib with oxaliplatin-based therapy in patients with metastatic colorectal cancer. ^ Using a panel of colon cancer cell lines and murine models, I demonstrate that administration of oxaliplatin, a commonly utilized chemotherapy for colorectal cancer, results in an increased activation of Src. The activation occurs acutely in some, but not all, colorectal carcinoma cell lines. Cell lines selected for oxaliplatin resistance are further increased in Src activity. Treatment of cell lines with dasatinib, a non-selective pharmacologic inhibitor of the Src family kinases synergistically killed some, but not all cell lines. Cell lines with the highest acute activation of Src after oxaliplatin administration were the most sensitive to the combination therapy. Previous work demonstrated that siRNA to Src increased sensitivity to oxaliplatin, suggesting that the effects of dasatinib are primarily due to its ability to inhibit Src in these cell lines. ^ To examine the mechanism underlying these results, I examined the effects of reactive oxygen species (ROS), as previous studies have demonstrated that platinum chemotherapeutics result in intracellular oxidative stress. I demonstrated that oxaliplatin-induced reactive oxygen species were higher in the cell lines with Src activation, relative to those in which Src was not activated. This oxaliplatin-induced Src activation was blocked by the administration of anti-oxidants, thereby demonstrating that synergistic killing between dasatinib and oxaliplatin was associated with the ability of the latter to generate ROS. ^ In a murine model of colorectal cancer metastasis to the liver, the combination of dasatinib and oxaliplatin was more effective in reducing tumor volume than either agent alone. However, when oxaliplatin resistant cell lines were treated with a combination of oxaliplatin and AZD0530, an inhibitor in the clinic with increased specificity for Src, no additional benefit was seen, although Src was activated by oxaliplatin and Src substrates were inhibited. The indolent growth of oxaliplatin-resistant cells, unlike the growth of oxaliplatin resistant tumors in patients, precludes definitive interpretation of these results. ^ To further explore Src activation in patients with oxaliplatin exposure and resistance, an immunohistochemistry analysis of tumor tissue from resected liver metastases of colorectal cancer was performed. Utilizing a tissue microarray, staining for phosphorylated Src and FAK demonstrated strong staining of tumor relative to stromal and normal liver. In patients recently exposed to oxaliplatin, there was increased FAK activation, supporting the clinical relevance of the prior preclinical studies. ^ To pursue the potential clinical benefit of the combination of Src inhibition with oxaliplatin, a phase IB clinical trial was completed. Thirty patients with refractory metastatic colorectal cancer were treated with a combination of 5-FU, oxaliplatin, an epidermal-growth factor receptor monoclonal antibody, and dasatinib. The recommended phase II dose of dasatinib was established, and toxicities were quantified. Pharmacodynamic studies demonstrated increased phosphorylation of the Src substrate paxillin after dasatinib therapy. Tumor biopsies were obtained and Src expression levels were quantitated. Clinical benefit was seen with the combination, including a response rate of 20% and disease control rate of 56%, prompting a larger clinical study. ^ In summary, although Src is constitutively activated in metastatic colorectal cancer, administration of oxaliplatin chemotherapy can further increase its activity, through a reactive oxygen species dependent manner. Inhibition of Src in combination with oxaliplatin provides additional benefit in vitro, in preclinical animal models, and in the clinic. Further study of Src inhibition in the clinic and identification of predictive biomarkers of response will be required to further advance this promising therapeutic target. ^

Graph based study of allergen cross-reactivity of plant lipid transfer proteins (LTPs) using microarray in a multicenter study.

The study of cross-reactivity in allergy is key to both understanding. the allergic response of many patients and providing them with a rational treatment In the present study, protein microarrays and a co-sensitization graph approach were used in conjunction with an allergen microarray immunoassay. This enabled us to include a wide number of proteins and a large number of patients, and to study sensitization profiles among members of the LTP family. Fourteen LTPs from the most frequent plant food-induced allergies in the geographical area studied were printed into a microarray specifically designed for this research. 212 patients with fruit allergy and 117 food-tolerant pollen allergic subjects were recruited from seven regions of Spain with different pollen profiles, and their sera were tested with allergen microarray. This approach has proven itself to be a good tool to study cross-reactivity between members of LTP family, and could become a useful strategy to analyze other families of allergens.

Identification of eukaryotic mRNAs that are translated at reduced cap binding complex eIF4F concentrations using a cDNA microarray

Although most eukaryotic mRNAs need a functional cap binding complex eIF4F for efficient 5′ end- dependent scanning to initiate translation, picornaviral, hepatitis C viral, and a few cellular RNAs have been shown to be translated by internal ribosome entry, a mechanism that can operate in the presence of low levels of functional eIF4F. To identify cellular mRNAs that can be translated when eIF4F is depleted or in low abundance and that, therefore, may contain internal ribosome entry sites, mRNAs that remained associated with polysomes were isolated from human cells after infection with poliovirus and were identified by using a cDNA microarray. Approximately 200 of the 7000 mRNAs analyzed remained associated with polysomes under these conditions. Among the gene products encoded by these polysome-associated mRNAs were immediate-early transcription factors, kinases, and phosphatases of the mitogen-activated protein kinase pathways and several protooncogenes, including c-myc and Pim-1. In addition, the mRNA encoding Cyr61, a secreted factor that can promote angiogenesis and tumor growth, was selectively mobilized into polysomes when eIF4F concentrations were reduced, although its overall abundance changed only slightly. Subsequent tests confirmed the presence of internal ribosome entry sites in the 5′ noncoding regions of both Cyr61 and Pim-1 mRNAs. Overall, this study suggests that diverse mRNAs whose gene products have been implicated in a variety of stress responses, including inflammation, angiogenesis, and the response to serum, can use translational initiation mechanisms that require little or no intact cap binding protein complex eIF4F.

Analysis of transforming activity of human synovial sarcoma-associated chimeric protein SYT-SSX1 bound to chromatin remodeling factor hBRM/hSNF2α

Human synovial sarcoma has been shown to exclusively harbor the chromosomal translocation t(X;18) that produces the chimeric gene SYT-SSX. However, the role of SYT-SSX in cellular transformation remains unclear. In this study, we have established 3Y1 rat fibroblast cell lines that constitutively express SYT, SSX1, and SYT-SSX1 and found that SYT-SSX1 promoted growth rate in culture, anchorage-independent growth in soft agar, and tumor formation in nude mice. Deletion of the N-terminal 181 amino acids of SYT-SSX1 caused loss of its transforming activity. Furthermore, association of SYT-SSX1 with the chromatin remodeling factor hBRM/hSNF2α, which regulates transcription, was demonstrated in both SYT-SSX1-expressing 3Y1 cells and in the human synovial sarcoma cell line HS-SY-II. The binding region between the two molecules was shown to reside within the N-terminal 181 amino acids stretch (aa 1–181) of SYT-SSX1 and 50 amino acids (aa 156–205) of hBRM/hSNF2α and we found that the overexpression of this binding region of hBRM/hSNF2α significantly suppressed the anchorage-independent growth of SYT-SSX1-expressing 3Y1 cells. To analyze the transcriptional regulation by SYT-SSX1, we established conditional expression system of SYT-SSX1 and examined the gene expression profiles. The down-regulation of potential tumor suppressor DCC was observed among 1,176 genes analyzed by microarray analysis, and semi-quantitative reverse transcription–PCR confirmed this finding. These data clearly demonstrate transforming activity of human oncogene SYT-SSX1 and also involvement of chromatin remodeling factor hBRM/hSNF2α in human cancer.

Tissue-specific gene expression in soybean (Glycine max) detected by cDNA microarray analysis

We have constructed cDNA microarrays for soybean (Glycine max L. Merrill), containing approximately 4,100 Unigene ESTs derived from axenic roots, to evaluate their application and utility for functional genomics of organ differentiation in legumes. We assessed microarray technology by conducting studies to evaluate the accuracy of microarray data and have found them to be both reliable and reproducible in repeat hybridisations. Several ESTs showed high levels (>50 fold) of differential expression in either root or shoot tissue of soybean. A small number of physiologically interesting, and differentially expressed sequences found by microarray analysis were verified by both quantitative real-time RT-PCR and Northern blot analysis. There was a linear correlation (r(2) = 0.99, over 5 orders of magnitude) between microarray and quantitative real-time RT-PCR data. Microarray analysis of soybean has enormous potential not only for the discovery of new genes involved in tissue differentiation and function, but also to study the expression of previously characterised genes, gene networks and gene interactions in wild-type, mutant or transgenic; plants.