Predicting cancer involvement of genes from heterogeneous data

Background: Systematic approaches for identifying proteins involved in different types of cancer are needed. Experimental techniques such as microarrays are being used to characterize cancer, but validating their results can be a laborious task. Computational approaches are used to prioritize between genes putatively involved in cancer, usually based on further analyzing experimental data. Results: We implemented a systematic method using the PIANA software that predicts cancer involvement of genes by integrating heterogeneous datasets. Specifically, we produced lists of genes likely to be involved in cancer by relying on: (i) protein-protein interactions; (ii) differential expression data; and (iii) structural and functional properties of cancer genes. The integrative approach that combines multiple sources of data obtained positive predictive values ranging from 23% (on a list of 811 genes) to 73% (on a list of 22 genes), outperforming the use of any of the data sources alone. We analyze a list of 20 cancer gene predictions, finding that most of them have been recently linked to cancer in literature. Conclusion: Our approach to identifying and prioritizing candidate cancer genes can be used to produce lists of genes likely to be involved in cancer. Our results suggest that differential expression studies yielding high numbers of candidate cancer genes can be filtered using protein interaction networks.

Networking of differentially expressed genes in human cancer cells resistant to methotrexate

BACKGROUND: The need for an integrated view of data obtained from high-throughput technologies gave rise to network analyses. These are especially useful to rationalize how external perturbations propagate through the expression of genes. To address this issue in the case of drug resistance, we constructed biological association networks of genes differentially expressed in cell lines resistant to methotrexate (MTX). METHODS: Seven cell lines representative of different types of cancer, including colon cancer (HT29 and Caco2), breast cancer (MCF-7 and MDA-MB-468), pancreatic cancer (MIA PaCa-2), erythroblastic leukemia (K562) and osteosarcoma (Saos-2), were used. The differential expression pattern between sensitive and MTX-resistant cells was determined by whole human genome microarrays and analyzed with the GeneSpring GX software package. Genes deregulated in common between the different cancer cell lines served to generate biological association networks using the Pathway Architect software. RESULTS: Dikkopf homolog-1 (DKK1) is a highly interconnected node in the network generated with genes in common between the two colon cancer cell lines, and functional validations of this target using small interfering RNAs (siRNAs) showed a chemosensitization toward MTX. Members of the UDP-glucuronosyltransferase 1A (UGT1A) family formed a network of genes differentially expressed in the two breast cancer cell lines. siRNA treatment against UGT1A also showed an increase in MTX sensitivity. Eukaryotic translation elongation factor 1 alpha 1 (EEF1A1) was overexpressed among the pancreatic cancer, leukemia and osteosarcoma cell lines, and siRNA treatment against EEF1A1 produced a chemosensitization toward MTX. CONCLUSIONS: Biological association networks identified DKK1, UGT1As and EEF1A1 as important gene nodes in MTX-resistance. Treatments using siRNA technology against these three genes showed chemosensitization toward MTX.

Networking of differentially expressed genes in human cancer cells resistant to methotrexate

BACKGROUND: The need for an integrated view of data obtained from high-throughput technologies gave rise to network analyses. These are especially useful to rationalize how external perturbations propagate through the expression of genes. To address this issue in the case of drug resistance, we constructed biological association networks of genes differentially expressed in cell lines resistant to methotrexate (MTX). METHODS: Seven cell lines representative of different types of cancer, including colon cancer (HT29 and Caco2), breast cancer (MCF-7 and MDA-MB-468), pancreatic cancer (MIA PaCa-2), erythroblastic leukemia (K562) and osteosarcoma (Saos-2), were used. The differential expression pattern between sensitive and MTX-resistant cells was determined by whole human genome microarrays and analyzed with the GeneSpring GX software package. Genes deregulated in common between the different cancer cell lines served to generate biological association networks using the Pathway Architect software. RESULTS: Dikkopf homolog-1 (DKK1) is a highly interconnected node in the network generated with genes in common between the two colon cancer cell lines, and functional validations of this target using small interfering RNAs (siRNAs) showed a chemosensitization toward MTX. Members of the UDP-glucuronosyltransferase 1A (UGT1A) family formed a network of genes differentially expressed in the two breast cancer cell lines. siRNA treatment against UGT1A also showed an increase in MTX sensitivity. Eukaryotic translation elongation factor 1 alpha 1 (EEF1A1) was overexpressed among the pancreatic cancer, leukemia and osteosarcoma cell lines, and siRNA treatment against EEF1A1 produced a chemosensitization toward MTX. CONCLUSIONS: Biological association networks identified DKK1, UGT1As and EEF1A1 as important gene nodes in MTX-resistance. Treatments using siRNA technology against these three genes showed chemosensitization toward MTX.

