Short-term arginine deprivation results in large-scale modulation of hepatic gene expression in both normal and tumor cells: microarray bioinformatic analysis
Nutrition & Metabolism volume 3, Article number: 37 (2006)
We have reported arginine-sensitive regulation of LAT1 amino acid transporter (SLC 7A5) in normal rodent hepatic cells with loss of arginine sensitivity and high level constitutive expression in tumor cells. We hypothesized that liver cell gene expression is highly sensitive to alterations in the amino acid microenvironment and that tumor cells may differ substantially in gene sets sensitive to amino acid availability. To assess the potential number and classes of hepatic genes sensitive to arginine availability at the RNA level and compare these between normal and tumor cells, we used an Affymetrix microarray approach, a paired in vitro model of normal rat hepatic cells and a tumorigenic derivative with triplicate independent replicates. Cells were exposed to arginine-deficient or control conditions for 18 hours in medium formulated to maintain differentiated function.
Initial two-way analysis with a p-value of 0.05 identified 1419 genes in normal cells versus 2175 in tumor cells whose expression was altered in arginine-deficient conditions relative to controls, representing 9–14% of the rat genome. More stringent bioinformatic analysis with 9-way comparisons and a minimum of 2-fold variation narrowed this set to 56 arginine-responsive genes in normal liver cells and 162 in tumor cells. Approximately half the arginine-responsive genes in normal cells overlap with those in tumor cells. Of these, the majority was increased in expression and included multiple growth, survival, and stress-related genes. GADD45, TA1/LAT1, and caspases 11 and 12 were among this group. Previously known amino acid regulated genes were among the pool in both cell types. Available cDNA probes allowed independent validation of microarray data for multiple genes. Among genes downregulated under arginine-deficient conditions were multiple genes involved in cholesterol and fatty acid metabolism. Expression of low-density lipoprotein receptor was decreased in both normal and tumor cells.
Arginine-sensitive regulation appears to be an important homeostatic mechanism to coordinate cell response and nutrient availability in hepatic cells. Genes predicted as arginine-responsive in stringent microarray data analysis were confirmed by Northern blot and RT-PCR. Although the profile of arginine-responsive genes is altered and increased, a considerable portion of the "arginome" is maintained upon neoplastic transformation.
Cell growth is dependent upon availability of essential amino acids for protein synthesis and this relationship makes amino acid-dependent control of gene expression an important area of study [1–3]. We previously reported that levels of the tumor-associated glycoprotein amino acid transporter, TA1/LAT1/CD98 light chain/SLC7A5, increase in normal hepatic cells under low arginine conditions, while levels are constitutive and high in hepatic tumor cells [4, 5]. Upregulation of this gene is associated with multiple cancer types and we and others have hypothesized that increased expression may provide an adaptive advantage in the tumor microenvironment where nutrients are limiting [6, 7]. Loss of nutrient-sensitive regulation may comprise a subset of the loss of responsiveness to negative growth signals or autonomy from positive growth factors, characteristic of cancer cells in general . We thus manipulated arginine concentration in a culture medium that maintains hepatic differentiation as a means to investigate, using a paired normal and tumorigenic cell type, how many and what types of hepatic genes are responsive to a transient change in amino acid levels. Affymetrix gene chips and bioinformatic analyses were used in a nutrigenomics approach. The goals of the study were to: 1) assess the scope of arginine-responsive hepatic gene expression using a well-defined in vitro rat model of normal and tumorigenic cells, 2) determine to what extent amino acid responsive regulation is retained upon transformation, and 3) provide a microarray dataset predicting novel genes and pathways subject to amino acid (arginine) regulation.
Results and discussion
Triplicate RAE 230 arrays containing a total of 15,923 rat genes were hybridized for both normal and tumor cells grown under arginine-sufficient (+) and arginine-deficient (-) conditions for 18 hours in three independent pair-wise comparison experiments. Microarray Suite™ 5.0 software analysis of the resulting datasets revealed a mean of 8611 genes or 54% of the total rat array scored "present" (expressed) for normal cells versus 8355 genes or 53% in tumor cells (Table 1). Of the genes expressed, pair-wise comparisons of Arg + and - conditions revealed no change in expression for 7042 (82%) of genes in normal cells versus 6173 (74%) in tumor cells. Using a p-value of 0.05 or less and pair-wise comparisons for the three experiments in normal cells, expression of 811 genes increased, 608 decreased and 97 showed marginal change relative to arginine-sufficient controls. For tumor cells, transient arginine deficiency resulted in increased expression of 1249 genes, decreased expression of 926 genes, and marginal changes in 109. Thus, expression of as many as 2175 genes was scored as arginine-sensitive when the data were analyzed as independent pairs. This represented 8.9% of genes in normal cells and 13.6% of genes in tumor cells. At the time of the data collection, nearly half of these genes lacked complete functional annotation in the rat genomic database.
Because the number of genes in the preliminary analysis comprised a large portion of the rat genome, additional filters were applied to the dataset using GeneSpring™ 6.1 and 7.0 software analysis in order to identify the most reproducible and significant alterations resulting from transient arginine deficiency. For each cell type, all + samples were compared to all – samples, resulting in 9 comparisons. A concordance threshold of 66% or greater was set and 'arginine-responsive' increases or decreases were defined as a change of 2-fold and higher under arginine-deficient conditions. All other genes expressed were defined as "no change".
Using these more stringent parameters, the total number of arginine-responsive genes in normal cells and tumor cells was significantly reduced, to 56 and 162 respectively. A summary of the data comparing normal and tumor cells by functional gene class and direction of change is shown in Table 2. All 56 and 162 genes in normal and tumor cells respectively listed by accession numbers and corresponding direction of change relative to arginine-sufficient conditions (I = increase; D = decrease) are listed in supplemental Tables 4 and 5 [see additional files 1 and 2]. The entire 12 microarray dataset is available in MIAME format  in the NCBI Gene Expression Omnibus (GEO) database, confirmation #GSE2275. Figure 1 presents Venn diagrams of overlap between the arginine-responsive gene subsets of normal and tumor cells. Just over half the arginine-responsive genes in normal cells are also responsive in tumor cells. A surprising finding was the large number (133) of additional genes demonstrating arginine-responsive expression in tumor cells. More than one third of these were listed as transcribed sequences with no listed functional information or homology to known sequences at the time of analysis. Not surprisingly, genes with increased expression upon arginine depletion in both normal and tumor cells include a large number of genes involved in cell cycle and growth regulation, cell death and apoptosis, and stress response. Several of these were known amino acid-responsive genes such as asparagine synthetase . Among other genes scoring as increased in both cell types was growth-arrest DNA damage-inducible 45 alpha (GADD45) and caspase 12. Genes down-regulated in both normal and tumor cells include plasminogen activator urokinase and low density lipoprotein receptor. Genes that exhibited arginine-responsive regulation in normal cells but a loss in tumor cells (or vice versa) may represent 'interesting' subsets in which the potential contribution of amino-acid regulation to the malignant phenotype could be explored. Genes whose expression was altered in tumor cells but not normal cells included tumor suppressor retinoblastoma-like 2, stress-response protein 70 kDa heat shock protein precursor and cell-surface linked signal transducer MAP-kinase phosphatase, cell cycle regulator cyclin dependent kinase inhibitor 2C (p18) and Ras. Genes scoring as changed in normal but "no change" in tumor cells were mainly in cholesterol and steroid biosynthesis pathways. Related to these were farnesyl diphosphate synthase and insulin-induced gene 1/growth response protein CL-6 (Insig-1/CL-6), an ER-residing membrane protein that blocks the proteolytic activation of sterol-regulatory element binding protein (SREBPs) transcription factors that activate the synthesis of fatty acids and cholesterol [11–13].
Independent verification of expression was sought for a subset of genes via Northern blot analysis (Figure 2) and RT-PCR (Figure 3). Genes were selected in which cDNA probes were available with appropriate known controls to validate selected cases in which arginine-deprivation predicted an increase, decrease or no significant change in expression. Though the level of expression for some genes was found to differ substantially between normal and tumor samples, between arginine + and - conditions respectively, or between independent experimental replicates, trends predicted in microarray data analysis were generally validated by the confirmatory RNA expression analysis and densitometry. That is, LAT1, GADD45 and 4F2 increased in both cell types under conditions of arginine deprivation; LDLr decreased in both cell types upon arginine deprivation; LAT2 and GAPDH showed no significant change in either cell type, though LAT2 expression was very low in tumor cells; p21 decreased in normal cells but was below detection limits in tumor cells; farnesyl diphosphate synthase (FDPS) decreased in normal cells but not in tumor cells upon arginine deprivation; Insig 1 showed large variability with decreases in some replicates but not in others. The use of different probe sets between the microarray and Northern or PCR analyses, differences in quantitation and normalization methodology or undetermined other variables may account for differences in fold change values observed between these techniques.
Given its enormous potential for insight into regulatory physiology, nutrigenomics, the use of genomic tools to study nutrient-gene interactions, is an area of expanding interest in health and disease [14–16]. We found global changes in gene expression when hepatic cells were subjected to short term deficiency of arginine. Arginine, a precursor of proteins and other compounds critical to normal metabolism, is considered essential for growth and development of many mammals [17, 18]. Our experimental system was one of amino acid imbalance, but not of total arginine starvation, since the medium was supplemented with 5% non-dialyzed fetal bovine serum and the exposure to deficient conditions was less than 24 hr. In fact, arginine-free medium has been widely used as a means to deplete primary hepatocyte cultures of small numbers of contaminating nonparenchymal cells because all the enzymes necessary for arginine synthesis are expressed by the normal hepatocyte . Importantly, viability was well preserved over the experimental timeframe (data not shown). Thus, the observed changes in gene expression were not the result of impending cell death but are more likely a response to nutrient stress.
Stress caused by changes in the cellular environment including nutrients is recognized as an important physiological parameter to which eukaryotic cells have evolved complex cellular responses [19–22]. Although it is not possible to differentiate primary effects from secondary during this timeframe, overall, cell growth and maintenance, metabolism, stress-related and apoptotic gene classes appear most predominant. Genes associated with metabolism and stress response that were present in the rat chip and showed a significant p-value for altered expression upon arginine deprivation are presented in Table 3. Because many genes associated with stress response in human cells were not present or identifiable on the RAE 230 chip however, it is not possible to fully assess arginine response of these gene sets in this model. With the exception of a subset of genes involved in cholesterol and steroid biosynthesis, there was no overall loss of amino-acid responsive regulation in the tumor cells compared to normal cells. In fact, to our surprise, in assessing the entire rat transcriptome a larger proportion of genes was found to be arginine-responsive in tumor cells compared to the normal cell from which they were derived. The use of this paired cell set in the experimental design is an important consideration as it minimizes variation due to other factors such as genetic strain or viral exposure. Our data not only add to the nutrigenomics literature and database in general, but also extend the finding of others in which changes in amino acid composition or starvation have been associated with specific patterns of hepatic gene response both in vivo and in vitro [23–25]. The data suggest that arginine deprivation of hepatic normal and tumor cells may induce an initial stress response very similar to the well known endoplasmic reticulum (ER) stress-response described in eukaryotic cell types . The location of products of several arginine-responsive genes including caspase 12, LDLr, Insig-1 and cholesterol/steroid biosynthesis all in the ER are consistent with the induction of a potent luminal ER stress response in both normal and tumor cells. Interestingly, although both cell types appear able to modulate LDLr, a regulator of cholesterol and steroid biosynthesis, in response to arginine deprivation, at some point downstream of the stress pathway, tumor cells appear to deviate by losing regulation of Insig-1 and metabolic genes directly involved in cholesterol and steroid biosynthesis. While it is tempting to speculate that loss of regulation within this gene cluster may be adaptive in a tumor microenvironment of limited nutrient availability, permitting malignant cells to persist in membrane synthesis required for repeated cell divisions, endocytosis, fusion or other processes, the relevance of these observed changes in gene expression in vitro to potential alterations in hepatic function and tumor growth in vivo are currently unknown. This hypothesis is testable but beyond the limits of the current project. The generation and reporting of nutrigenomic data is undeniably descriptive. However, datasets such as this one are valuable in that they facilitate hypothesis generation and further research. Much more study is needed to define specific genes and mechanistic pathways linking amino acids and other dietary nutrients to cancer risk.
The rat hepatocyte normal diploid cell line WB344 ("normal") and GP7TB, a chemically-transformed derivative capable of generating hepatic tumors upon transplantation to syngeneic rats ("tumor"), were obtained from the laboratory of Dr. William Coleman, Dept. Pathology, University of North Carolina, Chapel Hill and maintained as previously described . To assess arginine-responsive gene expression, at no more than 80% confluence, media was removed and cells were rinsed briefly with a custom formulation of Chee's Essential Medium (CEM) without arginine (GIBCO Invitrogen Corp., Carlsbad, CA). Medium was then replaced with CEM with or without L-Arginine HCl at 0.168 gm/liter for an additional 18 hours at 37 C and 5% CO2. Media included 5% fetal bovine serum. The CEM media formulation was utilized because it has been demonstrated to maintain differentiated function and gene expression in cultured hepatic cells .
RNA isolation and microarray analysis
Cells were harvested by lysis-extraction of total RNA with TriReagent (Molecular Research Center, Inc., Cincinnati, OH). RNA quantity and quality were assessed by absorbance at 260 and 280 nm using a Genesys 5 spectrophotometer (Spectronic, Leeds, UK). Electrophoresis of the RNA samples followed by ethidium bromide staining was used to further evaluate RNA samples. Only non-degraded RNAs without DNA contamination were utilized. Five μg of total RNA were converted to the first strand of cDNA by using SuperScript II RT (Invitrogen Corporation, Carlsbad, CA), and gel-purified T7-oligo dT(24) (W.M. Keck Foundation, New Haven, CT) as the primer. Second strand synthesis was performed using E coli Polymerase I, DNA ligase, E coli RNA H and T4 DNA Polymerase according to the manufacturer's instructions and reagents in the SuperScript Double-Stranded cDNA Synthesis Kit (Invitrogen). The resulting cDNA was purified by extraction with (25:24:1) phenol: chloroform: isoamyl alcohol (Ambion, Inc., Austin, TX), separation of the aqueous phase using phase lock gels and precipitation in ethanol. The purified cDNA subsequently served as template for production of biotin labeled cRNA transcript using the BioArray High Yield RNA Transcript Labeling Kit (ENZO Biochemical, New York, NY) and biotin-labeled UTP and CTP. Labeled cRNA was isolated using the RNeasy Mini Kit columns (Qiagen, Valencia, CA), quantified for purity, concentration and yield and subsequently fragmented in 100 mM potassium acetate-30 mM potassium acetate-40 mM Tris-acetate (pH 8.1) for 35 minutes at 94 C to generate 35–200 bp fragments, suitable for hybridization. Eukaryotic Target Hybridization was accomplished with 10 μg of fragmented cRNA, eukaryotic hybridization controls, herring sperm DNA (Promega Biosciences, Inc., San Luis Obsipo, CA), acetylated BSA (Invitrogen), 2× hybridization buffer and DEPC treated water to make the cocktail final concentration 100 mM MES and 1 M [Na+]. Following clarification of the cocktail by heating and centrifugation, rat 230A expression arrays (Affymetrix, Santa Clara, CA) were loaded with the target cRNA cocktail and rotated at 60 rpm for 16 hours at 45 C. Following hybridization, the cocktail was removed and the microarrays were washed with Non-stringent buffer [6× SSPE, 0.01% Tween 20 (Pierce, Rockford, IL) 0.005% Antifoam (Sigma, St Louis, MO)] at 30 C and Stringent buffer (100 mM MES, 0.1 M [Na+], 0.01% Tween 20) at 50 C in a Affymetrix GeneChip Fluidics Station 400 according to the manufacturer's protocol. Microarrays were subsequently stained with Streptavidin-conjugated Phycoerythrin (SAPE, Invitrogen, Carlsbad, CA) and staining intensity was antibody amplified using a biotinylated anti-streptavidin antibody (Vector Laboratories, Burlingame, CA) according to Affymetrix protocols. The GeneChips were scanned at 570 nm using the Agilent Technologies G2500A GeneArray scanner with MicroArray Suite software. All hybridization steps were performed at the Brown University COBRE Center for Genomics and Proteomics Microarray Facility. Each hybridization was performed in triplicate from separate RNA prepared from independent cell cultures of both normal and tumor cells.
Each gene on the 230A array is represented by 20 different 25-base cDNA oligonucleotides complementary to a cRNA target transcript (perfect match) together with specificity control oligonucleotides containing a single base substitution (mismatch) for each perfect-match. The combination of perfect-match and mismatch cDNA oligonucleotides for each gene is termed a probe set. Affymetrix-defined absolute mathematical algorithms describing perfect-match and mismatch hybridization intensities were used to define each gene as "present" or "absent" and assign a value. Binding intensity values were scaled to evaluate differential expression +/- arginine supplementation. Affymetrix-defined comparison mathematical algorithms determined whether a transcript was classified as "increased", "decreased", "marginally increased", "marginally decreased" or "not changed" and a fold change in expression was calculated. Only genes that had signal present were selected from the normalized data for further analysis. This was accomplished using GeneSpring version 7.0 (Silicon Genetics, Redwood City, CA) importing data from Affymetrix MicroArray Suite (MAS 5.0) into GeneSpring as tab delimited text files. Additional criteria/filters were used to classify a gene as significantly up- or down-regulated upon short term arginine starvation: (i) expression of a gene must be classified as increased or marginally increased (or decreased or marginally decreased) in each replicate compared to the arginine sufficient condition. (ii) The mean fold change for arginine starvation-induced gene expression must be greater than 2. Mean transcriptional expression of a given gene was calculated as the sum of the fold change in gene expression for each arginine-minus condition compared to each arginine-sufficient (normal) condition divided by 6. Standard errors in mean transcriptional expression of a gene were also calculated. Comparison between gene expression levels among groups was done using the Two-Way ANOVA test.
All microarray data were submitted in compliance with the Minimal Information about Microarray Experiments (MIAME) via the Gene Expression Omnibus (GEO) data repository under the title "A Bioinformatic Analysis of Arginine-Sensitive Regulation of rat Hepatic Gene Expression", confirmation #GSE2275 . The MIAME standard, developed by the Microarray Gene Expression Data Society, specifies information to permit experimental reproducibility, standardization of data reporting and sharing, and allows global access to the original microarray data . Annotation of the differentially expressed genes and classification into functional groups was done using the Database for Annotation, Visualization, and Integrated Discovery (DAVID) program (NIH).
Northern blot confirmation of gene expression
Aliquots (10 μg) of total RNA were size fractionated on 1% agarose/formaldehyde gels as described previously . After electrophoresis, gels were equilibrated in 1 M ammonium acetate and RNA was transferred to nylon membranes. Blots were baked for 2 hr at 80 C. Hybridization with cDNA probes was carried out at 42 C in ULTRAhyb hybridization solution (Ambion, Austin, TX) sequentially to random primed, 32P-labeled cDNAs of interest. Probes for TA1/LAT1, 4F2/CD98 and GAPDH were described previously . The following partial cDNAs were generated by PCR from known positive template RNA to confirm expression patterns of additional genes of interest: rat GADD45, a 706 bp PCR fragment representing nucleotides 1 to 706 [GenBank NM_024127.1]; rat Insig1/CL6, a 1392 bp fragment representing nucleotides 339 to 1730 [GenBank L13619]; rat LAT2, a 990 bp fragment corresponding to nucleotides 580 to 1570 [GenBank NM_053442.1]. These partial cDNAs were PCR amplified under individually optimized conditions with template RNA from positive control rat tissues or cell lines and subcloned into pCR4 Blunt-TOPO plasmid vector. Rat p21 cDNA, a 316 bp fragment representing nucleotides 318 to 634 [GenBank U24174] was the kind gift of Dr. Phillip Gruppuso, RI Hospital. The identity of all cDNA probes was confirmed by DNA sequencing (Keck Biotechnology Laboratory, Yale University (New Haven, CT). Purified inserts were digested from vector with EcoR1. For LAT2, EcoR1 digestion resulted in two insert products and the larger, 660 bp fragment was used as probe. After hybridization and washing, blots were exposed to X-ray film in the presence of intensifying screens or a Cyclone Phosphorimager (Model A431201) (Downers Grove, IL).
RT-PCR confirmation of gene expression
Semi-quantitative Reverse Transcription Polymerase Chain Reaction was used to assess relative expression of the following genes by amplification of these products: rat LDLr, a 531 bp fragment representing nucleotides 106 to 636 [GenBank NM_175762]; rat farnesyl diphosphate synthase (FDPS), a 656 bp fragment representing nucleotides 229 to 884 [GenBank NM_031840.1]; rat GADD45, a 706 bp fragment representing nucleotides 1 to 706 [GenBank NM_024127.1]; Insig-1, a 447 bp fragment representing nucleotides 1284 to 1730 [GenBank L13619] and rat GAPDH, a 452 bp fragment representing nucleotides 1369 to 1820 [GenBank NM_017008]. RNAs used for Northern blot analysis were used to generate cDNA and PCR products. Rat adult liver total RNA served as a positive control. Three micrograms of total RNA in a final volume of 20 μl was DNase I (Invitrogen, Carlsbad, CA) treated at room temperature for 15 minutes and the reaction stopped by addition of 25 mM EDTA followed by incubation at 65 C for 10 minutes to inactivate the DNase I enzyme. DNase treated total RNA was reverse-transcribed with and without reverse transcriptase for 60 minutes at 50 C using Superscript III Reverse Transcriptase (Invitrogen, Carlsbad, CA), Oligo (dT)12–18 (Invitrogen, Carlsbad, CA), RNase OUT ribonuclease inhibitor (Invitrogen, Carlsbad, CA, and 10 mM DNTP set (Invitrogen, Carlsbad, CA. Platinum Taq DNA polymerase (Invitrogen, Carlsbad, CA) was used for the RT-PCR reactions, which included 100 ng cDNA per 20 μl containing a final concentration of 1.5–2.25 mM MgCl2 0.2–0.3 mM dNTP, 0.5 μM each forward and reverse primers, and 0.025 U Platinum Taq DNA polymerase. PCR reactions were run in a Px2 Thermal Cycler (Themo Electron Corp.). LDLr and FDPS reactions were 94 C for 60 s followed by 94 C for 30 s, 55 C for 30 s, 72 C for 30 s (30 cycles), and a final incubation at 72 C for 10 minutes. GADD45 reactions were cycled at 94 C for 60 s, 50 C for 45 s, 72 C for 90 s for (40 cycles) and a final incubation at 72 C for 10 minutes. Insig-1 reactions were cycled at 94 C for 30 s, 55 C for 30 s, 72 C for 90 s (35 cycles) and a final incubation at 72 C for 10 minutes. GAPDH was cycled at 94 C for 60 s, followed by 94 C for 30 s, 53 C for 30 s, 72 C for 30 s (30 cycles), and a final incubation at 72 C for 10 minutes. Electrophoresis of PCR products was performed in 2% agarose, 1× TAE gel and visualized by ethidium bromide staining. Densitometry was run on LabWorks Software and values were ratioed to GAPDH densitometry values.
L-type amino acid transporter, type 1
Kilberg MS, Pan YX, Chen H, Leung-Pineda V: Nutritional control of gene expression: how mammalian cells respond to amino acid limitation. Annu Rev Ntr. 2004
Kimball SR, Jefferson LS: Amino acids as regulators of gene expression. Nutr Metab (London). 2004, 1: 3-10.1186/1743-7075-1-3.
Cousins RJ: Nutritional regulation of gene expression. Am J Med. 1999, 106: 20S-23S. 10.1016/S0002-9343(98)00342-8.
Shultz VD, Campbell W, Karr S, Hixson DC, Thompson NL: TA1 oncofetal rat liver cDNA and putative amino acid permease: temporal correlation with c-myc during acute CCl liver injury and variation of RNA levels in response to amino acids in hepatocyte cultures. Toxicol Appl Pharmacol. 1999, 154: 84-96. 10.1006/taap.1998.8555.
Campbell WA, Sah DE, Medina MM, Albina JE, Coleman WB, Thompson NL: TA1/LAT-1/CD98 light chain and system L activity but not 4F2/CD98 heavy chain, respond to arginine availability in rat hepatic cells. J Biol Chem. 2000, 275: 5347-5354. 10.1074/jbc.275.8.5347.
Wolfe DA, Panzica MA, Wang S, Bassily NH, Thompson NL: Expression of a highly conserved oncofetal gene, TA1/E16, in human colon carcinoma and other primary tumors: homology to S. mansoni amino acid permease and C. elegans gene products. Cancer Res. 1996, 56: 5012-5022.
Fuchs BC, Bode BP: Amino acid transporters ASCT2 and LAT1 in cancer: partners in crime?. Semin Cancer Biol. 2005, 15: 254-266. 10.1016/j.semcancer.2005.04.005.
Hanahan D, Weinberg RA: The hallmarks of cancer. Cell. 2005, 100: 57-70. 10.1016/S0092-8674(00)81683-9.
Brazma A, Hingamp P, Quackenbush J, Sherlock G, Spellman P, Stoeckert C, Aach J, Ansorge W, Ball CA, Causton HC, Gaasterland T, Glenisson P, Holstege FCP, Kim IF, Markowitz V, Matese JC, Parkinson H, Robinson A, Sarkans U, Schulze-Kremer S, Steward J, Taylor R, Vilo J, Vingron M: Minimum information about a microarray experiment (MIAME) – toward standards for microarray data. Nat Genet. 2001, 29: 365-371. 10.1038/ng1201-365.
Farfournoux P, Bruhat A, Jousse C: Amino acid regulation of gene expression. Biochem J. 2000, 351: 1-12. 10.1042/0264-6021:3510001.
Attie AD: Insig: a significant integrator of nutrient and hormonal signals. J Clin Invest. 2004, 113: 1112-1114. 10.1172/JCI200421450.
Lee JN, Ye J: Proteolytic activation of sterol regulatory element-binding protein induced by cellular stress through depletion of Insig-1. J Biol Chem. 2004, 279: 45257-45265. 10.1074/jbc.M408235200.
Sun LP, Li L, Goldstein JL, Brown MS: Insig requirement for sterol-mediated inhibition of Scap/SREBP binding to COPII proteins in vitro. J Biol Chem. 2005, 280: 26483-26490. 10.1074/jbc.M504041200.
Corthesy-Theulaz I, den Dunnen JT, Ferre P, Geurts JM, Muller M, van Belzen N, van Ommen B: Nutrigenomics: the impact of biomics technology on nutrition. Ann Nutr Metab. 2005, 49: 355-365. 10.1159/000088315.
Mutch DM, Wahli W, Williamson G: Nutrigenomics and nutrigenetics: the emerging faces of nutrition. FASEB J. 2005, 19: 1602-1616. 10.1096/fj.05-3911rev.
Saito K, Arai S, Kato H: A nutrigenomics database – integrated repository for publications and associated microarray data in nutrigenomics research. Br J Nutr. 2005, 94: 493-495. 10.1079/BJN20051536.
Grillo MA, Colombatto S: Arginine revisited: Minireview article. Amino Acids. 2004, 26: 345-351.
Evoy D, Lieberman MD, Fahey TJ, Daly JM: Immunonutrition: the role of arginine. Nutr. 1998, 14: 611-617. 10.1016/S0899-9007(98)00005-7.
Pan Y, Chen H, Siu F, Kilberg MS: Amino acid deprivation and endoplasmic reticulum stress induce expression of multiple activating transcription factor-3 mRNA species that, when overexpressed in HepG2 cells, modulate transcription by the human asparagine synthetase promoter. J Biol Chem. 2003, 278: 38402-38412. 10.1074/jbc.M304574200.
Hinnebusch AG: Translational regulation of gcn4 and the general amino acid control of yeast. Annu Rev Microbiol. 2005, 59: 407-450. 10.1146/annurev.micro.59.031805.133833.
Huang Z, Tunnacliffe A: Gene induction by desiccation stress in human cell cultures. FEBS Lett. 2005, 579: 4973-4977. 10.1016/j.febslet.2005.07.084.
Durukawa-Hibi Y, Kobayashi Y, Chen C, Motoyama N: FOXO transcription factors in cell-cycle regulation and the response to oxidative stress. Antioxd Redox Signal. 2005, 7: 752-760. 10.1089/ars.2005.7.752.
Nelsen C, Rickheim D, Tucker M, McKenzie T, Hansen L, Pestell R, Albrecht J: Amino acids regulate hepatocyte proliferation through modulation of cyclin D1 expression. J Biol Chem. 2003, 278: 25853-25858. 10.1074/jbc.M302360200.
Zhang J, Underwood LE, D'ercole AJ: Hepatic mRNAs up-regulated by starvation: an expression profile determined by suppression subtractive hybridization. FASEB J. 2001, 15: 1261-1263. 10.1096/fj.01-0260com.
Spiegelberg T, Bishop JO: Tissue-specific gene expression in mouse hepatocytes cultured in growth-restricting medium. Mol Cell Biol. 1988, 8: 3338-3344.
Gene Expression Omnibus (GEO) data repository. http://www.ncbi.nlm.nih.gov/projects/geo/
Supported by American Institute for Cancer Research grant 02A127; NIH 1RO1 CA73611; P20 RR017695 and P20 RR015578.
Bonnie Lau is supported by a National Research Service Award F30 ES013639-01.
Declaration of competing interests
The author(s) declare that they have no competing interests.
HWL carried out portions of the bioinformatic analysis, submitted the MIAME dataset to GEO, prepared manuscript figures and tables for the original manuscript and participated in the initial text draft. CS prepared the cRNAs and conducted the microarray hybridization, scanning, raw data collection and supervised primary bioinformatic data analysis. ALB isolated total RNA from cell cultures, conducted Northern blot hybridization, generated PCR probes, RT-PCR quantitation and figure preparation, and participated in data analysis. BL participated in the initial microarray bioinformatic analysis, RT-PCR and data interpretation. CF carried out the cell culture and participated in figure preparation and unpublished confirmatory experiments. ES conducted additional bioinformatics analysis assessing specific classes of genes including stress response, cholesterol biosynthesis and other categories for a table in the revised manuscript; NLT designed experiments, supervised the overall project and prepared the final version of the manuscript for submission.
Electronic supplementary material
Additional file 1: Major functional classes of arginine-sensitive genes in normal hepatic cells. GeneSpring functional categories of arginine-responsive genes in normal hepatic cells based on 9-way comparison of three independent data sets. Numbers of genes in each class and subclass are shown in parentheses along with gene accession numbers. Increase (I) and Decrease (D) defined as a two fold or higher change in expression from + arg to - arg. (XLS 26 KB)
Additional file 2: Major functional classes of arginine-sensitive genes in hepatic tumor cells. GeneSpring functional categories of arginine-responsive genes in tumorigenic hepatic cells based on 9-way comparison of three independent data sets. Numbers of genes in each class and subclass are shown in parentheses along with gene accession numbers. Increase (I) and Decrease (D) defined as a two fold or higher change in expression from + arg to - arg. (XLS 60 KB)
About this article
Cite this article
Leong, H.X., Simkevich, C., Lesieur-Brooks, A. et al. Short-term arginine deprivation results in large-scale modulation of hepatic gene expression in both normal and tumor cells: microarray bioinformatic analysis. Nutr Metab (Lond) 3, 37 (2006). https://doi.org/10.1186/1743-7075-3-37