Human breast tumor slices: A model for identification of vitamin D regulated genes in the tumor microenvironment

While many studies have addressed the direct effects of 1 alpha,25(OH)(2)D(3) on breast cancer (BC) cells, stromal-epithelial interactions, which are important for the tumor development, have been largely ignored. In addition, high concentrations of the hormone, which cannot be attained in vivo, have been used. Our aim was to establish a more physiological breast cancer model, represented by BC tissue slices, which maintain epithelial-mesenchymal interactions, cultured with a relatively low 1 alpha,25(OH)(2)D(3) concentration, in order to evaluate the vitamin D pathway. Freshly excised human BC samples were sliced and cultured in complete culture media containing vehicle, 0.5 nM or 100 nM 1 alpha,25(OH)(2)D(3) for 24 h. BC slices remained viable for at least 24 h, as evaluated by preserved tissue morphology in hematoxylin and eosin (HE) stained sections and bromodeoxyuridine (BrdU) incorporation by 10% of tumor cells. VDR mRNA expression was detected in all samples and CYP24A1 mRNA expression was induced by 1 alpha,25(OH)(2)D(3) in both concentrations (but mainly with 100 nM). Our results indicate that the vitamin D signaling pathway is functional in BC slices, a model which preserves stromal-epithelial interactions and mimics in vivo conditions. (C) 2010 Elsevier Ltd. All rights reserved.

Arrow Plot: a new graphical tool for selecting up and down regulated genes and genes differentially expressed on samples subgroups

Background: A common task in analyzing microarray data is to determine which genes are differentially expressed across two (or more) kind of tissue samples or samples submitted under experimental conditions. Several statistical methods have been proposed to accomplish this goal, generally based on measures of distance between classes. It is well known that biological samples are heterogeneous because of factors such as molecular subtypes or genetic background that are often unknown to the experimenter. For instance, in experiments which involve molecular classification of tumors it is important to identify significant subtypes of cancer. Bimodal or multimodal distributions often reflect the presence of subsamples mixtures. Consequently, there can be genes differentially expressed on sample subgroups which are missed if usual statistical approaches are used. In this paper we propose a new graphical tool which not only identifies genes with up and down regulations, but also genes with differential expression in different subclasses, that are usually missed if current statistical methods are used. This tool is based on two measures of distance between samples, namely the overlapping coefficient (OVL) between two densities and the area under the receiver operating characteristic (ROC) curve. The methodology proposed here was implemented in the open-source R software. Results: This method was applied to a publicly available dataset, as well as to a simulated dataset. We compared our results with the ones obtained using some of the standard methods for detecting differentially expressed genes, namely Welch t-statistic, fold change (FC), rank products (RP), average difference (AD), weighted average difference (WAD), moderated t-statistic (modT), intensity-based moderated t-statistic (ibmT), significance analysis of microarrays (samT) and area under the ROC curve (AUC). On both datasets all differentially expressed genes with bimodal or multimodal distributions were not selected by all standard selection procedures. We also compared our results with (i) area between ROC curve and rising area (ABCR) and (ii) the test for not proper ROC curves (TNRC). We found our methodology more comprehensive, because it detects both bimodal and multimodal distributions and different variances can be considered on both samples. Another advantage of our method is that we can analyze graphically the behavior of different kinds of differentially expressed genes. Conclusion: Our results indicate that the arrow plot represents a new flexible and useful tool for the analysis of gene expression profiles from microarrays.

Study of the retinoic acid regulated genes in the rodent hypothalamus controlling growth and appetite and the potential for the epigenetic control

Dissertação de mestrado em Bioquímica (área de especialização em Biomedicina)

Identification of corticosteroid-regulated genes in cardiomyocytes by serial analysis of gene expression.

Corticosteroids (aldosterone, cortisol/corticosterone) exert direct functional effects on cardiomyocytes. However, gene networks activated by corticosteroids in cardiomyocytes, as well as the involvement of the mineralocorticoid receptor (MR) vs the glucocorticoid receptor (GR) in these effects, remain largely unknown. Here we characterized the corticosteroid-dependent transcriptome in primary culture of neonatal mouse cardiomyocytes treated with 10(-6) M aldosterone, a concentration predicted to occupy both MR and GR. Serial analysis of gene expression revealed 101 aldosterone-regulated genes. The MR/GR specificity was characterized for one regulated transcript, namely ecto-ADP-ribosyltransferase-3 (Art3). Using cardiomyocytes from GR(null/null) or MR(null/null) mice we demonstrate that in GR(null/null) cardiomyocytes the response is abrogated, but it is fully maintained in MR(null/null) cardiomyocytes. We conclude that Art3 expression is regulated exclusively via the GR. Our study identifies a new set of corticosteroid-regulated genes in cardiomyocytes and demonstrates a new approach to studying the selectivity of MR- vs GR-dependent effects.

