Gene expression databases and data mining.

The DNA microarray technology has arguably caught the attention of the worldwide life science community and is now systematically supporting major discoveries in many fields of study. The majority of the initial technical challenges of conducting experiments are being resolved, only to be replaced with new informatics hurdles, including statistical analysis, data visualization, interpretation, and storage. Two systems of databases, one containing expression data and one containing annotation data are quickly becoming essential knowledge repositories of the research community. This present paper surveys several databases, which are considered "pillars" of research and important nodes in the network. This paper focuses on a generalized workflow scheme typical for microarray experiments using two examples related to cancer research. The workflow is used to reference appropriate databases and tools for each step in the process of array experimentation. Additionally, benefits and drawbacks of current array databases are addressed, and suggestions are made for their improvement.

Classification of human astrocytic gliomas on the basis of gene expression: a correlated group of genes with angiogenic activity emerges as a strong predictor of subtypes.

The development of targeted treatment strategies adapted to individual patients requires identification of the different tumor classes according to their biology and prognosis. We focus here on the molecular aspects underlying these differences, in terms of sets of genes that control pathogenesis of the different subtypes of astrocytic glioma. By performing cDNA-array analysis of 53 patient biopsies, comprising low-grade astrocytoma, secondary glioblastoma (respective recurrent high-grade tumors), and newly diagnosed primary glioblastoma, we demonstrate that human gliomas can be differentiated according to their gene expression. We found that low-grade astrocytoma have the most specific and similar expression profiles, whereas primary glioblastoma exhibit much larger variation between tumors. Secondary glioblastoma display features of both other groups. We identified several sets of genes with relatively highly correlated expression within groups that: (a). can be associated with specific biological functions; and (b). effectively differentiate tumor class. One prominent gene cluster discriminating primary versus nonprimary glioblastoma comprises mostly genes involved in angiogenesis, including VEGF fms-related tyrosine kinase 1 but also IGFBP2, that has not yet been directly linked to angiogenesis. In situ hybridization demonstrating coexpression of IGFBP2 and VEGF in pseudopalisading cells surrounding tumor necrosis provided further evidence for a possible involvement of IGFBP2 in angiogenesis. The separating groups of genes were found by the unsupervised coupled two-way clustering method, and their classification power was validated by a supervised construction of a nearly perfect glioma classifier.

MYC regulation of a "poor-prognosis" metastatic cancer cell state.

Gene expression signatures are used in the clinic as prognostic tools to determine the risk of individual patients with localized breast tumors developing distant metastasis. We lack a clear understanding, however, of whether these correlative biomarkers link to a common biological network that regulates metastasis. We find that the c-MYC oncoprotein coordinately regulates the expression of 13 different "poor-outcome" cancer signatures. In addition, functional inactivation of MYC in human breast cancer cells specifically inhibits distant metastasis in vivo and invasive behavior in vitro of these cells. These results suggest that MYC oncogene activity (as marked by "poor-prognosis" signature expression) may be necessary for the translocation of poor-outcome human breast tumors to distant sites.

Selecting control genes for RT-QPCR using public microarray data.

Background: Gene expression analysis has emerged as a major biological research area, with real-time quantitative reverse transcription PCR (RT-QPCR) being one of the most accurate and widely used techniques for expression profiling of selected genes. In order to obtain results that are comparable across assays, a stable normalization strategy is required. In general, the normalization of PCR measurements between different samples uses one to several control genes (e. g. housekeeping genes), from which a baseline reference level is constructed. Thus, the choice of the control genes is of utmost importance, yet there is not a generally accepted standard technique for screening a large number of candidates and identifying the best ones. Results: We propose a novel approach for scoring and ranking candidate genes for their suitability as control genes. Our approach relies on publicly available microarray data and allows the combination of multiple data sets originating from different platforms and/or representing different pathologies. The use of microarray data allows the screening of tens of thousands of genes, producing very comprehensive lists of candidates. We also provide two lists of candidate control genes: one which is breast cancer-specific and one with more general applicability. Two genes from the breast cancer list which had not been previously used as control genes are identified and validated by RT-QPCR. Open source R functions are available at http://www.isrec.isb-sib.ch/similar to vpopovic/research/ Conclusion: We proposed a new method for identifying candidate control genes for RT-QPCR which was able to rank thousands of genes according to some predefined suitability criteria and we applied it to the case of breast cancer. We also empirically showed that translating the results from microarray to PCR platform was achievable.

Comparative transcriptomics of rice reveals an ancient pattern of response to microbial colonization.

