An expressed sequence tag (EST) analysis approach was undertaken to identify major genes involved in cold acclimation of Rhododendron, a broad-leaf, woody evergreen species. Two cDNA libraries were constructed, one from winter-collected (cold-acclimated, CA; leaf freezing tolerance -53 degrees C) leaves, and the other from summer-collected (non-acclimated, NA; leaf freezing tolerance -7 degrees C) leaves of field-grown Rhododendron catawbiense plants. A total of 862 5'-end high-quality ESTs were generated by sequencing cDNA clones from the two libraries (423 from CA and 439 from NA library). Only about 6.3% of assembled unique transcripts were shared between the libraries, suggesting remarkable differences in gene expression between CA and NA leaves. Analysis of the relative frequency at which specific cDNAs were picked from each library indicated that four genes or gene families were highly abundant in the CA library including early light-induced proteins (ELIP), dehydrins/late embryogenesis abundant proteins (LEA), cytochrome P450, and beta-amylase. Similarly, seven genes or gene families were highly abundant in the NA library and included chlorophyll a/b-binding protein, NADH dehydrogenase subunit I, plastidic aldolase, and serine:glyoxylate aminotransferase, among others. Northern blot analyses for seven selected abundant genes confirmed their preferential expression in either CA or NA leaf tissues. Our results suggest that osmotic regulation, desiccation tolerance, photoinhibition tolerance, and photosynthesis adjustment are some of the key components of cold adaptation in Rhododendron.