Epigenome-wide association studies of occupational exposure to benzene and formaldehyde

Sufficient evidence supports a relationship between certain myeloid neoplasms and exposure to benzene or formaldehyde. DNA methylation could underlie benzene- and formaldehyde-induced health outcomes, but data in exposed human populations are limited. We conducted two cross-sectional epigenome-wide association studies (EWAS), one in workers exposed to benzene and another in workers exposed to formaldehyde. Using HumanMethylation450 BeadChips, we investigated differences in blood cell DNA methylation among 50 benzene-exposed subjects and 48 controls, and among 31 formaldehyde-exposed subjects and 40 controls. We performed CpG-level and regional-level analyses. In the benzene EWAS, we found genome-wide significant alterations, i.e., FWER-controlled P-values <0.05, in the mean and variance of methylation at 22 and 318 CpG sites, respectively, and in mean methylation of a large genomic region. Pathway analysis of genes corresponding to benzene-associated differential methylation sites revealed an impact on the AMPK signalling pathway. In formaldehyde-exposed subjects compared to controls, 9 CpGs in the DUSP22 gene promoter had genome-wide significant decreased methylation variability and a large region of the HOXA5 promoter with 44 CpGs was hypomethylated. Our findings suggest that DNA methylation may contribute to the pathogenesis of diseases related to benzene and formaldehyde exposure. Aberrant expression and methylation of HOXA5 previously has been shown to be clinically significant in myeloid leukaemias. The tumour suppressor gene DUSP22 is a potential biomarker of exposure to formaldehyde, and irregularities have been associated with multiple exposures and diseases.

Consensus on the Key Characteristics of Immunotoxic Agents as a Basis for Hazard Identification

BACKGROUND

Key characteristics (KCs), properties of agents or exposures that confer potential hazard, have been developed for carcinogens and other toxicant classes. KCs have been used in the systematic assessment of hazards and to identify assay and data gaps that limit screening and risk assessment. Many of the mechanisms through which pharmaceuticals and occupational or environmental agents modulate immune function are well recognized. Thus KCs could be identified for immunoactive substances and applied to improve hazard assessment of immunodulatory agents.

OBJECTIVES

The goal was to generate a consensus-based synthesis of scientific evidence describing the KCs of agents known to cause immunotoxicity and potential applications, such as assays to measure the KCs.

METHODS

A committee of 18 experts with diverse specialties identified 10 KCs of immunotoxic agents, namely, 1) covalently binds to proteins to form novel antigens, 2) affects antigen processing and presentation, 3) alters immune cell signaling, 4) alters immune cell proliferation, 5) modifies cellular differentiation, 6) alters immune cell-cell communication, 7) alters effector function of specific cell types, 8) alters immune cell trafficking, 9) alters cell death processes, and 10) breaks down immune tolerance. The group considered how these KCs could influence immune processes and contribute to hypersensitivity, inappropriate enhancement, immunosuppression, or autoimmunity.

DISCUSSION

KCs can be used to improve efforts to identify agents that cause immunotoxicity via one or more mechanisms, to develop better testing and biomarker approaches to evaluate immunotoxicity, and to enable a more comprehensive and mechanistic understanding of adverse effects of exposures on the immune system. https://doi.org/10.1289/EHP10800.

Key Characteristics of Cardiovascular Toxicants

BACKGROUND

The concept of chemical agents having properties that confer potential hazard called key characteristics (KCs) was first developed to identify carcinogenic hazards. Identification of KCs of cardiovascular (CV) toxicants could facilitate the systematic assessment of CV hazards and understanding of assay and data gaps associated with current approaches.

OBJECTIVES

We sought to develop a consensus-based synthesis of scientific evidence on the KCs of chemical and nonchemical agents known to cause CV toxicity along with methods to measure them.

METHODS

An expert working group was convened to discuss mechanisms associated with CV toxicity.

RESULTS

The group identified 12 KCs of CV toxicants, defined as exogenous agents that adversely interfere with function of the CV system. The KCs were organized into those primarily affecting cardiac tissue (numbers 1-4 below), the vascular system (5-7), or both (8-12), as follows: 1) impairs regulation of cardiac excitability, 2) impairs cardiac contractility and relaxation, 3) induces cardiomyocyte injury and death, 4) induces proliferation of valve stroma, 5) impacts endothelial and vascular function, 6) alters hemostasis, 7) causes dyslipidemia, 8) impairs mitochondrial function, 9) modifies autonomic nervous system activity, 10) induces oxidative stress, 11) causes inflammation, and 12) alters hormone signaling.

DISCUSSION

These 12 KCs can be used to help identify pharmaceuticals and environmental pollutants as CV toxicants, as well as to better understand the mechanistic underpinnings of their toxicity. For example, evidence exists that fine particulate matter [PM ≤2.5μm in aerodynamic diameter (PM2.5)] air pollution, arsenic, anthracycline drugs, and other exogenous chemicals possess one or more of the described KCs. In conclusion, the KCs could be used to identify potential CV toxicants and to define a set of test methods to evaluate CV toxicity in a more comprehensive and standardized manner than current approaches. https://doi.org/10.1289/EHP9321.

Applying genome-wide CRISPR to identify known and novel genes and pathways that modulate formaldehyde toxicity

Formaldehyde (FA), a ubiquitous environmental pollutant, is classified as a Group I human carcinogen by the International Agency for Research on Cancer. Previously, we reported that FA induced hematotoxicity and chromosomal aneuploidy in exposed workers and toxicity in bone marrow and hematopoietic stem cells of experimental animals. Using functional toxicogenomic profiling in yeast, we identified genes and cellular processes modulating eukaryotic FA cytotoxicity. Although we validated some of these findings in yeast, many specific genes, pathways and mechanisms of action of FA in human cells are not known. In the current study, we applied genome-wide, loss-of-function CRISPR screening to identify modulators of FA toxicity in the human hematopoietic K562 cell line. We assessed the cellular genetic determinants of susceptibility and resistance to FA at 40, 100 and 150 μM (IC10, IC20 and IC60, respectively) at two time points, day 8 and day 20. We identified multiple candidate genes that increase sensitivity (e.g. ADH5, ESD and FANC family) or resistance (e.g. FASN and KDM6A) to FA when disrupted. Pathway analysis revealed a major role for the FA metabolism and Fanconi anemia pathway in FA tolerance, consistent with findings from previous studies. Additional network analyses revealed potential new roles for one-carbon metabolism, fatty acid synthesis and mTOR signaling in modulating FA toxicity. Validation of these novel findings will further enhance our understanding of FA toxicity in human cells. Our findings support the utility of CRISPR-based functional genomics screening of environmental chemicals.

Using the Key Characteristics of Carcinogens to Develop Research on Chemical Mixtures and Cancer

BACKGROUND

People are exposed to numerous chemicals throughout their lifetimes. Many of these chemicals display one or more of the key characteristics of carcinogens or interact with processes described in the hallmarks of cancer. Therefore, evaluating the effects of chemical mixtures on cancer development is an important pursuit. Challenges involved in designing research studies to evaluate the joint action of chemicals on cancer risk include the time taken to perform the experiments because of the long latency and choosing an appropriate experimental design.

OBJECTIVES

The objectives of this work are to present the case for developing a research program on mixtures of environmental chemicals and cancer risk and describe recommended approaches.

METHODS

A working group comprising the coauthors focused attention on the design of mixtures studies to inform cancer risk assessment as part of a larger effort to refine the key characteristics of carcinogens and explore their application. Working group members reviewed the key characteristics of carcinogens, hallmarks of cancer, and mixtures research for other disease end points. The group discussed options for developing tractable projects to evaluate the joint effects of environmental chemicals on cancer development.

RESULTS AND DISCUSSION

Three approaches for developing a research program to evaluate the effects of mixtures on cancer development were proposed: a chemical screening approach, a transgenic model-based approach, and a disease-centered approach. Advantages and disadvantages of each are discussed. https://doi.org/10.1289/EHP8525.

Formaldehyde-induced hematopoietic stem and progenitor cell toxicity in mouse lung and nose

