ABSTRACT: Cowpea (V. unguiculata L. Walp.) is an important tropical grain legume. Asparagus bean (V. unguiculata ssp. sesquipedialis) is a distinctive subspecies of cowpea, which is considered one of the top ten Asian vegetables. It can be adapted to a wide range of environmental stimuli such as drought and heat. Nevertheless, it is an extremely cold-sensitive tropical species. Improvement of chilling tolerance in asparagus bean may significantly increase its production and prolong its supply. However, gene regulation and signaling pathways related to cold response in this crop remain unknown. Using Illumina sequencing technology, modification of global gene expression in response to chilling stress in two asparagus bean cultivars-"Dubai bean" and "Ningjiang-3", which are tolerant and sensitive to chilling, respectively-were investigated. More than 1.8 million clean reads were obtained from each sample. After de novo assembly, 88,869 unigenes were finally generated with a mean length of 635 bp. Of these unigenes, 41,925 (47.18%) had functional annotations when aligned to public protein databases. Further, we identified 3,510 differentially expressed genes (DEGs) in Dubai bean, including 2,103 up-regulated genes and 1,407 down-regulated genes. While in Ningjiang-3, we found 2,868 DEGs, 1,786 of which were increasing and the others were decreasing. 1,744 DEGs were commonly regulated in two cultivars, suggesting that some genes play fundamental roles in asparagus bean during cold stress. Functional classification of the DEGs in two cultivars using Mercator pipeline indicated that RNA, protein, signaling, stress and hormone metabolism were five major groups. In RNA group, analysis of TFs in DREB subfamily showed that ICE1-CBF3-COR cold responsive cascade may also exist in asparagus bean. Our study is the first to provide the transcriptome sequence resource for asparagus bean, which will accelerate breeding cold resistant asparagus bean varieties through genetic engineering, and advance our knowledge of the genes involved in the complex regulatory networks of this plant under cold stress.
