TOXICOLOGICAL SCIENCES 100(1), 75–87 (2007) doi:10.1093/toxsci/kfm200 Advance Access publication August 6, 2007

Exposure to Arsenic at Levels Found in U.S. Drinking Water Modifies Expression in the Mouse Lung

Angeline S. Andrew,*,†,‡,1 Viviane Bernardo,§,{ Linda A. Warnke,k Jennifer C. Davey,‡,kj Thomas Hampton,‡,kj Rebecca A. Mason,*,kj Jessica E. Thorpe,k Michael A. Ihnat,k and Joshua W. Hamilton†,‡,kj *Department of Community and Family Medicine, Dartmouth Medical School, Lebanon, New Hampshire 03756; †Norris Cotton Cancer Center, Dartmouth- Hitchcock Medical Center, Lebanon, New Hampshire 03756; ‡Center for Environmental Health Sciences, Dartmouth Medical School, Hanover, New Hampshire Downloaded from https://academic.oup.com/toxsci/article/100/1/75/1624780 by guest on 29 September 2021 03755; §Thayer School of Engineering/Computer Sciences Department, Dartmouth College, Hanover, New Hampshire 03755; {Health Informatics Department, Federal University of Sao Paulo/Escola Paulista de Medicina-UNIFESP/EPM, Sao Paulo, SP, Brazil; kDepartment of Cell Biology, University of Oklahoma Health Sciences Center, Oklahoma City, Oklahoma 73190; and kjDepartment of Pharmacology and Toxicology, Dartmouth Medical School, Hanover, New Hampshire 03755

Received April 6, 2007; accepted July 30, 2007

smoking and other risk factors to induce mutations, DNA The mechanisms of action of drinking water arsenic in the lung adducts, and cancer risk (Ahsan and Thomas, 2004; Chen and the threshold for biologic effects remain controversial. Our et al., 2004; Evans et al., 2004; Rossman et al., 2002). The study utilizes Affymetrix 22,690 transcript oligonucleotide micro- mechanisms of action and the threshold for biologic effects and arrays to assess the long-term effects of increasing doses of drinking water arsenic on expression levels in the mouse lung. disease risk, however, remain controversial (Kitchin, 2001). Mice were exposed at levels commonly found in contaminated Many of the epidemiologic studies and in vivo / in vitro drinking water wells in the United States (0, 0.1, 1 ppb), as well experiments have been conducted at high doses with acute as the 50 ppb former maximum contaminant level, for 5 weeks. exposure. Thus far, microarray studies have explored arsenic- The expression profiles revealed modification of a number of induced expression changes in cell culture models as well important signaling pathways, many with corroborating evidence as several animal organs; however, the effects in the adult lung of arsenic responsiveness. We observed statistically significant have not been reported (Durham and Snow, 2006; Shi et al., expression changes for transcripts involved in angiogenesis, lipid 2004). In contrast to a strengthening of similar effects with metabolism, oxygen transport, apoptosis, cell cycle, and immune increasing dose that is observed with many compounds, arsenic response. Validation by reverse transcription–PCR and immuno- may activate completely different pathways at low versus high blot assays confirmed expression changes for a subset of tran- doses (Andrew et al., 2003; Barchowsky et al., 1999; Lau scripts. These data identify arsenic-modified signaling pathways that will help guide investigations into mechanisms of arsenic’s et al., 2004; Soucy et al., 2003). Relatively little is known health effects and clarify the threshold for biologic effects and about the effects of arsenic at levels of exposure that are potential disease risk. common to drinking water in the United States, particularly in Key Words: arsenic; apoptosis; cell cycle; drinking water; the lung. immune response; lung; microarray; oxygen transport. Our study utilizes oligonucleotide microarrays to assess the long-term effects of levels of drinking water arsenic commonly found in contaminated wells in the United States, as well as the 50 ppb former standard (0, 0.1, 1, and 50 ppb) on expression Drinking water arsenic exposure is an established human levels in the mouse lung. We present a functional approach health risk associated with cardiovascular disease, diabetes, and to microarray analysis in which we focus on statistically cancer at multiple organ sites including the lung, bladder, and significant expression modifications and investigate both the skin (IARC, 2004). Arsenic exposure has been associated with biologic function and how each regulates another using impaired lung function and bronchiectasis (De et al., 2004; updated bioinformatics analysis tools. Smith et al., 2006). In an Indian study, arsenic-associated skin lesions were independently associated with increased risk of respiratory illness (odds ratio 4.9; 95% confidence interval 3.2– 7.5) compared to individuals with no skin lesions (Ghosh et al., MATERIALS AND METHODS 2007). Arsenic has also been found to synergize with cigarette Four groups of adult (7–8 weeks old) male C57/BL6 mice were exposed to 1 To whom correspondence should be addressed at Dartmouth Medical increasing concentrations of inorganic arsenic (n ¼ 4 at 0.1 lg/l [ppb], n ¼ 3 School, 7927 Rubin 860, One Medical Center Drive, Lebanon, NH 03756. Fax: at 1 lg/l [ppb], and n ¼ 3at50lg/l [ppb]) (sodium arsenite; LabTech, (603) 653-9093. E-mail: [email protected]. Pittsburgh, PA) in their drinking water for a period of 5 weeks. A control group

Ó The Author 2007. Published by Oxford University Press on behalf of the Society of Toxicology. All rights reserved. For Permissions, please email: [email protected] 76 ANDREW ET AL. of mice (n ¼ 3) consumed uncontaminated drinking water. Lungs from We further characterized the functional effects of all the statistically additional groups of mice exposed to 0 or 10 lg/l (ppb) were used in the reverse significant arsenic-modified transcripts (includes all of the SAM two-class transcription (RT)–PCR and immunoblot analyses. All mice were fed AIN-76A selected transcripts) by implementing the Database for Annotation, Visualiza- arsenic-free chow (Harlan-Teklad, Madison, WI) to control for other dietary tion, and Integrated Discovery (DAVID) search engine (Table 2) sources of arsenic (total arsenic < 5 ppb). Levels of arsenic in food and water (Dennis et al., 2003). This bioinformatic tool identifies functional processes were confirmed by Inductively Coupled Plasma Mass Spectrometry analysis by that are significantly overrepresented ( p < 0.05) by the modified transcripts. the Dartmouth Trace Metal Facility. Lungs were removed and immediately Table 3 lists the biologic roles of all the significant arsenic-modified transcripts placed in RNAlater to stabilize the RNA levels after sacrificing the animals. (all those selected by SAM multiclass or two-class analysis at any exposure) Exposed and unexposed animals were sacrificed at the same time of day. All from the Pathway Studio ResNet 4.0 database (updated April 2007) (Ariadne protocols were approved by the University of Oklahoma Institutional Animal Genomics, Rockville, MD). This Pathway Studio ResNet database catalogs Care and Use Committee. relationships between biologic entities based on the published literature and RNA was isolated using Qiagen RNeasy columns (Qiagen Inc., Valencia, was also used to identify the direct interactions between the arsenic-modified CA), and DNAse treatment was performed using Ambion DNAfree reagents transcripts at each arsenic dose (Fig. 3 and Supplementary Data). (Austin, TX) according to the manufacturer’s instructions. The expression Histologic slices of formalin-fixed, paraffin-embedded lung tissue from Downloaded from https://academic.oup.com/toxsci/article/100/1/75/1624780 by guest on 29 September 2021 profiles were generated using the Affymetrix GeneChip Technology—chip arsenic-exposed and unexposed animals were stained with hematoxylin and GeneChip Murine Genome 430 oligonucleotide arrays (Affymetrix, Santa Clara, eosin and evaluated by a trained pathologist to assess the number of neutrophils CA) which simultaneously tested 22,690 transcripts using one chip for each (Supplementary Data). mouse on the integrated GeneChip Instrument System in the Dartmouth The levels of Nr4a1, Hsp70, Ahr, and Cyclin D1 protein were assessed by Microarray Core Facility. Our experiment was performed in compliance with the immunoblotting using sodium dodecyl sulfate–polyacrylamide gel electropho- Minimum Information About a Microarray Experiment checklist for standard- resis (SDS–PAGE) to resolve from mouse lung tissue (Fig. 4). Frozen ization guidelines for microarray experiments. mouse lung was weighed and homogenized with EBC lysis buffer (50mM Tris, The raw microarray data were preprocessed and normalized using pH 8.0, 100mM NaCl) containing 10 ll/ml PMSF, 5 ll/ml aprotinin, and 1 ll/ml GeneTraffic version 3.2 which is a microarray data management and analysis leupeptin. NP-40 Lysis buffer (10%) was added at a 5% vol/vol ratio. After client-server application (Stratagene, La Jolla, CA). The control (0 lg/l) group centrifugation for 15 min at 14,000 rpm at 4°C the lysates were boiled for was set as the baseline, and data were normalized using Robust Multi-Chip 5 min and clarified by centrifugation at 13,000 rpm for 10 min. Equal amounts Analysis (Irizarry et al., 2003). of cell lysates were resolved by electrophoresis on 7.5% and 10–20% SDS– The statistical significance of expression changes relative to the control polyacrylamide gels. Electrophoresis was performed at constant voltage (200 V), group was assessed using methods based on modified t- or F-tests that adjust then the resolved proteins were transferred from the polyacrylamide gel to for multiple comparisons. This adjustment bounds a false discovery rate polyvinylidene difluoride membrane (PVDF, Immobilon-P; Millipore, Bedford, probability (FDR), i.e., chance that a transcript regarded as significant is a false MA) by semi-dry transfer (Hoeffer Semiphor, San Fransisco, CA) for 1 h at positive, to select the transcripts (Benjamini and Hochberg, 1995). This method constant current (100 mA) using transfer buffer (25mM Tris, 192mM glycine, is implemented in the Significance Analysis of Microarrays (SAM) application 20% (vol/vol) methanol, 0.01% SDS). To eliminate nonspecific interactions of version 1.13 (Tusher et al., 2001) implemented in The Institute of Genomic antibodies with the membrane, the PVDF membrane was blocked with Tris- Research (TIGR) MultiExperiment Viewer (TIGR MeV) version 3.1 (TIGR, Tween Buffered Saline (TTBS) (10mM Tris–HCl, pH 8.0, 150mM NaCl, Rockville, MD). We performed SAM multiclass and two-class comparisons 0.05% Tween-20) containing 5% milk (7.5 g/150 ml) for 1 h at room using 1000 permutations and selected significant transcripts at a FDR of 5%. temperature. Membranes were incubated with the primary antibody: anti-Nur77/ We focused our attention on the group of differentially expressed transcripts Nr4a1 (Pharmingen, San Diego, CA), Hsp70/Hspa1b (Transduction Labora- tories, Lexington, KY) diluted 1:1000, anti-Ahr (BIOMOL, Plymouth Meeting, that were called significant by the two-class SAM or multiclass SAM (dose PA) diluted 1:5000, or Cyclin D1 antibody (Santa Cruz Biotechnology Inc., response) analysis (provided as Supplementary Data). The transcripts that were Santa Cruz, CA) diluted 1:500 in TTBS overnight at 4°C. Actin was used as significant by multiclass analysis or were modified at a minimum of two doses a loading control, and the antibody was diluted 1:50,000 in TTBS for 1 h of arsenic are shown in Table 1. Transcripts are grouped by the lowest arsenic (Calbiochem, San Diego, CA). The membranes were washed 3 times with dose with statistically significant transcript expression modification compared TTBS. The Nur77/Nr4a1 and Hsp70 membranes were incubated with with control. To identify subgroups of transcripts with similar patterns of horseradish peroxidase (HRP)–linked goat anti-rabbit IgG (Santa Cruz expression among the statistically significant group, the multiclass selected Biotechnology Inc.) 1:3000 in TTBS with 5% milk (1.5 g/30 ml) for 1 h at transcripts were clustered by TIGR MeV with a hierarchical clustering with the room temperature. The Ahr, Cyclin D1, and Actin membranes were incubated complete linkage algorithm and Pearson correlation metric (Fig. 1). with HRP-linked goat anti-mouse (Bio-Rad Laboratories, Inc., Hercules, CA) We selected transcripts for validation by real-time PCR using independent 1:2000, 1:3000, and 1:2000, respectively, in TTBS with 5% milk (1.5 g/30 ml). primer sets based on the microarray results and hypothesized involvement in After 3 washes with TTBS, protein bands were visualized by enhanced arsenic pathogenesis in the lung (Fig. 2). Taqman primer–probe sets for each chemiluminescence using the Amersham ECL Plus Western Blotting Detection selected transcript were obtained from Applied Biosystems Inc. (ABI, Foster system (GE Healthcare, Piscataway, NJ) and film (Lumi-Film; Roche City, CA). Real-time RT–PCR was performed using the ABI PRISM sequence Molecular Biochemicals, Indianapolis, IN). detection system and software. Briefly, total RNA (0.5 lg) was reverse transcribed using 100 U Moloney Murine Leukemia Virus reverse transcriptase in a mixture with oligo-dT and dNTPs according to the instructions provided with the Qiagen Omniscript kit (Qiagen Inc.). Samples were reverse transcribed RESULTS in a MJ Research PTC-100 thermocycler (MJ Research Inc., Watertown, MA) for 60 min at 44°C and the reaction terminated by heating to 95°C for 10 min. To address the controversy over the biological effects of Expression of specific was assessed by real-time PCR using 10 ng total arsenic at levels commonly found in U.S. drinking water, we RNA, 400nM primers, 200nM probe, and TaqMan Universal PCR Master Mix assessed expression patterns associated with arsenic exposure (ABI). Relative quantitation was performed using a standard curve consisting in the mouse lung. To evaluate the statistical significance of of serial dilutions of pooled sample cDNA from the same source as the test RNA with each plate. Relative expression levels of each gene were normalized these changes in expression between the treated and control to 18s rRNA. Statistical significance was assessed by one-way ANOVA with animals, we started with the hypothesis that arsenic would have Newman-Keuls post-test using GraphPad Prism software (San Diego, CA). a dose-responsive effect on expression. In Figure 1, we TABLE 1 Multiple Dose Significantly Modified Transcriptsa Organized by Arsenic Dose

Fold change Probes significantd (n) Trend Figure 1 Symbol Gene name Functionb gene ID (þ/)c 0.1 ppb 1 ppb 50 ppb 0.1 ppb 1 ppb 50 ppb Multi cluster Downloaded fromhttps://academic.oup.com/toxsci/article/100/1/75/1624780bygueston29September2021

Down 0.1 ppb Angptl4 Angiopoietin-like 4 Negative regulation of apoptosis, 57875 1.0 0.2 0.2 1 1 c1 inhibits angiogenesis, inhibits lipoprotein lipase Cd68 CD68 antigen Macrosialin (CD68), macrophage- 12514 0.6 0.6 0.2 1 1 c1 specific, increased by aging in brain Cpt1a Carnitine palmitoyltransferase 1a, Fatty acid metabolism 12894 0.6 0.6 0.3 1 1 c1 liver Cxcl7 Chemokine (C-X-C motif) ligand 7 Iimmune response, chemokine, 57349 1.1 0.8 0.3 1 1 c1 cytokine MICROARRAY ARSENIC WATER DRINKING D4Wsu53e DNA segment, Chr 4, Wayne State NA 27981 0.9 0.4 0.7 2 1 c1 University 53 Eif4b Eukaryotic translation initiation Regulation of translational initiation 75705 0.4 0.5 0.4 2 1 c3 factor 4B Hbb-b1 Hemoglobin, beta adult major chain Hemoglobin complex member, 15129 0.6 0.6 0.3 1 hemopoiesis, oxygen transport, iron ion binding Hbp1 High mobility group box Transcriptional repressor, substrate 73389 0.5 0.6 0.1 1 1 c1 transcription factor 1 for p38 MAP kinase Ms4a8a Membrane-spanning 4-domains, A, Receptor activity 64381 0.4 0.2 0.2 1 1 c1 8A Pabpc1 Poly A binding protein, RNA binding 18458 0.5 0.6 0.6 1 1 c3 cytoplasmic 1 S100a9 S100 calcium-binding protein A9 Inflammatory response, neutrophil 20202 2.4 2.7 2.0 1 (calgranulin B) release Siat8d ST8 alpha-N-acetyl-neuraminide Protein amino acid glycosylation 20452 0.5 1.5 1.2 1 1 c3 alpha-2,8-sialyltransferase 4 Down 0.1 and 1 ppb Alas2 Aminolevulinic acid synthase 2, Heme biosynthesis, metabolism 5- 11656 1.5 1.4 1.2 1 1 1 c3 erythroid aminolevulinate synthase activity Bpgm 2,3-Bisphosphoglycerate mutase Glycolysis, hemoglobin activity 12183 1.2 1.2 1.0 2 1 1 c3 Cd53 CD53 antigen macrophage protection against 12508 0.5 0.7 0.6 1 1 1 c3 LPS-induced oxidative stress and UVB irradiation Ear1 Eosinophil-ssociated, ribonuclease A Hydrolase activity, endonuclease 13586 1.2 1.3 0.4 1 1 1 c1 family, member 1 activity Ear2 Eosinophil-associated, ribonuclease Endonuclease activity 13587 1.0 1.2 0.3 2 1 1 c1 A family, 2 Hba-a1 Hemoglobin alpha, adult chain 1 Hemoglobin complex member, 15122 1.4 1.2 1.1 1 1 oxygen transport, iron ion binding 77 TABLE 1—Continued 78

Fold change Probes significantd (n) Entrez Trend Figure 1 Symbol Gene name Functionb gene ID (þ/)c 0.1 ppb 1 ppb 50 ppb 0.1 ppb 1 ppb 50 ppb Multi cluster

Iigp1 Interferon inducible GTPase 1 Induced by various gram-negative 60440 0.6 0.8 0.6 1 1 1 c3 Downloaded fromhttps://academic.oup.com/toxsci/article/100/1/75/1624780bygueston29September2021 lipopolysaccharides, aids development of intracellular resistance during the interferon response to infection Ltgb2 Integrin beta 2 Activated T cell proliferation, cell- 16414 0.6 0.8 0.5 1 1 matrix adhesion Mkrn1 Makorin, ring finger protein, 1 Nucleic acid binding 54484 0.8 0.9 0.7 2 1 1 c3 Snca Synuclein, alpha Protects against oxidative stress via 20617 0.8 0.7 0.6 2 1 1 c1 inactivating the c-Jun-N-terminal kinase Down 0.1 and 50 ppb Cte1 Cytosolic acyl-CoA thioesterase 1 Acyl-CoA thioesterase activity, long- 26897 0.6 0.8 0.8 2 1 1 c1 chain fatty acid metabolism Eif3s5 Eukaryotic translation initiation Expression of a truncated eIF3e 66085 0.6 0.7 0.6 1 1 factor 3, subunit 5 (epsilon) causes malignant transformation

Fech Ferrochelatase Exon 10-deleted ferrochelatase 14151 0.5 0.3 0.5 1 1 AL. ET ANDREW heterozygous mice exhibited skin photosensitivity Gltscr2 Glioma tumor suppressor candidate NA 68077 0.5 0.7 0.5 2 1 region gene 2 Hlx H2.0-like homeo box gene NA 15284 0.3 0.5 0.6 1 1 Rgl1 Ral guanine nucleotide dissociation Guanyl-nucleotide exchange factor 19731 0.4 0.4 0.6 1 1 1 c3 stimulator-like 1 activity small GTPase mediated signal transduction Down 0.1, 1, and 50 ppb AI448196 Expressed sequence AI448196 NA 102910 0.6 0.8 0.7 1 1 1 1 c3 Stk17b Serine/threonine kinase 17b Plays critical roles in T cell apoptosis 98267 0.6 0.8 0.9 2 1 2 1 c3 (apoptosis-inducing) and memory T cell development Down 1 ppb Ahr Aryl-hydrocarbon receptor Cell cycle, xenobiotic metabolism, 11622 0.5 1.0 0.6 1 1 c3 regulation of transcription, DNA- dependent Cd79a CD79A antigen (immunoglobulin- Cell surface receptor linked signal 12518 0.4 1.0 0.8 1 1 c3 associated alpha) transduction, defense response, couples the B-cell antigen receptor to distal signaling pathways Ian6 Immune-associated nucleotide 6 Belongs immuno-associated 231931 0.5 1.0 0.8 1 1 c3 nucleotide (IAN) subfamily of nucleotide-binding proteins Lmo2 LIM domain only 2 involved in T-cell tumorigenesis due 16909 0.4 0.9 0.8 1 1 c3 to reprogramming of gene expression after enforced expression in T-cell precursors Down 1, 50 ppb Coro1a Coronin, actin-binding protein 1A Actin binding 12721 0.4 0.9 0.8 2 2 c3 Dbp D site albumin promoter-binding Regulation of transcription, DNA 13170 0.5 1.0 1.2 1 1 2 c3 protein binding, circadian rhythm

Plk2 Polo-like kinase 2 (Drosophila) ATP binding, protein serine/threonine 20620 0.2 0.8 0.7 1 1 1 c3 Downloaded fromhttps://academic.oup.com/toxsci/article/100/1/75/1624780bygueston29September2021 kinase activity cell cycle, protein amino acid phosphorylation Ptprc Protein tyrosine phosphatase, Cellular defense response, protein 19264 0.4 1.0 0.7 1 1 c3 receptor type, C amino acid dephosphorylation Down 50 ppb Agtrl1 Angiotensin receptor–like 1 Inhibits glucose-stimulated insulin 23796 0.7 1.5 1.0 1 1 c3 secretion Fus Fusion, derived from t(12;16) Positive regulation of transcription 233908 0.5 0.3 0.5 2 2 c3 malignant liposarcoma (human) from Pol II promoter Igh-VJ558 Immunoglobulin heavy chain (J558 B-cell antigen recognition. 16061 0.2 0.2 1.2 3 1 c4 family) Igj Immunoglobulin joining chain Humoral immune response, antigen 16069 0.3 0.0 1.1 1 1 c4 MICROARRAY ARSENIC WATER DRINKING binding. Igk Immunoglobulin kappa chain Immunoglobulin light chain variable 243469 0.1 0.3 1.1 1 1 c4 complex region Igk-V28 Immunoglobulin kappa chain Humoral immune response, antigen 16114 0.1 0.8 1.8 3 3 c4 variable 28 (V28) binding. Igk-V8 Immunoglobulin kappa chain Humoral immune response, antigen 16123 0.1 0.3 1.1 1 1 c4 variable 8 (V8) binding. Ltb Lymphotoxin B Tumor necrosis factor receptor 16994 0.5 0.6 0.9 1 1 c3 binding, immune response Nr1d1 Nuclear receptor subfamily 1, Regulation of transcription, DNA- 217166 0.6 0.7 0.9 1 1 c3 group D, member 1 dependent Nr4a1 Nuclear receptor subfamily 4, Regulation of transcription, steroid 15370 0.1 0.9 1.3 1 1 c3 group A, member 1 hormone receptor, apoptosis Down multiclass only Cd79b CD79B antigen Humoral immune response, maintains 15985 0.4 1.0 1.1 1 c3 preB cell and immature B cell survival and to mediate B cell differentiation Cirbp Cold inducible RNA binding protein RNA binding 12696 0.4 0.8 0.8 1 c3 Elf5 E74-like factor 5 Regulation of transcription 13711 0.1 0.6 0.4 1 c3 Myo10 Myosin X ATP binding, cytoskeleton 17909 0.0 0.6 0.6 1 c3 organization and biogenesis Myo6 Myosin VI ATP binding, cytoskeleton 17920 0.1 0.5 0.3 1 c3 organization and biogenesis Up multiclass only Fkbp5, Fkbp51 FK506 binding protein 5 Protein folding, glucocorticoid 14229 þ 0.2 0.6 0.8 1 c2 signaling Hspa8 Heat shock protein 8 ATPase activity, chaperone activity, 15481 þ 0.2 1.0 0.3 1 c2 regulation of cell cycle Mosc1 MOCO sulphurase C-terminal NA 66112 þ 0.1 0.8 0.6 1 c2 domain containing 1 MT2 Metallothionein 2 Metal ion binding 17750 þ 0.3 0.6 0.7 79 80 ANDREW ET AL.

hierarchically clustered transcripts that were selected by the multiclass SAM analysis (59 transcripts) and observed that the cluster Figure 1 data broke into four main branches, indicating subgroups of transcripts (labeled Clusters 1–4). Visual inspection of cluster 1 indicates decreased levels of transcripts involved in controlling )