Formaldehyde (FA), an economically important and ubiquitous chemical, has been classified as a human carcinogen and myeloid leukemogen. However, the underlying mechanisms of leukemogenesis remain unclear. Unlike many classical leukemogens that damage hematopoietic stem/progenitor cells (HSC/HPC) directly in the bone marrow, FA-as the smallest, most reactive aldehyde-is thought to be incapable of reaching the bone marrow through inhalation exposure. A recent breakthrough study discovered that mouse lung contains functional HSC/HPC that can produce blood cells and travel bi-directionally between the lung and bone marrow, while another early study reported the presence of HSC/HPC in rat nose. Based on these findings, we hypothesized that FA inhalation could induce toxicity in HSC/HPC present in mouse lung and/or nose rather than in the bone marrow. To test this hypothesis, we adapted a commercially available protocol for culturing burst-forming unit-erythroid (BFU-E) and colony-forming unit-granulocyte, macrophage (CFU-GM) colonies from bone marrow and spleen to also enable culture of these colonies from mouse lung and nose, a novel application of this assay. We reported that in vivo exposure to FA at 3 mg/m3 or ex vivo exposure up to 400 µM FA decreased the formation of both colony types from mouse lung and nose as well as from bone marrow and spleen. These findings, to the best of our knowledge, are the first empirically to show that FA exposure can damage mouse pulmonary and olfactory HSC/HPC and provide potential biological plausibility for the induction of leukemia at the sites of entry rather than the bone marrow.

The Key Characteristics of Carcinogens: Relationship to the Hallmarks of Cancer, Relevant Biomarkers, and Assays to Measure Them

The key characteristics (KC) of human carcinogens provide a uniform approach to evaluating mechanistic evidence in cancer hazard identification. Refinements to the approach were requested by organizations and individuals applying the KCs. We assembled an expert committee with knowledge of carcinogenesis and experience in applying the KCs in cancer hazard identification. We leveraged this expertise and examined the literature to more clearly describe each KC, identify current and emerging assays and in vivo biomarkers that can be used to measure them, and make recommendations for future assay development. We found that the KCs are clearly distinct from the Hallmarks of Cancer, that interrelationships among the KCs can be leveraged to strengthen the KC approach (and an understanding of environmental carcinogenesis), and that the KC approach is applicable to the systematic evaluation of a broad range of potential cancer hazards in vivo and in vitro We identified gaps in coverage of the KCs by current assays. Future efforts should expand the breadth, specificity, and sensitivity of validated assays and biomarkers that can measure the 10 KCs. Refinement of the KC approach will enhance and accelerate carcinogen identification, a first step in cancer prevention.See all articles in this CEBP Focus section, "Environmental Carcinogenesis: Pathways to Prevention."

In utero exposure to endocrine disrupting chemicals, micro-RNA profiles, and fetal growth: a pilot study protocol

Background: The developing fetus is particularly vulnerable to the effects of endocrine disrupting chemicals (EDCs). Molecular fingerprints of EDCs can be identified via microRNA (miRNA) expression profiles and may be etiologically implicated in the developmental origin of disease (DOHaD).

Methods/design: This pilot study includes pregnant women at high risk (smoking at conception), and low risk (non-smoking at conception) for SGA birth (birthweight<10^th percentile for gestational age). We have randomly selected 12 mothers (3 high-risk SGA birth, 3 low-risk SGA birth, 3 high-risk non-SGA birth, 3 low-risk non-SGA birth), with EDC measurements from gestational week 17. All offspring are female. We aim to test the stability of our samples (maternal serum, cord blood, placenta tissue), observe the differential expression of miRNA profiles over time (gestational weeks 17, 25, 33, 37, birth), and study the consistency between maternal EDC measures and miRNA expression profiles across our repeated measures.

Expected impact of the study for Public Health: Results from this pilot study will inform the development of a larger cohort wide analysis, and will impact the current state of knowledge in the fields of public health, epigenetics, and the DOHaD.

Proposed Key Characteristics of Female Reproductive Toxicants as an Approach for Organizing and Evaluating Mechanistic Data in Hazard Assessment

BACKGROUND

Identification of female reproductive toxicants is currently based largely on integrated epidemiological and in vivo toxicology data and, to a lesser degree, on mechanistic data. A uniform approach to systematically search, organize, integrate, and evaluate mechanistic evidence of female reproductive toxicity from various data types is lacking.

OBJECTIVE

We sought to apply a key characteristics approach similar to that pioneered for carcinogen hazard identification to female reproductive toxicant hazard identification.

METHODS

A working group of international experts was convened to discuss mechanisms associated with chemical-induced female reproductive toxicity and identified 10 key characteristics of chemicals that cause female reproductive toxicity: 1) alters hormone receptor signaling; alters reproductive hormone production, secretion, or metabolism; 2) chemical or metabolite is genotoxic; 3) induces epigenetic alterations; 4) causes mitochondrial dysfunction; 5) induces oxidative stress; 6) alters immune function; 7) alters cell signal transduction; 8) alters direct cell–cell interactions; 9) alters survival, proliferation, cell death, or metabolic pathways; and 10) alters microtubules and associated structures. As proof of principle, cyclophosphamide and diethylstilbestrol (DES), for which both human and animal studies have demonstrated female reproductive toxicity, display at least 5 and 3 key characteristics, respectively. 2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD), for which the epidemiological evidence is mixed, exhibits 5 key characteristics.

DISCUSSION

Future efforts should focus on evaluating the proposed key characteristics against additional known and suspected female reproductive toxicants. Chemicals that exhibit one or more of the key characteristics could be prioritized for additional evaluation and testing. A key characteristics approach has the potential to integrate with pathway-based toxicity testing to improve prediction of female reproductive toxicity in chemicals and potentially prevent some toxicants from entering common use. https://doi.org/10.1289/EHP4971.

Identification of gene expression predictors of occupational benzene exposure

BACKGROUND

Previously, using microarrays and mRNA-Sequencing (mRNA-Seq) we found that occupational exposure to a range of benzene levels perturbed gene expression in peripheral blood mononuclear cells.

OBJECTIVES

In the current study, we sought to identify gene expression biomarkers predictive of benzene exposure below 1 part per million (ppm), the occupational standard in the U.S.

METHODS

First, we used the nCounter platform to validate altered expression of 30 genes in 33 unexposed controls and 57 subjects exposed to benzene (<1 to ≥5 ppm). Second, we used SuperLearner (SL) to identify a minimal number of genes for which altered expression could predict <1 ppm benzene exposure, in 44 subjects with a mean air benzene level of 0.55±0.248 ppm (minimum 0.203ppm).

RESULTS

nCounter and microarray expression levels were highly correlated (coefficients >0.7, p<0.05) for 26 microarray-selected genes. nCounter and mRNA-Seq levels were poorly correlated for 4 mRNA-Seq-selected genes. Using negative binomial regression with adjustment for covariates and multiple testing, we confirmed differential expression of 23 microarray-selected genes in the entire benzene-exposed group, and 27 genes in the <1 ppm-exposed subgroup, compared with the control group. Using SL, we identified 3 pairs of genes that could predict <1 ppm benzene exposure with cross-validated AUC estimates >0.9 (p<0.0001) and were not predictive of other exposures (nickel, arsenic, smoking, stress). The predictive gene pairs are PRG2/CLEC5A, NFKBI/CLEC5A, and ACSL1/CLEC5A. They play roles in innate immunity and inflammatory responses.

CONCLUSIONS

Using nCounter and SL, we validated the altered expression of multiple mRNAs by benzene and identified gene pairs predictive of exposure to benzene at levels below the US occupational standard of 1ppm.

Assessing health risks from multiple environmental stressors: Moving from G×E to I×E

Research on disease causation often attempts to isolate the effects of individual factors, including individual genes or environmental factors. This reductionist approach has generated many discoveries, but misses important interactive and cumulative effects that may help explain the broad range of variability in disease occurrence observed across studies and individuals. A disease rarely results from a single factor, and instead results from a broader combination of factors, characterized here as intrinsic (I) and extrinsic (E) factors. Intrinsic vulnerability or resilience emanates from a variety of both fixed and shifting biological factors including genetic traits, while extrinsic factors comprise all biologically-relevant external stressors encountered across the lifespan. The I×E concept incorporates the multi-factorial and dynamic nature of health and disease and provides a unified, conceptual basis for integrating results from multiple areas of research, including genomics, G×E, developmental origins of health and disease, and the exposome. We describe the utility of the I×E concept to better understand and characterize the cumulative impact of multiple extrinsic and intrinsic factors on individual and population health. New research methods increasingly facilitate the measurement of multifactorial and interactive effects in epidemiological and toxicological studies. Tiered or indicator-based approaches can guide the selection of potentially relevant I and E factors for study and quantification, and exposomics methods may eventually produce results that can be used to generate a response function over the life course. Quantitative data on I×E interactive effects should generate a better understanding of the variability in human response to environmental factors. The proposed I×E concept highlights the role for broader study design in order to identify extrinsic and intrinsic factors amenable to interventions at the individual and population levels in order to enhance resilience, reduce vulnerability and improve health.