Gene expression following induction of regeneration in Drosophila wing imaginal discs

Background: Regeneration is the ability of an organism to rebuild a body part that has been damaged or amputated, and can be studied at the molecular level using model organisms. Drosophila imaginal discs, which are the larval primordia of adult cuticular structures, are capable of undergoing regenerative growth after transplantation and in vivo culture into the adult abdomen. Results: Using expression profile analyses, we studied the regenerative behaviour of wing discs at 0, 24 and 72 hours after fragmentation and implantation into adult females. Based on expression level, we generated a catalogue of genes with putative role in wing disc regeneration, identifying four classes: 1) genes with differential expression within the first 24 hours; 2) genes with differential expression between 24 and 72 hours; 3) genes that changed significantly in expression levels between the two time periods; 4) genes with a sustained increase or decrease in their expression levels throughout regeneration. Among these genes, we identified members of the JNK and Notch signalling pathways and chromatin regulators. Through computational analysis, we recognized putative binding sites for transcription factors downstream of these pathways that are conserved in multiple Drosophilids, indicating a potential relationship between members of the different gene classes. Experimental data from genetic mutants provide evidence of a requirement of selected genes in wing disc regeneration. Conclusions: We have been able to distinguish various classes of genes involved in early and late steps of the regeneration process. Our data suggests the integration of signalling pathways in the promoters of regulated genes.

Gene expression following induction of regeneration in Drosophila wing imaginal discs

Background: Regeneration is the ability of an organism to rebuild a body part that has been damaged or amputated, and can be studied at the molecular level using model organisms. Drosophila imaginal discs, which are the larval primordia of adult cuticular structures, are capable of undergoing regenerative growth after transplantation and in vivo culture into the adult abdomen. Results: Using expression profile analyses, we studied the regenerative behaviour of wing discs at 0, 24 and 72 hours after fragmentation and implantation into adult females. Based on expression level, we generated a catalogue of genes with putative role in wing disc regeneration, identifying four classes: 1) genes with differential expression within the first 24 hours; 2) genes with differential expression between 24 and 72 hours; 3) genes that changed significantly in expression levels between the two time periods; 4) genes with a sustained increase or decrease in their expression levels throughout regeneration. Among these genes, we identified members of the JNK and Notch signalling pathways and chromatin regulators. Through computational analysis, we recognized putative binding sites for transcription factors downstream of these pathways that are conserved in multiple Drosophilids, indicating a potential relationship between members of the different gene classes. Experimental data from genetic mutants provide evidence of a requirement of selected genes in wing disc regeneration. Conclusions: We have been able to distinguish various classes of genes involved in early and late steps of the regeneration process. Our data suggests the integration of signalling pathways in the promoters of regulated genes.

Genetic and genomic analysis modeling of germline c-MYC overexpression and cancer susceptibility

