“Unmixing” Tissue Gene Expression Signatures from Tumor Biopsies

Emergent molecular measurement methods, such as DNA microarray, qRTPCR, andmany others, offer tremendous promise for the personalized treatment of cancer. Thesetechnologies measure the amount of specific proteins, RNA, DNA or other moleculartargets from tumor specimens with the goal of “fingerprinting” individual cancers. Tumorspecimens are heterogeneous; an individual specimen typically contains unknownamounts of multiple tissues types. Thus, the measured molecular concentrations resultfrom an unknown mixture of tissue types, and must be normalized to account for thecomposition of the mixture.For example, a breast tumor biopsy may contain normal, dysplastic and cancerousepithelial cells, as well as stromal components (fatty and connective tissue) and bloodand lymphatic vessels. Our diagnostic interest focuses solely on the dysplastic andcancerous epithelial cells. The remaining tissue components serve to “contaminate”the signal of interest. The proportion of each of the tissue components changes asa function of patient characteristics (e.g., age), and varies spatially across the tumorregion. Because each of the tissue components produces a different molecular signature,and the amount of each tissue type is specimen dependent, we must estimate the tissuecomposition of the specimen, and adjust the molecular signal for this composition.Using the idea of a chemical mass balance, we consider the total measured concentrationsto be a weighted sum of the individual tissue signatures, where weightsare determined by the relative amounts of the different tissue types. We develop acompositional source apportionment model to estimate the relative amounts of tissuecomponents in a tumor specimen. We then use these estimates to infer the tissuespecificconcentrations of key molecular targets for sub-typing individual tumors. Weanticipate these specific measurements will greatly improve our ability to discriminatebetween different classes of tumors, and allow more precise matching of each patient tothe appropriate treatment

Whole genome analysis of p38 SAPK-mediated gene expression upon stress

Background: Cells have the ability to respond and adapt to environmental changes through activation of stress-activated protein kinases (SAPKs). Although p38 SAPK signalling is known to participate in the regulation of gene expression little is known on the molecular mechanisms used by this SAPK to regulate stress-responsive genes and the overall set of genes regulated by p38 in response to different stimuli.Results: Here, we report a whole genome expression analyses on mouse embryonic fibroblasts (MEFs) treated with three different p38 SAPK activating-stimuli, namely osmostress, the cytokine TNFα and the protein synthesis inhibitor anisomycin. We have found that the activation kinetics of p38α SAPK in response to these insults is different and also leads to a complex gene pattern response specific for a given stress with a restricted set of overlapping genes. In addition, we have analysed the contribution of p38α the major p38 family member present in MEFs, to the overall stress-induced transcriptional response by using both a chemical inhibitor (SB203580) and p38α deficient (p38α-/-) MEFs. We show here that p38 SAPK dependency ranged between 60% and 88% depending on the treatments and that there is a very good overlap between the inhibitor treatment and the ko cells. Furthermore, we have found that the dependency of SAPK varies depending on the time the cells are subjected to osmostress. Conclusions: Our genome-wide transcriptional analyses shows a selective response to specific stimuli and a restricted common response of up to 20% of the stress up-regulated early genes that involves an important set of transcription factors, which might be critical for either cell adaptation or preparation for continuous extra-cellular changes. Interestingly, up to 85% of the up-regulated genes are under the transcriptional control of p38 SAPK. Thus, activation of p38 SAPK is critical to elicit the early gene expression program required for cell adaptation to stress.

Characterization of human gene expression changes after olive oil ingestion: an exploratory approach

Olive oil consumption is protective against risk factors for cardiovascular and cancer diseases. A nutrigenomic approach was performed to assess whether changes in gene expression could occur in human peripheral blood mononuclear cells after oli ve oil ingestion at postprandial state. Six healthy male volunteers ingested, at fasting state, 50 ml of olive oil. Prior to intervention a 1-week washout period with a controlled diet and sunflower oil as the only source of fat was followed. During the 3 days before and on the intervention day, a very low-phenolic compound diet was followed. At baseline (0 h) and at post-ingestion (6 h), total RNA was isolated and gene expression (29,082 genes) was evaluated by microarray. From microarray data, nutrient-gene interactions were observed in genes related to metabolism, cellular processes, cancer, and atherosclerosis (e.g. USP48 by 2.16; OGT by 1.68-fold change) and associated processes such as inflammation (e.g. AKAP13 by 2.30; IL-10 by 1.66-fold change) and DNA damage (e.g. DCLRE1C by 1.47; POLK by 1.44- fold change). When results obtained by microarray were verified by qRT-PCR in nine genes, full concordance was achieved only in the case of up-regulated genes. Changes were observed at a real-life dose of olive oil, as it is daily consumed in some Mediterranean areas. Our results support the hypothesis that postprandial protective changes related to olive oil consumption could be mediated through gene expression changes.

