Winter rapeseed (Brassica rapa L.) is a major oilseed crop in Northern China, where its production was severely affected by chilling and freezing stress. However, not much is known about the role of differentially accumulated proteins (DAPs) during the chilling and freezing stress. In this study, isobaric tag for relative and absolute quantification (iTRAQ) technology was performed to identify DAPs under freezing stress. To explore the molecular mechanisms of cold stress tolerance at the cellular and protein levels, the morphological and physiological differences in the shoot apical meristem (SAM) of two winter rapeseed varieties, Longyou 7 (cold-tolerant) and Lenox (cold-sensitive), were explored in field-grown plants. Compared to Lenox, Longyou 7 had a lower SAM height and higher collar diameter. The level of malondialdehyde (MDA) and indole-3-acetic acid (IAA) content was also decreased. Simultaneously, the soluble sugars (SS) content, superoxide dismutase (SOD) activity, peroxidase (POD) activity, soluble protein (SP) content, and collar diameter were increased in Longyou 7 as compared to Lenox. A total of 6330 proteins were identified. Among this, 98, 107, 183 and 111 DAPs were expressed in L7 CK/Le CK, L7 d/Le d, Le d/Le CK and L7 d/L7 CK, respectively. Quantitative real-time PCR (RT-qPCR) analysis of the coding genes for seventeen randomly selected DAPs was performed for validation. These DAPs were identified based on gene ontology enrichment analysis, which revealed that glutathione transferase activity, carbohydrate-binding, glutathione binding, metabolic process, and IAA response were closely associated with the cold stress response. In addition, some cold-induced proteins, such as glutathione S-transferase phi 2(GSTF2), might play an essential role during cold acclimation in the SAM of Brassica rapa. The present study provides valuable information on the involvement of DAPs during cold stress responses in Brassica rapa L, and hence could be used for breeding experiments.
