Diatoms are unicellular, photosynthetic, eukaryotic algae with a ubiquitous distribution in water environments and they play an important role in the carbon cycle. Molecular or morphological changes in these species under ecological stress conditions are expected to serve as early indicators of toxicity and can point to a global impact on the entire ecosystem. Thalassiosira pseudonana, a marine diatom and the first with a fully sequenced genome has been selected as an aquatic model organism for ecotoxicological studies using molecular tools. A customized DNA microarray containing probes for the available gene sequences has been developed and tested to analyze the effects of a common pollutant, benzo(a)pyrene (BaP), at a sub-lethal concentration. This approach in diatoms has helped to elucidate pathway/metabolic processes involved in the mode of action of this pollutant, including lipid metabolism, silicon metabolism and stress response. A dose-response of BaP on diatoms has been made and the effect of this compound on the expression of selected genes was assessed by quantitative real time-PCR. Up-regulation of the long-chain acyl-CoA synthetase and the anti-apoptotic transmembrane Bax inhibitor, as well as down-regulation of silicon transporter 1 and a heat shock factor was confirmed at lower concentrations of BaP, but not the heat-shock protein 20. The study has allowed the identification of molecular biomarkers to BaP to be later on integrated into environmental monitoring for water quality assessment.