Pro-inflammatory gene expression and neurotoxic effects of activated microglia are attenuated by absence of CCAAT/enhancer binding protein beta

Background. Microglia and astrocytes respond to homeostatic disturbances with profound changes of gene expression. This response, known as glial activation or neuroinflammation, can be detrimental to the surrounding tissue. The transcription factor CCAAT/enhancer binding protein ß (C/EBPß) is an important regulator of gene expression in inflammation but little is known about its involvement in glial activation. To explore the functional role of C/EBPß in glial activation we have analyzed pro-inflammatory gene expression and neurotoxicity in murine wild type and C/EBPß-null glial cultures. Methods. Due to fertility and mortality problems associated with the C/EBPß-null genotype we developed a protocol to prepare mixed glial cultures from cerebral cortex of a single mouse embryo with high yield. Wild-type and C/EBPß-null glial cultures were compared in terms of total cell density by Hoechst-33258 staining; microglial content by CD11b immunocytochemistry; astroglial content by GFAP western blot; gene expression by quantitative real-time PCR, western blot, immunocytochemistry and Griess reaction; and microglial neurotoxicity by estimating MAP2 content in neuronal/microglial cocultures. C/EBPß DNA binding activity was evaluated by electrophoretic mobility shift assay and quantitative chromatin immunoprecipitation. Results. C/EBPß mRNA and protein levels, as well as DNA binding, were increased in glial cultures by treatment with lipopolysaccharide (LPS) or LPS + interferon ¿ (IFN¿). Quantitative chromatin immunoprecipitation showed binding of C/EBPß to pro-inflammatory gene promoters in glial activation in a stimulus- and gene-dependent manner. In agreement with these results, LPS and LPS+IFN¿ induced different transcriptional patterns between pro-inflammatory cytokines and NO synthase-2 genes. Furthermore, the expressions of IL-1ß and NO synthase-2, and consequent NO production, were reduced in the absence of C/EBPß. In addition, neurotoxicity elicited by LPS+IFN¿-treated microglia co-cultured with neurons was completely abolished by the absence of C/EBPß in microglia.

Nutritional and insulin regulation of fatty acid synthetase and leptin gene expression through ADD1/SREBP1.

The ability to regulate specific genes of energy metabolism in response to fasting and feeding is an important adaptation allowing survival of intermittent food supplies. However, little is known about transcription factors involved in such responses in higher organisms. We show here that gene expression in adipose tissue for adipocyte determination differentiation dependent factor (ADD) 1/sterol regulatory element binding protein (SREBP) 1, a basic-helix-loop-helix protein that has a dual DNA-binding specificity, is reduced dramatically upon fasting and elevated upon refeeding; this parallels closely the regulation of two adipose cell genes that are crucial in energy homeostasis, fatty acid synthetase (FAS) and leptin. This elevation of ADD1/SREBP1, leptin, and FAS that is induced by feeding in vivo is mimicked by exposure of cultured adipocytes to insulin, the classic hormone of the fed state. We also show that the promoters for both leptin and FAS are transactivated by ADD1/SREBP1. A mutation in the basic domain of ADD1/SREBP1 that allows E-box binding but destroys sterol regulatory element-1 binding prevents leptin gene transactivation but has no effect on the increase in FAS promoter function. Molecular dissection of the FAS promoter shows that most if not all of this action of ADD1/SREBP1 is through an E-box motif at -64 to -59, contained with a sequence identified previously as the major insulin response element of this gene. These results indicate that ADD1/SREBP1 is a key transcription factor linking changes in nutritional status and insulin levels to the expression of certain genes that regulate systemic energy metabolism.

Characterization of Arabidopsis FPS isozymes and FPS gene expression analysis provide insight into the biosynthesis of isoprenoid precursors in seeds