n apoptosis, fatty acid metabolism, and chemokines at the 0.1 (

d and 1 ppb doses, but not at 50 ppb. Cluster 2 shows increased levels of transcripts involved in embryonic limb morphogen- esis, glucocorticoid signaling, and antiapoptosis, particularly at the 1 and 50 ppb doses. The immune response transcripts in cluster 4 are decreased mainly at the 50 ppb dose (Igj, Downloaded from https://academic.oup.com/toxsci/article/100/1/75/1624780 by guest on 29 September 2021 Igh-VJ558, Igk-V28). We observed a consistent pattern of st two different arsenic doses. decreased expression at all arsenic doses in the largest group, 0.01 ppb, control vs. 1 ppb, control vs. 50 ppb) or cluster 3. These genes are involved in a variety of processes including signal transduction, hemoglobin activity, glycolysis, 0.3 1 c4 transcription, apoptosis, and glycosylation. Arsenic exposure at 1 ppb strongly decreased Alas2, Siat8d, Agtrl, and Ear1 and increased expression of Zfp145, Fkbp5, and Hspa1b. 0.6

Two-class SAM analysis of each arsenic exposure dose

Fold change Probes significant compared with control (0 ppb) revealed 94 statistically signifi- 0.3 1.3 1.5 1 c2 0.7 0.4 1.60.2 1.2 1.2 1.4cantly modified 2 c2 1 transcripts c2 for 0.1 versus control. We chose 0.1 ppb 1 ppb 50 ppb 0.1 ppb 1 ppb 50 ppb Multi

, 2005). this lowest dose (0.1 ppb) as a low level that has been associated with adverse health effects in a Finnish study and is c ) et al. þ þ þ / / / /