Background: Germline genetic variation is associated with the differential expression of many human genes. The phenotypic effects of this type of variation may be important when considering susceptibility to common genetic diseases. Three regions at 8q24 have recently been identified to independently confer risk of prostate cancer. Variation at 8q24 has also recently been associated with risk of breast and colorectal cancer. However, none of the risk variants map at or relatively close to known genes, with c-MYC mapping a few hundred kilobases distally. Results: This study identifies cis-regulators of germline c-MYC expression in immortalized lymphocytes of HapMap individuals. Quantitative analysis of c-MYC expression in normal prostate tissues suggests an association between overexpression and variants in Region 1 of prostate cancer risk. Somatic c-MYC overexpression correlates with prostate cancer progression and more aggressive tumor forms, which was also a pathological variable associated with Region 1. Expression profiling analysis and modeling of transcriptional regulatory networks predicts a functional association between MYC and the prostate tumor suppressor KLF6. Analysis of MYC/Myc-driven cell transformation and tumorigenesis substantiates a model in which MYC overexpression promotes transformation by down-regulating KLF6. In this model, a feedback loop through E-cadherin down-regulation causes further transactivation of c-MYC.Conclusion: This study proposes that variation at putative 8q24 cis-regulator(s) of transcription can significantly alter germline c-MYC expression levels and, thus, contribute to prostate cancer susceptibility by down-regulating the prostate tumor suppressor KLF6 gene.

Transcriptome analyses of primitively eusocial wasps reveal novel insights into the evolution of sociality and the origin of alternative phenotypes

BACKGROUND: Understanding how alternative phenotypes arise from the same genome is a major challenge in modern biology. Eusociality in insects requires the evolution of two alternative phenotypes - workers, who sacrifice personal reproduction, and queens, who realize that reproduction. Extensive work on honeybees and ants has revealed the molecular basis of derived queen and worker phenotypes in highly eusocial lineages, but we lack equivalent deep-level analyses of wasps and of primitively eusocial species, the latter of which can reveal how phenotypic decoupling first occurs in the early stages of eusocial evolution. RESULTS: We sequenced 20 Gbp of transcriptomes derived from brains of different behavioral castes of the primitively eusocial tropical paper wasp Polistes canadensis. Surprisingly, 75% of the 2,442 genes differentially expressed between phenotypes were novel, having no significant homology with described sequences. Moreover, 90% of these novel genes were significantly upregulated in workers relative to queens. Differential expression of novel genes in the early stages of sociality may be important in facilitating the evolution of worker behavioral complexity in eusocial evolution. We also found surprisingly low correlation in the identity and direction of expression of differentially expressed genes across similar phenotypes in different social lineages, supporting the idea that social evolution in different lineages requires substantial de novo rewiring of molecular pathways. CONCLUSIONS: These genomic resources for aculeate wasps and first transcriptome-wide insights into the origin of castes bring us closer to a more general understanding of eusocial evolution and how phenotypic diversity arises from the same genome.

Differential regulation of horizontally acquired and core genome genes by the bacterial modulator H-NS

Horizontal acquisition of DNA by bacteria dramatically increases genetic diversity and hence successful bacterial colonization of several niches, including the human host. A relevant issue is how this newly acquired DNA interacts and integrates in the regulatory networks of the bacterial cell. The global modulator H-NS targets both core genome and HGT genes and silences gene expression in response to external stimuli such as osmolarity and temperature. Here we provide evidence that H-NS discriminates and differentially modulates core and HGT DNA. As an example of this, plasmid R27-encoded H-NS protein has evolved to selectively silence HGT genes and does not interfere with core genome regulation. In turn, differential regulation of both gene lineages by resident chromosomal H-NS requires a helper protein: the Hha protein. Tight silencing of HGT DNA is accomplished by H-NS-Hha complexes. In contrast, core genes are modulated by H-NS homoligomers. Remarkably, the presence of Hha-like proteins is restricted to the Enterobacteriaceae. In addition, conjugative plasmids encoding H-NS variants have hitherto been isolated only from members of the family. Thus, the H-NS system in enteric bacteria presents unique evolutionary features. The capacity to selectively discriminate between core and HGT DNA may help to maintain horizontally transmitted DNA in silent form and may give these bacteria a competitive advantage in adapting to new environments, including host colonization.

