Gene expression regulation and lineage evolution: the North and South tale of the hybrid polyploid Squalius alburnoides complex

The evolution of hybrid polyploid vertebrates, their viability and their perpetuation over evolutionary time have always been questions of great interest. However, little is known about the impact of hybridization and polyploidization on the regulatory networks that guarantee the appropriate quantitative and qualitative gene expression programme. The Squalius alburnoides complex of hybrid fish is an attractive system to address these questions, as it includes a wide variety of diploid and polyploid forms, and intricate systems of genetic exchange. Through the study of genome-specific allele expression of seven housekeeping and tissue-specific genes, we found that a gene copy silencing mechanism of dosage compensation exists throughout the distribution range of the complex. Here we show that the allele-specific patterns of silencing vary within the complex, according to the geographical origin and the type of genome involved in the hybridization process. In southern populations, triploids of S. alburnoides show an overall tendency for silencing the allele from the minority genome, while northern population polyploids exhibit preferential biallelic gene expression patterns, irrespective of genomic composition. The present findings further suggest that gene copy silencing and variable expression of specific allele combinations may be important processes in vertebrate polyploid evolution.

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.

Angiogenesis-related gene expression profile with independent prognostic value in advanced ovarian carcinoma.

BACKGROUND Ovarian carcinoma is the most important cause of gynecological cancer-related mortality in Western societies. Despite the improved median overall survival in patients receiving chemotherapy regimens such as paclitaxel and carboplatin combination, relapse still occurs in most advanced diseased patients. Increased angiogenesis is associated with rapid recurrence and decreased survival in ovarian cancer. This study was planned to identify an angiogenesis-related gene expression profile with prognostic value in advanced ovarian carcinoma patients. METHODOLOGY/PRINCIPAL FINDINGS RNAs were collected from formalin-fixed paraffin-embedded samples of 61 patients with III/IV FIGO stage ovarian cancer who underwent surgical cytoreduction and received a carboplatin plus paclitaxel regimen. Expression levels of 82 angiogenesis related genes were measured by quantitative real-time polymerase chain reaction using TaqMan low-density arrays. A 34-gene-profile which was able to predict the overall survival of ovarian carcinoma patients was identified. After a leave-one-out cross validation, the profile distinguished two groups of patients with different outcomes. Median overall survival and progression-free survival for the high risk group was 28.3 and 15.0 months, respectively, and was not reached by patients in the low risk group at the end of follow-up. Moreover, the profile maintained an independent prognostic value in the multivariate analysis. The hazard ratio for death was 2.3 (95% CI, 1.5 to 3.2; p<0.001). CONCLUSIONS/SIGNIFICANCE It is possible to generate a prognostic model for advanced ovarian carcinoma based on angiogenesis-related genes using formalin-fixed paraffin-embedded samples. The present results are consistent with the increasing weight of angiogenesis genes in the prognosis of ovarian carcinoma.

Gene expression profiling in breast cancer: understanding the molecular basis of histologic grade to improve prognosis.

BACKGROUND: Histologic grade in breast cancer provides clinically important prognostic information. However, 30%-60% of tumors are classified as histologic grade 2. This grade is associated with an intermediate risk of recurrence and is thus not informative for clinical decision making. We examined whether histologic grade was associated with gene expression profiles of breast cancers and whether such profiles could be used to improve histologic grading. METHODS: We analyzed microarray data from 189 invasive breast carcinomas and from three published gene expression datasets from breast carcinomas. We identified differentially expressed genes in a training set of 64 estrogen receptor (ER)-positive tumor samples by comparing expression profiles between histologic grade 3 tumors and histologic grade 1 tumors and used the expression of these genes to define the gene expression grade index. Data from 597 independent tumors were used to evaluate the association between relapse-free survival and the gene expression grade index in a Kaplan-Meier analysis. All statistical tests were two-sided. RESULTS: We identified 97 genes in our training set that were associated with histologic grade; most of these genes were involved in cell cycle regulation and proliferation. In validation datasets, the gene expression grade index was strongly associated with histologic grade 1 and 3 status; however, among histologic grade 2 tumors, the index spanned the values for histologic grade 1-3 tumors. Among patients with histologic grade 2 tumors, a high gene expression grade index was associated with a higher risk of recurrence than a low gene expression grade index (hazard ratio = 3.61, 95% confidence interval = 2.25 to 5.78; P < .001, log-rank test). CONCLUSIONS: Gene expression grade index appeared to reclassify patients with histologic grade 2 tumors into two groups with high versus low risks of recurrence. This approach may improve the accuracy of tumor grading and thus its prognostic value.