/ found in U.S. drinking water sources (Kurttio et al., 1999). We þ þ Trend ( detected 26 modified transcripts for 1 ppb arsenic versus con-

.’’ trol and 37 modified transcripts for the 50 ppb arsenic versus

Continued control analysis (FDR < 0.05). Table 1 lists the subset of these 13653 15511 Entrez

gene ID transcripts that were selected as statistically significant by the multiclass SAM or by two-class SAM analyses at a minimum of two arsenic doses, organized by the lowest dose with signi- TABLE 1— ficant expression modification (all significantly modified tran-

ww.informatics.jax.org) (Eppig scripts are provided in the Supplementary Table). The direction

b of expression modification for each transcript that was significantly modified is marked with a ‘‘’’ (indicating decreased expression) or a ‘‘þ’’ (indicating increased expression). The ratio of each treated group versus control and the number of modified probe sets are shown in the adjoining columns. The ’’ and decreased expression designated ‘‘

þ modified transcripts were involved functional processes including T-cell differentiation. inhibition of caspase activation, telomere maintenance

Embryonic limb morphogenesis 235320 transcription factors, immune response, oxygen transport, cell

Informatics database (http://w cycle, oxidative stress, or fatty acid metabolism. While many transcripts showed consistent responses at all doses of arsenic, several were most dramatically decreased at the 0.1 ppb dose, including Angptl4 and Cxcl7 (Table 1, Down 0.1 ppb). Others such as Ear1 and Ear2 were decreased strongly at the 0.1 and 1 ppb doses, but not at 50 ppb (Table 1, Down 0.1, 1 ppb). A few of the transcripts had a differential response to the 0.1 ppb compared to 1 or 50 ppb arsenic (Table 1, containing 16 Expressed sequence AI467657 NA 102538 Early growth response 1 Transcription factor, zinc ion binding, Heat shock protein 1B Anti-apoptosis, DNA repair, Zinc finger and BTB domain Mixed multiclass only). Specifically, Egr1 levels were in- creased at 0.1 ppb, but decreased at the higher doses. In contrast, Hsp70 and Zfp145 levels were increased only above 1 ppb. In Table 1, the ‘‘Down 1, 50 ppb’’ section shows transcripts that were most strongly decreased at the higher Table 1 includes genesFunctional that annotation were from selected the as Mouse statistically Genome significant by eitherDenotes the the number SAM of multiclass probes for analysis each or gene that by were the selected two-class as statistically SAM significant analysis (FDR 0.05) for by at either the lea two-class SAM analysis (control vs. Increased expression with arsenic exposure designated ‘‘ doses, including Coro1a, Dbp, Plk2, and Ptprc. Likewise, the a b c d AI467657 Egr1 Hspa1b, Hsp70 Zfp145, Zfbt16 Symbol Gene name Function the multiclass (dose-response trend control, 0.01, 1, and 50 ppb). Mixed multiclass only immune response transcripts (Igh-VJ558, Igj, Igk, Igk-V28, DRINKING WATER ARSENIC MICROARRAY 81 Downloaded from https://academic.oup.com/toxsci/article/100/1/75/1624780 by guest on 29 September 2021

FIG. 1. Hierarchical clustering of transcripts significantly modified by arsenic. Transcripts were selected by SAM multiclass analysis using an FDR threshold of 0.05. Each block represents the expression level of a probe set in a single animal (green and red coloring indicate decreased and increased expression, respectively). Arrays are grouped by arsenic dose: control, 0.1, 1, and 50 ppb. Transcripts were clustered by a hierarchical clustering using the complete linkage algorithm and Pearson correlation metric. 82 ANDREW ET AL. Downloaded from https://academic.oup.com/toxsci/article/100/1/75/1624780 by guest on 29 September 2021

FIG. 2. Consistency between expression changes assessed by microarray and RT–PCR. Microarray-derived expression data are represented by blue solid lines/squares, and RT–PCR results are shown as black dashed lines/triangles. The graphs show the mean log2 ratio of arsenic exposed to unexposed group ± SD (* indicates p < 0.05 for RT–PCR by one-way ANOVA with Newman-Keuls posttest).

Ltb) were decreased almost exclusively with exposure to 50 arsenic treatment by implementing the DAVID search engine. ppb arsenic (Table 1, Down 50 ppb). We identified biologic functional processes that are over- Since the number of modifications that we could validate by represented by the significantly modified genes from the two- RT–PCR was limited by the quantity of sample remaining, we class analysis for each arsenic dose. Expression of transcripts chose to validate changes in genes with known biologic func- involved in response to abiotic stimulus was decreased at the tions that we hypothesized could be involved in arsenic lower doses (Table 2, 0.1 ppb row 6, 1 ppb row 3) and pathogenesis in the lung. These functions included oxygen increased at 50 ppb (row 5). The 0.1 and 1 ppb doses both transport and control of cell proliferation/apoptosis. The vali- decreased expression of hydrolases (Table 2, 0.1 ppb row 3, 1 dation graphs in Figure 2 show that for most of the transcripts, ppb row 1). Apoptosis and immune response were notable the expression profile derived from the array (blue line) closely categories modified by the 50 ppb dose, with potential impli- mirrors the expression observed by real-time PCR (black line). cations for disease risk (Table 2, 50 ppb rows 1, 3, 4). Analysis For example, the decrease in aminolevulinic acid synthase 2 of the biologic roles of all statistically significantly modified (Alas2) expression validated by real-time PCR was observed at transcripts (Table 3) indicated that arsenic modified a number all arsenic doses (Fig. 2). Likewise, decreases in expression of of transcripts that encode receptors (row 1), ligands (row 3), the oxygen transporter, hemoglobin alpha, adult chain 1 (Hba-a1) metabolic enzymes (row 11), transcription factors (row 13), were also validated by real-time PCR. and transporters (row 16). We further characterized the functional processes repre- The pathway map in Figure 3 was generated by Pathway sented by the transcripts that were significantly modified by Studio and illustrates direct interactions between arsenic-modified DRINKING WATER ARSENIC MICROARRAY 83

TABLE 2 Significantly Overrepresented Functional Processes of Arsenic-Modified Genesa

System Term Trend (þ/–)b Count Percent p Value

0.1 ppb GOTERM_MF_ALL Endoribonuclease activity 4 4.7 0.001 GOTERM_CC_ALL Lytic vacuole 4 4.7 0.015 SP_PIR_KEYWORDS Hydrolase 12 14.1 0.003 SP_PIR_KEYWORDS Signal 17 20.0 0.002 GOTERM_BP_ALL Lipid metabolism 8 9.4 0.007 GOTERM_BP_ALL Response to abiotic 7 8.2 0.016

stimulus Downloaded from https://academic.oup.com/toxsci/article/100/1/75/1624780 by guest on 29 September 2021 GOTERM_MF_ALL Translation initiation 3 3.5 0.039 factor activity 1 ppb SP_PIR_KEYWORDS Hydrolase 6 22.2 0.009 SP_PIR_KEYWORDS Phosphorylation 7 25.9 0.005 GOTERM_BP_ALL Response to abiotic 5 18.5 0.002 stimulus SP_PIR_KEYWORDS Transferase 5 18.5 0.044 50 ppb INTERPRO_NAME IPR007110: 6 17.1 0.006 Immunoglobulin-like GOTERM_MF_ALL ATP binding 6 17.1 0.031 GOTERM_BP_ALL Immune response 5 14.3 0.042 GOTERM_BP_ALL Apoptosis þ 3 25.0 0.007 GOTERM_BP_ALL Response to abiotic þ 3 25.0 0.014 stimulus

Note. GOTERM, Gene Ontology database; BP, biologic process; MF, molecular function; INTERPRO, European Bioinformatics Institute protein database; SP_PIR, Protein Information Resource. aAnalysis included all two-class SAM statistically modified transcripts. bIncreased expression with arsenic exposure designated ‘‘þ’’ and decreased expression designated ‘‘.’’

transcripts by each dose compared with control and identifies TABLE 3 common regulators. For example, Figure 3 shows that arsenic Biologic Role of Arsenic-Modified Genes decreases expression of the Cyclin D1/CCND1 regulator HBP1 at all three doses. This pathway map also shows common regulators, Arsenic-modified Total for example, PPARA modifies the expression of three arsenic Biologic role entities (n) entities (n) responsive genes—NR1D1, CPT1A,andCd36. Other central Nuclear receptors 3 58 nodes include MAPK1, SP1, and TP53. Receptors 11 1043 Histologic slices of formalin-fixed, paraffin-embedded lung Ligands 6 538 tissue from the mice were examined to look for structural Acetylases 1 39 differences between the arsenic unexposed versus exposed Ph 4 439 Glycosyltransferases 2 199 lungs (Supplementary Data). Pathologic examination did not Gtpase regulators 1 99 reveal any visible abnormalities in the structure of the lungs. Extracellular proteins 2 237 The pathologist did not detect differences between lung sam- Phosphatases 3 412 ples in the types of inflammatory cells present. Particularly, she Gtp binding proteins 2 273 did not observe any infiltration of neutrophils with arsenic Metabolic enzymes 10 1587 Kinases 4 874 exposure, indicating that the expression modification occurred Transcription factors 12 2133 within a static set of cells. Proteases 2 563 Lastly, we used immunoblot analysis to assess the protein Ubiquitin ligases 1 365 levels associated with several arsenic-modified transcripts. As Transporters 5 1163 shown in Figure 4, mice exposed to arsenic had increased Gpcr 2 2813 levels of Nr4a1/Nur77 and Hsp70 compared with unexposed Note. Analysis included all multiclass and two-class SAM statistically animals, while Ahr and cyclin D1 levels were reduced in modified transcripts at any exposure. arsenic-exposed animals. 84 ANDREW ET AL. Downloaded from https://academic.oup.com/toxsci/article/100/1/75/1624780 by guest on 29 September 2021

FIG. 3. Direct pathway relations between arsenic-modified genes. Genes are represented by shapes; those increased with arsenic exposure are colored red, those with decreased expression in response to arsenic are colored green (ratio: control vs. 50 ppb arsenic, see Supplementary Data for 0.1 and 1 ppb). Other genes in the pathway that were not significantly modified by arsenic are colored gray. Relationship lines were added by Pathway Studio between symbols based on evidence from the literature: depicted by colored lines and arrows in Figure 3 (red, positive regulation; green, negative regulation; gray, not specified).

DISCUSSION a number of important signaling pathways, many with corrob- orating evidence of arsenic responsiveness. Drinking water arsenic exposure is associated with increased Many of the modified transcripts showed consistently risk of lung diseases, particularly at levels exceeding 100 ppb decreased expression in arsenic exposed animals, while the found internationally (Smith and Smith, 2004). The effects of few increases were strongest at the higher doses. The immune lower doses and the biologic mechanisms remain unclear. In response transcripts showed clear evidence of a dose-re- particular, few in vivo studies have examined the biologic sponsive decrease at 1 and 50 ppb, but were unaffected, or even effects of drinking water arsenic at levels 50 ppb and below increased, at the 0.1 ppb dose. We also observed effects on that are highly relevant to U.S. drinking water levels. To several transcripts (e.g., S100a9, Ms4a8a, Angptl4, Cxcl7, address the need for data, our study evaluated differential Ear1, Ear2) that appear to be strongest at the lower doses of expression patterns in the lungs of mice exposed to drinking arsenic. In some cases, the responses even suggest a biphasic water arsenic (0.1, 1, 50 ppb) for a period of 5 weeks using effect that could be interpreted as being beneficial at 0.1 ppb microarrays. The expression profiles revealed modification of and adverse at higher levels (e.g., Egr1, Hsp70). DRINKING WATER ARSENIC MICROARRAY 85

thioesterases (Cte1), enzymes that regulate the levels of free fatty acid and coenzyme A by catalyzing the hydrolysis of acyl- CoAs (Hunt et al., 2006). Likewise, expression of carnitine palmitoyltransferase 1A (Cpt1a), the enzyme that catalyzes the primary rate controlling step of fatty acid oxidation, was decreased in arsenic-exposed animals (Cook and Park, 1999). Changes in the oxidative modification of lipoproteins may be related to the pathogenesis of peripheral artery disease (Steinberg et al., 1989). Arsenic-exposed animals had decreased expression of a number of transcripts involved in oxygen transport (Alas2, Downloaded from https://academic.oup.com/toxsci/article/100/1/75/1624780 by guest on 29 September 2021 Hbb-b1, Hba-a1, Bpgm). The effects on multiple hemoglobin complex and oxygen transporter transcripts were particularly FIG. 4. Arsenic exposure is associated with modified protein levels in notable at the lower doses. Alas2 catalyzes the first step in the mouse lung. Immunoblot analysis using antibodies specific to Nr4a1, Ahr, heme biosynthetic pathway (Kramer et al., 2000). The hemo- Cyclin D1, Hsp70, and Actin was performed on protein extracted from globin complex members Hbb-b1 and Hba-a1 actually bind and unexposed and exposed (10 ppb arsenic) animals (n ¼ 2). transport oxygen. Bpgm is an enzyme that controls the levels of an allosteric effector of hemoglobin and the dissociation of oxygen (Garel et al., 1990). We observed significant modification of pathways that We also observed differential expression of several apoptotic involve the central regulators MAPK1, CCND1, SP1, and transcripts on the microarray, including decreases in Ahr, TP53 (Fig. 3). Previous studies have reported that arsenic Stk17B, Nr4a1, Egr1, and Angptl4 and increases in Hspa1b/ activates MAPK1/ERK signaling in lung cells (Samet et al., Hsp70 and ZFP145/ZBTB16 (above 1 ppb). Nr4a1 is a nuclear 1998). MAPK1/ERK regulates Egr1 expression in response to orphan steroid receptor that induces apoptosis (Rajpal et al., a number of stimuli (Guha et al., 2001). As observed 2003). Interestingly, the protein levels of Nr4a1 were higher in previously in the bladder, arsenic had dose-dependent effects the arsenic-exposed mice, indicating a possible role for post- on Egr1 (increased at 0.1 ppb, decreased at 1 and 50 ppb), translational modification. We also observed increased levels possibly via the epidermal growth factor receptor, MAPK1/ of Hspa1b/Hsp70 with arsenic exposure levels 1 ppb and ERK, and SP1 signaling pathways (Adnane et al., 1999; Luster higher, as observed previously in the lymphocytes of arsenic- and Simeonova, 2004; Simeonova et al., 2000). We observed exposed individuals in Bangladesh (Argos et al., 2006). the same effect of arsenic on Egr1 expression in vitro by RT– Immunoblot analysis confirmed increased protein levels of PCR following treatment of cultured epithelial cells (data not Hsp70 with arsenic exposure. Hsp70 is induced in response shown). Egr1 is a zinc finger transcription factor that regulates to stress and is important for DNA repair and maintaining many processes, including differentiation, tumor suppression, genomic stability (Hunt et al., 2004). cell cycle arrest, and apoptosis, and its expression level is The most dramatic effect of arsenic exposure in this related to lung cancer survival (Ferraro et al., 2005). experiment was a decrease in transcripts involved in immune MAPK1/ERK signaling also plays a role in regulation of the response, including Igh-VJ558, Igj, Igk-V28, Igk-V8, Cd79b, cell cycle control protein cyclin D1 (CCND1) (Du et al., 2006). Cxcl7, Ian6, s100a9, and Ahr (see Fig. 1, cluster 4). Igh-VJ558 The decreased expression of the transcriptional repressor HBP1 is involved in B-cell antigen recognition and Igj, Igk-V28, and that we observed with arsenic is consistent with previous Igk-V8 in humoral immune response and antigen binding. reports that arsenic alters cyclin D1 levels (Sampson et al., Decreases in these transcripts were evident mainly at 50 ppb 2001; Vogt and Rossman, 2001). Arsenic exposure increased (Table 1, Down 50 ppb). Expression of the precursor plasma Hspa8/Hsc70, which binds to cyclin D1, as part of the chemokine Cxcl7 that may affect tumor development by catalytically active cyclin D1-cdk4 complex (Diehl et al., 2003). attracting immunocompetent cells was also decreased in the Although we did not observe statistically significant changes in arsenic-exposed groups, particularly at 0.1 and 1 ppb (Van CCND1 transcript levels by microarray, our immunoblot results Damme et al., 2004). Similarly, the chemokines Cxcl2, Ccl3, show decreased cyclin D1 protein levels in the lungs of arsenic- Ccrl2, Cxcl3, Ccl4, and Ccl20 were decreased in the exposed animals. lymphocytes of arsenic-exposed individuals from Bangladesh The decreased expression of transcripts involved in lipid (Argos et al., 2006). We also observed decreased expression of metabolic processes was strongest in the animals exposed to Calgranulin B/s100a9 with arsenic exposure, particularly at the the lowest dose, 0.1 ppb. In addition to other roles, Angptl4 lowest dose, 0.1 ppb. Calgranulin B is involved in the inhibits lipoprotein lipase, an enzyme that regulates triglyceride inflammatory response to stimuli, including the release of clearance and lipid homeostasis (Koster et al., 2005). This low neutrophils and is regulated by Ahr signaling (Temchura et al., dose of arsenic was also associated with decreased Acyl-CoA 2005). Ahr, which was also decreased at the transcript and 86 ANDREW ET AL. protein level in arsenic-exposed mice, is involved in the REFERENCES generation of regulatory T cells (Funatake et al., 2005). SP1 cooperatively regulates expression of AHR, which modifies Adnane, J., Shao, Z., and Robbins, P. D. (1999). Cyclin D1 associates with the TBP-associated factor TAF(II)250 to regulate Sp1-mediated transcription. cell cycle progression. The direction of this effect is dependent Oncogene 18, 239–247. on the presence or absence of exogenous ligands, such as Ahsan, H., and Thomas, D. C. (2004). Lung cancer etiology: Independent and Polycyclic Aromatic Hydrocarbons PAHs (Marlowe and Puga, joint effects of genetics, tobacco, and arsenic. JAMA 292, 3026–3029. 2005; Wang et al., 1999). Andrew, A. S., Warren, A. J., Barchowsky, A., Temple, K. A., Klei, L., We also observed differential expression of several apriori Soucy, N. V., O’Hara, K. A., and Hamilton, J. W. (2003). Genomic and hypothesized transcripts based on previous work. Induction of proteomic profiling of responses to toxic metals in human lung cells. metallothionein by other toxic metals including cadmium, mer- Environ. Health Perspect. 111, 825–835. cury, and copper is well documented; however, its response to Argos, M., Kibriya, M. G., Parvez, F., Jasmine, F., Rakibuz-Zaman, M., and Ahsan, H. (2006). Gene expression profiles in peripheral lymphocytes by Downloaded from https://academic.oup.com/toxsci/article/100/1/75/1624780 by guest on 29 September 2021 arsenic is less clear. Consistent with previous reports in other arsenic exposure and skin lesion status in a Bangladeshi population. Cancer organs and cells, the classic metal ion binding transcripts Mt1 and Epidemiol. Biomarkers Prev. 15, 1367–1375. Mt2 were induced by low doses of drinking water arsenic exposure Barchowsky, A., Roussel, R. R., Klei, L. R., James, P. E., Ganju, N., in the mouse lung. In utero exposure to arsenic also induces Smith, K. R., and Dudek, E. J. (1999). Low levels of arsenic trioxide metallothionein levels in the lung (Shen et al., 2007). Mt knockout stimulate proliferative signals in primary vascular cells without activating studies suggest that metallothionein protects mice against arsenic- stress effector pathways. Toxicol. Appl. Pharmacol. 159, 65–75. induced toxicity (Liu et al., 2000; Zheng et al., 2003). Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate—A practical and powerful approach to multiple testing. J. Roy. Stat. In summary, this study demonstrates that exposure to Soc. B Met. 57, 289–300. drinking water arsenic at levels commonly found in U.S. Chen, C. L., Hsu, L. I., Chiou, H. Y., Hsueh, Y. M., Chen, S. Y., Wu, M. M., drinking water is associated with altered expression profiles in and Chen, C. J. (2004). Ingested arsenic, cigarette smoking, and lung cancer the mouse lung. Pathways involved in immune response, angio- risk: A follow-up study in arseniasis-endemic areas in Taiwan. JAMA 292, genesis, lipid metabolism, and oxygen transport were partic- 2984–2990. ularly affected. With further validation, some of these genes Cook, G. A., and Park, E. A. (1999). Expression and regulation of carnitine may be useful as intermediate biomarkers in studies in- palmitoyltransferase-Ialpha and -Ibeta genes. Am. J. Med. Sci. 318, 43–48. vestigating arsenic-associated diseases. Future investigation of De, B. K., Majumdar, D., Sen, S., Guru, S., and Kundu, S. (2004). Pulmonary the effects of arsenic on the function of these pathways is involvement in chronic arsenic poisoning from drinking contaminated ground-water. J. Assoc. Physicians India 52, 395–400. warranted. In particular, little is known about the responses and Dennis, G., Jr, Sherman, B. T., Hosack, D. A., Yang, J., Gao, W., Lane, H. C., health effects of arsenic exposure at doses below 0.1 ppb, and and Lempicki, R. A. (2003). DAVID: Database for annotation, visualization, the possibility that arsenic may be an essential element at some and integrated discovery. Genome Biol. 4, P3. very low trace level remains to be explored. These data identify Diehl, J. A., Yang, W., Rimerman, R. A., Xiao, H., and Emili, A. (2003). pathways that will help guide investigations into mechanisms Hsc70 regulates accumulation of cyclin D1 and cyclin D1-dependent protein of arsenic’s health effects and help clarify the threshold for kinase. Mol. Cell. Biol. 23, 1764–1774. biologic effects and potential disease risk. Du, H. J., Tang, N., Liu, B. C., You, B. R., Shen, F. H., Ye, M., Gao, A., and Huang, C. (2006). Benzo[a]pyrene-induced cell cycle progression is through ERKs/cyclin D1 pathway and requires the activation of JNKs and p38 mapk in human diploid lung fibroblasts. Mol. Cell. Biochem. 287, 79–89. FUNDING Durham, T. R., and Snow, E. T. (2006). Metal ions and carcinogenesis. EXS 96, The National Cancer Institute, National Institutes of Health 97–130. (NIH); the National Institute of Environmental Health Eppig, J. T., Bult, C. J., Kadin, J. A., Richardson, J. E., Blake, J. A., Anagnostopoulos, A., Baldarelli, R. M., Baya, M., Beal, J. S., Bello, S. M., Sciences, NIH; the National Center for Research Resources et al. (2005). The mouse genome database (MGD): From genes to mice—A (NCRR), NIH (grant numbers CA099500, CA102327, P42 community resource for mouse biology. Nucleic Acids Res. 33, D471–D4475. ES007373, P20RR018787). Evans, C. D., LaDow, K., Schumann, B. L., Savage, R. E., Jr., Caruso, J., Vonderheide, A., Succop, P., and Talaska, G. (2004). Effect of arsenic on benzo[a]pyrene DNA adduct levels in mouse skin and lung. Carcinogenesis SUPPLEMENTARY DATA (Oxf) 25, 493–497. Ferraro, B., Bepler, G., Sharma, S., Cantor, A., and Haura, E. B. (2005). EGR1 Supplementary Table and other data are available online at predicts PTEN and survival in patients with non-small-cell lung cancer. http://toxsci.oxfordjournals.org/. J. Clin. Oncol. 23, 1921–1926. Funatake, C. J., Marshall, N. B., Steppan, L. B., Mourich, D. V., and Kerkvliet, N. I. (2005). Cutting edge: Activation of the aryl hydrocarbon ACKNOWLEDGMENTS receptor by 2,3,7,8-tetrachlorodibenzo-p-dioxin generates a population of CD4þ CD25þ cells with characteristics of regulatory T cells. J. Immunol. Its contents are solely the responsibility of the authors and do not necessarily 175, 4184–4188. represent the official views of the NIEHS or NIH. The authors do not have any Garel, M. C., Lemarchandel, V., Prehu, M. O., Calvin, M. C., Arous, N., competing financial interests. Rosa, R., Rosa, J., and Cohen-Solal, M. (1990). Natural and artificial mutants DRINKING WATER ARSENIC MICROARRAY 87

of the human 2,3-bisphosphoglycerate as a tool for the evaluation of activation of known and novel apoptotic pathways by Nur77 orphan steroid structure-function relationships. Biomed. Biochim. Acta 49, S166–S1171. receptor. EMBO J. 22, 6526–6536. Ghosh, P., Banerjee, M., De Chaudhuri, S., Chowdhury, R., Das, J. K., Rossman, T. G., Uddin, A. N., Burns, F. J., and Bosland, M. C. (2002). Mukherjee, A., Sarkar, A. K., Mondal, L., Baidya, K., Sau, T. J., et al. Arsenite cocarcinogenesis: An animal model derived from genetic (2007). Comparison of health effects between individuals with and without toxicology studies. Environ. Health Perspect. 110((Suppl. 5)), 749–752. skin lesions in the population exposed to arsenic through drinking water in Samet, J. M., Graves, L. M., Quay, J., Dailey, L. A., Devlin, R. B., Ghio, A. J., West Bengal, India. J. Expo. Sci. Environ. Epidemiol. 17, 215–223. Wu, W., Bromberg, P. A., and Reed, W. (1998). Activation of MAPKs in Guha, M., O’Connell, M. A., Pawlinski, R., Hollis, A., McGovern, P., human bronchial epithelial cells exposed to metals 275, L551–L558. Yan, S. F., Stern, D., and Mackman, N. (2001). Lipopolysaccharide Sampson, E. M., Haque, Z. K., Ku, M. C., Tevosian, S. G., Albanese, C., activation of the MEK-ERK1/2 pathway in human monocytic cells mediates Pestell, R. G., Paulson, K. E., and Yee, A. S. (2001). Negative regulation of tissue factor and tumor necrosis factor alpha expression by inducing Elk-1 the Wnt-beta-catenin pathway by the transcriptional repressor HBP1. EMBO phosphorylation and Egr-1 expression. Blood 98, 1429–1439. J. 20, 4500–4511. Hunt, C. R., Dix, D. J., Sharma, G. G., Pandita, R. K., Gupta, A., Funk, M., and Shen, J., Liu, J., Xie, Y., Diwan, B. A., and Waalkes, M. P. (2007). Fetal onset Downloaded from https://academic.oup.com/toxsci/article/100/1/75/1624780 by guest on 29 September 2021 Pandita, T. K. (2004). Genomic instability and enhanced radiosensitivity in of aberrant gene expression relevant to pulmonary carcinogenesis in lung Hsp70.1- and Hsp70.3-deficient mice. Mol. Cell. Biol. 24, 899–911. adenocarcinoma development induced by in utero arsenic exposure. Toxicol. Hunt, M. C., Rautanen, A., Westin, M. A., Svensson, L. T., and Alexson, S. E. Sci. 95, 313–320. (2006). Analysis of the mouse and human acyl-CoA thioesterase (ACOT) Shi, H., Shi, X., and Liu, K. J. (2004). Oxidative mechanism of arsenic toxicity gene clusters shows that convergent, functional evolution results in a reduced and carcinogenesis. Mol. Cell. Biochem. 255, 67–78. number of human peroxisomal ACOTs. FASEB J. 20, 1855–1864. Simeonova, P. P., Wang, S., Toriuma, W., Kommineni, V., Matheson, J., IARC. (2004). Some Drinking-water Disinfectants and Contaminants, Includ- Unimye, N., Kayama, F., Harki, D., Ding, M., Vallyathan, V., et al. (2000). ing Arsenic. International Agency for Research on Cancer (IARC) Monographs Arsenic mediates cell proliferation and gene expression in the bladder on the Evaluation of Carcinogenic Risks to Humans, pp. 39–267. Lyon France. epithelium: Association with activating protein-1 transactivation. Cancer Irizarry, R. A., Hobbs, B., Collin, F., Beazer-Barclay, Y. D., Antonellis, K. J., Res. 60, 3445–3453. Scherf, U., and Speed, T. P. (2003). Exploration, normalization, and Smith, A. H., Marshall, G., Yuan, Y., Ferreccio, C., Liaw, J., von summaries of high density oligonucleotide array probe level data. Bio- Ehrenstein, O., Steinmaus, C., Bates, M. N., and Selvin, S. (2006). Increased statistics 4, 249–264. mortality from lung cancer and bronchiectasis in young adults after exposure Kitchin, K. T. (2001). Recent advances in arsenic carcinogenesis: Modes of to arsenic in utero and in early childhood. Environ. Health Perspect. 114, action, animal model systems, and methylated arsenic metabolites. Toxicol. 1293–1296. Appl. Pharmacol. 172, 249–261. Smith, A. H., and Smith, M. M. (2004). Arsenic drinking water regulations in Koster, A., Chao, Y. B., Mosior, M., Ford, A., Gonzalez-Dewhitt, P. A., developing countries with extensive exposure. Toxicology 198, 39–44. Hale, J. E., Li, D., Qiu, Y., Fraser, C. C., Yang, D. D., et al. (2005). Transgenic Angptl4 overexpression and targeted disruption of Angptl4 and Soucy, N. V., Ihnat, M. A., Kamat, C. D., Hess, L., Post, M. J., Klei, L. R., Angptl3: Regulation of triglyceride metabolism. Endocrinology 146, 4943– Clark, C., and Barchowsky, A. (2003). Arsenic stimulates angiogenesis and 4950. tumorigenesis in vivo. Toxicol. Sci. 76, 271–279. Kramer, M. F., Gunaratne, P., and Ferreira, G. C. (2000). Transcriptional Steinberg, D., Carew, T. E., Fielding, C., Fogelman, A. M., Mahley, R. W., regulation of the murine erythroid-specific 5-aminolevulinate synthase gene. Sniderman, A. D., and Zilversmit, D. B. (1989). Lipoproteins and the Gene 247, 153–166. pathogenesis of atherosclerosis. Circulation 80, 719–723. Kurttio, P., Pukkala, E., Kahelin, H., Auvinen, A., and Pekkanen, J. (1999). Temchura, V. V., Frericks, M., Nacken, W., and Esser, C. (2005). Role of the Arsenic concentrations in well water and risk of bladder and kidney cancer aryl hydrocarbon receptor in thymocyte emigration in vivo. Eur. J. Immunol. in Finland. 107, 705–710. 35, 2738–2747. Lau, A. T., Li, M., Xie, R., He, Q. Y., and Chiu, J. F. (2004). Opposed arsenite- Tusher, V. G., Tibshirani, R., and Chu, G. (2001). Significance analysis of induced signaling pathways promote cell proliferation or apoptosis in microarrays applied to the ionizing radiation response. Proc. Natl. Acad. Sci. cultured lung cells. Carcinogenesis (Oxf) 25, 21–28. USA 98, 5116–5121. Liu, J., Liu, Y., Goyer, R. A., Achanzar, W., and Waalkes, M. P. (2000). Van Damme, J., Struyf, S., and Opdenakker, G. (2004). Chemokine-protease Metallothionein-I/II null mice are more sensitive than wild-type mice to the interactions in cancer. Semin. Cancer Biol. 14, 201–208. hepatotoxic and nephrotoxic effects of chronic oral or injected inorganic Vogt, B. L., and Rossman, T. G. (2001). Effects of arsenite on p53, p21 and arsenicals. Toxicol. Sci. 55, 460–467. cyclin D expression in normal human fibroblasts—A possible mechanism for Luster, M. I., and Simeonova, P. P. (2004). Arsenic and urinary bladder cell arsenite’s comutagenicity. Mutat. Res. 478, 159–168. proliferation. Toxicol. Appl. Pharmacol. 198, 419–423. Wang, F., Wang, W., and Safe, S. (1999). Regulation of constitutive gene Marlowe, J. L., and Puga, A. (2005). Aryl hydrocarbon receptor, cell cycle expression through interactions of Sp1 protein with the nuclear aryl regulation, toxicity, and tumorigenesis. J. Cell. Biochem. 96, 1174–1184. hydrocarbon receptor complex. Biochemistry 38, 11490–11500. Rajpal, A., Cho, Y. A., Yelent, B., Koza-Taylor, P. H., Li, D., Chen, E., Zheng, X. H., Watts, G. S., Vaught, S., and Gandolfi, A. J. (2003). Low-level Whang, M., Kang, C., Turi, T. G., and Winoto, A. (2003). Transcriptional arsenite induced gene expression in HEK293 cells. Toxicology 187, 39–48.