Discovery of cold-resistance genes in Vitis amurensis using bud-based quantitative trait locus mapping and RNA-seq
BMC Genomics volume 23, Article number: 551 (2022)
In cold regions, low temperature is the main limiting factor affecting grape production. As an important breeding resource, V. amurensis Rupr. has played a crucial role in the discovery of genes which confer cold resistance in grapes. Thus far, many cold-resistance genes have been reported based on the study of V. amurensis. In order to identify more candidate genes related to cold resistance in V. amurensis, QTL mapping and RNA-seq was conducted based on the hybrid population and different cold-resistance cultivars in this study.
In this study, highly cold-resistant grape cultivar ‘Shuangyou’ (SY) which belongs to V. amurensis, and cold-sensitive cultivar ‘Red Globe’ (RG) which belongs to Vitis vinifera L. were used to identify cold resistance genes. Cold-resistance quantitative trait locus (QTL) mapping was performed based on genetic population construction through interspecific crossing of ‘Shuangyou’ and ‘Red Globe’. Additionally, transcriptome analysis was conducted for the dormant buds of these two cultivars at different periods. Based on transcriptome analysis and QTL mapping, many new structural genes and transcription factors which relate to V. amurensis cold resistance were discovered, including CORs (VaCOR413IM), GSTs (VaGST-APIC, VaGST-PARB, VaGSTF9 and VaGSTF13), ARFs (VaIAA27 and VaSAUR71), ERFs (VaAIL1), MYBs (VaMYBR2, VaMYBLL and VaMYB3R-1) and bHLHs (VaICE1 and VabHLH30).
This discovery of candidate cold-resistance genes will provide an important theoretical reference for grape cold-resistance mechanisms, research, and cold-resistant grape cultivar breeding in the future.
Grape is one of the most important fruit species in the world and it is also the main cultivated fruit crop in China. China is located in a region with a typical continental monsoon climate. The grape cultivars in most grape production regions are sensitive to cold. The grapevines should be buried in soil when the temperature falls below -15 ℃ to survive the cold environment in winter. This not only increases the production cost of the vineyard, but also causes a series of problems such as tree body damage, soil erosion, and soil horizon destruction that considerably limited the development of grape production. Many grape species such as Vitis riparia Michx, Vitis rupestris Scheele and Vitis labrusca L. indigenous to North America, Vitis amurensis Rupr. and Vitis yeshanensis J.X.Chen from Eastern Asia possess strong resistance to cold conditions . The cold resistance of V. amurensis is the highest among resistant species, and its cane can survive in temperatures of -40 to -50 ℃. Therefore, it is an excellent germplasm resource for cold-resistance breeding as seed-parent.
Thus far, there have been many reports related to genetic and transcriptional regulation mechanisms of grape cold resistance based on V. amurensis. Interaction of VaPAT1 and VaIDD3 in V. amurensis could enhance the expression of VaLOX, which plays a crucial role in jasmonic acid synthesis and improves cold resistance . Genome-wide association analysis shows that phosphoglycerate kinase gene (PGK) in V. amurensis significantly influences cold resistance . A number of loci and candidate genes related to grapevine cane cold-hardiness was discovered based on the hybrid population constructed using ‘Cabernet Sauvignon’ and ‘Zuoyouhong’ . Ethylene response factor VaERF092 improves cold resistance in V. amurensis by binding to the GCC-box element in the promoter region of VaWRKY33, thus promoting its expression . Furthermore, candidate cold-resistance genes identified based on V. amurensis were reported, including VaDof17d , VaAQUILO , VaERD15 , and VaWRKY12 . Candidate cold-resistance genes in V. vinifera were discovered, including VvNAC17 , VV-circATS1 , and VvBAP  and the function of these genes has been preliminarily verified in Arabidopsis thaliana. In summary, V. amurensis plays an important role in the discovery of cold-resistance genes. However, plants exposed to chilling (0 – 15 ℃) and freezing (< 0 ℃) stress show differing gene expression patterns [13, 14], and most previous studies have focused on chilling stress.
In this study, based on the hybrid population created through interspecific crossing of cold-resistance cultivar ‘Shuangyou’ (SY) and cold-sensitive cultivar ‘Red Globe’ (RG) and the high-density genetic linkage map , we conducted QTL mapping for cold resistance of grape dormant buds firstly and then transcriptome analysis was conducted for ‘SY’ and ‘RG’ buds at different dormant stages. Finally, candidate genes related to cold resistance of grape buds were screened by QTL mapping and RNA-seq.
Identification of cold resistance
Cold resistance identification of ‘SY’ and ‘RG’ buds at four different dormant stages was evaluated by Differential Thermal Analysis (DTA) analysis. The cold resistance of the buds of both cultivars increased with time when measured over 60 days after dormancy (DAD) began. At 60DAD, the cold resistance of the grape cultivars ‘SY’ and ‘RG’ reached the peak value of the heat release curve, and the cold resistance of ‘SY’ was significantly higher than ‘RG’ at 90 days after dormancy. The cold resistance of ‘RG’ was significantly reduced (Fig. 1a). Furthermore, 149 hybrid progenies of ‘SY’ and ‘RG’ were identified for cold resistance after dormancy, and the results showed continuous variation (Fig. 1b). These results indicated that the cold resistance of grape cultivar ‘SY’ was significantly higher than that of grape variety ‘RG’ and the cold resistance in grape was a typical quantitative trait controlled by multiple genes. Finally, eight individuals which showed greater cold resistance were selected.
Transcriptome and differentially expressed gene analysis
In order to discover the candidate genes for cold resistance in V. amurensis, we conducted RNA-seq analysis for the dormant buds of grape varieties ‘SY’ and ‘RG’ at 0 days and 60 days after dormancy. After removing low-quality reads and adapters, a total of 77.43 Gb clean data was obtained for subsequent research, and the average clean data of each sample was 6.45 Gb. The clean data obtained by RNA-seq has been uploaded to NCBI Sequence Read Archive (SRA) with the Accession Number PRJNA787359. In total 43,047 unigenes were obtained, including 42,413 annotated genes and 634 unannotated unigenes (Supplementary file 1). These 643 unannotated unigenes were then aligned with Swiss-Prot, GO, KEGG and Pfam databases by using BLAST and HMMER software and finally annotated in Supplementary file 2. The FPKM value which was calculated by the comparison of sequenced reads with obtained RNA-seq database represents the expression of each unigene (Supplementary file 3). To confirm the reliability and rationality of the experiment, we calculated the Pearson’s correlation coefficient for all genes expression levels between each sample and refraction based on their FPKM value distribution. These coefficients are presented in the form of a correlation matrix map (Fig. 2a). A total of 6413 genes were differentially expressed in SY0d vs. RG0d (|[log2 (fold change]|) > 1 and adjusted P < 0.05) after differential expression analysis, among which 3644 were upregulated and 2769 were downregulated; 6312 genes were differentially expressed in SY60d vs. RG60d, among which 4002 were upregulated and 2310 were downregulated. A total of 3136 differentially expressed genes were found between SY0d vs. RG0d and SY60d vs. RG60d (Fig. 2b, c, and Supplementary file 4).
Cold resistance genes discovery based on RNA-seq
In order to cope with low temperature stress, plants produce a series of complex physiological and biochemical reactions, among which cold regulated proteins (CORs) play a crucial role in the plant response to cold stress. In this study, transcriptome analysis was conducted on the buds of ‘SY’ and ‘RG’ at different dormant periods and a cold-resistance gene, COR413IM (Vitvi10g00719) was discovered. The qRT-PCR analysis showed that after cold treatment, expression of VaCOR413IM was significantly induced and upregulated both in ‘SY’ and ‘RG’, the expression level in ‘SY’ was significantly higher than in ‘RG’ at 6 h after cold treated (Fig. 3).
To further explore candidate genes related to cold resistance in grapes, we conducted KEGG and GO enrichment analysis for 6312 differentially expressed genes in SY60d vs. RG60d. KEGG enrichment analysis showed that 225 genes were significantly enriched in the following pathways: flavonoid biosynthesis, glutathione metabolism, circadian rhythms-plant, phenylpropanoid biosynthesis, flavone and flavanol biosynthesis, linoleic acid metabolism, and glycosaminoglycan degradation (adjusted P < 0.05) (Fig. 4a). Among these, the expression levels of 28 genes in SY60d were significantly higher than those in RG60d (Fig. 4b and Supplementary file 5). Many studies have shown that hormone signal transduction plays an important role in plant cold resistance. KEGG analysis revealed that the expression of 16 genes was significantly higher in SY60d than in RG60d (Fig. 4c and Supplementary file 6). After analyzing the differential expression of selected genes in different comparison groups (RG60d vs RG0d, SY60d vs SY0d and SY0d vs RG0d) with |[log2FC|> 1 and adjusted Pvalue < 0.05, we preliminarily screened six genes for differential expression verification (Table 1). The results showed that the expression of GST-APIC (Vitvi07g00283), VaGST-PARB (Vitvi07g02190), VaGSTF9 (Vitvi12g00080) and VaGSTF13(Vitvi07g02189) in ‘SY’ were significantly upregulatedafter cold treatment, while there was almost no significant difference in ‘RG’ (P < 0.05); IAA27 gene in ‘SY’ was significantly upregulated after cold treatment while downregulated in ‘RG’. SAU71 was significantly induced both in ‘SY’ and ‘RG’ after cold treatment and the expression level in ‘SY’ was significantly higher than that in ‘RG’ at 4 h (Fig. 5).
According to GO enrichment analysis, 2012 genes were significantly enriched in 30 catalogues, which were classified into categories of ‘Molecular function’, ‘Cellular component’, and ‘Biological process’ (Fig. 6a). In order to further explore candidate genes involved in transcriptional regulation mechanisms of cold resistance in grape, we focused on the analysis of genes in ‘DNA binding transcription factor activity’ catalogued under ‘Molecular function’. Finally, we identified 26 genes in which the expression level in SY60d was higher than in RG60d (Fig. 6b and Supplementary file 7). After analyzing the differential expression of selected genes in different comparison groups (RG60d vs RG0d, SY60d vs SY0d and SY0d vs RG0d), we preliminarily screened 1 gene for differential expression verification (Table 1). The results showed that ethylene response factor AIL1 was significantly upregulated in ‘SY’ and ‘RG’, and the expression level in ‘SY’ was significantly higher than in ‘RG’ (Fig. 7).
Cold resistance genes discovery based on QTL mapping
Based on the cold identification of hybrid offspring in 2021 and our constructed genetic linkage map, we conducted QTL mapping to further discover the candidate genes related to grape bud cold resistance. Four QTLs related to grape cold resistance were identified on LG2, LG3, LG6 and LG17, the phenotypic variation they explained ranged from 0.56% to 8.18% and the confidence intervals were 670 172–710 587, 488 381–524 863, 251 863–295 208, and 345 640–383 091, respectively (Table 2, Fig. 8a). Based on these confidence intervals, 20 candidate genes were discovered (Fig. 8b and Supplementary file 8). After analyzing the differential expression of selected genes in different comparison groups (RG60d vs. RG0d, SY60d vs. SY0d, and SY0d vs. RG0d), we preliminarily screened three genes for differential expression verification (Table 1). The results showed that the MYB transcription factor MYBR2 was significantly upregulated in ‘SY’ while there was no significant difference in ‘RG’. The expression level in ‘SY’ was significantly higher than in ‘RG’ after cold treatment. Auxin response factor ARF2A was downregulated in ‘SY’ while no significant difference was seen in ‘RG’. Transcription factor ICE1 was upregulated both in ‘SY’ and ‘RG’, the expression level in ‘SY’ was significantly higher than that in ‘RG’ (Fig. 9).
Transcriptional regulation factors involved in grape buds cold resistance discovery
Additionally, we conducted analysis of ERF, bHLH, MYB, and WRKY transcription factors for their roles in resistance to environment stress. A total of 33 (ERF), 34 (bHLH), 42 (MYB) and 43 (WRKY) genes were found (Fig. 10). After analyzing the differential expression of selected genes in different comparison groups (RG60d vs. RG0d, SY60d vs. SY0d and SY0d vs. RG0d), we preliminarily screened seven genes for differential expression verification (Table 1). The results showed that the expression of WRKY72 was not different in ‘SY’ and ‘RG’ at different cold treatment stages. bHLH transcription factor bHLH68 was downregulated in ‘SY’ and upregulated in ‘RG’, and the expression level in ‘SY’ was significantly higher than in ‘RG’. bHLH156 was downregulated in ‘SY’ and there was no significant difference in the expression of ‘RG’, the expression level in ‘SY’ was significantly higher than that in ‘RG’. bHLH059 was upregulated in ‘SY’ and ‘RG’, the expression level in ‘RG’ was significantly higher than in ‘SY’. bHLH30 was upregulated in ‘SY’ and there was no significant difference in expression in ‘RG’, the expression level in ‘SY’ was significantly higher than that in ‘RG’. MYB transcription factor MYBLL was upregulated both in ‘SY’ and ‘RG’, the expression level in ‘SY’ was significantly higher than in ‘RG’. MYB3R-1 was upregulated in ‘SY’ and there was no significant difference in ‘RG’, the expression level in ‘SY’ was significantly higher than in ‘RG’ (Fig. 11).
In this study, some individuals showed higher cold resistance than the male parent ‘SY’. The formation of heterobeltiosis may be the additive effects of several desired dominant alleles or the combined effect of different alleles at the same gene locus, or a combination of both  and the genetic differences between parents are the primary cause of it. Heterobeltiosis can help breeders efficiently screen for superior parents and predict the heterosis of parental combinations [17, 18]. Transgressive offspring in our study provided important grape cold resistance resources and they can also be used as material for underlying genetic and molecular mechanisms of grape cold resistance.
Thus far, several COR genes have been identified in plants, including COR27, COR47, COR78, COR15A and COR6.6 [19,20,21,22,23,24,25,26,27,28]. COR413, which was discovered in Arabidopsis, belongs to the COR subfamily and contains three proteins including COR413-plasma membrane (COR413PM), COR413-thylakoid membrane (COR413TM), and COR413-inner membrane (COR413IM) protein [29,30,31]. At present, COR413 proteins have been identified in rice, wheat, sorghum, and other species [29, 30, 32, 33]. As an important member of the COR gene family, no reports have been reported on the cold resistance properties of COR413 in fruit trees and the functional research of COR413 has mainly been focused on tomato [20, 34]. In our study, we first reported the potential cold resistance properties of VaCOR413IM in V. amurensis and have provided an important reference for future research of COR413 in grape cold resistance. Furthermore, as a transcription factor that positively regulates COR gene expression, INDUCER OF CBF EXPRESSION (ICE1) plays an important role in plant cold resistance . Previously, VaICE1 in V. amurensis has been successfully cloned and its function has been verified in Arabidopsis thaliana and tobacco [36, 37]. In our study, we discovered a new VaICE1 gene on LG17 in V. amurensis through cold-resistance QTL mapping and RNA-seq, providing a new reference for further research on the regulation of ICE1 in grape cold resistance.
Under cold stress, plants will produce a large number of reactive oxygen species (ROS) which can cause serious injury to their growth and development . Glutathione S transferase (GST) can scavenge the ROS in plant tissue , therefore, playing an important role in improving cold resistance in plants. Multiple GST genes in pumpkin are induced by cold environments . PtrERF9 improves cold resistance by promoting the expression of PtrGSTU17 and reducing ROS content in trifoliate orange under cold stress . Until now, GST genes involved in grape cold resistance have not been reported. In our study, expression of four GST genes including VaGST-APIC, VaGST-PARB, VaGSTF9, and VaGSTF13 were significantly induced and upregulated in the cold-resistant grape cultivar ‘SY’, indicating their potential role in grape cold resistance. Future studies should aim to verify their function.
Auxin response factor (ARF) and ethylene responsive factor (ERF) play important roles in plant resistance to abiotic stress [42,43,44,45,46,47,48,49]. In wheat, TaARF8, TaARF9, and TaARF21 are involved in the response to cold stress . IAA14 in Arabidopsis thaliana plays an important role in mediating auxin synthesis and cold resistance . Additionally, several ARF genes in sorghum are upregulated under cold stress . In this study, three ARF genes IAA27, SAUR71, and ARF2A were discovered after QTL mapping and KEGG enrichment analysis. Expression of IAA27 and SAUR71 was significantly upregulated in ‘SY’ under cold stress, indicating that they could be used as candidate genes for further research into the mechanism of auxin in regulating grape cold resistance. In recent years, an increasing number of studies have reported the important regulation effect of ERF in the cold resistance of fruit trees [41, 48, 52, 53]. ERF genes involved in grape cold resistance have also been reported, including VaERF057  and VaERF092 . In this study, an ethylene response factor VaAIL1 which may positively regulate grape cold resistance was discovered, and the specific mechanism remains to be further explored.
As the largest transcription factor family, many reports have showed the role of MYB in plant cold resistance such as: MdMYB88/124 , MdMYB23 , and MdMYB308L  in apple; PbrMYB5  in pear; and DgMYB2  in Chrysanthemum. Studies of describing the role of MYB transcription factors in grape cold resistance have been reported . In this study, three MYB transcription factors VaMYBR2, VaMYBLL, and VaMYB3R-1, which may be involved in V. amurensis cold resistance were preliminarily screened through QTL mapping and transcriptome analysis, providing an important reference for further studies on MYB transcription factors involved in regulating grape cold resistance. Moreover, bHLH transcription factors MdbHLH3 , MdbHLH33 , and VabHLH1 ; and the WRKY transcription factors CsWRKY46 , SmWRKY26, and SmWRKY32  play important roles in plant cold resistance. In this study, several bHLH transcription factors including bHLH68, bHLH156, bHLH059, bHLH30, and a WRKY transcription factor WRKY72 were discovered based on the transcriptome analysis of different grape cultivars under cold treatment. After qRT-PCR verification, the expression level of VabHLH30 was significantly induced and upregulated in ‘SY’, providing the possibility of bHLH30 participating in the cold resistance of V. amurensis. The mechanisms related to its regulations require further study.
Until now, there have been no reports of these candidate genes in relation to cold resistance. Interestingly, previous studies have shown potential interactions of these candidate genes in cold resistance. For example, the interaction between MdMYB308L and MdbHLH33 could positively regulate the cold resistance in apple  and PtrERF9 acts downstream of ethylene signaling and functions positively in cold tolerance via regulating PtrGSTU17 in Citrus .
Based on the QTL mapping and transcriptome analysis, we discovered five structural genes and eight transcription factors including VaCOR413IM (Vitvi10g00719), VaGST-APIC (Vitvi07g00283), VaGST-PARB (Vitvi07g02190), VaGSTF9 (Vitvi12g00080), VaGSTF13 (Vitvi07g02189), VaIAA27 (Vitvi11g00301), VaSAUR71 (Vitvi08g02133), VaAIL1 (Vitvi09g00108), VaMYBR2 (Vitvi06g00024), VaMYBLL (Vitvi17g00238), VaMYB3R-1 (Vitvi11g00320), VaICE1 (Vitvi17g00037) and VabHLH30 (Vitvi01g00232). The candidate genes identified in our study will provide an important reference for research into cold-resistance mechanisms and breeding in grape species.
Grape cultivar ‘Shuangyou’(SY) (V. amurensis Rupr.) and ‘Red Globe’(RG) (V. vinifera L.) and their hybrid population were cultivated in the Grape Experimental Garden of Shenyang Agricultural University (23°24’N, 41°50’E), China. In our study, October 20 in 2021 was defined as the beginning of bud dormancy, Dormant buds of ‘SY’ and ‘RG’ at 0, 30, 60 and 90 days after dormancy (DAD) were used for cold resistance identification. In order to discover candidate genes related to cold resistance, bud samples of ‘SY’ and ‘RG’ were collected at 0 DAD and 60 DAD in 2021. Three biological replicates were collected for each cultivar with at least 20 buds per replicate, the selected buds were frozen in liquid nitrogen and stored until preparation for RNA-seq analysis. The hybrid population was created in 2009 with ‘SY’ as the male parent and ‘RG’ as the female parent. Canes with buds from 149 individuals of the hybrid population were randomly collected, five canes (at least 15 buds) in each replicate were collected from three biological replicates of each individual in October 20, 2021. The canes were then placed in a refrigerator at 0 ℃ for two months to induce dormancy in the buds on the canes. After the buds were completely dormant, they were then used as materials for cold-resistance gene identification.
Bud cold resistance identification
Differential Thermal Analysis (DTA) was used to determine the cold resistance of dormant buds at 0DAD, 30DAD, 60DAD and 90DAD in a program-controlled cooling climate chamber (GDJS-100, Hefei Anke Environmental Testing Equipment Co., LTD., Hefei City). Full buds were removed from canes using flat cuts and then placed in a dish with moist blotting paper. The buds were placed with the flat cut downward into the thermoelectric modules (TEM) of the Tenney environmental test chamber and each module contained at least 15 buds. The methodology was as follows: temperature was decreased from room temperature to 4 °C, held at this temperature for 1 h (h); the temperature was decreased from 4 °C to -40 °C at a rate of 4 °C/h, then held at -40 °C for 1 h, and returned to room temperature over a 1 h period. The Data Loggers (CR300, Cambpell Scientific, Logan,UT) with Expansion Peripherals (AM15/32B, Cambpell Scientific, Logan, UT) is connected with the temperature sensor, scanning and collecting data every 10 s, the accuracy of which is ± 0.5℃. Meanwhile, the collected data is displayed synchronously on the computer monitor. Through the LoggerNet4.0 data collector management software (Campbell Scientific, Logan, UT), the process of the test can be observed in real time. The data is amplified and transmitted by signal amplification system and stored in Excel form to SD card. Origin2021 was used to draw the cooling curve. Cold treatment leads to an instantaneous rise in temperature and the peak value of the heat release curve is the lethal temperature of dormant buds (LTDB) .
Dormant bud transcriptome analysis
RNA integrity was assessed using the RNA Nano 6000 Assay Kit and the Bioanalyzer 2100 system (Agilent Technologies, Santa Clara, CA, USA). The input material for the RNA sample preparation was 1 μg RNA per sample. Sequencing libraries were generated using the NEBNext® Ultra™ RNA Library Prep Kit (New England Biolabs, Ipswich, MA, USA) and then sequenced on an Illumina Novaseq platform, finally, 150 bp paired-end reads were generated. Clean reads were obtained by removing reads containing adapter, ploy-N and low-quality reads from the raw data. The high quality and paired-end clean reads were aligned to the reference genome (https://www.ncbi.nlm.nih.gov/genome/401) using HISAT 2v2.0.5 software and the mapped reads of each sample were assembled by StringTie. The fragments per kilobase per million (FPKM) of each gene was calculated based on the length of the gene and the number of reads mapped to this gene. Differential expression analysis was performed using the DESeq2 R package (1.20.0) and genes with an adjusted P-value < 0.05 found by DESeq2 were assigned as differentially expressed. Gene Ontology (GO) enrichment analysis of differentially expressed genes was implemented by the cluster Profiler R package. Enrichment analysis of pathways in KEGG database [65,66,67] were performed using a hyper geometric test.
QTL mapping of dormant bud
The LTDB (mean value of three replicates) of each genotype collected in 2021 was used for QTL mapping. A multiple QTL mapping (MQM) method was used to find significant QTLs after a 1000-permutation test (α = 0.05) based on the R/qtl package . The max.qtl was set to 10 for forward selection. A 1-LOD confidence interval corresponding to the 95% confidence interval was calculated by using the “lodint” function. The explained phenotypic variation of each QTL (PEV) was estimated using the “fitqtl” function. The candidate gene search was conducted based on the confidence interval of each QTL on the integrated map. Genes were selected according to the 12X.v3 version of the Grape Genome database (https://urgi.versailles.inra. fr/Species/Vitis/Annotations).
qRT-PCR validation of candidate genes
Dormant bud (60 DAD) of grape cultivars ‘SY’ and ‘RG’ which treated at -16 °C for 2 h, 4 h, and 6 h in program-controlled cooling climate chamber and then these samples were used for total RNA extraction according to the manufacturer’s instructions of Plant Total RNA Isolation Kit (SK8631; Sangon Biotech, Shanghai, China). The PrimeScript™ RT-PCR Kit (RR047A; TaKaRa Bio, Kusatsu, Japan) was used to conduct cDNA synthesis and the cDNA was diluted five times. Quantitative real-time PCR (qRT-PCR) was conducted in ABI QuantStudio™ 6 Flex System (Applied Biosystems). The relative expression level of selected genes was normalized to grapevine β-actin  and calculated using the 2−ΔΔCT method. All reactions were performed using three biological replicates. The primers used in this study are listed in Supplementary file 9.
Availability of data and materials
The data supporting the results presented in this article are included as additional files. The raw sequencing data were deposited at the NCBI Sequence Read Archive (SRA) with the accession Number PRJNA787359.
Quantitative trait loci
Logarithm of odds
Phosphoglycerate kinase gene
Sequence Read Archive
Cold regulated proteins
Reactive oxygen species
Glutathione S transferase
Auxin response factor
Ethylene responsive factor
Fennell A. Freezing tolerance and injury in grapevines. J Crop Improv. 2004;10(1–2):201–35.
Wang Z, Wong DCJ, Wang Y, Xu G, Ren C, Liu Y, Kuang Y, Fan P, Li S, Xin H, et al. GRAS-domain transcription factor PAT1 regulates jasmonic acid biosynthesis in grape cold stress response. Plant Physiol. 2021;186(3):1660–78.
Wang Y, Xin H, Fan P, Zhang J, Liu Y, Dong Y, Wang Z, Yang Y, Zhang Q, Ming R, et al. The genome of Shanputao (Vitis amurensis) provides a new insight into cold tolerance of grapevine. Plant J. 2021;105(6):1495–506.
Su K, Xing H, Guo Y, Zhao F, Liu Z, Li K, Li Y, Guo X. High-density genetic linkage map construction and cane cold hardiness QTL mapping for Vitis based on restriction site-associated DNA sequencing. BMC Genomics. 2020;21(1):419.
Sun X, Zhang L, Wong DCJ, Wang Y, Zhu Z, Xu G, Wang Q, Li S, Liang Z, Xin H. The ethylene response factor VaERF092 from Amur grape regulates the transcription factor VaWRKY33, improving cold tolerance. Plant J. 2019;99(5):988–1002.
Wang Z, Wang Y, Tong Q, Xu G, Xu M, Li H, Fan P, Li S, Liang Z. Transcriptomic analysis of grapevine Dof transcription factor gene family in response to cold stress and functional analyses of the VaDof17d gene. Planta. 2021;253(2):55.
Sun X, Matus JT, Wong DCJ, Wang Z, Chai F, Zhang L, Fang T, Zhao L, Wang Y, Han Y, et al. The GARP/MYB-related grape transcription factor AQUILO improves cold tolerance and promotes the accumulation of raffinose family oligosaccharides. J Exp Bot. 2018;69(7):1749–64.
Yu D, Zhang L, Zhao K, Niu R, Zhai H, Zhang J. VaERD15, a transcription factor gene associated with cold-tolerance in Chinese Wild Vitis amurensis. Front Plant Sci. 2017;8:297.
Zhang L, Zhao T, Sun X, Wang Y, Du C, Zhu Z, Gichuki DK, Wang Q, Li S, Xin H. Overexpression of VaWRKY12, a transcription factor from Vitis amurensis with increased nuclear localization under low temperature, enhances cold tolerance of plants. Plant Mol Biol. 2019;100(1–2):95–110.
Ju YL, Yue XF, Min Z, Wang XH, Fang YL, Zhang JX. VvNAC17, a novel stress-responsive grapevine (Vitis vinifera L.) NAC transcription factor, increases sensitivity to abscisic acid and enhances salinity, freezing, and drought tolerance in transgenic Arabidopsis. Plant Physiol Biochem. 2020;146:98–111.
Gao Z, Li J, Luo M, Li H, Chen Q, Wang L, Song S, Zhao L, Xu W, Zhang C, et al. Characterization and cloning of grape circular RNAs identified the cold resistance-related Vv-circATS1. Plant Physiol. 2019;180(2):966–85.
Hou L, Zhang G, Zhao F, Zhu D, Fan X, Zhang Z, Liu X. VvBAP1 is involved in cold tolerance in Vitis vinifera L. Front Plant Sci. 2018;9:726.
Liu J, Shi Y, Yang S. Insights into the regulation of C-repeat binding factors in plant cold signaling. J Integr Plant Biol. 2018;60(9):780–95.
Ding Y, Shi Y, Yang S. Molecular regulation of plant responses to environmental temperatures. Mol Plant. 2020;13(4):544–64.
Guo Y, Shi G, Liu Z, Zhao Y, Yang X, Zhu J, Li K, Guo X. Using specific length amplified fragment sequencing to construct the high-density genetic map for Vitis (Vitis vinifera L. x Vitis amurensis Rupr.). Front Plant Sci. 2015;6:393.
Jones DF. Gene action in heterosis. Genetics. 1957;42:93–103.
Gowda M, Longin CFH, Lein V, Reif JC. Relevance of specific versus general combining ability in winter wheat. Crop Sci. 2012;52(6):2494–500.
Tao Z, Xian LN, Kai FJ, Hua FD, Qing H, Qian-h-Hua Y, Li Y, XianQi W, Ying-Jiang C, Jia-Kui Z. Relationship between heterosis and parental genetic distance based on molecular markers for functional genes related to yield traits in rice. Rice Sci. 2010;17(4):288–95.
Hajela RK, Horvath DP, Gilmour SJ, Thomashow MF. Molecular cloning and expression of cor (Cold-Regulated) genes in Arabidopsis thaliana. Plant Physiol. 1990;93(3):1246–52.
Zhang L, Guo X, Zhang Z, Wang A, Zhu J. Cold-regulated gene LeCOR413PM2 confers cold stress tolerance in tomato plants. Gene. 2021;764:145097.
Hwarari D, Guan Y, Ahmad B, Movahedi A, Min T, Hao Z, Lu Y, Chen J, Yang L. ICE-CBF-COR signaling cascade and its regulation in plants responding to cold stress. Int J Mol Sci. 2022;23(3):1549.
Guo W, Ward RW, Thomashow MF. Characterization of a cold-regulated wheat gene related to Arabidopsis cor47. Plant Physiol. 1992;100(2):915–22.
Horvath DP, McLarney BK, Thomashow MF. Regulation of Arabidopsis thaliana L. (Heyn) cor78 in response to low temperature. Plant Physiol. 1993;103(4):1047–53.
Uemura M, Gilmour SJ, Thomashow MF, Steponkus PL. Effects of COR6.6 and COR15am polypeptides encoded by COR (cold-regulated) genes of Arabidopsis thaliana on the freeze-induced fusion and leakage of liposomes. Plant Physiol. 1996;111(1):313–27.
Zhou M, Wu L, Liang J, Shen C, Lin J. Expression analysis and functional characterization of a novel cold-responsive gene CbCOR15a from Capsella bursa-pastoris. Mol Biol Rep. 2012;39(5):5169–79.
Thalhammer A, Bryant G, Sulpice R, Hincha DK. Disordered cold regulated15 proteins protect chloroplast membranes during freezing through binding and folding, but do not stabilize chloroplast enzymes in vivo. Plant Physiol. 2014;166(1):190–201.
Jia Y, Ding Y, Shi Y, Zhang X, Gong Z, Yang S. The cbfs triple mutants reveal the essential functions of CBFs in cold acclimation and allow the definition of CBF regulons in Arabidopsis. New Phytol. 2016;212(2):345–53.
Zhao C, Zhang Z, Xie S, Si T, Li Y, Zhu JK. Mutational evidence for the critical role of CBF transcription factors in cold acclimation in Arabidopsis. Plant Physiol. 2016;171(4):2744–59.
Breton G, Danyluk J, Charron JB, Sarhan F. Expression profiling and bioinformatic analyses of a novel stress-regulated multispanning transmembrane protein family from cereals and Arabidopsis. Plant Physiol. 2003;132(1):64–74.
Okawa K, Nakayama K, Kakizaki T, Yamashita T, Inaba T. Identification and characterization of Cor413im proteins as novel components of the chloroplast inner envelope. Plant Cell Environ. 2008;31(10):1470–83.
Ma X, Wang G, Zhao W, Yang M, Ma N, Kong F, Dong X, Meng Q. SlCOR413IM1: a novel cold-regulation gene from tomato, enhances drought stress tolerance in tobacco. J Plant Physiol. 2017;216:88–99.
Wang J, Zuo K-J, Qin J, Zhang L, Su L, Liu J, Ling H, Zhao J-Y, Cao Y-F, Tang K-X. Isolation and bioinformatics analyses of a COR413-like gene from Gossypium barbadense. Acta Physiol Plant. 2006;29(1):1–9.
Chen Y, Jiang J, Chang Q, Gu C, Song A, Chen S, Dong B, Chen F. Cold acclimation induces freezing tolerance via antioxidative enzymes, proline metabolism and gene expression changes in two chrysanthemum species. Mol Biol Rep. 2014;41(2):815–22.
Ma X, Chen C, Yang M, Dong X, Lv W, Meng Q. Cold-regulated protein (SlCOR413IM1) confers chilling stress tolerance in tomato plants. Plant Physiol Biochem. 2018;124:29–39.
Shi Y, Ding Y, Yang S. Molecular regulation of CBF signaling in cold acclimation. Trends Plant Sci. 2018;23(7):623–37.
Dong C, Zhang Z, Ren J, Qin Y, Huang J, Wang Y, Cai B, Wang B, Tao J. Stress-responsive gene ICE1 from Vitis amurensis increases cold tolerance in tobacco. Plant Physiol Biochem. 2013;71:212–7.
Xu W, Jiao Y, Li R, Zhang N, Xiao D, Ding X, Wang Z. Chinese wild-growing Vitis amurensis ICE1 and ICE2 encode MYC-type bHLH transcription activators that regulate cold tolerance in Arabidopsis. PLoS One. 2014;9(7):e102303.
Considine MJ, Foyer CH. Redox regulation of plant development. Antioxid Redox Signal. 2014;21(9):1305–26.
Cui XY, Gao Y, Guo J, Yu TF, Zheng WJ, Liu YW, Chen J, Xu ZS, Ma YZ. BES/BZR transcription factor TaBZR2 positively regulates drought responses by activation of TaGST1. Plant Physiol. 2019;180(1):605–20.
Abdul Kayum M, Nath UK, Park JI, Biswas MK, Choi EK, Song JY, Kim HT, Nou IS. Genome-wide identification, characterization, and expression profiling of Glutathione S-Transferase (GST) Family in pumpkin reveals likely role in cold-stress tolerance. Genes (Basel). 2018;9(2):84.
Zhang Y, Ming R, Khan M, Wang Y, Dahro B, Xiao W, Li C, Liu JH. ERF9 of Poncirus trifoliata (L.) Raf. undergoes feedback regulation by ethylene and modulates cold tolerance via regulating a glutathione S-transferase U17 gene. Plant Biotechnol J. 2022;20(1):183–200.
Schmidt R, Mieulet D, Hubberten HM, Obata T, Hoefgen R, Fernie AR, Fisahn J, San Segundo B, Guiderdoni E, Schippers JH, et al. Salt-responsive ERF1 regulates reactive oxygen species-dependent signaling during the initial response to salt stress in rice. Plant Cell. 2013;25(6):2115–31.
Lin T, Yang W, Lu W, Wang Y, Qi X. Transcription factors PvERF15 and PvMTF-1 form a cadmium stress transcriptional pathway. Plant Physiol. 2017;173(3):1565–73.
Zhuo C, Liang L, Zhao Y, Guo Z, Lu S. A cold responsive ethylene responsive factor from Medicago falcata confers cold tolerance by up-regulation of polyamine turnover, antioxidant protection, and proline accumulation. Plant Cell Environ. 2018;41(9):2021–32.
Bouzroud S, Gasparini K, Hu G, Barbosa MAM, Rosa BL, Fahr M, Bendaou N, Bouzayen M, Zsogon A, Smouni A, et al. Down regulation and loss of auxin response factor 4 function using CRISPR/Cas9 alters plant growth, stomatal function and improves tomato tolerance to salinity and osmotic stress. Genes (Basel). 2020;11(3):272.
Xu L, Wang D, Liu S, Fang Z, Su S, Guo C, Zhao C, Tang Y. Comprehensive atlas of wheat (Triticum aestivum L.) auxin response factor expression during male reproductive development and abiotic stress. Front Plant Sci. 2020;11:586144.
Su L, Xu M, Zhang J, Wang Y, Lei Y, Li Q. Genome-wide identification of auxin response factor (ARF) family in kiwifruit (Actinidia chinensis) and analysis of their inducible involvements in abiotic stresses. Physiol Mol Biol Plants. 2021;27(6):1261–76.
Wang Y, Jiang H, Mao Z, Liu W, Jiang S, Xu H, Su M, Zhang J, Wang N, Zhang Z, et al. Ethylene increases the cold tolerance of apple via the MdERF1B-MdCIbHLH1 regulatory module. Plant J. 2021;106(2):379–93.
Tang Y, Du G, Xiang J, Hu C, Li X, Wang W, Zhu H, Qiao L, Zhao C, Wang J, et al. Genome-wide identification of auxin response factor (ARF) gene family and the miR160-ARF18-mediated response to salt stress in peanut (Arachis hypogaea L.). Genomics. 2022;114(1):171–84.
Aslam M, Sugita K, Qin Y, Rahman A. Aux/IAA14 regulates microRNA-mediated cold stress response in Arabidopsis roots. Int J Mol Sci. 2020;21(22):8441.
Chen D, Wang W, Wu Y, Xie H, Zhao L, Zeng Q, Zhan Y. Expression and distribution of the auxin response factors in Sorghum bicolor during development and temperature stress. Int J Mol Sci. 2019;20(19):4816.
Wang M, Dai W, Du J, Ming R, Dahro B, Liu JH. ERF109 of trifoliate orange (Poncirus trifoliata (L.) Raf.) contributes to cold tolerance by directly regulating expression of Prx1 involved in antioxidative process. Plant Biotechnol J. 2019;17(7):1316–32.
Khan M, Hu J, Dahro B, Ming R, Zhang Y, Wang Y, Alhag A, Li C, Liu JH. ERF108 from Poncirus trifoliata (L.) Raf. functions in cold tolerance by modulating raffinose synthesis through transcriptional regulation of PtrRafS. Plant J. 2021;108(3):705–24.
Sun X, Zhao T, Gan S, Ren X, Fang L, Karungo SK, Wang Y, Chen L, Li S, Xin H. Ethylene positively regulates cold tolerance in grapevine by modulating the expression of ethylene response factor 057. Sci Rep. 2016;6:24066.
Xie Y, Chen P, Yan Y, Bao C, Li X, Wang L, Shen X, Li H, Liu X, Niu C, et al. An atypical R2R3 MYB transcription factor increases cold hardiness by CBF-dependent and CBF-independent pathways in apple. New Phytol. 2018;218(1):201–18.
An JP, Li R, Qu FJ, You CX, Wang XF, Hao YJ. R2R3-MYB transcription factor MdMYB23 is involved in the cold tolerance and proanthocyanidin accumulation in apple. Plant J. 2018;96(3):562–77.
An JP, Wang XF, Zhang XW, Xu HF, Bi SQ, You CX, Hao YJ. An apple MYB transcription factor regulates cold tolerance and anthocyanin accumulation and undergoes MIEL1-mediated degradation. Plant Biotechnol J. 2020;18(2):337–53.
Xing C, Liu Y, Zhao L, Zhang S, Huang X. A novel MYB transcription factor regulates ascorbic acid synthesis and affects cold tolerance. Plant Cell Environ. 2019;42(3):832–45.
Yang X, Luo Y, Bai H, Li X, Tang S, Liao X, Zhang L, Liu Q. DgMYB2 improves cold resistance in chrysanthemum by directly targeting DgGPX1. Hortic Res.2022;9:uhab028.
Xie XB, Li S, Zhang RF, Zhao J, Chen YC, Zhao Q, Yao YX, You CX, Zhang XS, Hao YJ. The bHLH transcription factor MdbHLH3 promotes anthocyanin accumulation and fruit colouration in response to low temperature in apples. Plant Cell Environ. 2012;35(11):1884–97.
Xu W, Zhang N, Jiao Y, Li R, Xiao D, Wang Z. The grapevine basic helix-loop-helix (bHLH) transcription factor positively modulates CBF-pathway and confers tolerance to cold-stress in Arabidopsis. Mol Biol Rep. 2014;41(8):5329–42.
Zhang Y, Yu H, Yang X, Li Q, Ling J, Wang H, Gu X, Huang S, Jiang W. CsWRKY46, a WRKY transcription factor from cucumber, confers cold resistance in transgenic-plant by regulating a set of cold-stress responsive genes in an ABA-dependent manner. Plant Physiol Biochem. 2016;108:478–87.
Yang Y, Liu J, Zhou X, Liu S, Zhuang Y. Identification of WRKY gene family and characterization of cold stress-responsive WRKY genes in eggplant. PeerJ. 2020;8:e8777.
Mills LJ, Ferguson JC, Keller M. Cold-hardiness evaluation of grapevine buds and cane tissues. Am J Enol Vitcult. 2006;57(2):194-200.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28(11):1947–51.
Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49(D1):D545–51.
Broman KW, Wu H, Sen S, Churchill GA. R/qtl: QTL mapping in experimental crosses. Bioinformatics (Oxford, England). 2003;19(7):889–90.
Fujimori N, Enoki S, Suzuki A, Naznin HA, Shimizu M, Suzuki S. Grape apoplasmic β-1,3-glucanase confers fungal disease resistance in Arabidopsis. Sci Hortic. 2016;200:105–10.
The English language editing in this article was supported by editage.
The collection of plant material comply with relevant institutional, national, and international guidelines and legislation. All of the material is owned by the authors and/or no permissions are required.
The research was supported by the National Natural Science Foundation of China (Grant No.31972368), the China Agriculture Research System (Grant No. CARS-29-yc-6) and the Liaoning key R&D Program (Grant No. 2020JH2/ 10200032). Funding bodies were not involved in the design of the study and collection, analysis, interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Unigenes obtained based on RNA-seq.
Annotation analysis for unannotated unigenes.
FPKM value of each unigene.
Differentially expressed genes between SY60d vs. RG60d.
The expression levels of 28 selected genes.
Gene expression in hormone signal transduction pathway.
Candidate genes related to cold resistance in‘DNA binding transcription factor activity’ catalogue.
Candidate cold resistance genes selected in QTL mapping.
The primers used in this study.
About this article
Cite this article
Ma, X., Zhao, F., Su, K. et al. Discovery of cold-resistance genes in Vitis amurensis using bud-based quantitative trait locus mapping and RNA-seq. BMC Genomics 23, 551 (2022). https://doi.org/10.1186/s12864-022-08788-y
- Vitus amurensis
- Grape bud
- Cold resistance breeding
- Transcriptome analysis
- QTL mapping