Steady-state expression of self-regulated genes.

MOTIVATION: Regulatory gene networks contain generic modules such as feedback loops that are essential for the regulation of many biological functions. The study of the stochastic mechanisms of gene regulation is instrumental for the understanding of how cells maintain their expression at levels commensurate with their biological role, as well as to engineer gene expression switches of appropriate behavior. The lack of precise knowledge on the steady-state distribution of gene expression requires the use of Gillespie algorithms and Monte-Carlo approximations. METHODOLOGY: In this study, we provide new exact formulas and efficient numerical algorithms for computing/modeling the steady-state of a class of self-regulated genes, and we use it to model/compute the stochastic expression of a gene of interest in an engineered network introduced in mammalian cells. The behavior of the genetic network is then analyzed experimentally in living cells. RESULTS: Stochastic models often reveal counter-intuitive experimental behaviors, and we find that this genetic architecture displays a unimodal behavior in mammalian cells, which was unexpected given its known bimodal response in unicellular organisms. We provide a molecular rationale for this behavior, and we implement it in the mathematical picture to explain the experimental results obtained from this network.

Contribution of CgPDR1-Regulated Genes in Enhanced Virulence of Azole-Resistant Candida glabrata.

In Candida glabrata, the transcription factor CgPdr1 is involved in resistance to azole antifungals via upregulation of ATP binding cassette (ABC)-transporter genes including at least CgCDR1, CgCDR2 and CgSNQ2. A high diversity of GOF (gain-of-function) mutations in CgPDR1 exists for the upregulation of ABC-transporters. These mutations enhance C. glabrata virulence in animal models, thus indicating that CgPDR1 might regulate the expression of yet unidentified virulence factors. We hypothesized that CgPdr1-dependent virulence factor(s) should be commonly regulated by all GOF mutations in CgPDR1. As deduced from transcript profiling with microarrays, a high number of genes (up to 385) were differentially regulated by a selected number (7) of GOF mutations expressed in the same genetic background. Surprisingly, the transcriptional profiles resulting from expression of GOF mutations showed minimal overlap in co-regulated genes. Only two genes, CgCDR1 and PUP1 (for PDR1upregulated and encoding a mitochondrial protein), were commonly upregulated by all tested GOFs. While both genes mediated azole resistance, although to different extents, their deletions in an azole-resistant isolate led to a reduction of virulence and decreased tissue burden as compared to clinical parents. As expected from their role in C. glabrata virulence, the two genes were expressed as well in vitro and in vivo. The individual overexpression of these two genes in a CgPDR1-independent manner could partially restore phenotypes obtained in clinical isolates. These data therefore demonstrate that at least these two CgPDR1-dependent and -upregulated genes contribute to the enhanced virulence of C. glabrata that acquired azole resistance.

Cloning of glucocorticoid-regulated genes in C6/ST1 rat glioma phenotypic reversion

The C6 rat glioma cell line is responsive to glucocorticoid hormones. C6 variants that are hyper-responsive (ST1) and resistant (P7) to hormone treatment have been derived previously. Glucocorticoid treatment of ST1 cells leads to complete reversion of the transformed phenotype and loss of tumorigenic potential. Production of C type retrovirus particles is also induced by glucocorticoids in ST1 cells. Cloning of the genes regulated by glucocorticoids in this cell system was used here as a strategy to uncover the gene products involved in the transformed-to-normal phenotypic change. Construction of a cDNA library from glucocorticoid-treated ST1 cells and screening by differential hybridization resulted in the isolation of three cellular sequences that code for rat metallothioneins (C27 and C41) and α1-acid glycoprotein (C36). Northern blot analysis revealed that expression of these genes was dramatically induced by hydrocortisone in ST1 but not in P7 cells. Viral genomic RNA was used to isolate and characterize retrovirus-related sequences that could also be responsible for the phenotypic reversion phenomenon.

Identification and classification of p53-regulated genes

Sequence-specific transactivation by p53 is essential to its role as a tumor suppressor. A modified tetracycline-inducible system was established to search for transcripts that were activated soon after p53 induction. Among 9,954 unique transcripts identified by serial analysis of gene expression, 34 were increased more than 10-fold; 31 of these had not previously been known to be regulated by p53. The transcription patterns of these genes, as well as previously described p53-regulated genes, were evaluated and classified in a panel of widely studied colorectal cancer cell lines. “Class I” genes were uniformly induced by p53 in all cell lines; “class II” genes were induced in a subset of the lines; and “class III” genes were not induced in any of the lines. These genes were also distinguished by the timing of their induction, their induction by clinically relevant chemotherapeutic agents, the absolute requirement for p53 in this induction, and their inducibility by p73, a p53 homolog. The results revealed substantial heterogeneity in the transcriptional responses to p53, even in cells derived from a single epithelial cell type, and pave the way to a deeper understanding of p53 tumor suppressor action.