Glomalean fungi induce and colonize symbiotic tissue called arbuscular mycorrhiza on the roots of most land plants. Other fungi also colonize plants but cause disease not symbiosis. Whole-transcriptome analysis using a custom-designed Affymetrix Gene-Chip and confirmation with real-time RT-PCR revealed 224 genes affected during arbuscular mycorrhizal symbiosis. We compared these transcription profiles with those from rice roots that were colonized by pathogens (Magnaporthe grisea and Fusarium moniliforme). Over 40% of genes showed differential regulation caused by both the symbiotic and at least one of the pathogenic interactions. A set of genes was similarly expressed in all three associations, revealing a conserved response to fungal colonization. The responses that were shared between pathogen and symbiont infection may play a role in compatibility. Likewise, the responses that are different may cause disease. Some of the genes that respond to mycorrhizal colonization may be involved in the uptake of phosphate. Indeed, phosphate addition mimicked the effect of mycorrhiza on 8% of the tested genes. We found that 34% of the mycorrhiza-associated rice genes were also associated with mycorrhiza in dicots, revealing a conserved pattern of response between the two angiosperm classes.

Host and invader impact of transfer of the clc genomic island into Pseudomonas aeruginosa PAO1

Genomic islands, large potentially mobile regions of bacterial chromosomes, are a major contributor to bacteria evolution. Here, we investigated the fitness cost and phenotypic differences between the bacterium Pseudomonas aeruginosa PAO1 and a derivative carrying one integrated copy of the clc element, a 103-kb genomic island [and integrative and conjugative element (ICE)] originating in Pseudomonas sp. strain B13 and a close relative of genomic islands found in clinical and environmental isolates of P. aeruginosa. By using a combination of whole genome transcriptome profiling, phenotypic arrays, competition experiments, and biofilm formation studies, only few differences became apparent, such as reduced biofilm growth and fourfold stationary phase repression of genes involved in acetoin metabolism in PAO1 containing the clc element. In contrast, PAO1 carrying the clc element acquired the capacity to grow on 3-chlorobenzoate and 2-aminophenol as sole carbon and energy substrates. No fitness loss >1% was detectable in competition experiments between PAO1 and PAO1 carrying the clc element. The genes from the clc element were not silent in PAO1, and excision was observed, although transfer of clc from PAO1 to other recipient bacteria was reduced by two orders of magnitude. Our results indicate that newly acquired mobile DNA not necessarily invoke an important fitness cost on their host. Absence of immediate detriment to the host may have contributed to the wide distribution of genomic islands like clc in bacterial genomes

Defining the role of common variation in the genomic and biological architecture of adult human height.

Using genome-wide data from 253,288 individuals, we identified 697 variants at genome-wide significance that together explained one-fifth of the heritability for adult height. By testing different numbers of variants in independent studies, we show that the most strongly associated ∼2,000, ∼3,700 and ∼9,500 SNPs explained ∼21%, ∼24% and ∼29% of phenotypic variance. Furthermore, all common variants together captured 60% of heritability. The 697 variants clustered in 423 loci were enriched for genes, pathways and tissue types known to be involved in growth and together implicated genes and pathways not highlighted in earlier efforts, such as signaling by fibroblast growth factors, WNT/β-catenin and chondroitin sulfate-related genes. We identified several genes and pathways not previously connected with human skeletal growth, including mTOR, osteoglycin and binding of hyaluronic acid. Our results indicate a genetic architecture for human height that is characterized by a very large but finite number (thousands) of causal variants.

STAT3 links IL-22 signaling in intestinal epithelial cells to mucosal wound healing.

Signal transducer and activator of transcription (STAT) 3 is a pleiotropic transcription factor with important functions in cytokine signaling in a variety of tissues. However, the role of STAT3 in the intestinal epithelium is not well understood. We demonstrate that development of colonic inflammation is associated with the induction of STAT3 activity in intestinal epithelial cells (IECs). Studies in genetically engineered mice showed that epithelial STAT3 activation in dextran sodium sulfate colitis is dependent on interleukin (IL)-22 rather than IL-6. IL-22 was secreted by colonic CD11c(+) cells in response to Toll-like receptor stimulation. Conditional knockout mice with an IEC-specific deletion of STAT3 activity were highly susceptible to experimental colitis, indicating that epithelial STAT3 regulates gut homeostasis. STAT3(IEC-KO) mice, upon induction of colitis, showed a striking defect of epithelial restitution. Gene chip analysis indicated that STAT3 regulates the cellular stress response, apoptosis, and pathways associated with wound healing in IECs. Consistently, both IL-22 and epithelial STAT3 were found to be important in wound-healing experiments in vivo. In summary, our data suggest that intestinal epithelial STAT3 activation regulates immune homeostasis in the gut by promoting IL-22-dependent mucosal wound healing.

Notch signaling is required for normal prostatic epithelial cell proliferation and differentiation.

