Biogeochemical mapping of selenium in Italian agricultural soils was accomplished by measuring the Se concentration of representative samples of wheat grains from 71 provinces. The range of the concentration values averaged on a provincial basis was 7-245 ng Se g(-1). A multiple regression model based on six geochemical and pedoclimatic variables was developed to interpret the observed data and to predict Se concentration of wheat in areas where analytical data were missing and in the different Italian soil regions. The statistical model explained only part of the observed variance, but succeeded in identifying Se-enriched as well as Se-depleted areas with an acceptable level of agreement with the biogeochemical map based on measured Se in wheat. Furthermore, the model showed that within the range of concentrations measured in Italian soils, Se-bioaccessibility is controlled not only by the Se content of the soil parent rocks, but also by their overall geochemical nature (carbonatic vs. silicatic) and by pedoclimatic variables (temperature and rain intensity excursions) related to fluctuations of soil moisture and pH. Overall, several Se-marginal and Se-deficient areas were identified on the Italian territory. The implications of these findings for public health are discussed briefly.