The effects of heavy metals in agricultural soils have received special attention due to their potential for accumulation in crops, which can affect species at all trophic levels. Therefore, there is a critical need for reliable bioassays for assessing risk levels due to heavy metals in agricultural soil. In the present study, we used microarrays to investigate changes in gene expression of Lycopersicon esculentum in response to Cd-, Cr-, Hg-, or Pb-spiked soil. Exposure to (1)/10 median lethal concentrations (LC50) of Cd, Cr, Hg, or Pb for 7 days resulted in expression changes in 29 Cd-specific, 58 Cr-specific, 192 Hg-specific and 864 Pb-specific genes as determined by microarray analysis, whereas conventional morphological and physiological bioassays did not reveal any toxicant stresses. Hierarchical clustering analysis showed that the characteristic gene expression profiles induced by Cd, Cr, Hg, and Pb were distinct from not only the control but also one another. Furthermore, a total of three genes related to "ion transport" for Cd, 14 genes related to "external encapsulating structure organization", "reproductive developmental process", "lipid metabolic process" and "response to stimulus" for Cr, 11 genes related to "cellular metabolic process" and "cellular response to stimulus" for Hg, 78 genes related to 20 biological processes (e.g., DNA metabolic process, monosaccharide catabolic process, cell division) for Pb were identified and selected as their potential biomarkers. These findings demonstrated that microarray-based analysis of Lycopersicon esculentum was a sensitive tool for the early detection of potential toxicity of heavy metals in agricultural soil, as well as an effective tool for identifying the heavy metal-specific genes, which should be useful for assessing risk levels due to heavy metals in agricultural soil.