Transcriptional ontogeny of the developing liver.

BACKGROUND: During embryogenesis the liver is derived from endodermal cells lining the digestive tract. These endodermal progenitor cells contribute to forming the parenchyma of a number of organs including the liver and pancreas. Early in organogenesis the fetal liver is populated by hematopoietic stem cells, the source for a number of blood cells including nucleated erythrocytes. A comprehensive analysis of the transcriptional changes that occur during the early stages of development to adulthood in the liver was carried out. RESULTS: We characterized gene expression changes in the developing mouse liver at gestational days (GD) 11.5, 12.5, 13.5, 14.5, 16.5, and 19 and in the neonate (postnatal day (PND) 7 and 32) compared to that in the adult liver (PND67) using full-genome microarrays. The fetal liver, and to a lesser extent the neonatal liver, exhibited dramatic differences in gene expression compared to adults. Canonical pathway analysis of the fetal liver signature demonstrated increases in functions important in cell replication and DNA fidelity whereas most metabolic pathways of intermediary metabolism were under expressed. Comparison of the dataset to a number of previously published microarray datasets revealed 1) a striking similarity between the fetal liver and that of the pancreas in both mice and humans, 2) a nucleated erythrocyte signature in the fetus and 3) under expression of most xenobiotic metabolism genes throughout development, with the exception of a number of transporters associated with either hematopoietic cells or cell proliferation in hepatocytes. CONCLUSIONS: Overall, these findings reveal the complexity of gene expression changes during liver development and maturation, and provide a foundation to predict responses to chemical and drug exposure as a function of early life-stages.

Gene expression of normal human epidermal keratinocytes modulated by trivalent arsenicals.

Chronic exposure to inorganic arsenic (iAs) is associated with the development of benign and malignant human skin lesions including nonmelanoma skin cancers. The precise arsenical form(s) responsible for this carcinogenic effect are unknown, although trivalent inorganic arsenic (iAs(III)) and two of its toxic metabolites, monomethylarsonous acid (MMA(III)) and methylarsinous acid (DMA(III)), are attractive candidates. In an effort to better understand and compare their toxic effects in the skin, we compared the global gene expression profiles of normal human epidermal keratinocytes (NHEKs) exposed to varying noncytotoxic/slightly cytotoxic concentrations of iAs(III), MMA(III), and DMA(III) for 24 h. Exposure to each arsenical treatment group exhibited a dose effect in the number of altered genes and the magnitude of expression change in NHEKs. The most significant gene expression changes associated with iAs(III) and MMA(III) exposure were consistent with several key events believed to be important to As-driven skin carcinogenesis, namely induction of oxidative stress, increased transcript levels of keratinocyte growth factors, and modulation of MAPK and NF-κB pathways. At both comparable arsenical concentrations and comparable NHEK toxicity, greater potential carcinogenic effects were observed in MMA(III)-exposed NHEKs than those exposed to iAs(III), including involvement of more proinflammatory signals and increased transcript levels of more growth factor genes. In contrast, none of these above-mentioned transcriptional trends were among the most significantly altered functions in the DMA(III) treatment group. This study suggests the relative capacity of each of the tested arsenicals to drive suspected key events in As-mediated skin carcinogenesis is MMA(III) > iAs(III) with little contribution from DMA(III).

Splice variant specific increase in Ca2+/calmodulin-dependent protein kinase 1-gamma mRNA expression in response to acute pyrethroid exposure.

In mammals, pyrethroids are neurotoxicants that interfere with ion channel function in excitable neuronal membranes. Previous work demonstrated increases in the expression of Ca(2+)/calmodulin-dependent protein kinase 1-gamma (Camk1g) mRNA following acute deltamethrin and permethrin exposure. In the rat, this gene is expressed as two distinct splice variants, Camk1g1 and Camk1g2. The present study tests the hypothesis that changes in Camk1g mRNA expression in the rat following acute pyrethroid exposure are due to a specific increase in the Camk1g1 splice variant and not the Camk1g2 splice variant. Long-Evans rats were acutely exposed to permethrin, deltamethrin, or corn oil vehicle. Frontal cortex was collected at 6 h postdosing. In addition, rats were exposed to permethrin (100 mg/kg) or deltamethrin (3 mg/kg), and frontal cortex was collected at 1, 3, 6, 9, 12, or 24 h along with time-matched vehicle controls. Expression of Camk1g1 and Camk1g2 mRNA was measured by quantitative real-time RT-PCR and quantified using the 2(-Delta Delta C)T method. Dose-dependent increases in Camk1g1 mRNA expression were observed for both pyrethroids at 6 h. In addition, a dose-dependent increase in Camk1g2 was observed at 6 h although it was very small in magnitude. The increases in Camk1g1 expression for deltamethrin and permethrin peak between 3 and 6 h postexposure and returns to control levels by 9 h. There was no increase in CAMK1G1 protein as measured with Western blots. The present data demonstrate that pyrethroid-induced changes in Camk1g are driven mainly by increased expression of the Camk1g1 splice variant.