Formaldehyde induces toxicity in mouse bone marrow and hematopoietic stem/progenitor cells and enhances benzene-induced adverse effects

FA in air for 2 weeks, mimicking occupational exposure, then measured complete blood counts, nucleated BM cell count, and myeloid progenitor colony formation. We also investigated potential mechanisms of FA toxicity, including reactive oxygen species (ROS) generation, apoptosis, and hematopoietic growth factor and receptor levels. FA exposure significantly reduced nucleated BM cells and BM-derived colony-forming unit-granulocyte-macrophage (CFU-GM) and burst-forming unit-erythroid (BFU-E); down-regulated GM-CSFRα and EPOR expression; increased ROS in nucleated BM, spleen and CFU-GM cells; and increased apoptosis in nucleated spleen and CFU-GM cells. FA and BZ each similarly altered BM mature cells and stem/progenitor counts, BM and CFU-GM ROS, and apoptosis in spleen and CFU-GM but had differential effects on other end points. Co-exposure was more potent for several end points. Thus, FA is toxic to the mouse hematopoietic system, including BM stem/progenitor cells, and it enhances BZ-induced toxic effects. Our findings suggest that FA may induce BM toxicity by affecting myeloid progenitor growth and survival through oxidative damage and reduced expression levels of GM-CSFRα and EPOR.

Identification of Genes That Modulate Susceptibility to Formaldehyde and Imatinib by Functional Genomic Screening in Human Haploid KBM7 Cells

Though current functional genomic screening systems are useful for investigating human susceptibility to chemical toxicity, they have limitations. Well-established, high-throughput yeast mutant screens identify only evolutionarily conserved processes. RNA interference can be applied in human cells but is limited by incomplete gene knockout and off-target effects. Human haploid cell screening is advantageous as it requires knockdown of only a single copy of each gene. A human haploid cell mutant library (KBM7-Mu), derived from a chronic myeloid leukemia (CML) patient, was recently developed and has been used to identify genes that modulate sensitivity to infectious agents and pharmaceutical drugs. Here, we sought to improve the KBM7-Mu screening process to enable efficient screening of environmental chemicals. We developed a semi-solid medium based screening approach that cultures individual mutant colonies from chemically resistant cells, faster (by 2-3 weeks) and with less labor than the original liquid medium-based approach. As proof of principle, we identified genetic mutants that confer resistance to the carcinogen formaldehyde (FA, 12 genes, 18 hits) and the CML chemotherapeutic agent imatinib (6 genes, 13 hits). Validation experiments conducted on KBM7 mutants lacking each of the 18 genes confirmed resistance of 6 FA mutants (CTC1, FCRLA, GOT1, LPR5, M1AP, and MAP2K5) and 1 imatinib-resistant mutant (LYRM9). Despite the improvements to the method, it remains technically challenging to limit false positive findings. Nonetheless, our findings demonstrate the broad applicability of this optimized haploid approach to screen toxic chemicals to identify novel susceptibility genes and gain insight into potential mechanisms of toxicity.

Data on megakaryocytes in the bone marrow of mice exposed to formaldehyde

Previously, we reported that occupational exposure to formaldehyde (FA) exposure in factory workers reduced platelet counts, http://dx.doi.org/10.1158/1055-9965.EPI-09-0762[1], while exposure in mice increased platelet counts http://dx.doi.org/10.1371/journal.pone.0074974[2]. Bone marrow megakaryocyte (MK) numbers were also increased in exposed mice, as determined qualitatively. The data presented here are from a quantitative evaluation of MK numbers in the bone marrow histopathological slides from the previous FA exposure experiments in mice. Bone marrow slides were prepared using a single 5 μm section of femur from 2 mice randomly selected from each exposure group (n=9) treated with 0, 0.5 and 3.0 mg/m(3) FA by nose-only inhalation. MKs were systemically counted and average MK frequency was calculated as the total MK per slide divided by the number of fields evaluated. Data are presented visually as microscopy views and graphically as MK frequency.

Induction of centrosome amplification by formaldehyde, but not hydroquinone, in human lymphoblastoid TK6 cells

Benzene and formaldehyde (FA) are important industrial chemicals and environmental pollutants that cause leukemia by inducing DNA damage and chromosome aberrations in hematopoietic stem cells (HSC), the target cells for leukemia. Our previous studies showed that workers exposed to benzene and FA exhibit increased levels of aneuploidy in their blood cells. As centrosome amplification is a common phenomenon in human cancers, including leukemia, and is associated with aneuploidy in carcinogenesis, we hypothesized that benzene and FA would induce centrosome amplification in vitro. We treated human lymphoblastoid TK6 cells with a range of concentrations of hydroquinone (HQ, a benzene metabolite) or FA for 24 h, allowed the cells to recover in fresh medium for 24 h, and examined centrosome amplification; chromosomal gain, loss, and breakage; and cytotoxicity. We included melphalan and etoposide, chemotherapeutic drugs that cause therapy-related acute myeloid leukemia and that have been shown to induce centrosome amplification as well as chromosomal aneuploidy and breakage, as positive controls. Melphalan and etoposide induced centrosome amplification and chromosome gain and breakage in a dose-dependent manner, at cytotoxic concentrations. HQ, though cytotoxic, did not induce centrosome amplification or any chromosomal aberration. FA-induced centrosome amplification and cytotoxicity, but did not induce chromosomal aberrations. Our data suggest, for the first time, that centrosome amplification is a potential mechanism underlying FA-induced leukemogenesis, but not benzene-induced leukemogenesis, as mediated through HQ. Future studies are needed to delineate the mechanisms of centrosome amplification and its association with DNA damage, chromosomal aneuploidy and carcinogenesis, following exposure to FA.

Functional genomic screening approaches in mechanistic toxicology and potential future applications of CRISPR-Cas9

Characterizing variability in the extent and nature of responses to environmental exposures is a critical aspect of human health risk assessment. Chemical toxicants act by many different mechanisms, however, and the genes involved in adverse outcome pathways (AOPs) and AOP networks are not yet characterized. Functional genomic approaches can reveal both toxicity pathways and susceptibility genes, through knockdown or knockout of all non-essential genes in a cell of interest, and identification of genes associated with a toxicity phenotype following toxicant exposure. Screening approaches in yeast and human near-haploid leukemic KBM7 cells have identified roles for genes and pathways involved in response to many toxicants but are limited by partial homology among yeast and human genes and limited relevance to normal diploid cells. RNA interference (RNAi) suppresses mRNA expression level but is limited by off-target effects (OTEs) and incomplete knockdown. The recently developed gene editing approach called clustered regularly interspaced short palindrome repeats-associated nuclease (CRISPR)-Cas9, can precisely knock-out most regions of the genome at the DNA level with fewer OTEs than RNAi, in multiple human cell types, thus overcoming the limitations of the other approaches. It has been used to identify genes involved in the response to chemical and microbial toxicants in several human cell types and could readily be extended to the systematic screening of large numbers of environmental chemicals. CRISPR-Cas9 can also repress and activate gene expression, including that of non-coding RNA, with near-saturation, thus offering the potential to more fully characterize AOPs and AOP networks. Finally, CRISPR-Cas9 can generate complex animal models in which to conduct preclinical toxicity testing at the level of individual genotypes or haplotypes. Therefore, CRISPR-Cas9 is a powerful and flexible functional genomic screening approach that can be harnessed to provide unprecedented mechanistic insight in the field of modern toxicology.

Chromosome-wide aneuploidy study of cultured circulating myeloid progenitor cells from workers occupationally exposed to formaldehyde

Formaldehyde (FA) is an economically important industrial chemical to which millions of people worldwide are exposed environmentally and occupationally. Recently, the International Agency for Cancer Research concluded that there is sufficient evidence that FA causes leukemia, particularly myeloid leukemia. To evaluate the biological plausibility of this association, we employed a chromosome-wide aneuploidy study approach, which allows the evaluation of aneuploidy and structural chromosome aberrations (SCAs) of all 24 chromosomes simultaneously, to analyze cultured myeloid progenitor cells from 29 workers exposed to relatively high levels of FA and 23 unexposed controls. We found statistically significant increases in the frequencies of monosomy, trisomy, tetrasomy and SCAs of multiple chromosomes in exposed workers compared with controls, with particularly notable effects for monosomy 1 [P = 6.02E-06, incidence rate ratio (IRR) = 2.31], monosomy 5 (P = 9.01E-06; IRR = 2.24), monosomy 7 (P = 1.57E-05; IRR = 2.17), trisomy 5 (P = 1.98E-05; IRR = 3.40) and SCAs of chromosome 5 (P = 0.024; IRR = 4.15). The detection of increased levels of monosomy 7 and SCAs of chromosome 5 is particularly relevant as they are frequently observed in acute myeloid leukemia. Our findings provide further evidence that leukemia-related cytogenetic changes can occur in the circulating myeloid progenitor cells of healthy workers exposed to FA, which may be a potential mechanism underlying FA-induced leukemogenesis.

