cDNA microarrays detect activation of a myogenic transcription program by the PAX3-FKHR fusion oncogene

Alveolar rhabdomyosarcoma is an aggressive pediatric cancer of striated muscle characterized in 60% of cases by a t(2;13)(q35;q14). This results in the fusion of PAX3, a developmental transcription factor required for limb myogenesis, with FKHR, a member of the forkhead family of transcription factors. The resultant PAX3-FKHR gene possesses transforming properties; however, the effects of this chimeric oncogene on gene expression are largely unknown. To investigate the actions of these transcription factors, both Pax3 and PAX3-FKHR were introduced into NIH 3T3 cells, and the resultant gene expression changes were analyzed with a murine cDNA microarray containing 2,225 elements. We found that PAX3-FKHR but not PAX3 activated a myogenic transcription program including the induction of transcription factors MyoD, Myogenin, Six1, and Slug as well as a battery of genes involved in several aspects of muscle function. Notable among this group were the growth factor gene Igf2 and its binding protein Igfbp5. Relevance of this model was suggested by verification that three of these genes (IGFBP5, HSIX1, and Slug) were also expressed in alveolar rhabdomyosarcoma cell lines. This study utilizes cDNA microarrays to elucidate the pattern of gene expression induced by an oncogenic transcription factor and demonstrates the profound myogenic properties of PAX3-FKHR in NIH 3T3 cells.

Assessment of clone identity and sequence fidelity for 1189 IMAGE cDNA clones

This report documents the error rate in a commercially distributed subset of the IMAGE Consortium mouse cDNA clone collection. After isolation of plasmid DNA from 1189 bacterial stock cultures, only 62.2% were uncontaminated and contained cDNA inserts that had significant sequence identity to published data for the ordered clones. An agarose gel electrophoresis pre-screening strategy identified 361 stock cultures that appeared to contain two or more plasmid species. Isolation of individual colonies from these stocks demonstrated that 7.1% of the original 1189 stocks contained both a correct and an incorrect plasmid. 5.9% of the original 1189 stocks contained multiple, distinct, incorrect plasmids, indicating the likelihood of multiple contaminating events. While only 739 of the stocks purchased contained the desired cDNA clone, agarose gel pre-screening, colony isolation and similarity searching of dbEST allowed for the identification of an additional 420 clones that would have otherwise been discarded. Considering the high error rate in this subset of the IMAGE cDNA clone set, the use of sequence verified clones for cDNA microarray construction is warranted. When this is not possible, pre-screening non-sequence verified clones with agarose gel electrophoresis provides an inexpensive and efficient method to eliminate contaminated clones from the probe set.

An evaluation of the performance of cDNA microarrays for detecting changes in global mRNA expression

The cDNA microarray is one technological approach that has the potential to accurately measure changes in global mRNA expression levels. We report an assessment of an optimized cDNA microarray platform to generate accurate, precise and reliable data consistent with the objective of using microarrays as an acquisition platform to populate gene expression databases. The study design consisted of two independent evaluations with 70 arrays from two different manufactured lots and used three human tissue sources as samples: placenta, brain and heart. Overall signal response was linear over three orders of magnitude and the sensitivity for any element was estimated to be 2 pg mRNA. The calculated coefficient of variation for differential expression for all non-differentiated elements was 12–14% across the entire signal range and did not vary with array batch or tissue source. The minimum detectable fold change for differential expression was 1.4. Accuracy, in terms of bias (observed minus expected differential expression ratio), was less than 1 part in 10 000 for all non-differentiated elements. The results presented in this report demonstrate the reproducible performance of the cDNA microarray technology platform and the methods provide a useful framework for evaluating other technologies that monitor changes in global mRNA expression.

Host microarray analysis reveals a role for the Salmonella response regulator phoP in human macrophage cell death