K-Ras mutant fraction in A/J mouse lung increases as a function of benzo[a]pyrene dose.

K-Ras mutant fraction (MF) was measured to examine the default assumption of low-dose linearity in the benzo[a]pyrene (B[a]P) mutational response. Groups of 10 male A/J mice (7- to 9-weeks old) received a single i.p. injection of 0, 0.05, 0.5, 5, or 50 mg/kg B[a]P and were sacrificed 28 days after treatment. K-Ras codon 12 TGT and GAT MFs in lung DNAs were measured using Allele-specific Competitive Blocker-PCR (ACB-PCR). The K-Ras codon 12 TGT geometric mean MF was 3.88 x 10(-4) in controls, indicating an average of 1 mutation in every approximately 1,288 lung cells. The K-Ras codon 12 TGT geometric mean MFs were as follows: 3.56 x 10(-4); 6.19 x 10(-4); 2.02 x 10(-3), and 3.50 x 10(-3) for the 0.05, 0.5, 5, and 50 mg/kg B[a]P treatment groups, respectively. The 5 and 50 mg/kg dose groups had TGT MFs significantly higher than did controls. Although 10(-5) is considered as the limit of accurate ACB-PCR quantitation, K-Ras codon 12 GAT geometric mean MFs were as follows: 8.38 x 10(-7), 1.47 x 10(-6), 2.19 x 10(-6), 5.71 x 10(-6), and 8.99 x 10(-6) for the 0, 0.05, 0.5, 5, and 50 mg/kg B[a]P treatment groups, respectively. The K-Ras TGT and GAT MFs increased in a B[a]P-dose-dependent manner, with response approximately linear over the 0.05 to 5 mg/kg dose range. K-Ras MF increased with B[a]P adduct burden measured for identical doses in a separate study. Thus, ACB-PCR may be useful in characterizing the shape of a dose-response curve at low doses and establishing relationships between DNA adducts and tumor-associated mutations.

Transcriptional changes associated with reduced spontaneous liver tumor incidence in mice chronically exposed to high dose arsenic.

Exposure of male C3H mice in utero (from gestational days 8-18) to 85ppm sodium arsenite via the dams' drinking water has previously been shown to increase liver tumor incidence by 2 years of age. However, in our companion study (Ahlborn et al., 2009), continuous exposure to 85ppm sodium arsenic (from gestational day 8 to postnatal day 365) did not result in increased tumor incidence, but rather in a significant reduction (0% tumor incidence). The purpose of the present study was to examine the gene expression responses that may lead to the apparent protective effect of continuous arsenic exposure. Genes in many functional categories including cellular growth and proliferation, gene expression, cell death, oxidative stress, protein ubiquitination, and mitochondrial dysfunction were altered by continuous arsenic treatment. Many of these genes are known to be involved in liver cancer. One such gene associated with rodent hepatocarcinogenesis, Scd1, encodes stearoyl-CoA desaturase and was down-regulated by continuous arsenic treatment. An overlap between the genes in our study affected by continuous arsenic exposure and those from the literature affected by long-term caloric restriction suggests that reduction in the spontaneous tumor incidence under both conditions may involve similar gene pathways such as fatty acid metabolism, apoptosis, and stress response.

Global gene expression profiling of hyperkeratotic skin lesions from inner Mongolians chronically exposed to arsenic.

The skin is an organ that is highly sensitive to chronic arsenic (As) exposure. Skin lesions such as hyperkeratoses (HKs) are common early manifestations of arsenicosis in humans. HKs can be precursor lesions of nonmelanoma skin cancers (NMSCs), but the driving forces behind their formation and how they may ultimately progress to NMSCs are unknown. The goal of this study was to examine the global gene expression profiles of As-related HKs in an effort to better understand gene expression changes that are potentially associated with early stages of As carcinogenesis. HK biopsies were removed from individuals living in an arsenicosis-endemic region in Inner Mongolia who had been exposed to high As levels in their drinking water for >20 years. Gene expression profiling was performed on RNA isolated from 7 individuals in this group and from 4 lesion-free skin samples from healthy individuals. Consistent with the pathological characteristics of the HK lesions, major functional categories and known canonical pathways represented by altered transcripts include those involved in development, differentiation, apoptosis, proliferation, and stress response. The results of this study may help define a signature profile of gene expression changes associated with long-term As exposure in the skin.