Differential regulation of horizontally acquired and core genome genes by the bacterial modulator H-NS

Horizontal acquisition of DNA by bacteria dramatically increases genetic diversity and hence successful bacterial colonization of several niches, including the human host. A relevant issue is how this newly acquired DNA interacts and integrates in the regulatory networks of the bacterial cell. The global modulator H-NS targets both core genome and HGT genes and silences gene expression in response to external stimuli such as osmolarity and temperature. Here we provide evidence that H-NS discriminates and differentially modulates core and HGT DNA. As an example of this, plasmid R27-encoded H-NS protein has evolved to selectively silence HGT genes and does not interfere with core genome regulation. In turn, differential regulation of both gene lineages by resident chromosomal H-NS requires a helper protein: the Hha protein. Tight silencing of HGT DNA is accomplished by H-NS-Hha complexes. In contrast, core genes are modulated by H-NS homoligomers. Remarkably, the presence of Hha-like proteins is restricted to the Enterobacteriaceae. In addition, conjugative plasmids encoding H-NS variants have hitherto been isolated only from members of the family. Thus, the H-NS system in enteric bacteria presents unique evolutionary features. The capacity to selectively discriminate between core and HGT DNA may help to maintain horizontally transmitted DNA in silent form and may give these bacteria a competitive advantage in adapting to new environments, including host colonization.

Teleost fish larvae adapt to dietary arachidonic acid supply through modulation of the expression of lipid metabolism and stress response genes

Dietary fatty acid supply can affect stress response in fish during early development. Although knowledge on the mechanisms involved in fatty acid regulation of stress tolerance is scarce, it has often been hypothesised that eicosanoid profiles can influence cortisol production. Genomic cortisol actions are mediated by cytosolic receptors which may respond to cellular fatty acid signalling. An experiment was designed to test the effects of feeding gilthead sea-bream larvae with four microdiets, containing graded arachidonic acid (ARA) levels (0·4, 0·8, 1·5 and 3·0 %), on the expression of genes involved in stress response (steroidogenic acute regulatory protein, glucocorticoid receptor and phosphoenolpyruvate carboxykinase), lipid and, particularly, eicosanoid metabolism (hormone-sensitive lipase, PPARα, phospholipase A2, cyclo-oxygenase-2 and 5-lipoxygenase), as determined by real-time quantitative PCR. Fish fatty acid phenotypes reflected dietary fatty acid profiles. Growth performance, survival after acute stress and similar whole-body basal cortisol levels suggested that sea-bream larvae could tolerate a wide range of dietary ARA levels. Transcription of all genes analysed was significantly reduced at dietary ARA levels above 0·4 %. Nonetheless, despite practical suppression of phospholipase A2 transcription, higher leukotriene B4 levels were detected in larvae fed 3·0 % ARA, whereas a similar trend was observed regarding PGE2 production. The present study demonstrates that adaptation to a wide range of dietary ARA levels in gilthead sea-bream larvae involves the modulation of the expression of genes related to eicosanoid synthesis, lipid metabolism and stress response. The roles of ARA, other polyunsaturates and eicosanoids as signals in this process are discussed.

Comparative effects of oleoyl-estrone and a specific b3-adrenergic agonist (CL316, 243) on the expression of genes involved in energy metabolism of rat white adipose tissue