Bacterial pathogens manipulate host cells to promote pathogen survival and dissemination. We used a 22,571 human cDNA microarray to identify host pathways that are affected by the Salmonella enterica subspecies typhimurium phoP gene, a transcription factor required for virulence, by comparing the expression profiles of human monocytic tissue culture cells infected with either the wild-type bacteria or a phoP∷Tn10 mutant strain. Both wild-type and phoP∷Tn10 bacteria induced a common set of genes, many of which are proinflammatory. Differentially expressed genes included those that affect host cell death, suggesting that the phoP regulatory system controls bacterial genes that alter macrophage survival. Subsequent experiments showed that the phoP∷Tn10 mutant strain is defective for killing both cultured and primary human macrophages but is able to replicate intracellularly. These experiments indicate that phoP plays a role in Salmonella-induced human macrophage cell death.

The full Bayesian significance test for mixture models: results in gene expression clustering

Gene clustering is a useful exploratory technique to group together genes with similar expression levels under distinct cell cycle phases or distinct conditions. It helps the biologist to identify potentially meaningful relationships between genes. In this study, we propose a clustering method based on multivariate normal mixture models, where the number of clusters is predicted via sequential hypothesis tests: at each step, the method considers a mixture model of m components (m = 2 in the first step) and tests if in fact it should be m - 1. If the hypothesis is rejected, m is increased and a new test is carried out. The method continues (increasing m) until the hypothesis is accepted. The theoretical core of the method is the full Bayesian significance test, an intuitive Bayesian approach, which needs no model complexity penalization nor positive probabilities for sharp hypotheses. Numerical experiments were based on a cDNA microarray dataset consisting of expression levels of 205 genes belonging to four functional categories, for 10 distinct strains of Saccharomyces cerevisiae. To analyze the method's sensitivity to data dimension, we performed principal components analysis on the original dataset and predicted the number of classes using 2 to 10 principal components. Compared to Mclust (model-based clustering), our method shows more consistent results.

Searching for molecular markers in head and neck squamous cell carcinomas (HNSCC) by statistical and bioinformatic analysis of larynx-derived SAGE libraries

Background: Head and neck squamous cell carcinoma (HNSCC) is one of the most common malignancies in humans. The average 5-year survival rate is one of the lowest among aggressive cancers, showing no significant improvement in recent years. When detected early, HNSCC has a good prognosis, but most patients present metastatic disease at the time of diagnosis, which significantly reduces survival rate. Despite extensive research, no molecular markers are currently available for diagnostic or prognostic purposes. Methods: Aiming to identify differentially-expressed genes involved in laryngeal squamous cell carcinoma (LSCC) development and progression, we generated individual Serial Analysis of Gene Expression (SAGE) libraries from a metastatic and non-metastatic larynx carcinoma, as well as from a normal larynx mucosa sample. Approximately 54,000 unique tags were sequenced in three libraries. Results: Statistical data analysis identified a subset of 1,216 differentially expressed tags between tumor and normal libraries, and 894 differentially expressed tags between metastatic and non-metastatic carcinomas. Three genes displaying differential regulation, one down-regulated (KRT31) and two up-regulated (BST2, MFAP2), as well as one with a non-significant differential expression pattern (GNA15) in our SAGE data were selected for real-time polymerase chain reaction (PCR) in a set of HNSCC samples. Consistent with our statistical analysis, quantitative PCR confirmed the upregulation of BST2 and MFAP2 and the downregulation of KRT31 when samples of HNSCC were compared to tumor-free surgical margins. As expected, GNA15 presented a non-significant differential expression pattern when tumor samples were compared to normal tissues. Conclusion: To the best of our knowledge, this is the first study reporting SAGE data in head and neck squamous cell tumors. Statistical analysis was effective in identifying differentially expressed genes reportedly involved in cancer development. The differential expression of a subset of genes was confirmed in additional larynx carcinoma samples and in carcinomas from a distinct head and neck subsite. This result suggests the existence of potential common biomarkers for prognosis and targeted-therapy development in this heterogeneous type of tumor.

Gene Expression Profiles in Radiation Workers Occupationally Exposed to Ionizing Radiation