Emerging approaches in predictive toxicology

Predictive toxicology plays an important role in the assessment of toxicity of chemicals and the drug development process. While there are several well-established in vitro and in vivo assays that are suitable for predictive toxicology, recent advances in high-throughput analytical technologies and model systems are expected to have a major impact on the field of predictive toxicology. This commentary provides an overview of the state of the current science and a brief discussion on future perspectives for the field of predictive toxicology for human toxicity. Computational models for predictive toxicology, needs for further refinement and obstacles to expand computational models to include additional classes of chemical compounds are highlighted. Functional and comparative genomics approaches in predictive toxicology are discussed with an emphasis on successful utilization of recently developed model systems for high-throughput analysis. The advantages of three-dimensional model systems and stem cells and their use in predictive toxicology testing are also described.

Characterization of changes in gene expression and biochemical pathways at low levels of benzene exposure

Benzene, a ubiquitous environmental pollutant, causes acute myeloid leukemia (AML). Recently, through transcriptome profiling of peripheral blood mononuclear cells (PBMC), we reported dose-dependent effects of benzene exposure on gene expression and biochemical pathways in 83 workers exposed across four airborne concentration ranges (from <1 ppm to >10 ppm) compared with 42 subjects with non-workplace ambient exposure levels. Here, we further characterize these dose-dependent effects with continuous benzene exposure in all 125 study subjects. We estimated air benzene exposure levels in the 42 environmentally-exposed subjects from their unmetabolized urinary benzene levels. We used a novel non-parametric, data-adaptive model selection method to estimate the change with dose in the expression of each gene. We describe non-parametric approaches to model pathway responses and used these to estimate the dose responses of the AML pathway and 4 other pathways of interest. The response patterns of majority of genes as captured by mean estimates of the first and second principal components of the dose-response for the five pathways and the profiles of 6 AML pathway response-representative genes (identified by clustering) exhibited similar apparent supra-linear responses. Responses at or below 0.1 ppm benzene were observed for altered expression of AML pathway genes and CYP2E1. Together, these data show that benzene alters disease-relevant pathways and genes in a dose-dependent manner, with effects apparent at doses as low as 100 ppb in air. Studies with extensive exposure assessment of subjects exposed in the low-dose range between 10 ppb and 1 ppm are needed to confirm these findings.

Application of toxicogenomic profiling to evaluate effects of benzene and formaldehyde: from yeast to human

Genetic variation underlies a significant proportion of the individual variation in human susceptibility to toxicants. The primary current approaches to identify gene-environment (GxE) associations, genome-wide association studies and candidate gene association studies, require large exposed and control populations and an understanding of toxicity genes and pathways, respectively. This limits their application in the study of GxE associations for the leukemogens benzene and formaldehyde, whose toxicity has long been a focus of our research. As an alternative approach, our published work has applied innovative in vitro functional genomics testing systems, including unbiased functional screening assays in yeast and a near-haploid human bone marrow cell line. Through comparative genomic and computational analyses of the resulting data, human genes and pathways that may modulate susceptibility to benzene and formaldehyde were identified, and the roles of several genes in mammalian cell models were validated. In populations occupationally exposed to low levels of benzene, we applied peripheral blood mononuclear cell transcriptomics and chromosome-wide aneuploidy studies in lymphocytes. In this review, we describe our comprehensive toxicogenomic approach and the potential mechanisms of toxicity and susceptibility genes identified for benzene and formaldehyde, as well as related studies conducted by other researchers.

Inhaled formaldehyde induces DNA-protein crosslinks and oxidative stress in bone marrow and other distant organs of exposed mice

Formaldehyde (FA), a major industrial chemical and ubiquitous environmental pollutant, has been classified as a leukemogen. The causal relationship remains unclear, however, due to limited evidence that FA induces toxicity in bone marrow, the site of leukemia induction, and in other distal organs. Although induction of DNA-protein crosslinks (DPC), a hallmark of FA toxicity, was not previously detected in the bone marrow of FA-exposed rats and monkeys in studies published in the 1980s, our recent studies showed increased DPC in the bone marrow, liver, kidney, and testes of exposed Kunming mice. To confirm these preliminary results, in the current study we exposed BALB/c mice to 0, 0.5, 1.0, and 3.0 mg m(-3) FA (8 hr per day, for 7 consecutive days) by nose-only inhalation and measured DPC levels in bone marrow and other organs of exposed mice. As oxidative stress is a potential mechanism of FA toxicity, we also measured glutathione (GSH), reactive oxygen species (ROS), and malondialdehyde (MDA), in the bone marrow, peripheral blood mononuclear cells, lung, liver, spleen, and testes of exposed mice. Significant dose-dependent increases in DPC, decreases in GSH, and increases in ROS and MDA were observed in all organs examined (except for DPC in lung). Bone marrow was among the organs with the strongest effects for DPC, GSH, and ROS. In conclusion, exposure of mice to FA by inhalation induced genotoxicity and oxidative stress in bone marrow and other organs. These findings strengthen the biological plausibility of FA-induced leukemogenesis and systemic toxicity.

Global gene expression response of a population exposed to benzene: a pilot study exploring the use of RNA-sequencing technology

The mechanism of toxicity of the leukemogen benzene is not entirely known. This pilot study used RNA-sequencing (RNA-seq) technology to examine the effect of benzene exposure on gene expression in peripheral blood mononuclear cells obtained from 10 workers occupationally exposed to high levels of benzene (≥5 ppm) in air and 10 matched unexposed control workers, from a large study (n = 125) in which gene expression was previously measured by microarray. RNA-seq is more sensitive and has a wider dynamic range for the quantification of gene expression. Further, it has the ability to detect novel transcripts and alternative splice variants. The main conclusions from our analysis of the 20 workers by RNA-seq are as follows: The Pearson correlation between the two technical replicates for the RNA-seq experiments was 0.98 and the correlation between RNA-seq and microarray signals for the 20 subjects was around 0.6. 60% of the transcripts with detected reads from the RNA-seq experiments did not have corresponding probes on the microarrays. Fifty-three percent of the transcripts detected by RNA-seq and 99% of those with probes on the microarray were protein-coding. There was a significant overlap (P < 0.05) in transcripts declared differentially expressed due to benzene exposure using the two technologies. About 20% of the transcripts declared differentially expressed using the RNA-seq data were non-coding transcripts. Six transcripts were determined (false-discovery rate < 0.05) to be alternatively spliced as a result of benzene exposure. Overall, this pilot study shows that RNA-seq can complement the information obtained by microarray in the analysis of changes in transcript expression from chemical exposures.

Analysis of the transcriptome in molecular epidemiology studies

The human transcriptome is complex, comprising multiple transcript types, mostly in the form of non-coding RNA (ncRNA). The majority of ncRNA is of the long form (lncRNA, ≥ 200 bp), which plays an important role in gene regulation through multiple mechanisms including epigenetics, chromatin modification, control of transcription factor binding, and regulation of alternative splicing. Both mRNA and ncRNA exhibit additional variability in the form of alternative splicing and RNA editing. All aspects of the human transcriptome can potentially be dysregulated by environmental exposures. Next-generation RNA sequencing (RNA-Seq) is the best available methodology to measure this although it has limitations, including experimental bias. The third phase of the MicroArray Quality Control Consortium project (MAQC-III), also called Sequencing Quality Control (SeQC), aims to address these limitations through standardization of experimental and bioinformatic methodologies. A limited number of toxicogenomic studies have been conducted to date using RNA-Seq. This review describes the complexity of the human transcriptome, the application of transcriptomics by RNA-Seq or microarray in molecular epidemiology studies, and limitations of these approaches including the type of cell or tissue analyzed, experimental variation, and confounding. By using good study designs with precise, individual exposure measurements, sufficient power and incorporation of phenotypic anchors, studies in human populations can identify biomarkers of exposure and/or early effect and elucidate mechanisms of action underlying associated diseases, even at low doses. Analysis of datasets at the pathway level can compensate for some of the limitations of RNA-Seq and, as more datasets become available, will increasingly elucidate the exposure-disease continuum.