Background: The combination of oleoyl-estrone (OE) and a selective b3-adrenergic agonist (B3A; CL316,243) treatment in rats results in a profound and rapid wasting of body reserves (lipid). Methods: In the present study we investigated the effect of OE (oral gavage) and/or B3A (subcutaneous constant infusion) administration for 10 days to overweight male rats, compared with controls, on three distinct white adipose tissue (WAT) sites: subcutaneous inguinal, retroperitoneal and epididymal. Tissue weight, DNA (and, from these values cellularity), cAMP content and the expression of several key energy handling metabolism and control genes were analyzed and computed in relation to the whole site mass. Results: Both OE and B3A significantly decreased WAT mass, with no loss of DNA (cell numbers). OE decreased and B3A increased cAMP. Gene expression patterns were markedly different for OE and B3A. OE tended to decrease expression of most genes studied, with no changes (versus controls) of lipolytic but decrease of lipogenic enzyme genes. The effects of B3A were widely different, with a generalized increase in the expression of most genes, including the adrenergic receptors, and, especially the uncoupling protein UCP1. Discussion: OE and B3A, elicit widely different responses in WAT gene expression, end producing similar effects, such as shrinking of WAT, loss of fat, maintenance of cell numbers. OE acted essentially on the balance of lipolysislipogenesis and the blocking of the uptake of substrates; its decrease of synthesis favouring lipolysis. B3A induced a shotgun increase in the expression of most regulatory systems in the adipocyte, an effect that in the end favoured again the loss of lipid; this barely selective increase probably produces inefficiency, which coupled with the increase in UCP1 expression may help WAT to waste energy through thermogenesis. Conclusions: There were considerable differences in the responses of the three WAT sites. OE in general lowered gene expression and stealthily induced a substrate imbalance. B3A increasing the expression of most genes enhanced energy waste through inefficiency rather than through specific pathway activation. There was not a synergistic effect between OE and B3A in WAT, but their combined action increased WAT energy waste.

Enhanced Botrytis cinerea resistance of Arabidopsis plants grown in compost may be explained by increased expression of defense-related genes, as revealed by microarray analysis

Composts are the products obtained after the aerobic degradation of different types of organic matter waste and can be used as substrates or substrate/soil amendments for plant cultivation. There is a small but increasing number of reports that suggest that foliar diseases may be reduced when using compost, rather than standard substrates, as growing medium. The purpose of this study was to examine the gene expression alteration produced by the compost to gain knowledge of the mechanisms involved in compost-induced systemic resistance. A compost from olive marc and olive tree leaves was able to induce resistance against Botrytis cinerea in Arabidopsis, unlike the standard substrate, perlite. Microarray analyses revealed that 178 genes were differently expressed, with a fold change cut-off of 1, of which 155 were up-regulated and 23 were down-regulated in compost-grown, as against perlite-grown plants. A functional enrichment study of up-regulated genes revealed that 38 Gene Ontology terms were significantly enriched. Response to stress, biotic stimulus, other organism, bacterium, fungus, chemical and abiotic stimulus, SA and ABA stimulus, oxidative stress, water, temperature and cold were significantly enriched, as were immune and defense responses, systemic acquired resistance, secondary metabolic process and oxireductase activity. Interestingly, PR1 expression, which was equally enhanced by growing the plants in compost and by B. cinerea inoculation, was further boosted in compost-grown pathogen-inoculated plants. Compost triggered a plant response that shares similarities with both systemic acquired resistance and ABA-dependent/independent abiotic stress responses.

Differential effects of polyphenols and alcohol of red wine on the expression of adhesion molecules and inflammatory cytokines related to atherosclerosis: a randomized clinical trial

Background: Few clinical studies have focused on the alcoholindependent cardiovascular effects of the phenolic compounds of red wine (RW). Objective: We aimed to evaluate the effects of ethanol and phenolic compounds of RW on the expression of inflammatory biomarkers related to atherosclerosis in subjects at high risk of cardiovascular disease. Design: Sixty-seven high-risk, male volunteers were included in a randomized, crossover consumption trial. After a washout period, all subjects received RW (30 g alcohol/d), the equivalent amount of dealcoholized red wine (DRW), or gin (30 g alcohol/d) for 4 wk. Before and after each intervention period, 7 cellular and 18 serum inflammatory biomarkers were evaluated. Results: Alcohol increased IL-10 and decreased macrophage-derived chemokine concentrations, whereas the phenolic compounds of RW decreased serum concentrations of intercellular adhesion molecule- 1, E-selectin, and IL-6 and inhibited the expression of lymphocyte function-associated antigen 1 in T lymphocytes and macrophage-1 receptor, Sialil-Lewis X, and C-C chemokine receptor type 2 expression in monocytes. Both ethanol and phenolic compounds of RW downregulated serum concentrations of CD40 antigen, CD40 ligand, IL-16, monocyte chemotactic protein-1, and vascular cell adhesion molecule-1. Conclusion: The results suggest that the phenolic content of RW may modulate leukocyte adhesion molecules, whereas both ethanol and polyphenols of RW may modulate soluble inflammatory mediators in high-risk patients. The trial was registered in the International Standard Randomized Controlled Trial Number Register at http://www. isrctn.org/ as ISRCTN88720134

