Even though nicotine has been shown to modulate mRNA expression of a variety of genes, a comprehensive high-throughput study of the effects of nicotine on the tissue-specific gene expression profiles has been lacking in the literature. In this study, cDNA microarrays containing 1117 genes and ESTs were used to assess the transcriptional response to chronic nicotine treatment in rat, based on four brain regions, i.e. prefrontal cortex (PFC), nucleus accumbens (NAs), ventral tegmental area (VTA), and amygdala (AMYG). On the basis of a non-parametric resampling method, an index (called jackknifed reliability index, JRI) was proposed, and employed to determine the inherent measurement error across multiple arrays used in this study. Upon removal of the outliers, the mean correlation coefficient between duplicate measurements increased to 0.978+/-0.0035 from 0.941+/-0.045. Results from principal component analysis and pairwise correlations suggested that brain regions studied were highly similar in terms of their absolute expression levels, but exhibited divergent transcriptional responses to chronic nicotine administration. For example, PFC and NAs were significantly more similar to each other (r=0.7; P<10(-14)) than to either VTA or AMYG. Furthermore, we confirmed our microarray results for two representative genes, i.e. the weak inward rectifier K(+) channel (TWIK-1), and phosphate and tensin homolog (PTEN) by using real-time quantitative RT-PCR technique. Finally, a number of genes, involved in MAPK, phosphatidylinositol, and EGFR signaling pathways, were identified and proposed as possible targets in response to nicotine administration.