Elucidation of the complex processes involved in plant growth requires analysis of the spatial gene expression patterns in all affected tissues. This protocol extension is an adaptation of a protocol that describes how to use barcoded oligo-dT microarrays to evaluate spatial global gene expression profiles in mammalian tissue to enable it to be applied to plant material. Here, we explain the required adjustments for preparing and treating plant tissue sections on the array surface, specifically in regard to how to permeabilize and remove the tissue. Once the tissue has been removed, the cDNA-mRNA hybrid that is left on the slide is processed in the same way as cDNA obtained during experiments on mammalian tissue; thus the later stages of the protocol are not included here, and readers should follow the accompanying protocol for those. We have previously used our protocol to generate high-quality sequencing libraries for Arabidopsis thaliana inflorescence, Populus tremula developing and dormant leaf buds, and Picea abies female cones. However, we anticipate that the protocol can be adapted to other tissue types and species. The entire protocol for preparing samples and processing libraries can be completed in 3-4 d.