Gene expression profiling of isolated tumour cells from anaplastic large cell lymphomas: insights into its cellular origin, pathogenesis and relation to Hodgkin lymphoma.

Anaplastic large cell lymphoma (ALCL) is a main type of T-cell lymphomas and comprises three distinct entities: systemic anaplastic lymphoma kinase (ALK) positive, systemic ALK(-) and cutaneous ALK(-) ALCL (cALCL). Little is known about their pathogenesis and their cellular origin, and morphological and immunophenotypical overlap exists between ALK(-) ALCL and classical Hodgkin lymphoma (cHL). We conducted gene expression profiling of microdissected lymphoma cells of five ALK(+) and four ALK(-) systemic ALCL, seven cALCL and sixteen cHL, and of eight subsets of normal T and NK cells. The analysis supports a derivation of ALCL from activated T cells, but the lymphoma cells acquired a gene expression pattern hampering an assignment to a CD4(+), CD8(+) or CD30(+) T-cell origin. Indeed, ALCL display a down-modulation of many T-cell characteristic molecules. All ALCL types show significant expression of NFkappaB target genes and upregulation of genes involved in oncogenesis (e.g. EZH2). Surprisingly, few genes are differentially expressed between systemic and cALCL despite their different clinical behaviour, and between ALK(-) ALCL and cHL despite their different cellular origin. ALK(+) ALCL are characterized by expression of genes regulated by pathways constitutively activated by ALK. This study provides multiple novel insights into the molecular biology and pathogenesis of ALCL.

Test of four colon cancer risk-scores in formalin fixed paraffin embedded microarray gene expression data.

BACKGROUND: Prognosis prediction for resected primary colon cancer is based on the T-stage Node Metastasis (TNM) staging system. We investigated if four well-documented gene expression risk scores can improve patient stratification. METHODS: Microarray-based versions of risk-scores were applied to a large independent cohort of 688 stage II/III tumors from the PETACC-3 trial. Prognostic value for relapse-free survival (RFS), survival after relapse (SAR), and overall survival (OS) was assessed by regression analysis. To assess improvement over a reference, prognostic model was assessed with the area under curve (AUC) of receiver operating characteristic (ROC) curves. All statistical tests were two-sided, except the AUC increase. RESULTS: All four risk scores (RSs) showed a statistically significant association (single-test, P < .0167) with OS or RFS in univariate models, but with HRs below 1.38 per interquartile range. Three scores were predictors of shorter RFS, one of shorter SAR. Each RS could only marginally improve an RFS or OS model with the known factors T-stage, N-stage, and microsatellite instability (MSI) status (AUC gains < 0.025 units). The pairwise interscore discordance was never high (maximal Spearman correlation = 0.563) A combined score showed a trend to higher prognostic value and higher AUC increase for OS (HR = 1.74, 95% confidence interval [CI] = 1.44 to 2.10, P < .001, AUC from 0.6918 to 0.7321) and RFS (HR = 1.56, 95% CI = 1.33 to 1.84, P < .001, AUC from 0.6723 to 0.6945) than any single score. CONCLUSIONS: The four tested gene expression-based risk scores provide prognostic information but contribute only marginally to improving models based on established risk factors. A combination of the risk scores might provide more robust information. Predictors of RFS and SAR might need to be different.

ET-1 and NOS III gene expression regulation by plaque-free and plaque-prone hemodynamic conditions

