That benzene causes acute myeloid leukaemia is settled. What is not settled is the shape of the curve at low exposure, and that shape determines every occupational limit set from it. Applying two different historical exposure reconstructions to the same cohort of workers produces a linear relationship under one and a quadratic one under the other. Objectives: To review the published evidence on benzene exposure and adult acute myeloid leukaemia: what risk is reported at high and low cumulative exposure, how the dose-response relationship has been modelled, why estimates differ between cohorts, and what follows for exposure limits.
Material and Methods: Narrative review of occupational cohort studies, case-cohort analyses, meta-analyses, risk assessments and agency evaluations indexed in PubMed and the major occupational health and epidemiology journals. Studies were eligible if they reported the association between quantified benzene exposure and myeloid malignancy in a defined adult population. No new patient data were generated and no participants were recruited; every estimate quoted is one published by the original investigators.
Results: In a case-cohort analysis following 110,631 Chinese workers between 1972 and 1999, risk of myelodysplastic syndrome or acute myeloid leukaemia rose with cumulative exposure, with an excess relative risk per 100 parts per million-years of 0.30 (95% CI −0.02 to 1.09, trend p = 0.08), corresponding to a rate ratio of 1.30. Pooled analysis gave an odds ratio of 1.94 (95% CI 0.95 to 3.95) for cumulative exposure below 40 parts per million-years. Among petroleum distribution workers, myelodysplastic syndrome risk was elevated above 2.93 against below 0.348 parts per million-years, at an odds ratio of 4.33 (95% CI 1.31 to 14.3), an exposure range more than ten times lower than the pooled analysis low category. Significant dose-response for acute myeloid leukaemia has been reported from cumulative exposures just below 1 part per million-year. In the historically pivotal cohort, one exposure reconstruction predicted a linear dose-response and another a quadratic one, with the departure from linearity borderline at p = 0.08. Conclusion: The causal relationship is not in question and the quantitative basis for regulation is weaker than its use implies. Effects are being detected at
exposures well below limits derived from cohorts whose lowest category exceeds them by an order of magnitude. Abbreviations: AML – Acute Myeloid Leukaemia, MDS – Myelodysplastic Syndrome, AMML – Acute Myelomonocytic Leukaemia, ANLL – Acute Non-Lymphocytic Leukaemia, ppm – Parts Per Million, ERR – Excess Relative Risk, RR – Rate Ratio, OR – Odds Ratio, OEL – Occupational Exposure Limit, TWA – Time-Weighted Average, TLV – Threshold Limit Value, LMS – Linearized Multistage, IARC – International Agency for Research on Cancer.
Keywords: Acute Myeloid Leukaemia; Benzene; Occupational Exposure; Dose-Response; Exposure Limits
time-weighted average limit stands at 1 part per million in several jurisdictions, with one advisory body lowering its threshold value to 0.02 parts per million [5,7].
Those limits are derived quantitatively from cohort data by extrapolating a dose-response relationship down from the high exposures historically experienced [8,9]. The extrapolation depends on the assumed shape of the curve, and that shape is where the epidemiology becomes contested rather than settled [10,11].
This review sets out what risk is reported at high and low exposure, how the relationship has been modelled, why cohorts disagree, and what follows for the limits derived from them.
Search strategy: PubMed and the archives of the major occupational health, epidemiology and haematology journals were searched for studies of benzene and myeloid malignancy. Search terms combined benzene with acute myeloid leukaemia, myelodysplastic syndrome, occupational exposure, cumulative exposure, parts per million-years, dose-response and exposure limits. Reference lists of retrieved cohort analyses and agency evaluations were screened by hand.
Eligibility: Studies were included if they reported the association between quantified benzene exposure and myeloid malignancy in a defined adult population, or modelled that relationship for risk assessment. Studies of haematological effects short of malignancy are cited where they bear on the plausibility of low-dose effects and are identified as such.
Data extraction and handling: Cohort, follow-up period, exposure assessment method, exposure metric, lag applied and reported effect estimates were extracted as published. The exposure range covered is stated alongside every estimate, because the central finding of this review is that the word low denotes ranges differing by more than an order of magnitude between literatures. No pooling or re-analysis was performed and no participant-level data were obtained. Figures 1 to 4 and Tables 1 to 3 reproduce published values only. This is a narrative synthesis and was not registered as a systematic review.
The largest prospective evidence: A multicentre case-cohort study followed 110,631 Chinese workers between 1972 and 1999, estimating exposures by hierarchical modelling of occupational factors calibrated against historical routine measurements [3,12]. For the exposure window of two to under ten years, risk of myelodysplastic syndrome or acute myeloid leukaemia rose with cumulative exposure, giving an excess relative risk per 100 parts per million-years of 0.30 (95% CI −0.02 to 1.09, trend p = 0.08) and a corresponding rate ratio of 1.30 (Figure 1) [12].
Cumulative exposure and average intensity were highly correlated at 0.94 in that window, so the two metrics could not be separated and the cumulative results were reported [12]. Across the studies compared, a positive exposure response was observed for all leukaemias or acute myeloid leukaemia, with risks highest at the highest exposures [12].
Results below 400 parts per million-years were generally compatible across studies, given wide confidence intervals reflecting small numbers of leukaemia cases in most cohorts [12]. Above that level, risks in the historically pivotal rubber hydrochloride cohort exceeded those in the Chinese study, probably reflecting higher exposures within that category [12,13].
Earlier analysis of the same Chinese population reported cause-specific mortality and incidence of haematopoietic malignancy across a retrospective cohort design [12], and the original dose-related incidence analysis established the exposure-response that later work refined [20]. Exposure assessment in that cohort has itself been described in detail, which is unusual and valuable [21].
Figure 1: Rate ratio at 100 parts per million-years in a cohort of 110,631 workers [12] What is reported at lower exposure: A systematic review and meta-analysis synthesising the evidence by cumulative exposure category reported a pooled odds ratio of 1.94 (95% CI 0.95 to 3.95) for exposures below 40 parts per million-years (Figure 2) [14]. The interval crosses unity, so the pooled low-exposure estimate is compatible with no effect as well as with a near doubling.
Studies conducted at genuinely low exposure report larger effects. Among petroleum distribution workers, myelodysplastic syndrome risk was elevated in those exposed above 2.93 parts per million-years compared with below 0.348, at an odds ratio of 4.33 (95% CI 1.31 to 14.3) [15,16]. A significant dose response for acute myeloid leukaemia has been reported from cumulative exposures just below 1 part per million-year [17], and an elevenfold risk from exposures well below that level [16,18].
Effects short of malignancy appear lower still. Total white cell and platelet counts were significantly reduced in workers exposed at 1 part per million compared with unexposed controls [16,19], and increased aneuploidy has been reported in germ cells at the same concentration [16].
Whether these lower-exposure findings should carry equal weight is contested. They rest on few cases and wide intervals, and the populations differ from the historical cohorts in co-exposures and in smoking prevalence [14,15]. They are nonetheless the only evidence generated at concentrations workers actually encounter today [5,7].
Exposure Principal reported
range finding Curve shape Rubber depends on Crump 1994 [10] High hydrochloride exposure estimates used Benzene and the Hayes et al. 1997 Chinese dose-related Wide [20] workers incidence of leukaemia Pooled odds ratio Khalade et al. 2010 Under 40 Pooled 1.94, interval [14] ppm-years crossing unity Myelodysplasia Schnatter et al. Petroleum Under 3 elevated at very low 2013 [15] distribution ppm-years exposure Significant dose Stenehjem et al. Offshore Under 1 response below 1 2015 [17] petroleum ppm-year ppm-year Linet et al. 2019 110,631 Rate ratio 1.30 at Wide [12] workers 100 ppm-years Table 1: Principal studies reviewed, with exposure range and the finding as reported by the original investigators What low means: The comparison above is not like for like, and the mismatch is the point. The pooled analysis defines its lowest category as under 40 parts per million-years [14]. The petroleum worker study divides its entire cohort into tertiles at 0.348 and 2.93 parts per million-years, so its highest exposure group sits more than tenfold below the pooled analysis's lowest (Figure 3) [15].
An effect detected in the second study occurs at exposures that the first would classify as effectively zero. Investigators have noted this explicitly, observing that their tertiles are much lower than the low cumulative exposure category used in pooled work and than those in the historically pivotal cohorts [15].
This matters directly for regulation. Occupational limits set at around 1 part per million as an eight-hour average [5,7] correspond, over a working lifetime, to cumulative exposures in the range where effects are now being reported rather than the range from which the
The curve shape is not determined by the data: The historically pivotal cohort showed a significant excess of acute myelocytic or monocytic leukaemia, with 8 to 10 cases observed against 1.61 expected and a strong dose-response trend [10,13]. Quantitative risk estimates proved robust to whether the myelomonocytic endpoint or all leukaemia was modelled, and to which of two exposure reconstructions was used, differing by no more than twofold, provided a linear model was applied [10].
The qualification is decisive. One set of exposure estimates predicted a linear dose response and the other predicted a quadratic one, with the departure from linearity borderline significant at p = 0.08 (Figure 4) [10]. Same workers, same leukaemias, different reconstruction of what they breathed decades earlier, different curve.
A linear curve implies risk at any exposure and no safe threshold; a quadratic one implies risk falling away steeply at low doses and a defensible limit. The two lead to different regulatory conclusions from identical outcome data, and the choice between them rests on retrospective exposure estimation rather than on the leukaemias observed [8,10,11].
Does benzene cause Settled Agency evaluation [1,2] myeloid leukaemia?
Is there a Yes, at high Multiple cohorts [12,20] dose-response? exposure What is the curve Depends on reconstruction Contested shape at low dose? [10,11] Is there a safe Effects reported below 1 Not established threshold? ppm-year [17,18] Do lymphoid No exposure-response for Debated neoplasms occur? lymphoma [3] Table 2: What is settled and what remains contested How benzene causes the disease: The mechanism supports rather than settles the low-dose question. Benzene metabolites damage haematopoietic progenitors through oxidative injury, topoisomerase inhibition and disruption of chromosome segregation, producing the clastogenic effects characteristic of the exposure [27,28]. Chromosome-wide analysis in exposed workers demonstrated aneuploidy across the genome rather than at isolated loci [29].
Susceptibility varies between individuals. Polymorphism in the gene encoding NAD(P)H quinone oxidoreductase, which detoxifies benzene metabolites, associates with haematological toxicity in exposed workers [30], so a population-average dose-response necessarily averages across people at different risk from the same air concentration.
Where the injury produces a clone that persists, the resulting disease follows the pathway established for secondary myeloid neoplasms generally [31], and would be classified and managed under current criteria for myeloid malignancy [32,33]. Benzene sits alongside prior chemotherapy and radiation among the established environmental causes of an otherwise largely idiopathic disease [34]. Other cohorts and other endpoints: Refinery and petrochemical cohorts have been followed for mortality over several decades [25], and hospital-based case-control work has examined personal, lifestyle and environmental risk factors by disease subtype [24]. Analyses of B-lymphocyte malignancy in industrial cohorts address the contested lymphoid question separately [26], where no exposure-response has been demonstrated despite raised counts in exposed workers [3,4].
High-exposure Consistent Exposures no longer cohorts dose-response encountered [12,13] Pooled low-exposure Combines small Interval crosses unity [14] estimate studies Contemporary Relevant to current Few cases, wide intervals low-exposure studies limits [15,17] Mechanistic and Effects at 1 ppm Not malignancy biomarker work demonstrated endpoints [19,29]
Two things are true simultaneously and are frequently confused. Benzene causes acute myeloid leukaemia, and the evidence for that is as strong as occupational epidemiology produces [1,2,20]. The quantitative relationship used to set exposure limits is considerably less secure than its regulatory application suggests [8,10,11].
The insecurity is specific rather than general. At high cumulative exposure the cohorts broadly agree [12,13]. Below roughly 40 parts per million-years the confidence intervals widen enough to include no effect [14], and the extrapolation into that region is where limits are actually derived.
Historical exposure in these cohorts was not measured contemporaneously for most workers. It has been reconstructed decades later from job histories, process descriptions and sparse industrial hygiene samples [10,12,21]. Two competent reconstructions of the same cohort yielded linear and quadratic dose-responses respectively [10].
The outcome data did not change between those analyses; only the assumption about what workers breathed. That places the weight of a regulatory decision on the least well-measured element of the study, and it is why later work has invested heavily in hierarchical exposure modelling calibrated against routine measurements [12,21].
Modern occupational exposures are far below those in the historical cohorts, so the practical question is what happens in a range those cohorts barely covered [5,15]. Studies conducted at that range report effects: myelodysplasia at an odds ratio of 4.33 above 2.93 parts per million-years [15], dose-response for acute myeloid leukaemia below 1 part per million-year [17], and haematological changes at 1 part per million [16,19].
These findings are individually imprecise, with wide intervals reflecting small case numbers [14,15]. They are nonetheless consistent in direction, and they arise in exactly the range that limits were extrapolated into rather than derived from, which is the situation in which extrapolation is least reliable [8,11].
The limitations belong to the underlying literature. Historical exposure was reconstructed rather than measured for most workers in the formative cohorts, and different reconstructions produce different dose-response shapes from identical outcomes [10,21]. Leukaemia case numbers are small in most cohorts, producing wide confidence intervals that limit what low-exposure categories can establish [12,14]. Diagnostic classification changed across the periods covered, and some cases originally recorded as acute myeloid leukaemia were reclassified as myelodysplasia or myeloproliferative disease on review [15,22]. Co-exposure to other solvents is common in the industries studied and is variably adjusted for [15,20]. And smoking is both a benzene source and a leukaemia risk factor, which confounds occupational estimates in cohorts where it is incompletely recorded [6,14].
Two gaps deserve naming. No study reviewed here resolves the curve shape at low exposure, and the question cannot be settled by enlarging cohorts whose exposures were reconstructed by the disputed methods [10,12]. And no analysis reviewed reports myeloid malignancy risk against biological exposure markers measured prospectively, which would remove the reconstruction problem entirely [16,23].
Benzene causes acute myeloid leukaemia, and cohorts of exposed workers agree at high cumulative exposure, with a rate ratio of 1.30 at 100 parts per million-years in the largest prospective study. Below about 40 parts per million-years the pooled estimate is 1.94 with an interval crossing unity, while studies conducted at genuinely low exposure report myelodysplasia at an odds ratio of 4.33 above 2.93 parts per million-years and dose-response for acute myeloid leukaemia below 1 part per million-year, ranges that pooled work classifies as effectively unexposed. The curve shape at low dose, which determines every limit extrapolated from these data, is not fixed by the outcomes: applying two historical exposure reconstructions to the same cohort produced a linear relationship under one and a quadratic one under the other, with the departure from linearity borderline at p = 0.08. The regulatory question turns on the least well-measured element of the evidence.
The authors thank the occupational cohorts and their investigators whose published analyses form the substance of this article.
No dedicated funding from public, commercial or charitable sources supported this review.
The authors declare no competing interests.
This article reviews previously published studies only. No new patient data were collected and no participants were recruited, so ethical approval was not required.