Alterations in serum immunoglobulin levels in workers occupationally exposed to trichloroethylene

Trichloroethylene (TCE) has been associated with a variety of immunotoxic effects and may be associated with an increased risk of non-Hodgkin lymphoma (NHL). Altered serum immunoglobulin (Ig) levels have been reported in NHL patients and in animals exposed to TCE. Recently, we reported that occupational exposure to TCE is associated with immunosuppressive effects and immune dysfunction, including suppression of B-cell counts and activation, even at relatively low levels. We hypothesized that TCE exposure would also affect Ig levels in humans. We measured serum levels of IgG, IgM and IgE, by enzyme-linked immunosorbent assay, in TCE-exposed workers (n = 80) and unexposed controls (n = 45), matched by age and gender, in a cross-sectional, molecular epidemiology study of occupational exposure to TCE in Guangdong, China. Exposed workers had about a 17.5% decline in serum levels of IgG compared with unexposed controls (P = 0.0002). Similarly, serum levels of IgM were reduced by about 38% in workers exposed to TCE compared with unexposed controls (P < 0.0001). Serum levels of both IgG and IgM were significantly decreased in workers exposed to TCE levels below 12 p.p.m., the median exposure level. Adjustment for B-cell counts had minimal impact on our findings. IgE levels were not significantly different between exposed and control subjects. These results provide further evidence that TCE is immunotoxic at relatively low exposure levels and provide additional biologic plausibility for the reported association of TCE with NHL.

Using bioinformatic approaches to identify pathways targeted by human leukemogens

We have applied bioinformatic approaches to identify pathways common to chemical leukemogens and to determine whether leukemogens could be distinguished from non-leukemogenic carcinogens. From all known and probable carcinogens classified by IARC and NTP, we identified 35 carcinogens that were associated with leukemia risk in human studies and 16 non-leukemogenic carcinogens. Using data on gene/protein targets available in the Comparative Toxicogenomics Database (CTD) for 29 of the leukemogens and 11 of the non-leukemogenic carcinogens, we analyzed for enrichment of all 250 human biochemical pathways in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. The top pathways targeted by the leukemogens included metabolism of xenobiotics by cytochrome P450, glutathione metabolism, neurotrophin signaling pathway, apoptosis, MAPK signaling, Toll-like receptor signaling and various cancer pathways. The 29 leukemogens formed 18 distinct clusters comprising 1 to 3 chemicals that did not correlate with known mechanism of action or with structural similarity as determined by 2D Tanimoto coefficients in the PubChem database. Unsupervised clustering and one-class support vector machines, based on the pathway data, were unable to distinguish the 29 leukemogens from 11 non-leukemogenic known and probable IARC carcinogens. However, using two-class random forests to estimate leukemogen and non-leukemogen patterns, we estimated a 76% chance of distinguishing a random leukemogen/non-leukemogen pair from each other.

Leukemia-related chromosomal loss detected in hematopoietic progenitor cells of benzene-exposed workers

Benzene exposure causes acute myeloid leukemia and hematotoxicity, shown as suppression of mature blood and myeloid progenitor cell numbers. As the leukemia-related aneuploidies monosomy 7 and trisomy 8 previously had been detected in the mature peripheral blood cells of exposed workers, we hypothesized that benzene could cause leukemia through the induction of these aneuploidies in hematopoietic stem and progenitor cells. We measured loss and gain of chromosomes 7 and 8 by fluorescence in situ hybridization in interphase colony-forming unit-granulocyte-macrophage (CFU-GM) cells cultured from otherwise healthy benzene-exposed (n=28) and unexposed (n=14) workers. CFU-GM monosomy 7 and 8 levels (but not trisomy) were significantly increased in subjects exposed to benzene overall, compared with levels in the control subjects (P=0.0055 and P=0.0034, respectively). Levels of monosomy 7 and 8 were significantly increased in subjects exposed to <10 p.p.m. (20%, P=0.0419 and 28%, P=0.0056, respectively) and ≥ 10 p.p.m. (48%, P=0.0045 and 32%, 0.0354) benzene, compared with controls, and significant exposure-response trends were detected (P(trend)=0.0033 and 0.0057). These data show that monosomies 7 and 8 are produced in a dose-dependent manner in the blood progenitor cells of workers exposed to benzene, and may be mechanistically relevant biomarkers of early effect for benzene and other leukemogens.

Current understanding of the mechanism of benzene-induced leukemia in humans: implications for risk assessment

Benzene causes acute myeloid leukemia and probably other hematological malignancies. As benzene also causes hematotoxicity even in workers exposed to levels below the US permissible occupational exposure limit of 1 part per million, further assessment of the health risks associated with its exposure, particularly at low levels, is needed. Here, we describe the probable mechanism by which benzene induces leukemia involving the targeting of critical genes and pathways through the induction of genetic, chromosomal or epigenetic abnormalities and genomic instability, in a hematopoietic stem cell (HSC); stromal cell dysregulation; apoptosis of HSCs and stromal cells and altered proliferation and differentiation of HSCs. These effects modulated by benzene-induced oxidative stress, aryl hydrocarbon receptor dysregulation and reduced immunosurveillance, lead to the generation of leukemic stem cells and subsequent clonal evolution to leukemia. A mode of action (MOA) approach to the risk assessment of benzene was recently proposed. This approach is limited, however, by the challenges of defining a simple stochastic MOA of benzene-induced leukemogenesis and of identifying relevant and quantifiable parameters associated with potential key events. An alternative risk assessment approach is the application of toxicogenomics and systems biology in human populations, animals and in vitro models of the HSC stem cell niche, exposed to a range of levels of benzene. These approaches will inform our understanding of the mechanisms of benzene toxicity and identify additional biomarkers of exposure, early effect and susceptibility useful for risk assessment.

Reproductive and developmental toxicity of formaldehyde: a systematic review

Formaldehyde, the recently classified carcinogen and ubiquitous environmental contaminant, has long been suspected of causing adverse reproductive and developmental effects, but previous reviews were inconclusive, due in part, to limitations in the design of many of the human population studies. In the current review, we systematically evaluated evidence of an association between formaldehyde exposure and adverse reproductive and developmental effects, in human populations and in vivo animal studies, in the peer-reviewed literature. The mostly retrospective human studies provided evidence of an association of maternal exposure with adverse reproductive and developmental effects. Further assessment of this association by meta-analysis revealed an increased risk of spontaneous abortion (1.76, 95% CI 1.20-2.59, p=0.002) and of all adverse pregnancy outcomes combined (1.54, 95% CI 1.27-1.88, p<0.001), in formaldehyde-exposed women, although differential recall, selection bias, or confounding cannot be ruled out. Evaluation of the animal studies including all routes of exposure, doses and dosing regimens studied, suggested positive associations between formaldehyde exposure and reproductive toxicity, mostly in males. Potential mechanisms underlying formaldehyde-induced reproductive and developmental toxicities, including chromosome and DNA damage (genotoxicity), oxidative stress, altered level and/or function of enzymes, hormones and proteins, apoptosis, toxicogenomic and epigenomic effects (such as DNA methylation), were identified. To clarify these associations, well-designed molecular epidemiologic studies, that include quantitative exposure assessment and diminish confounding factors, should examine both reproductive and developmental outcomes associated with exposure in males and females. Together with mechanistic and animal studies, this will allow us to better understand the systemic effect of formaldehyde exposure.

Benzene, the exposome and future investigations of leukemia etiology

Benzene exposure is associated with acute myeloid leukemia (AML), myelodysplastic syndromes (MDS), and probably lymphoma and childhood leukemia. Biological plausibility for a causal role of benzene in these diseases comes from its toxicity to hematopoietic stem cells (HSC) or progenitor cells, from which all leukemias and related disorders arise. The effect of this toxicity is manifest as lowered blood counts (hematotoxicity), even in individuals occupationally exposed to low levels of benzene. Benzene can induce AML/MDS via several well-characterized pathways associated with these diseases. Through its metabolites, benzene induces multiple alterations that likely contribute to the leukemogenic process, and appears to operate via multiple modes of action. To improve mechanistic understanding and for risk assessment purposes, it may be possible to measure several of the key events in these modes of action in an in vitro model of the bone marrow stem cell niche. Even though benzene is leukemogenic at relatively low occupational levels of exposure, it seems unlikely that it is a major cause of leukemia in the general population exposed to benzene in the ppb range. Other established non-genetic causes of AML, e.g. smoking, ionizing radiation and cancer chemotherapy, also only explain about 20% of AML incidence, leaving ∼80% unexplained. The question arises as to how to find the causes of the majority of de novo AMLs that remain unexplained. We propose that we should attempt to characterize the 'exposome' of human leukemia by using unbiased laboratory-based methods to find the unknown 'environmental' factors that contribute to leukemia etiology.