Résumé: Les environnements hémodynamiques, favorisant ou protégeant contre la formation de la plaque, induisent tout deux une augmentation de la production d'anion superoxide dans les cellules endothéliales (ECs). Par ailleurs, une régulation différente de l'expression des gènes a été décrite dans les cellules exposées à ces différentes conditions. Dans le but d'investiguer le rôle de l'augmentation du stress oxydatif dans l'expression des gènes régulée par le flux, nous avons d'abord exposé les EC à un flux unidirectionnel, non pulsé. Dans ces conditions, l'état oxydatif des cellules endothéliales est augmenté de façon transitoire. L'expression du gène de l'endothéline 1 (ET-1) est aussi induite de façon transitoire par un tel flux, alors que l'expression du gène de la nitiric oxyde synthase endothéliale (NOS III) est stimulé de façon durable. Au contraire, un flux unidirectionnel pulsé, qui induit une augmentation durable de la production d'anion superoxide, augmente aussi de façon durable l'expression des gènes de ET-1 comme de NOS III. Un flux oscillatoire (favorisant la plaque), qui lui aussi ,a des effets à long terme sur la production d'anion superoxide, a uniquement augmenté l'expression de ET-1. De plus, l'utilisation d'un antioxydant, a seulement partiellement inhibé la stimulation de l'expression du gène NOS III par le flux unidirectionnel pulsé, alors qu'il a complètement abrogé la stimulation de l'expression du gène ET-1 par le flux unidirectionnel pulsé et oscillatoire. Ceci suggère que les forces mécaniques régulent l'expression des gènes dans les EC par un double mécanisme dépendant et indépendant du stress oxidatif des cellules. Par ailleurs, ces résultats supportent ultérieurement l'hypothèse que la balance entre la réponse oxidative et anti-oxidante dans les cellules endothéliales exposées à un environnement hémodynamique est une des clés de la prédisposition à un dysfonctionnement endothélial observé dans des régions exposées à des flux perturbés. Abstract: Both plaque-free and plaque-prone hemodynamic environments induce an increase in the oxidative state of endothelial cells (ECs), whereas differential gene expression regulation was described in cells exposed to these conditions. In order to investigate the role of the increased oxidative state in flow-regulation of gene expression, we first exposed EC to non-pulsed unidirectional shear stress. These conditions only slightly increases ECs oxidative state and endothelin-1 (ET-1) mRNA expression, whereas endothelial nitric oxide synthase (NOS III) mRNA level were significantly up-regulated. On the contrary, both ET-1 and NOS III gene expression were significantly induced in EC exposed to pulsed-unidirectional flow (plaque-free). Only ET-1 gene expression was up-regulated by oscillatory flow (plaque-prone). Moreover, use of an antioxidant only partially inhibited NOS III gene up-regulation by unidirectional flow, whereas it completely abrogated ET-1 gene up-regulation by unidirectional and oscillatory flows. Thus suggesting that mechanical forces regulate gene expression in ECs both via oxidative stress-dependent and -independent mechanisms.

eVOC: a controlled vocabulary for unifying gene expression data.

Expression data contribute significantly to the biological value of the sequenced human genome, providing extensive information about gene structure and the pattern of gene expression. ESTs, together with SAGE libraries and microarray experiment information, provide a broad and rich view of the transcriptome. However, it is difficult to perform large-scale expression mining of the data generated by these diverse experimental approaches. Not only is the data stored in disparate locations, but there is frequent ambiguity in the meaning of terms used to describe the source of the material used in the experiment. Untangling semantic differences between the data provided by different resources is therefore largely reliant on the domain knowledge of a human expert. We present here eVOC, a system which associates labelled target cDNAs for microarray experiments, or cDNA libraries and their associated transcripts with controlled terms in a set of hierarchical vocabularies. eVOC consists of four orthogonal controlled vocabularies suitable for describing the domains of human gene expression data including Anatomical System, Cell Type, Pathology and Developmental Stage. We have curated and annotated 7016 cDNA libraries represented in dbEST, as well as 104 SAGE libraries,with expression information,and provide this as an integrated, public resource that allows the linking of transcripts and libraries with expression terms. Both the vocabularies and the vocabulary-annotated libraries can be retrieved from http://www.sanbi.ac.za/evoc/. Several groups are involved in developing this resource with the aim of unifying transcript expression information.