Notch pathway is crucial for stem/progenitor cell maintenance, growth and differentiation in a variety of tissues. Using a transgenic cell ablation approach, we found in our previous study that cells expressing Notch1 are crucial for prostate early development and re-growth. Here, we further define the role of Notch signaling in regulating prostatic epithelial cell growth and differentiation using biochemical and genetic approaches in ex vivo or in vivo systems. Treatment of developing prostate grown in culture with inhibitors of gamma-secretase/presenilin, which is required for Notch cleavage and activation, caused a robust increase in proliferation of epithelial cells co-expressing cytokeratin 8 and 14, lack of luminal/basal layer segregation and dramatically reduced branching morphogenesis. Using conditional Notch1 gene deletion mouse models, we found that inactivation of Notch1 signaling resulted in profound prostatic alterations, including increased tufting, bridging and enhanced epithelial proliferation. Cells within these lesions co-expressed both luminal and basal cell markers, a feature of prostatic epithelial cells in predifferentiation developmental stages. Microarray analysis revealed that the gene expression in a number of genetic networks was altered following Notch1 gene deletion in prostate. Furthermore, expression of Notch1 and its effector Hey-1 gene in human prostate adenocarcinomas were found significantly down-regulated compared to normal control tissues. Taken together, these data suggest that Notch signaling is critical for normal cell proliferation and differentiation in the prostate, and deregulation of this pathway may facilitate prostatic tumorigenesis.

Changes in the transcriptional profile of transporters in the intestine along the anterior-posterior and crypt-villus axes.

BACKGROUND: The purpose of this work was to characterize the expression of drug and nutrient carriers along the anterior-posterior and crypt-villus axes of the intestinal epithelium and to study the validity of utilizing whole gut tissue rather than purified epithelial cells to examine regional variations in gene expression. RESULTS: We have characterized the mRNA expression profiles of 76 % of all currently known transporters along the anterior-posterior axis of the gut. This is the first study to describe the expression profiles of the majority of all known transporters in the intestine. The expression profiles of transporters, as defined according to the Gene Ontology consortium, were measured in whole tissue of the murine duodenum, jejunum, ileum and colon using high-density microarrays. For nine transporters (Abca1, Abcc1, Abcc3, Abcg8, Slc10a2, Slc28a2, Slc2a1, Slc34a2 and Slc5a8), the mRNA profiles were further measured by RT-PCR in laser micro-dissected crypt and villus epithelial cells corresponding to the aforementioned intestinal regions. With respect to differentially regulated transporters, the colon had a distinct expression profile from small intestinal segments. The majority (59 % for p cutoff < or = 0.05) of transporter mRNA levels were constant across the intestinal sections studied. For the transporter subclass "carrier activity", which contains the majority of known carriers for biologically active compounds, a significant change (p < or = 0.05) along the anterior-posterior axis was observed. CONCLUSION: All nine transporters examined in laser-dissected material demonstrated good replication of the region-specific profiles revealed by microarray. Furthermore, we suggest that the distribution characteristics of Slc5a8 along the intestinal tract render it a suitable candidate carrier for monocarboxylate drugs in the posterior portion of the intestine. Our findings also predict that there is a significant difference in the absorption of carrier-mediated compounds in the different intestinal segments. The most pronounced differences can be expected between the adjoining segments ileum and colon, but the differences between the other adjoining segments are not negligible. Finally, for the examined genes, profiles measured in whole intestinal tissue extracts are representative of epithelial cell-only gene expression.

Histological and histochemical evaluation of human oral mucosa constructs developed by tissue engineering

Reconstruction of large oral mucosa defects is often challenging, since the shortage of healthy oral mucosa to replace the excised tissues is very common. In this context, tissue engineering techniques may provide a source of autologous tissues available for transplant in these patients. In this work, we developed a new model of artificial oral mucosa generated by tissue engineering using a fibrin-agarose scaffold. For that purpose, we generated primary cultures of human oral mucosa fibroblasts and keratinocytes from small biopsies of normal oral mucosa using enzymatic treatments. Then we determined the viability of the cultured cells by electron probe quantitative X-ray microanalysis, and we demonstrated that most of the cells in the primary cultures were alive and had high K/Na ratios. Once cell viability was determined, we used the cultured fibroblasts and keratinocytes to develop an artificial oral mucosa construct by using a fibrin-agarose extracellular matrix and a sequential culture technique using porous culture inserts. Histological analysis of the artificial tissues showed high similarities with normal oral mucosa controls. The epithelium of the oral substitutes had several layers, with desmosomes and apical microvilli and microplicae. Both the controls and the oral mucosa substitutes showed high suprabasal expression of cytokeratin 13 and low expression of cytokeratin 10. All these results suggest that our model of oral mucosa using fibrin-agarose scaffolds show several similarities with native human oral mucosa.

Fas ligand-induced proinflammatory transcriptional responses in reconstructed human epidermis. Recruitment of the epidermal growth factor receptor and activation of MAP kinases.