Differential effects of polyphenols and alcohol of red wine on the expression of adhesion molecules and inflammatory cytokines related to atherosclerosis: a randomized clinical trial

Background: Few clinical studies have focused on the alcoholindependent cardiovascular effects of the phenolic compounds of red wine (RW). Objective: We aimed to evaluate the effects of ethanol and phenolic compounds of RW on the expression of inflammatory biomarkers related to atherosclerosis in subjects at high risk of cardiovascular disease. Design: Sixty-seven high-risk, male volunteers were included in a randomized, crossover consumption trial. After a washout period, all subjects received RW (30 g alcohol/d), the equivalent amount of dealcoholized red wine (DRW), or gin (30 g alcohol/d) for 4 wk. Before and after each intervention period, 7 cellular and 18 serum inflammatory biomarkers were evaluated. Results: Alcohol increased IL-10 and decreased macrophage-derived chemokine concentrations, whereas the phenolic compounds of RW decreased serum concentrations of intercellular adhesion molecule- 1, E-selectin, and IL-6 and inhibited the expression of lymphocyte function-associated antigen 1 in T lymphocytes and macrophage-1 receptor, Sialil-Lewis X, and C-C chemokine receptor type 2 expression in monocytes. Both ethanol and phenolic compounds of RW downregulated serum concentrations of CD40 antigen, CD40 ligand, IL-16, monocyte chemotactic protein-1, and vascular cell adhesion molecule-1. Conclusion: The results suggest that the phenolic content of RW may modulate leukocyte adhesion molecules, whereas both ethanol and polyphenols of RW may modulate soluble inflammatory mediators in high-risk patients. The trial was registered in the International Standard Randomized Controlled Trial Number Register at http://www. isrctn.org/ as ISRCTN88720134

Cell fusion reprogramming leads to a specific hepatic expression pattern during mouse bone marrow derived hepatocyte formation In Vivo

The fusion of bone marrow (BM) hematopoietic cells with hepatocytes to generate BM derived hepatocytes (BMDH) is a natural process, which is enhanced in damaged tissues. However, the reprogramming needed to generate BMDH and the identity of the resultant cells is essentially unknown. In a mouse model of chronic liver damage, here we identify a modification in the chromatin structure of the hematopoietic nucleus during BMDH formation, accompanied by the loss of the key hematopoietic transcription factor PU.1/Sfpi1 (SFFV proviral integration 1) and gain of the key hepatic transcriptional regulator HNF-1A homeobox A (HNF-1A/Hnf1a). Through genome-wide expression analysis of laser captured BMDH, a differential gene expression pattern was detected and the chromatin changes observed were confirmed at the level of chromatin regulator genes. Similarly, Tranforming Growth Factor-β1 (TGF-β1) and neurotransmitter (e.g. Prostaglandin E Receptor 4 [Ptger4]) pathway genes were over-expressed. In summary, in vivo BMDH generation is a process in which the hematopoietic cell nucleus changes its identity and acquires hepatic features. These BMDHs have their own cell identity characterized by an expression pattern different from hematopoietic cells or hepatocytes. The role of these BMDHs in the liver requires further investigation.