Meta-analysis of gene expression profiles in breast cancer: toward a unified understanding of breast cancer subtyping and prognosis signatures.

INTRODUCTION: Breast cancer subtyping and prognosis have been studied extensively by gene expression profiling, resulting in disparate signatures with little overlap in their constituent genes. Although a previous study demonstrated a prognostic concordance among gene expression signatures, it was limited to only one dataset and did not fully elucidate how the different genes were related to one another nor did it examine the contribution of well-known biological processes of breast cancer tumorigenesis to their prognostic performance. METHOD: To address the above issues and to further validate these initial findings, we performed the largest meta-analysis of publicly available breast cancer gene expression and clinical data, which are comprised of 2,833 breast tumors. Gene coexpression modules of three key biological processes in breast cancer (namely, proliferation, estrogen receptor [ER], and HER2 signaling) were used to dissect the role of constituent genes of nine prognostic signatures. RESULTS: Using a meta-analytical approach, we consolidated the signatures associated with ER signaling, ERBB2 amplification, and proliferation. Previously published expression-based nomenclature of breast cancer 'intrinsic' subtypes can be mapped to the three modules, namely, the ER-/HER2- (basal-like), the HER2+ (HER2-like), and the low- and high-proliferation ER+/HER2- subtypes (luminal A and B). We showed that all nine prognostic signatures exhibited a similar prognostic performance in the entire dataset. Their prognostic abilities are due mostly to the detection of proliferation activity. Although ER- status (basal-like) and ERBB2+ expression status correspond to bad outcome, they seem to act through elevated expression of proliferation genes and thus contain only indirect information about prognosis. Clinical variables measuring the extent of tumor progression, such as tumor size and nodal status, still add independent prognostic information to proliferation genes. CONCLUSION: This meta-analysis unifies various results of previous gene expression studies in breast cancer. It reveals connections between traditional prognostic factors, expression-based subtyping, and prognostic signatures, highlighting the important role of proliferation in breast cancer prognosis.

Different Effects of BORIS/CTCFL on Stemness Gene Expression, Sphere Formation and Cell Survival in Epithelial Cancer Stem Cells.

Cancer stem cells are cancer cells characterized by stem cell properties and represent a small population of tumor cells that drives tumor development, progression, metastasis and drug resistance. To date, the molecular mechanisms that generate and regulate cancer stem cells are not well defined. BORIS (Brother of Regulator of Imprinted Sites) or CTCFL (CTCF-like) is a DNA-binding protein that is expressed in normal tissues only in germ cells and is re-activated in tumors. Recent evidences have highlighted the correlation of BORIS/CTCFL expression with poor overall survival of different cancer patients. We have previously shown an association of BORIS-expressing cells with stemness gene expression in embryonic cancer cells. Here, we studied the role of BORIS in epithelial tumor cells. Using BORIS-molecular beacon that was already validated, we were able to show the presence of BORIS mRNA in cancer stem cell-enriched populations (side population and spheres) of cervical, colon and breast tumor cells. BORIS silencing studies showed a decrease of sphere formation capacity in breast and colon tumor cells. Importantly, BORIS-silencing led to down-regulation of hTERT, stem cell (NANOG, OCT4, SOX2 and BMI1) and cancer stem cell markers (ABCG2, CD44 and ALDH1) genes. Conversely, BORIS-induction led to up-regulation of the same genes. These phenotypes were observed in cervical, colon and invasive breast tumor cells. However, a completely different behavior was observed in the non-invasive breast tumor cells (MCF7). Indeed, these cells acquired an epithelial mesenchymal transition phenotype after BORIS silencing. Our results demonstrate that BORIS is associated with cancer stem cell-enriched populations of several epithelial tumor cells and the different phenotypes depend on the origin of tumor cells.

