Blueberry (Vaccinium spp.) is a small berry with high economic value. Although cold storage can extend the storage time of blueberry to more than 60 days, it leads to chilling injury (CI) displaying as pedicle pits; and the samples of 0 °C-30 days was the critical point of CI. However, little is known about the mechanism and the molecular basis response to cold stress in blueberry have not been explained definitely. To comprehensively reveal the CI mechanisms in response to cold stress, we performed high-throughput RNA Seq analysis to investigate the gene regulation network in 0d (control) and 30d chilled blueberry. At the same time, the pitting and decay rate, electrolyte leakage (EL), malondialdehyde (MDA) proline content and GSH content were measured.Two cDNA libraries from 0d (control) and 30d chilled samples were constructed and sequenced, generating a total of 35,060 unigenes with an N50 length of 1348 bp. Of these, 1852 were differentially expressed, with 1167 upregulated and 685 downregulated. Forty-five cold-induced transcription factor (TF) families containing 1023 TFs were identified. The DEGs indicated biological processes such as stress responses; cell wall metabolism; abscisic acid, gibberellin, membrane lipid, energy metabolism, cellular components, and molecular functions were significantly responsed to cold storage. The transcriptional level of 40 DEGs were verified by qRT-PCR.The postharvest cold storage leads serious CI in blueberry, which substantially decreases the quality, storability and consumer acceptance. The MDA content, proline content, EL increased and the GSH content decreased in this chilled process. The biological processes such as stress responses, hormone metabolic processes were significantly affected by CI. Overall, the results obtained here are valuable for preventing CI under cold storage and could help to perfect the lack of the genetic information of non-model plant species.