Paternal smoking and risk of childhood acute lymphoblastic leukemia: systematic review and meta-analysis

Objective. To investigate the association between paternal smoking and childhood acute lymphoblastic leukemia (ALL). Method. We identified 18 published epidemiologic studies that reported data on both paternal smoking and childhood ALL risk. We performed a meta-analysis and analyzed dose-response relationships on ALL risk for smoking during preconception, during pregnancy, after birth, and ever smoking. Results. The summary odds ratio (OR) of childhood ALL associated with paternal smoking was 1.11 (95% Confidence Interval (CI): 1.05-1.18, I(2) = 18%) during any time period, 1.25 (95% CI: 1.08-1.46, I(2) = 53%) preconception; 1.24 (95% CI: 1.07-1.43, I(2) = 54%) during pregnancy, and 1.24 (95% CI: 0.96-1.60, I(2) = 64%) after birth, with a dose-response relationship between childhood ALL and paternal smoking preconception or after birth. Conclusion. The evidence supports a positive association between childhood ALL and paternal ever smoking and at each exposure time period examined. Future epidemiologic studies should assess paternal smoking during well-defined exposure windows and should include biomarkers to assess smoking exposure and toxicological mechanisms.

Global gene expression profiling of a population exposed to a range of benzene levels

BACKGROUND

Benzene, an established cause of acute myeloid leukemia (AML), may also cause one or more lymphoid malignancies in humans. Previously, we identified genes and pathways associated with exposure to high (> 10 ppm) levels of benzene through transcriptomic analyses of blood cells from a small number of occupationally exposed workers.

OBJECTIVES

The goals of this study were to identify potential biomarkers of benzene exposure and/or early effects and to elucidate mechanisms relevant to risk of hematotoxicity, leukemia, and lymphoid malignancy in occupationally exposed individuals, many of whom were exposed to benzene levels < 1 ppm, the current U.S. occupational standard.

METHODS

We analyzed global gene expression in the peripheral blood mononuclear cells of 125 workers exposed to benzene levels ranging from < 1 ppm to > 10 ppm. Study design and analysis with a mixed-effects model minimized potential confounding and experimental variability.

RESULTS

We observed highly significant widespread perturbation of gene expression at all exposure levels. The AML pathway was among the pathways most significantly associated with benzene exposure. Immune response pathways were associated with most exposure levels, potentially providing biological plausibility for an association between lymphoma and benzene exposure. We identified a 16-gene expression signature associated with all levels of benzene exposure.

CONCLUSIONS

Our findings suggest that chronic benzene exposure, even at levels below the current U.S. occupational standard, perturbs many genes, biological processes, and pathways. These findings expand our understanding of the mechanisms by which benzene may induce hematotoxicity, leukemia, and lymphoma and reveal relevant potential biomarkers associated with a range of exposures.

Chromosome-wide aneuploidy study (CWAS) in workers exposed to an established leukemogen, benzene

Evidence suggests that de novo, therapy-related and benzene-induced acute myeloid leukemias (AML) occur via similar cytogenetic and genetic pathways, several of which involve aneuploidy, the loss or gain of chromosomes. Aneuploidy of specific chromosomes has been detected in benzene-related leukemia patients as well as in healthy benzene-exposed workers, suggesting that aneuploidy precedes and may be a potential mechanism underlying benzene-induced leukemia. Here, we analyzed the peripheral blood lymphocytes of 47 exposed workers and 27 unexposed controls using a novel OctoChrome fluorescence in situ hybridization (FISH) technique that simultaneously detects aneuploidy in all 24 chromosomes. Through this chromosome-wide aneuploidy study (CWAS) approach, we found heterogeneity in the monosomy and trisomy rates of the 22 autosomes when plotted against continuous benzene exposure. In addition, statistically significant, chromosome-specific increases in the rates of monosomy [5, 6, 7, 10, 16 and 19] and trisomy [5, 6, 7, 8, 10, 14, 16, 21 and 22] were found to be dose dependently associated with benzene exposure. Furthermore, significantly higher rates of monosomy and trisomy were observed in a priori defined 'susceptible' chromosome sets compared with all other chromosomes. Together, these findings confirm that benzene exposure is associated with specific chromosomal aneuploidies in hematopoietic cells, which suggests that such aneuploidies may play roles in benzene-induced leukemogenesis.

An emerging role for epigenetic dysregulation in arsenic toxicity and carcinogenesis

BACKGROUND

Exposure to arsenic, an established human carcinogen, through consumption of highly contaminated drinking water is a worldwide public health concern. Several mechanisms by which arsenical compounds induce tumorigenesis have been proposed, including oxidative stress, genotoxic damage, and chromosomal abnormalities. Recent studies have suggested that epigenetic mechanisms may also mediate toxicity and carcinogenicity resulting from arsenic exposure.

OBJECTIVE

We examined the evidence supporting the roles of the three major epigenetic mechanisms-DNA methylation, histone modification, and microRNA (miRNA) expression-in arsenic toxicity and, in particular, carcinogenicity. We also investigated future research directions necessary to clarify epigenetic and other mechanisms in humans.

DATA SOURCES AND SYNTHESIS

We conducted a PubMed search of arsenic exposure and epigenetic modification through April 2010 and summarized the in vitro and in vivo research findings, from both our group and others, on arsenic-associated epigenetic alteration and its potential role in toxicity and carcinogenicity.

CONCLUSIONS

Arsenic exposure has been shown to alter methylation levels of both global DNA and gene promoters; histone acetylation, methylation, and phosphorylation; and miRNA expression, in studies analyzing mainly a limited number of epigenetic end points. Systematic epigenomic studies in human populations exposed to arsenic or in patients with arsenic-associated cancer have not yet been performed. Such studies would help to elucidate the relationship between arsenic exposure, epigenetic dysregulation, and carcinogenesis and are becoming feasible because of recent technological advancements.

Toxicogenomic profiling of chemically exposed humans in risk assessment

Gene-environment interactions contribute to complex disease development. The environmental contribution, in particular low-level and prevalent environmental exposures, may constitute much of the risk and contribute substantially to disease. Systematic risk evaluation of the majority of human chemical exposures, has not been conducted and is a goal of regulatory agencies in the U.S. and worldwide. With the recent recognition that toxicological approaches more predictive of effects in humans are required for risk assessment, in vitro human cell line data as well as animal data are being used to identify toxicity mechanisms that can be translated into biomarkers relevant to human exposure studies. In this review, we discuss how data from toxicogenomic studies of exposed human populations can inform risk assessment, by generating biomarkers of exposure, early effect, and/or susceptibility, elucidating mechanisms of action underlying exposure-related disease, and detecting response at low doses. Good experimental design incorporating precise, individual exposure measurements, phenotypic anchors (pre-disease or traditional toxicological markers), and a range of relevant exposure levels, is necessary. Further, toxicogenomic studies need to be designed with sufficient power to detect true effects of the exposure. As more studies are performed and incorporated into databases such as the Comparative Toxicogenomics Database (CTD) and Chemical Effects in Biological Systems (CEBS), data can be mined for classification of newly tested chemicals (hazard identification), and, for investigating the dose-response, and inter-relationship among genes, environment and disease in a systems biology approach (risk characterization).

Systems biology of human benzene exposure

Toxicogenomic studies, including genome-wide analyses of susceptibility genes (genomics), gene expression (transcriptomics), protein expression (proteomics), and epigenetic modifications (epigenomics), of human populations exposed to benzene are crucial to understanding gene-environment interactions, providing the ability to develop biomarkers of exposure, early effect and susceptibility. Comprehensive analysis of these toxicogenomic and epigenomic profiles by bioinformatics in the context of phenotypic endpoints, comprises systems biology, which has the potential to comprehensively define the mechanisms by which benzene causes leukemia. We have applied this approach to a molecular epidemiology study of workers exposed to benzene. Hematotoxicity, a significant decrease in almost all blood cell counts, was identified as a phenotypic effect of benzene that occurred even below 1 ppm benzene exposure. We found a significant decrease in the formation of progenitor colonies arising from bone marrow stem cells with increasing benzene exposure, showing that progenitor cells are more sensitive to the effects of benzene than mature blood cells, likely leading to the observed hematotoxicity. Analysis of transcriptomics by microarray in the peripheral blood mononuclear cells of exposed workers, identified genes and pathways (apoptosis, immune response, and inflammatory response) altered at high (>10 ppm) and low (<1 ppm) benzene levels. Serum proteomics by SELDI-TOF-MS revealed proteins consistently down-regulated in exposed workers. Preliminary epigenomics data showed effects of benzene on the DNA methylation of specific genes. Genomic screens for candidate genes involved in susceptibility to benzene toxicity are being undertaken in yeast, with subsequent confirmation by RNAi in human cells, to expand upon the findings from candidate gene analyses. Data on these and future biomarkers will be used to populate a large toxicogenomics database, to which we will apply bioinformatic approaches to understand the interactions among benzene toxicity, susceptibility genes, mRNA, and DNA methylation through a systems biology approach.