Arabidopsis thaliana contains two genes encoding farnesyl diphosphate (FPP) synthase (FPS), the prenyl diphoshate synthase that catalyzes the synthesis of FPP from isopentenyl diphosphate (IPP) and dimethylallyl diphosphate (DMAPP). In this study, we provide evidence that the two Arabidopsis short FPS isozymes FPS1S and FPS2 localize to the cytosol. Both enzymes were expressed in E. coli, purified and biochemically characterized. Despite FPS1S and FPS2 share more than 90% amino acid sequence identity, FPS2 was found to be more efficient as a catalyst, more sensitive to the inhibitory effect of NaCl, and more resistant to thermal inactivation than FPS1S. Homology modelling for FPS1S and FPS2 and analysis of the amino acid differences between the two enzymes revealed an increase in surface polarity and a greater capacity to form surface salt bridges of FPS2 compared to FPS1S. These factors most likely account for the enhanced thermostability of FPS2. Expression analysis of FPS::GUS genes in seeds showed that FPS1 and FPS2 display complementary patterns of expression particularly at late stages of seed development, which suggests that Arabidopsis seeds have two spatially segregated sources of FPP. Functional complementation studies of the Arabidopsis fps2 knockout mutant seed phenotypes demonstrated that under normal conditions FPS1S and FPS2 are functionally interchangeable. A putative role for FPS2 in maintaining seed germination capacity under adverse environmental conditions is discussed.

Gene expression profiles in rat mesenteric lymph nodes upon supplementation with Conjugated Linoleic Acid during gestation and suckling

Background Diet plays a role on the development of the immune system, and polyunsaturated fatty acids can modulate the expression of a variety of genes. Human milk contains conjugated linoleic acid (CLA), a fatty acid that seems to contribute to immune development. Indeed, recent studies carried out in our group in suckling animals have shown that the immune function is enhanced after feeding them with an 80:20 isomer mix composed of c9,t11 and t10,c12 CLA. However, little work has been done on the effects of CLA on gene expression, and even less regarding immune system development in early life. Results The expression profile of mesenteric lymph nodes from animals supplemented with CLA during gestation and suckling through dam's milk (Group A) or by oral gavage (Group B), supplemented just during suckling (Group C) and control animals (Group D) was determined with the aid of the specific GeneChip® Rat Genome 230 2.0 (Affymettrix). Bioinformatics analyses were performed using the GeneSpring GX software package v10.0.2 and lead to the identification of 89 genes differentially expressed in all three dietary approaches. Generation of a biological association network evidenced several genes, such as connective tissue growth factor (Ctgf), tissue inhibitor of metalloproteinase 1 (Timp1), galanin (Gal), synaptotagmin 1 (Syt1), growth factor receptor bound protein 2 (Grb2), actin gamma 2 (Actg2) and smooth muscle alpha actin (Acta2), as highly interconnected nodes of the resulting network. Gene underexpression was confirmed by Real-Time RT-PCR. Conclusions Ctgf, Timp1, Gal and Syt1, among others, are genes modulated by CLA supplementation that may have a role on mucosal immune responses in early life.

Gene expression profiles in rat mesenteric lymph nodes upon supplementation with Conjugated Linoleic Acid during gestation and suckling

Background Diet plays a role on the development of the immune system, and polyunsaturated fatty acids can modulate the expression of a variety of genes. Human milk contains conjugated linoleic acid (CLA), a fatty acid that seems to contribute to immune development. Indeed, recent studies carried out in our group in suckling animals have shown that the immune function is enhanced after feeding them with an 80:20 isomer mix composed of c9,t11 and t10,c12 CLA. However, little work has been done on the effects of CLA on gene expression, and even less regarding immune system development in early life. Results The expression profile of mesenteric lymph nodes from animals supplemented with CLA during gestation and suckling through dam's milk (Group A) or by oral gavage (Group B), supplemented just during suckling (Group C) and control animals (Group D) was determined with the aid of the specific GeneChip® Rat Genome 230 2.0 (Affymettrix). Bioinformatics analyses were performed using the GeneSpring GX software package v10.0.2 and lead to the identification of 89 genes differentially expressed in all three dietary approaches. Generation of a biological association network evidenced several genes, such as connective tissue growth factor (Ctgf), tissue inhibitor of metalloproteinase 1 (Timp1), galanin (Gal), synaptotagmin 1 (Syt1), growth factor receptor bound protein 2 (Grb2), actin gamma 2 (Actg2) and smooth muscle alpha actin (Acta2), as highly interconnected nodes of the resulting network. Gene underexpression was confirmed by Real-Time RT-PCR. Conclusions Ctgf, Timp1, Gal and Syt1, among others, are genes modulated by CLA supplementation that may have a role on mucosal immune responses in early life.

Thermal evolution of gene expression profiles in Drosophila subobscura