Ionizing radiation OR) imposes risks to human health and the environment. IR at low doses and low (lose rates has the potency to initiate carcinogenesis. Genotoxic environmental agents such as IR trigger a cascade of signal transduction pathways for cellular protection. In this study, using cDNA microarray technique, we monitored the gene expression profiles in lymphocytes derived from radiation-ex posed individuals (radiation workers). Physical dosimetry records on these patients indicated that the absorbed dose ranged from 0.696 to 39.088 mSv. Gene expression analysis revealed statistically significant transcriptional changes in a total of 78 genes (21 up-regulated and 57 clown-regulated) involved in several biological processes such as ubiquitin cycle (UHRF2 and PIAS1), DNA repair (LIG3, XPA, ERCC5, RAD52, DCLRE1C), cell cycle regulation/proliferation (RHOA, CABLES2, TGFB2, IL16), and stress response (GSTP1, PPP2R5A, DUSP22). Some of the genes that showed altered expression profiles in this study call be used as biomarkers for monitoring the chronic low level exposure in humans. Additionally, alterations in gene expression patterns observed in chronically exposed radiation workers reinforces the need for defining the effective radiation dose that causes immediate genetic damage as well as the long-term effects on genomic instability, including cancer.

A role for mammalian target of rapamycin (mTOR) pathway in non alcoholic steatohepatitis related-cirrhosis

Non-alcoholic fatty liver disease (NAFLD) encompasses the whole spectrum of steatosis, nonalcoholic steatohepatitis (NASH), and NASH-related cirrhosis (NASH/Cir). Although molecular advances have been made in this field, the pathogenesis of NAFLD is not completely understood. The gene expression profiling associated to NASH/Cir was assessed, in an attempt to better characterize the pathways involved in its etiopathogenesis. Methods: In the first step, we used cDNA microarray to evaluate the gene expression profiles in normal liver (n=3) and NASH/Cir samples (n=3) by GeneSifter (TM) analysis to identify differentially expressed genes and biological pathways. Second, tissue microarray was used to determine immunohistochemical expression of phosphorylated mTOR and 4E-BP1 in 11 normal liver samples, 10 NASH/Cir samples and in 37 samples of cirrhosis of other etiologies to further explore the involvement of the mTOR pathway evidenced by the gene expression analysis. Results: 138 and 106 genes were, respectively, up and down regulated in NASH/Cir in comparison to normal liver. Among the 9 pathways identified as significantly modulated in NASH/Cir, the participation of the mTOR pathway was confirmed, since expression of cytoplasmic and membrane phospho-mTOR were higher in NASH/Cir in comparison to cirrhosis of other etiologies and to normal liver. Conclusions: Recent findings have suggested a role for the cellular ""nutrient sensor"" mTOR in NAFLD and the present study corroborates the participation of this pathway in NASH/Cir. Phospho-mTOR evaluation might be of clinical utility as a potential marker for identification of NASH/Cir in cases mistakenly considered as cryptogenic cirrhosis owing to paucity of clinical data.

No-match ORESTES explored as tumor markers

Sequencing technologies and new bioinformatics tools have led to the complete sequencing of various genomes. However, information regarding the human transcriptome and its annotation is yet to be completed. The Human Cancer Genome Project, using ORESTES (open reading frame EST sequences) methodology, contributed to this objective by generating data from about 1.2 million expressed sequence tags. Approximately 30 of these sequences did not align to ESTs in the public databases and were considered no-match ORESTES. On the basis that a set of these ESTs could represent new transcripts, we constructed a cDNA microarray. This platform was used to hybridize against 12 different normal or tumor tissues. We identified 3421 transcribed regions not associated with annotated transcripts, representing 83.3 of the platform. The total number of differentially expressed sequences was 1007. Also, 28 of analyzed sequences could represent noncoding RNAs. Our data reinforces the knowledge of the human genome being pervasively transcribed, and point out molecular marker candidates for different cancers. To reinforce our data, we confirmed, by real-time PCR, the differential expression of three out of eight potentially tumor markers in prostate tissues. Lists of 1007 differentially expressed sequences, and the 291 potentially noncoding tumor markers were provided.

Differential gene expression between the biotrophic-like and saprotrophic mycelia of the witches` broom pathogen Moniliophthora perniciosa

Moniliophthora perniciosa is a hemibiotrophic fungus that causes witches` broom disease (WBD) in cacao. Marked dimorphism characterizes this fungus, showing a monokaryotic or biotrophic phase that causes disease symptoms and a later dikaryotic or saprotrophic phase. A combined strategy of DNA microarray, expressed sequence tag, and real-time reverse-transcriptase polymerase chain reaction analyses was employed to analyze differences between these two fungal stages in vitro. In all, 1,131 putative genes were hybridized with cDNA from different phases, resulting in 189 differentially expressed genes, and 4,595 reads were clusterized, producing 1,534 unigenes. The analysis of these genes, which represent approximately 21% of the total genes, indicates that the biotrophic-like phase undergoes carbon and nitrogen catabollite repression that correlates to the expression of phytopathogenicity genes. Moreover, downregulation of mitochondrial oxidative phosphorylation and the presence of a putative ngr1 of Saccharomyces cerevisiae could help explain its lower growth rate. In contrast, the saprotrophic mycelium expresses genes related to the metabolism of hexoses, ammonia, and oxidative phosphorylation, which could explain its faster growth. Antifungal toxins were upregulated and could prevent the colonization by competing fungi. This work significantly contributes to our understanding of the molecular mechanisms of WBD and, to our knowledge, is the first to analyze differential gene expression of the different phases of a hemibiotrophic fungus.

Molecular profiling of isolated histological components of Wilms tumor implicates a common role for the Wnt signaling pathway in kidney and tumor development

Wilms tumor (WT), a tumor composed of three histological components - blastema (BL), epithelia and stroma - is considered an appropriate model system to study the biological relationship between differentiation and tumorigenesis. To investigate molecular associations between nephrogenesis and WT, the gene expression pattern of individual cellular components was analyzed, using a customized platform containing 4,608 genes. WT gene expression patterns were compared to genes regulated during kidney differentiation. BL had a closer gene expression pattern to the earliest stage of normal renal development. The BL gene expression pattern was compared to that of fetal kidney (FK) and also between FK and mature kidney, identifying 25 common de-regulated genes supposedly involved in the earliest events of WT onset. Quantitative RT-PCR was performed, confirming the difference in expression levels for 13 of 16 genes (81.2%) in the initial set and 8 of 13 (61.5%) in an independent set of samples. An overrepresentation of genes belonging to the Wnt signaling pathway was identified, namely PLCG2, ROCK2 and adenomatous polyposis coli (APC). Activation of the Wnt pathway was confirmed in WT, using APC at protein level and PLCG2 at mRNA and protein level. APC showed positive nuclear immunostaining for an independent set of WT samples, similarly to the FK in week 11. Lack of PLCG2 expression was confirmed in WT and in FK until week 18. Taken together, these results provided molecular evidence of the recapitulation of the embryonic kidney by WT as well as involvement of the Wnt pathway in the earliest events of WT onset. Copyright (C) 2008 S. Karger AG, Basel.

Influence of the interaction between nodal fibroblast and breast cancer cells on gene expression

Our aim was to evaluate the interaction between breast cancer cells and nodal fibroblasts, by means of their gene expression profile. Fibroblast primary cultures were established from negative and positive lymph nodes from breast cancer patients and a similar gene expression pattern was identified, following cell culture. Fibroblasts and breast cancer cells (MDA-MB231, MDA-MB435, and MCF7) were cultured alone or co-cultured separated by a porous membrane (which allows passage of soluble factors) for comparison. Each breast cancer lineage exerted a particular effect on fibroblasts viability and transcriptional profile. However, fibroblasts from positive and negative nodes had a parallel transcriptional behavior when co-cultured with a specific breast cancer cell line. The effects of nodal fibroblasts on breast cancer cells were also investigated. MDA MB-231 cells viability and migration were enhanced by the presence of fibroblasts and accordingly, MDA-MB435 and MCF7 cells viability followed a similar pattern. MDA-MB231 gene expression profile, as evaluated by cDNA microarray, was influenced by the fibroblasts presence, and HNMT, COMT, FN3K, and SOD2 were confirmed downregulated in MDA-MB231 co-cultured cells with fibroblasts from both negative and positive nodes, in a new series of RT-PCR assays. In summary, transcriptional changes induced in breast cancer cells by fibroblasts from positive as well as negative nodes are very much alike in a specific lineage. However, fibroblasts effects are distinct in each one of the breast cancer lineages, suggesting that the inter-relationships between stromal and malignant cells are dependent on the intrinsic subtype of the tumor.

Gene expression profile of residual breast cancer after doxorubicin and cyclophosphamide neoadjuvant chemotherapy

In breast cancer patients, primary chemotherapy is associated with the same survival benefits as adjuvant chemotherapy. Residual tumors represent a clinical challenge, Lis they may be resistant to additional cycles of the same drugs. Our aim was to identify differential transcripts expressed in residual tumors, after neoadjuvant chemotherapy, that might be related with tumor resistance. Hence, 16 patients with paired tumor samples, collected before and after treatment (4 cycles doxorubicin/cyclophosphamide, AC) had their gene expression evaluated on cDNA microarray slides containing 4,608 genes. Three hundred and eighty-nine genes were differentially expressed (paired Student`s t-test, pFDR<0.01) between pre- and post-chemotherapy samples and among the regulated functions were the JNK cascade and cell death. Unsupervised hierarchical clustering identified one branch comprising exclusively, eight pre-chemotherapy samples and another branch, including the former correspondent eight post-chemotherapy samples and other 16 paired pre/post-chemotherapy samples. No differences in clinical and tumor parameters could explain this clustering. Another group of I I patients with paired samples had expression of selected genes determined by real-time RT-PCR and CTGF and DUSP1 were confirmed more expressed in post- as compared to pre-chemotherapy samples. After neoadjuvant chemotherapy some residual samples may retain their molecular signature while others present significant changes in their gene expression, probably induced by the treatment. CTGF and DUSP1 overexpression in residual samples may be a reflection of resistance to further administration of AC regimen.

Maternal embryonic leucine zipper kinase transcript abundance correlates with malignancy grade in human astrocytomas

We have performed cDNA microarray analyses to identify gene expression differences between highly invasive glioblastoma multiforme (GBM) and typically benign pilocytic astrocytomas (PA). Despite the significant clinical and pathological differences between the 2 tumor types, only 63 genes were found to exhibit 2-fold or greater overexpression in GBM as compared to PA. Forty percent of these genes are related to the regulation of the cell cycle and mitosis. QT-PCR validation of 6 overexpressed genes: MELK, AUKB, ASPM, PRC1, IL13RA2 and KIAA0101 confirmed at least a 5-fold increase in the average expression levels in GBM. Maternal embryonic leucine zipper kinase (MELK) exhibited the most statistically significant difference. A more detailed investigation of MELK expression was undertaken to study its oncogenic relevance. In the examination of more than 100 tumors of the central nervous system, we found progressively higher expression of MELK with astrocytoma grade and a noteworthy uniformity of high level expression in GBM. Similar level of overexpression was also observed in medulloblastoma. We found neither gene promoter hypomethylation nor amplification to be a factor in MELK expression, but were able to demonstrate that MELK knockdown in malignant astrocytoma cell lines caused a reduction in proliferation and anchorage-independent growth in in vitro assays. Our results indicate that GBM and PA differ by the expression of surprisingly few genes. Among them, MELK correlated with malignancy grade in astrocytomas and represents a therapeutic target for the management of the most frequent brain tumors in adult and children. (C) 2007 Wiley-Liss, Inc.