The benzene metabolite, hydroquinone and etoposide both induce endoreduplication in human lymphoblastoid TK6 cells

Both occupational exposure to the leukemogen benzene and in vitro exposure to its metabolite hydroquinone (HQ) lead to the induction of numerical and structural chromosome changes. Several studies have shown that HQ can form DNA adducts, disrupt microtubule assembly and inhibit DNA topoisomerase II (topo II) activity. As these are potential mechanisms underlying endoreduplication (END), a phenomenon that involves DNA amplification without corresponding cell division, we hypothesized that HQ could cause END. We measured END in the human lymphoblastoid cell line, TK6, treated with HQ (0-20 microM) and etoposide (0-0.2 microM) for 48 h. Etoposide was used as a positive control as it is a topo II poison and established human leukemogen that has previously been shown to induce END in Chinese hamster ovary cells. Both HQ and etoposide significantly induced END in a dose-dependent manner (P(trend) < 0.0001 and P(trend) = 0.0003, respectively). Since END may underlie the acquisition of high chromosome numbers by tumour cells, it may play a role in inducing genomic instability and subsequent carcinogenesis from HQ and etoposide. In order to further explore the cytogenetic effects of HQ and etoposide, we also examined specific structural changes. HQ did not induce translocations of chromosome 11 [t(11;?)] but significantly induced translocations of chromosome 21 [t(21;?)] and structural chromosome aberrations (SCA) (P(trend) = 0.0415 and P(trend) < 0.0001, respectively). Etoposide potently induced all these structural changes (P(trend) < 0.0001). The lack of an effect of HQ on t(11;?) and the reduced ability of HQ to induce t(21;?) and SCA, compared with etoposide, further suggests that HQ acts primarily as a topo II catalytic inhibitor rather than as a topo II poison in intact human cells.

Microarray analysis of gene expression in peripheral blood mononuclear cells from dioxin-exposed human subjects

Tetrachlorodibenzo-p-dioxin (TCDD) is classified as a human carcinogen and exerts toxic effects on the skin (chloracne). Effects on reproductive, immunological, and endocrine systems have also been observed in animal models. TCDD acts through the aryl hydrocarbon receptor (AhR) pathway influencing largely unknown gene networks. An industrial accident in Seveso, Italy in 1976 exposed thousands of people to substantial quantities of TCDD. Twenty years after the exposure, this study examines global gene expression in the mononuclear cells of 26 Seveso female never smokers, with similar age, alcohol consumption, use of medications, and background plasma levels of 22 dioxin congeners unrelated to the Seveso accident. Plasma dioxin levels were still elevated in the exposed subjects. We performed analyses in two different comparison groups. The first included high-exposed study subjects compared with individuals with background TCDD levels (average plasma levels 99.4 and 6.7ppt, respectively); the second compared subjects who developed chloracne after the accident, and those who did not develop this disease. Overall, we observed a modest alteration of gene expression based on dioxin levels or on chloracne status. In the comparison between high levels and background levels of TCDD, four histone genes were up-regulated and modified expression of HIST1H3H was confirmed by real-time PCR. In the comparison between chloracne case-control subjects, five hemoglobin genes were up-regulated. Pathway analysis revealed two major networks for each comparison, involving cell proliferation, apoptosis, immunological and hematological disease, and other pathways. Further examination of the role of these genes in dioxin induced-toxicity is warranted.

Molecular biomarkers for the study of childhood leukemia

Various specific chromosome rearrangements, including t(8;21), t(15;17), and inv(16), are found in acute myeloid leukemia (AML) and in childhood acute lymphocytic leukemia (ALL), t(12;21) and t(1;19) are common. We sequenced the translocation breakpoints of 56 patients with childhood ALL or AML harboring t(12;21), t(8;21), t(15;17), inv(16), and t(1;19), and demonstrated, with the notable exception of t(1;19), that these rearrangements are commonly detected in the neonatal blood spots (Guthrie cards) of the cases. These findings show that most childhood leukemias begin before birth and that maternal and perinatal exposures such as chemical and infectious agents are likely to be critical. Indeed, we have reported that exposure to indoor pesticides during pregnancy and the first year of life raises leukemia risk, but that later exposures do not. We have also examined aberrant gene methylation in different cytogenetic subgroups and have found striking differences between them, suggesting that epigenetic events are also important in the development of some forms of childhood leukemia. Further, at least two studies now show that the inactivating NAD(P)H:quinone acceptor oxidoreductase (NQO1) C609T polymorphism is positively associated with leukemias arising in the first 1-2 years of life and polymorphisms in the 5,10-methylenetetrahydrofolate reductase (MTHFR) gene have been associated with adult and childhood ALL. Thus, low folate intake and compounds that are detoxified by NQO1 may be important in elevating leukemia risk in children. Finally, we are exploring the use of proteomics to subclassify leukemia, because cytogenetic analysis is costly and time-consuming. Several proteins have been identified that may serve as useful biomarkers for rapidly identifying different forms of childhood leukemia.

Lack of increased genetic damage in 1,3-butadiene-exposed Chinese workers studied in relation to EPHX1 and GST genotypes

1,3-Butadiene (BD) is an important industrial chemical and pollutant. Its ability to induce genetic damage and cause hematological malignancies in humans is controversial. We have examined chromosome damage by fluorescence in situ hybridization (FISH) and mutations in the HPRT gene in the blood of Chinese workers exposed to BD. Peripheral blood samples were collected and cultured from 39 workers exposed to BD (median level 2 ppm, 6 h time-weighted average) and 38 matched controls in Yanshan, China. No difference in the level of aneuploidy or structural changes in chromosomes 1, 7, 8, and 12 was detected in metaphase cells from exposed subjects in comparison with matched controls, nor was there an increase in the frequency of HPRT mutations in the BD-exposed workers. Because genetic polymorphisms in glutathione S-transferase (GST) enzymes and microsomal epoxide hydrolase (EPHX1) may affect the genotoxic effects of BD and its metabolites, we also related chromosome alterations and gene mutations to GSTT1, GSTM1 and EPHX1 genotypes. Overall, there was no effect of variants in these genotypes on numerical or structural changes in chromosomes 1, 7, 8 and 12 or on HPRT mutant frequency in relation to BD exposure, but the GST genotypes did influence background levels of both hyperdiploidy and HPRT mutant frequency. In conclusion, our data show no increase in chromosomal aberrations or HPRT mutations among workers exposed to BD, even in potentially susceptible genetic subgroups. The study is, however, quite small and the levels of BD exposure are not extremely high, but our findings in China do support those from a similar study conducted in the Czech Republic. Together, these studies suggest that low levels of occupational BD exposure do not pose a significant risk of genetic damage.

Prenatal origin of chromosomal translocations in acute childhood leukemia: implications and future directions

We, and others, have demonstrated an in utero origin for translocations associated with childhood leukemia, with latency periods in some cases exceeding 10 years. The mechanism of generation of most of the translocations is thought to be aberrant repair following abortive apoptosis, rather than V(D)J recombination or exposure to topoisomerase II inhibitors. Folate supplementation may prevent some of the chromosome breakage leading to translocation formation. Translocations t(8;21) and t(12;21) have been shown to occur in the normal population (before birth) at a frequency that is 100-fold greater than the risk of developing the corresponding leukemia. In most instances, additional genetic changes are required for progression to leukemia. Tyrosine kinase receptor (RTK) mutations, which give cells a survival/proliferative advantage, are proposed to act cooperatively with fusion genes, leading to transformation. However, translocations and cooperating RTK mutations have not been identified for all leukemia subtypes, particularly in acute myeloid leukemia. The core binding transcriptional pathway is frequently targeted by translocation in utero. We propose that this pathway is highly sensitive during fetal hematopoiesis and may be targeted by mechanisms other than translocation. For each leukemia subtype it is important to characterize the corresponding leukemic stem cell, which is thought to be the initial target for translocation. This would help to elucidate the molecular pathways involved in the progression from preleukemic clone harboring a translocation to fully disseminated leukemia.