Identification of transforming growth factor-β- regulated genes in Caenorhabditis elegans by differential hybridization of arrayed cDNAs

Members of the transforming growth factor-β family play critical roles in body patterning, in both vertebrates and invertebrates. One transforming growth factor-β-related gene, dbl-1, has been shown to regulate body length and male ray patterning in Caenorhabditis elegans. We screened arrayed cDNAs to identify downstream target genes for the DBL-1 signaling by using differential hybridization. C. elegans cDNAs representing 7,584 independent genes were arrayed on a nylon membrane at high density and hybridized with 33P-labeled DNA probes synthesized from the mRNAs of wild-type, dbl-1, sma-2, and lon-2 worms. Signals for all the spots representing hybridized DNA were quantified and compared among strains. The screening identified 22 and 2 clones, which were positively and negatively regulated, respectively, by the DBL-1 signal. Northern hybridization confirmed the expression profiles of most of the clones, indicating good reliability of the differential hybridization using arrayed cDNAs. In situ hybridization analysis revealed the spatial and temporal expression patterns of each clone and showed that at least four genes, including the gene for the type I receptor for DBL-1, sma-6, were transcriptionally regulated by the DBL-1 signal.

Identification of diverse nerve growth factor-regulated genes by serial analysis of gene expression (SAGE) profiling

Neurotrophic factors such as nerve growth factor (NGF) promote a wide variety of responses in neurons, including differentiation, survival, plasticity, and repair. Such actions often require changes in gene expression. To identify the regulated genes and thereby to more fully understand the NGF mechanism, we carried out serial analysis of gene expression (SAGE) profiling of transcripts derived from rat PC12 cells before and after NGF-promoted neuronal differentiation. Multiple criteria supported the reliability of the profile. Approximately 157,000 SAGE tags were analyzed, representing at least 21,000 unique transcripts. Of these, nearly 800 were regulated by 6-fold or more in response to NGF. Approximately 150 of the regulated transcripts have been matched to named genes, the majority of which were not previously known to be NGF-responsive. Functional categorization of the regulated genes provides insight into the complex, integrated mechanism by which NGF promotes its multiple actions. It is anticipated that as genomic sequence information accrues the data derived here will continue to provide information about neurotrophic factor mechanisms.

The First Step of Gibberellin Biosynthesis in Pumpkin Is Catalyzed by at Least Two Copalyl Diphosphate Synthases Encoded by Differentially Regulated Genes

The first step in gibberellin biosynthesis is catalyzed by copalyl diphosphate synthase (CPS) and ent-kaurene synthase. We have cloned from pumpkin (Cucurbita maxima L.) two cDNAs, CmCPS1 and CmCPS2, that each encode a CPS. Both recombinant fusion CmCPS proteins were active in vitro. CPS are translocated into plastids and processed by cleavage of transit peptides. For CmCPS1 and CmCPS2, the putative transit peptides cannot exceed the first 99 and 107 amino acids, respectively, because longer N-terminal deletions abolished activity. Levels of both CmCPS transcripts were strictly regulated in an organ-specific and developmental manner. Both transcripts were almost undetectable in leaves and were abundant in petioles. CmCPS1 transcript levels were high in young cotyledons and low in roots. In contrast, CmCPS2 transcripts were undetectable in cotyledons but present at significant levels in roots. In hypocotyls, apices, and petioles, CmCPS1 transcript levels decreased with age much more rapidly than those of CmCPS2. We speculate that CmCPS1 expression is correlated with the early stages of organ development, whereas CmCPS2 expression is correlated with subsequent growth. In contrast, C. maxima ent-kaurene synthase transcripts were detected in every organ at almost constant levels. Thus, ent-kaurene biosynthesis may be regulated through control of CPS expression.

The Expression of Light-Regulated Genes in the High-Pigment-1 Mutant of Tomato1