Fas ligand (FasL) exerts potent proapoptotic and proinflammatory actions on epidermal keratinocytes and has been implicated in the pathogenesis of eczema, toxic epidermal necrolysis, and drug-induced skin eruptions. We used reconstructed human epidermis to investigate the mechanisms of FasL-induced inflammatory responses and their relationships with FasL-triggered caspase activity. Caspase activity was a potent antagonist of the pro-inflammatory gene expression triggered by FasL prior to the onset of cell death. Furthermore, we found that FasL-stimulated autocrine production of epidermal growth factor receptor (EGFR) ligands, and the subsequent activation of EGFR and ERK1 and ERK2 mitogen-activated protein kinases, were obligatory extracellular steps for the FasL-induced expression of a subset of inflammatory mediators, including CXCL8/interleukin (IL)-8, ICAM-1, IL-1alpha, IL-1beta, CCL20/MIP-3alpha, and thymic stromal lymphopoietin. These results expand the known physiological role of EGFR and its ligands from promoting keratinocyte mitogenesis and survival to mediating FasL-induced epidermal inflammation.

Increased expression of urokinase-type plasminogen activator mRNA determines adverse prognosis in ErbB2-positive primary breast cancer.

PURPOSE: To evaluate and validate mRNA expression markers capable of identifying patients with ErbB2-positive breast cancer associated with distant metastasis and reduced survival. PATIENTS AND METHODS: Expression of 60 genes involved in breast cancer biology was assessed by quantitative real-time PCR (qrt-PCR) in 317 primary breast cancer patients and correlated with clinical outcome data. Results were validated subsequently using two previously published and publicly available microarray data sets with different patient populations comprising 295 and 286 breast cancer samples, respectively. RESULTS: Of the 60 genes measured by qrt-PCR, urokinase-type plasminogen activator (uPA or PLAU) mRNA expression was the most significant marker associated with distant metastasis-free survival (MFS) by univariate Cox analysis in patients with ErbB2-positive tumors and an independent factor in multivariate analysis. Subsequent validation in two microarray data sets confirmed the prognostic value of uPA in ErbB2-positive tumors by both univariate and multivariate analysis. uPA mRNA expression was not significantly associated with MFS in ErbB2-negative tumors. Kaplan-Meier analysis showed in all three study populations that patients with ErbB2-positive/uPA-positive tumors exhibited significantly reduced MFS (hazard ratios [HR], 4.3; 95% CI, 1.6 to 11.8; HR, 2.7; 95% CI, 1.2 to 6.2; and, HR, 2.8; 95% CI, 1.1 to 7.1; all P < .02) as compared with the group with ErbB2-positive/uPA-negative tumors who exhibited similar outcome to those with ErbB2-negative tumors, irrespective of uPA status. CONCLUSION: After evaluation of 898 breast cancer patients, uPA mRNA expression emerged as a powerful prognostic indicator in ErbB2-positive tumors. These results were consistent among three independent study populations assayed by different techniques, including qrt-PCR and two microarray platforms.

Kallikrein genes are associated with lupus and glomerular basement membrane-specific antibody-induced nephritis in mice and humans.

Immune-mediated nephritis contributes to disease in systemic lupus erythematosus, Goodpasture syndrome (caused by antibodies specific for glomerular basement membrane [anti-GBM antibodies]), and spontaneous lupus nephritis. Inbred mouse strains differ in susceptibility to anti-GBM antibody-induced and spontaneous lupus nephritis. This study sought to clarify the genetic and molecular factors that maybe responsible for enhanced immune-mediated renal disease in these models. When the kidneys of 3 mouse strains sensitive to anti-GBM antibody-induced nephritis were compared with those of 2 control strains using microarray analysis, one-fifth of the underexpressed genes belonged to the kallikrein gene family,which encodes serine esterases. Mouse strains that upregulated renal and urinary kallikreins exhibited less evidence of disease. Antagonizing the kallikrein pathway augmented disease, while agonists dampened the severity of anti-GBM antibody-induced nephritis. In addition, nephritis-sensitive mouse strains had kallikrein haplotypes that were distinct from those of control strains, including several regulatory polymorphisms,some of which were associated with functional consequences. Indeed, increased susceptibility to anti-GBM antibody-induced nephritis and spontaneous lupus nephritis was achieved by breeding mice with a genetic interval harboring the kallikrein genes onto a disease-resistant background. Finally, both human SLE and spontaneous lupus nephritis were found to be associated with kallikrein genes, particularly KLK1 and the KLK3 promoter, when DNA SNPs from independent cohorts of SLE patients and controls were compared. Collectively, these studies suggest that kallikreins are protective disease-associated genes in anti-GBM antibody-induced nephritis and lupus.