Background: Despite its pervasiveness, the genetic basis of adaptation resulting in variation directly or indirectly related to temperature (climatic) gradients is poorly understood. By using 3-fold replicated laboratory thermal stocks covering much of the physiologically tolerable temperature range for the temperate (i.e., cold tolerant) species Drosophila subobscura we have assessed whole-genome transcriptional responses after three years of thermal adaptation, when the populations had already diverged for inversion frequencies, pre-adult life history components, and morphological traits. Total mRNA from each population was compared to a reference pool mRNA in a standard, highly replicated two-colour competitive hybridization experiment using cDNA microarrays.Results: A total of 306 (6.6%) cDNA clones were identified as 'differentially expressed' (following a false discovery rate correction) after contrasting the two furthest apart thermal selection regimes (i.e., 13°C vs . 22°C), also including four previously reported candidate genes for thermotolerance in Drosophila (Hsp26, Hsp68, Fst, and Treh). On the other hand, correlated patterns of gene expression were similar in cold- and warm-adapted populations. Analysis of functional categories defined by the Gene Ontology project point to an overrepresentation of genes involved in carbohydrate metabolism, nucleic acids metabolism and regulation of transcription among other categories. Although the location of differently expressed genes was approximately at random with respect to chromosomes, a physical mapping of 88 probes to the polytene chromosomes of D. subobscura has shown that a larger than expected number mapped inside inverted chromosomal segments.Conclusion: Our data suggest that a sizeable number of genes appear to be involved in thermal adaptation in Drosophila, with a substantial fraction implicated in metabolism. This apparently illustrates the formidable challenge to understanding the adaptive evolution of complex trait variation. Furthermore, some clustering of genes within inverted chromosomal sections was detected. Disentangling the effects of inversions will be obviously required in any future approach if we want to identify the relevant candidate genes.

Semiquantitative RT-PCR measurement of gene expression in rat tissues including a correction for varying cell size and number

Background: Current methodology of gene expression analysis limits the possibilities of comparison between cells/tissues of organs in which cell size and/or number changes as a consequence of the study (e.g. starvation). A method relating the abundance of specific mRNA copies per cell may allow direct comparison or different organs and/or changing physiological conditions. Methods: With a number of selected genes, we analysed the relationship of the number of bases and the fluorescence recorded at a present level using cDNA standards. A lineal relationship was found between the final number of bases and the length of the transcript. The constants of this equation and those of the relationship between fluorescence and number of bases in cDNA were determined and a general equation linking the length of the transcript and the initial number of copies of mRNA was deduced for a given pre-established fluorescence setting. This allowed the calculation of the concentration of the corresponding mRNAs per g of tissue. The inclusion of tissue RNA and the DNA content per cell, allowed the calculation of the mRNA copies per cell. Results: The application of this procedure to six genes: Arbp, cyclophilin, ChREBP, T4 deiodinase 2, acetyl-CoA carboxylase 1 and IRS-1, in liver and retroperitoneal adipose tissue of food-restricted rats allowed precise measures of their changes irrespective of the shrinking of the tissue, the loss of cells or changes in cell size, factors that deeply complicate the comparison between changing tissue conditions. The percentage results obtained with the present methods were essentially the same obtained with the delta-delta procedure and with individual cDNA standard curve quantitative RT-PCR estimation. Conclusion: The method presented allows the comparison (i.e. as copies of mRNA per cell) between different genes and tissues, establishing the degree of abundance of the different molecular species tested.

Minimization of biosynthetic costs in adaptive gene expression responses of yeast to environmental changes

Yeast successfully adapts to an environmental stress by altering physiology and fine-tuning metabolism. This fine-tuning is achieved through regulation of both gene expression and protein activity, and it is shaped by various physiological requirements. Such requirements impose a sustained evolutionary pressure that ultimately selects a specific gene expression profile, generating a suitable adaptive response to each environmental change. Although some of the requirements are stress specific, it is likely that others are common to various situations. We hypothesize that an evolutionary pressure for minimizing biosynthetic costs might have left signatures in the physicochemical properties of proteins whose gene expression is fine-tuned during adaptive responses. To test this hypothesis we analyze existing yeast transcriptomic data for such responses and investigate how several properties of proteins correlate to changes in gene expression. Our results reveal signatures that are consistent with a selective pressure for economy in protein synthesis during adaptive response of yeast to various types of stress. These signatures differentiate two groups of adaptive responses with respect to how cells manage expenditure in protein biosynthesis. In one group, significant trends towards downregulation of large proteins and upregulation of small ones are observed. In the other group we find no such trends. These results are consistent with resource limitation being important in the evolution of the first group of stress responses.