Three light-regulated genes, chlorophyll a/b-binding protein (CAB), ribulose-1,5-bisphosphate carboxylase/oxygenase small subunit, and chalcone synthase (CHS), are demonstrated to be up-regulated in the high-pigment-1 (hp-1) mutant of tomato (Lycopersicon esculentum Mill.) compared with wild type (WT). However, the pattern of up-regulation of the three genes depends on the light conditions, stage of development, and tissue studied. Compared with WT, the hp-1 mutant showed higher CAB gene expression in the dark after a single red-light pulse and in the pericarp of immature fruits. However, in vegetative tissues of light-grown seedlings and adult plants, CAB mRNA accumulation did not differ between WT and the hp-1 mutant. The ribulose-1,5-bisphosphate carboxylase/oxygenase small subunit mRNA accumulated to a higher level in the hp-1 mutant than WT under all light conditions and tissues studied, whereas CHS gene expression was up-regulated in de-etiolated vegetative hp-1-mutant tissues only. The CAB and CHS genes were shown to be phytochrome regulated and both phytochrome A and B1 play a role in CAB gene expression. These observations support the hypothesis that the HP-1 protein plays a general repressive role in phytochrome signal transduction.

Products of Proline Catabolism Can Induce Osmotically Regulated Genes in Rice1

Many plants accumulate high levels of free proline (Pro) in response to osmotic stress. This imino acid is widely believed to function as a protector or stabilizer of enzymes or membrane structures that are sensitive to dehydration or ionically induced damage. The present study provides evidence that the synthesis of Pro may have an additional effect. We found that intermediates in Pro biosynthesis and catabolism such as glutamine and Δ1-pyrroline-5-carboxylic acid (P5C) can increase the expression of several osmotically regulated genes in rice (Oryza sativa L.), including salT and dhn4. One millimolar P5C or its analog, 3,4-dehydroproline, produced a greater effect on gene expression than 1 mm l-Pro or 75 mm NaCl. These chemicals did not induce hsp70, S-adenosylmethionine synthetase, or another osmotically induced gene, Em, to any significant extent. Unlike NaCl, gene induction by P5C did not depend on the normal levels of either de novo protein synthesis or respiration, and did not raise abscisic acid levels significantly. P5C- and 3,4-dehydroproline-treated plants consumed less O2, had reduced NADPH levels, had increased NADH levels, and accumulated many osmolytes associated with osmotically stressed rice. These experiments indicate that osmotically induced increases in the concentrations of one or more intermediates in Pro metabolism could be influencing some of the characteristic responses to osmotic stress.

Gene repression by the ferric uptake regulator in Pseudomonas aeruginosa: cycle selection of iron-regulated genes.

The expression of at least 24 distinct genes of Pseudomonas aeruginosa PAO1 is under direct control of the "ferric uptake regulator" (Fur). Novel targets of the Fur protein were isolated in a powerful SELEX (systematic evolution of ligands by exponential enrichment)-like cycle selection consisting of in vitro DNA-Fur interaction, binding to anti-Fur antibody, purification on protein G, and PCR amplification. DNA fragments obtained after at least three exponential enrichment cycles were cloned and subjected to DNA mobility-shift assays and DNase I footprint analyses to verify the specific interaction with the Fur protein in vitro. Iron-dependent expression of the corresponding genes in vivo was monitored by RNase protection analysis. In total, 20 different DNA fragments were identified which represent actual Pseudomonas iron-regulated genes (PIGs). While four PIGs are identical to already known genes (pfeR, pvdS, tonB, and fumC, respectively), 16 PIGs represent previously unknown genes. Homology studies of the putative proteins encoded by the PIGs allowed us to speculate about their possible function. Two PIG products were highly similar to siderophore receptors from various species, and three PIG products were significantly homologous to alternative sigma factors. Furthermore, homologs of the Escherichia coli ORF1-tolQ, nuoA, stringent starvation protein Ssp, and of a two-component regulatory system similar to the Pseudomonas syringae LemA sensor kinase were identified. The putative gene products of seven additional PIGs did not show significant homologies to any known proteins. The PIGs were mapped on the P.aeruginosa chromosome. Their possible role in iron metabolism and virulence of P. aeruginosa is discussed.

Role of the GATA factors Gln3p and Nil1p of Saccharomyces cerevisiae in the expression of nitrogen-regulated genes.

We have isolated the NIL1 gene, whose product is an activator of the transcription of nitrogen-regulated genes, by virtue of the homology of its zinc-finger domain to that of the previously identified activator, the product of GLN3. Disruption of the chromosomal NIL1 gene enabled us to compare the effects of Gln3p and of Nil1p on the expression of the nitrogen-regulated genes GLN1, GDH2, and GAP1, coding respectively for glutamine synthetase, NAD-linked glutamate dehydrogenase, and general amino acid permease. Our results show that the nature of GATAAG sequence that serve as the upstream activation sequence elements for these genes determines their abilities to respond to Gln3p and Nil1p. The results further indicate that Gln3p is inactivated by an increase in the intracellular concentration of glutamine and that Nil1p is inactivated by an increase in intracellular glutamate.