Tamoxifen inhibits transforming growth factor-α gene expression in human breast carcinoma samples treated with triiodothyronine

Objectives: To examine the effects of triiodothyronine (T3), 17β-estradiol (E2), and tamoxifen (TAM) on transforming growth factor (TGF)-α gene expression in primary breast cancer cell cultures and interactions between the different treatments. Methods and results: Patients included in the study (no.=12) had been newly diagnosed with breast cancer. Fresh human breast carcinoma tissue was cut into 0.3-mm slices. These slices were placed in six 35-mm dishes on 2-ml organ culture medium. Dishes received the following treatments: dish 1: ethanol; dish 2: T3; dish 3: T3+TAM; dish 4: TAM; dish 5: E2; dish 6: E2+TAM. TGF-α mRNA content was normalized to glyceraldehyde-3-phosphate dehydrogenase mRNA levels. All tissues included in this study were positive for estrogen receptor (ER) and thyroid hormone receptor expression. Treatment with T3 for 48 h significantly increased TGF-α mRNA levels compared to controls (15-fold), and concomitant treatment with TAM reduced expression to 3.4-fold compared to controls. When only TAM was added to the culture medium, TGF-α mRNA expression increased 5.3-fold, significantly higher than with all other treatment modalities. Conclusion: We demonstrate that TGF-α mRNA expression is more efficiently upregulated by T3 than E2. Concomitant treatment with TAM had a mitigating effect on the T3 effect, while E2 induced TGF-α upregulation. Our findings show some similarities between primary culture and breast cancer cell lines, but also some important differences: a) induction of TGF-α, a mitogenic protein, by TAM; b) a differential effect of TAM that may depend on relative expression of ER α and β; and c) supraphysiological doses of T3 may induce mitogenic signals in breast cancer tissue under conditions of low circulating E2. ©2008, Editrice Kurtis.

Homeobox gene expression profile indicates HOXA5 as a candidate prognostic marker in oral squamous cell carcinoma

The search for molecular markers to improve diagnosis, individualize treatment and predict behavior of tumors has been the focus of several studies. This study aimed to analyze homeobox gene expression profile in oral squamous cell carcinoma (OSCC) as well as to investigate whether some of these genes are relevant molecular markers of prognosis and/or tumor aggressiveness. Homeobox gene expression levels were assessed by microarrays and qRT-PCR in OSCC tissues and adjacent non-cancerous matched tissues (margin), as well as in OSCC cell lines. Analysis of microarray data revealed the expression of 147 homeobox genes, including one set of six at least 2-fold up-regulated, and another set of 34 at least 2-fold down-regulated homeobox genes in OSCC. After qRT-PCR assays, the three most up-regulated homeobox genes (HOXA5, HOXD10 and HOXD11) revealed higher and statistically significant expression levels in OSCC samples when compared to margins. Patients presenting lower expression of HOXA5 had poorer prognosis compared to those with higher expression (P=0.03). Additionally, the status of HOXA5, HOXD10 and HOXD11 expression levels in OSCC cell lines also showed a significant up-regulation when compared to normal oral keratinocytes. Results confirm the presence of three significantly upregulated (>4-fold) homeobox genes (HOXA5, HOXD10 and HOXD11) in OSCC that may play a significant role in the pathogenesis of these tumors. Moreover, since lower levels of HOXA5 predict poor prognosis, this gene may be a novel candidate for development of therapeutic strategies in OSCC.

The TRANSFAC system on gene expression regulation

The TRANSFAC database on transcription factors and their DNA-binding sites and profiles (http://www.gene-regulation.de/) has been quantitatively extended and supplemented by a number of modules. These modules give information about pathologically relevant mutations in regulatory regions and transcription factor genes (PathoDB), scaffold/matrix attached regions (S/MARt DB), signal transduction (TRANSPATH) and gene expression sources (CYTOMER). Altogether, these distinct database modules constitute the TRANSFAC system. They are accompanied by a number of program routines for identifying potential transcription factor binding sites or for localizing individual components in the regulatory network of a cell.