Use of physiological constraints to identify quantitative design principles for gene expression in yeast adaptation to heat shock

Background: Understanding the relationship between gene expression changes, enzyme activity shifts, and the corresponding physiological adaptive response of organisms to environmental cues is crucial in explaining how cells cope with stress. For example, adaptation of yeast to heat shock involves a characteristic profile of changes to the expression levels of genes coding for enzymes of the glycolytic pathway and some of its branches. The experimental determination of changes in gene expression profiles provides a descriptive picture of the adaptive response to stress. However, it does not explain why a particular profile is selected for any given response. Results: We used mathematical models and analysis of in silico gene expression profiles (GEPs) to understand how changes in gene expression correlate to an efficient response of yeast cells to heat shock. An exhaustive set of GEPs, matched with the corresponding set of enzyme activities, was simulated and analyzed. The effectiveness of each profile in the response to heat shock was evaluated according to relevant physiological and functional criteria. The small subset of GEPs that lead to effective physiological responses after heat shock was identified as the result of the tuning of several evolutionary criteria. The experimentally observed transcriptional changes in response to heat shock belong to this set and can be explained by quantitative design principles at the physiological level that ultimately constrain changes in gene expression. Conclusion: Our theoretical approach suggests a method for understanding the combined effect of changes in the expression of multiple genes on the activity of metabolic pathways, and consequently on the adaptation of cellular metabolism to heat shock. This method identifies quantitative design principles that facilitate understating the response of the cell to stress.

Forage preservation (grazing vs. hay) fed to ewes affects the fatty acid profile of milk and CPT1B gene expression in the sheep mammary gland

Background: Alterations in lipid metabolism occur when animals are exposed to different feeding systems. In the last few decades, the characterisation of genes involved in fat metabolism and technological advances have enabled the study of the effect of diet on the milk fatty acid (FA) profile in the mammary gland and aided in the elucidation of the mechanisms of the response to diet. The aim of this study was to evaluate the effect of different forage diets (grazing vs. hay) near the time of ewe parturition on the relationship between the fatty acid profile and gene expression in the mammary gland of the Churra Tensina sheep breed. Results: In this study, the forage type affected the C18:2 cis-9 trans-11 (CLA) and long-chain saturated fatty acid (LCFA) content, with higher percentages during grazing than during hay feeding. This may suggest that these FAs act as regulatory factors for the transcriptional control of the carnitine palmitoyltransferase 1B (CPT1B) gene, which was more highly expressed in the grazing group (GRE). The most highly expressed gene in the mammary gland at the fifth week of lactation is CAAT/ enhancer- binding protein beta (CEBPB), possibly due to its role in milk fat synthesis in the mammary gland. More stable housekeeping genes in the ovine mammary gland that would be appropriate for use in gene expression studies were ribosomal protein L19 (RPL19) and glyceraldehyde- 3- phosphate dehydrogenase (GAPDH). Conclusions: Small changes in diet, such as the forage preservation (grazing vs. hay), can affect the milk fatty acid profile and the expression of the CPT1B gene, which is associated with the oxidation of fatty acids. When compared to hay fed indoors, grazing fresh low mountain pastures stimulates the milk content of CLA and LCFA via mammary uptake. In this sense, LCFA in milk may be acting as a regulatory factor for transcriptional control of the CPT1B gene, which was more highly expressed in the grazing group.

An activator/repressor dual system allows tight tetracycline-regulated gene expression in budding yeast

We have developed an activator/repressor expression system for budding yeast in which tetracyclines control in opposite ways the ability of tetR-based activator and repressor molecules to bind tetO promoters. This combination allows tight expression of tetO-driven genes, both in a direct (tetracycline-repressible) and reverse (tetracycline-inducible) dual system. Ssn6 and Tup1, that are components of a general repressor complex in yeast, have been tested for their repressing properties in the dual system, using lacZ and CLN2 as reporter genes. Ssn6 gives better results and allows complete switching-off of the regulated genes, although increasing the levels of the Tup1-based repressor by expressing it from a stronger promoter improves repressing efficiency of the latter. Effector-mediated shifts between expression and non-expression conditions are rapid. The dual system here described may be useful for the functional analysis of essential genes whose conditional expression can be tightly controlled by tetracyclines.

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.