Dimethylarsinic acid in drinking water changed the morphology of urinary bladder but not the expression of DNA repair genes of bladder transitional epithelium in F344 rats.

Inorganic arsenic increases urinary bladder transitional cell carcinoma in humans. In F344 rats, dimethylarsinic acid (DMA[V]) increases transitional cell carcinoma. Arsenic-induced inhibition of DNA repair has been reported in cultured cell lines and in lymphocytes of arsenic-exposed humans, but it has not been studied in urinary bladder. Should inhibition of DNA damage repair in transitional epithelium occur, it may contribute to carcinogenesis or cocarcinogenesis. We investigated morphology and expression of DNA repair genes in F344 rat transitional cells following up to 100 ppm DMA(V) in drinking water for four weeks. Mitochondria were very sensitive to DMA(V), and swollen mitochondria appeared to be the main source of vacuoles in the transitional epithelium. Real-time reverse transcriptase polymerase chain reaction (Real-Time RT PCR) showed the mRNA levels of tested DNA repair genes, ataxia telangectasia mutant (ATM), X-ray repair cross-complementing group 1 (XRCC1), excision repair cross-complementing group 3/xeroderma pigmentosum B (ERCC3/XPB), and DNA polymerase beta (Polbeta), were not altered by DMA(V). These data suggested that either DMA(V) does not affect DNA repair in the bladder or DMA(V) affects DNA repair without affecting baseline mRNA levels of repair genes. The possibility remains that DMA(V) may lower damage-induced increases in repair gene expression or cause post-translational modification of repair enzymes.

Dose response evaluation of gene expression profiles in the skin of K6/ODC mice exposed to sodium arsenite.

Chronic drinking water exposure to inorganic arsenic and its metabolites increases tumor frequency in the skin of K6/ODC transgenic mice. To identify potential biomarkers and modes of action for this skin tumorigenicity, we characterized gene expression profiles from analysis of K6/ODC mice administered 0, 0.05, 0.25, 1.0 and 10 ppm sodium arsenite in their drinking water for 4 weeks. Following exposure, total RNA was isolated from mouse skin and processed to biotin-labeled cRNA for microarray analyses. Skin gene expression was analyzed with Affymetrix Mouse Genome 430A 2.0 GeneChips, and pathway analysis was conducted with DAVID (NIH), Ingenuity Systems and MetaCore's GeneGo. Differential expression of several key genes was verified through qPCR. Only the highest dose (10 ppm) resulted in significantly altered KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways, including MAPK, regulation of actin cytoskeleton, Wnt, Jak-Stat, Tight junction, Toll-like, phosphatidylinositol and insulin signaling pathways. Approximately 20 genes exhibited a dose response, including several genes known to be associated with carcinogenesis or tumor progression including cyclin D1, CLIC4, Ephrin A1, STAT3 and DNA methyltransferase 3a. Although transcription changes in all identified genes have not previously been linked to arsenic carcinogenesis, their association with carcinogenesis in other systems suggests that these genes may play a role in the early stages of arsenic-induced skin carcinogenesis and can be considered potential biomarkers.

Toxicogenomic analysis incorporating operon-transcriptional coupling and toxicant concentration-expression response: analysis of MX-treated Salmonella.

BACKGROUND: Deficiencies in microarray technology cause unwanted variation in the hybridization signal, obscuring the true measurements of intracellular transcript levels. Here we describe a general method that can improve microarray analysis of toxicant-exposed cells that uses the intrinsic power of transcriptional coupling and toxicant concentration-expression response data. To illustrate this approach, we characterized changes in global gene expression induced in Salmonella typhimurium TA100 by 3-chloro-4-(dichloromethyl)-5-hydroxy-2(5H)-furanone (MX), the primary mutagen in chlorinated drinking water. We used the co-expression of genes within an operon and the monotonic increases or decreases in gene expression relative to increasing toxicant concentration to augment our identification of differentially expressed genes beyond Bayesian-t analysis. RESULTS: Operon analysis increased the number of altered genes by 95% from the list identified by a Bayesian t-test of control to the highest concentration of MX. Monotonic analysis added 46% more genes. A functional analysis of the resulting 448 differentially expressed genes yielded functional changes beyond what would be expected from only the mutagenic properties of MX. In addition to gene-expression changes in DNA-damage response, MX induced changes in expression of genes involved in membrane transport and porphyrin metabolism, among other biological processes. The disruption of porphyrin metabolism might be attributable to the structural similarity of MX, which is a chlorinated furanone, to ligands indigenous to the porphyrin metabolism pathway. Interestingly, our results indicate that the lexA regulon in Salmonella, which partially mediates the response to DNA damage, may contain only 60% of the genes present in this regulon in E. coli. In addition, nanH was found to be highly induced by MX and contains a putative lexA regulatory motif in its regulatory region, suggesting that it may be regulated by lexA. CONCLUSION: Operon and monotonic analyses improved the determination of differentially expressed genes beyond that of Bayesian-t analysis, showing that MX alters cellular metabolism involving pathways other than DNA damage. Because co-expression of similarly functioning genes also occurs in eukaryotes, this method has general applicability for improving analysis of toxicogenomic data.