Prenatal origin of TEL-AML1-positive acute lymphoblastic leukemia in children born in California

Acute lymphoblastic leukemia (ALL) is the most common form of childhood cancer. The peak incidence of ALL between ages 2 and 5 is accounted for by one subtype, referred to as common acute lymphoblastic leukemia (cALL). About 25% of cALL patients have the TEL-AML1 gene fusion derived from the t(12;21) chromosomal translocation. Recent evidence from retrospective analysis of neonatal blood spots (Guthrie cards) in Europe has demonstrated that this chromosome translocation may arise prenatally. The aim of our study was to determine whether TEL-AML1 fusions arise prenatally in a U.S. population of cALL patients. TEL-AML1-positive cALL cases (n = 14) were identified by fluorescence in situ hybridization, and the genomic breakpoints were identified by a streamlined long-distance PCR approach and sequenced. Clonotypic primers were designed for each patient breakpoint, and a nested PCR assay was used to determine the presence of the TEL-AML1 fusion sequence in neonatal Guthrie cards. Seven of 14 cases demonstrated clonotypic sequences on the archival Guthrie cards. The oldest patient that was positive was 6.7 years old at the time of diagnosis of leukemia. These results confirm previously published findings of a prenatal origin of TEL-AML1 in Europe by demonstrating its occurrence in a California-born population. Secondary changes were also similar to those described previously, with deletion of the second TEL allele being the most common. Other secondary changes included duplication of the fusion gene, trisomy 21, and monosomy X.

Regulation of erythropoietin gene expression depends on two different oxygen-sensing mechanisms

Erythropoietin (Epo), a glycoprotein hormone produced principally in the fetal kidney and in the adult liver in response to hypoxia, is the prime regulator of growth and differentiation in erythroid progenitor cells. The regulation of Epo gene expression is not fully understood, but two mechanisms have been proposed. One involves the participation of a heme protein capable of reversible oxygenation and the other depends on the intracellular concentration of reactive oxygen species (ROS), assumed to be a function of pO2. We have investigated the production of Epo in response to three stimuli, hypoxia, cobalt chloride, and the iron chelator desferrioxamine, in Hep3B cells. As expected, hypoxia caused a marked rise in Epo production. When the cells were exposed to the paired stimuli of hypoxia and cobalt no further increase was found. In contrast, chelation of iron under hypoxic conditions markedly enhanced Epo production, suggesting that the two stimuli act by separate pathways. The addition of carbon monoxide inhibited hypoxia-induced Epo production, independent of desferrioxamine concentration. Taken together these data support the concept that pO2 and ROS are sensed independently.

Erythroid gene expression is differentially regulated by erythropoietin, haemin and delta-aminolaevulinic acid in UT-7 cells

Erythropoietin (Epo) is essential for the later stages of erythropoiesis, acting to promote cell survival and proliferation, but its role in differentiation remains to be defined. The UT-7 cell line exhibits both erythroid and megakaryocytic characteristics and can be induced to differentiate along the erythroid pathway by Epo or the megakaryocytic pathway by phorbol myristic acetate. We have compared the effects of Epo and the chemical inducers, delta-aminolaevulinic acid (delta-ALA) and haemin on the differentiation capacity of UT-7 cells. Epo alone promoted relatively early events in erythroid maturation, without significant changes in haemoglobin production or morphology. GATA-2 and c-myb were down-regulated by Epo, and GATA-2 was further down-modulated by the inducers. Conversely, SCL expression was up-regulated by Epo and further increased by haemin and delta-ALA. Epo caused an increase in the proportion of cells expressing cell surface glycophorin A (GPA) and up-regulated beta- and gamma-globin by several fold. Both haemin and delta-ALA caused a de novo increase in alpha-globin expression as well as enhancing Epo-induced beta-globin expression, leading to a marked increase in haemoglobin production. These results suggest that haemoglobin production in UT-7 cells is limited by a deficiency of erythroid-specific aminolaevulinic acid synthase (ALAS-E) activity or globin synthesis as a consequence of their immaturity as a multipotential cell line.

Expression and mutagenesis of recombinant human and murine erythropoietins in Escherichia coli

Expression of the polypeptide hormone erythropoietin (EPO) in Escherichia coli by four bacterial expression vectors was examined. Complementary DNAs encoding human and murine EPO were amplified by polymerase chain reaction (PCR) and cloned into the glutathione-S-transferase (GST) fusion vector, pGEX-2T. Human EPO DNA was also cloned into the vectors, pET14b, pIN III-Omp A2 and pT7/7. Expression of human and murine EPO was obtained using constructs based on pGEX-2T. For constructs based on the other vectors, expression of EPO was absent or occurred at low levels, despite attempts to optimise conditions. Human and murine EPO, expressed as fusion proteins with GST, were partially soluble and displayed EPO bioactivity. Soluble GST-EPO fusion proteins were affinity purified on immobilised glutathione. Insoluble protein could also be purified by elution from gel slices following SDS-PAGE to yield either fusion protein or, after treatment with thrombin, unmodified EPO which was both soluble and bioactive. The pGEX expression system was evaluated as a means of analysing the structure-function relationships of EPO by in vitro mutagenesis. Three human and three murine EPO mutants were constructed and expressed as GST fusion proteins. Following purification, biological activity was evaluated using assays for bioactivity, immunoactivity and GST activity. The pGEX expression system complements eukaryotic systems described previously for expression of EPO and should provide much useful information about the structure-function relationships of the hormone.

Structure-function relationships of the erythropoietin molecule

The tertiary structure of erythropoietin (EPO) remains to be elucidated by X-ray crystallography. Although the amino acid sequence of EPO is known, the specific features that confer its biological activity are not well understood. In order to study the structure-function relationships of EPO by in vitro mutagenesis, we have used the vector pGEX-2T to express human and murine EPO fused to the carboxyl terminus of glutathione S-transferase (GST) in E. coli. The fusion proteins were the predicted size (46 kDa) by SDS-PAGE. GST-huEPO eluted from glutathione-agarose using reduced glutathione (GSH) was tested by radioimmunoassay and in a mouse spleen cell assay (MSCA). Dose-response curves parallel to recombinant human EPO (rHuEPO) were obtained in both assays. The ratio of immuno- to bioactivity was 4.7:1. Thus the presence of the 26 kDa GST protein at the end terminus of EPO does not abrogate biological activity. GST-mEPO also gave dose-response curves parallel to rHuEPO in the MSCA but not in the RIA. The wild-type murine and three mutant GST-EPO fusion proteins (166 Des-Arg, Glu 159-->Val, and Arg 163-->Glu) were tested in the MSCA and assayed for GST activity. The ratio of bioactivity to enzyme activity for the Arg 163-->Glu mutant was approximately one third of the value obtained for each of the other fusion proteins, indicating that arginine at 163 is functionally important for EPO activity. The availability of these human and murine gene constructs in pGEX should facilitate site-directed mutagenesis and permit detailed studies of the structure-function relationships for the two erythropoietins.

Auditory pathways in the budgerigar. II. Intratelencephalic pathways

The projections of two telencephalic areas in receipt of projections from the auditory relay nucleus of the thalamus (nucleus ovoidalis) were studied in the budgerigar (Melopsittacus undulatus) with autoradiographic methods. These nuclei are called field 'L' and neostriatum intermedium pars dorsolateralis (NIDL). The results show that neurons in both fields project laterally to a portion of the neostriatum intermedium pars ventrolateralis (NIVL) and rostrally to the rostromedial archistriatum. Horseradish peroxidase experiments confirm these projections and indicate that field 'L', NIDL and NIVL also receive projections from neurons in the hyperstriatum ventrale (HV). The projections of field 'L' and NIDL neurons to the rostromedial archistriatum may act as pathways subserving auditory feedback. Neurons in this portion of the archistriatum project to the contralateral field 'L', NIDL and NIVL. Furthermore, the medial archistriatal projection field of neurons within field 'L', NIDL and NIVL (i.e. rostromedial archistriatum) is located adjacent to a large archistriatal neuronal field projecting to the medulla, including the lateral reticular formation of the medulla. This large archistriatal field includes the nucleus archistriatalis robustus, the telencephalic nucleus identified as the source of projections to the lower motoneurons of the syrinx. Thus, projections from auditory telencephalic areas to the rostromedial archistriatum may serve functions related to processes associated with learning and vocal motor control.