Kidney toxicogenomics of chronic potassium bromate exposure in f344 male rats.

BACKGROUND: Potassium bromate (KBrO3), used in both the food and cosmetics industry, and a drinking water disinfection by-product, is a nephrotoxic compound and rodent carcinogen. To gain insight into the carcinogenic mechanism of action and provide possible biomarkers of KBrO3 exposure, the gene expression in kidneys from chronically exposed male F344 rats was investigated. METHODS: Male F344 rats were exposed to KBrO3 in drinking water for 52 and 100 wk. Kidneys were removed, frozen, and stored at -80°C, then used for Affymetrix microarray analysis. Gene expression patterns were examined using a non-carcinogenic (20 ppm) and carcinogenic dose (400 ppm) at 52 wk, and compared to 100 wk high dose (400 ppm) and adenoma gene expression. RESULTS: Statistical analysis revealed 144, 224, 43, and 994 genes out of 15866 from the 52 wk low, 52 wk high, 100 wk high, and adenomas respectively, were differentially expressed when compared to control kidneys. Gene ontology classification of the 52 wk high dose showed alterations of gene transcripts involved in oxidative stress, lipid metabolism, kidney function/ion transport, and cellular function. In a comparison of kidney development gene expression, alterations were seen in the adenomas but not in the 52 wk bromate-treated kidneys. However, the normal kidney from the high dose group resembled the adenoma expression pattern with early kidney development genes being up-regulated and adult phase genes being down-regulated. Moreover, eight genes were identified which could serve as biomarkers of carcinogenic exposure to bromate. The most promising of these was Pendrin, or Slc26a4, a solute carrier of chloride and iodide active in the kidney, thyroid, and inner ear. All these tissues are targets of KBrO3 toxicity. Expression array results were verified with quantitative real-time rtPCR. CONCLUSIONS: These data demonstrate that the 400 ppm carcinogenic dose of KBrO3 showed marked gene expression differences from the 20 ppm non-carcinogenic dose. Comparison of kidney development gene expression showed that the adenoma patterns were more characteristic of embryonic than adult kidneys, and that the normal kidney from the high dose group resembled the adenoma-like gene expression pattern. Taken together, the analysis from this study identifies potential biomarkers of exposure and illuminates a possible carcinogenic mode of action for KBrO3.

Quantitation of aberrant interlocus T-cell receptor rearrangements in mouse thymocytes and the effect of the herbicide 2,4-dichlorophenoxyacetic acid.

Small studies in human populations have suggested a correlation between the frequency of errors in antigen receptor gene assembly and lymphoid malignancy risk. In particular, agricultural workers exposed to pesticides have both an increased risk for lymphoma and an increased frequency of errors in antigen receptor gene assembly. In order to further investigate the potential of such errors to serve as a mechanistically based biomarker of lymphoid cancer risk, we have developed a sensitive PCR assay for quantifying errors of V(D)J recombination in the thymocytes of mice. This assay measures interlocus rearrangements between two T-cell receptor loci, V-gamma and J-beta, located on chromosomes 13 and 6, respectively. The baseline frequency in four strains of mice was determined at several ages (2-8 weeks of age) and was found to be stable at approximately 1.5 x 10(-5) per thymocyte. Strain AKR, which has a high susceptibility to T-cell lymphomas, did not show an elevated frequency of aberrant V(D)J events. We used this assay to examine the effects of the herbicide 2,4-dichlorophenoxyacetic acid (2,4-D) on the frequency of these events. Female B6C3F1 mice, 27 days of age, were exposed to 2,4-D by gavage at doses of 0, 3, 10, 30, and 100 mg/kg/day for 4 successive days and sacrificed on day 5. Thymus DNA was isolated and examined for illegitimate V(D)J recombination-mediated gene rearrangements. In addition, pregnant mice were exposed to 2,4-D and thymocytes from the offspring examined at 2 weeks of age. No significant increase in aberrant V(D)J rearrangements was found, indicating that under these conditions 2,4-D does not appear to effect this important mechanism of carcinogenesis.