Abstract
Asparagus (Asparagus officinalis L.) is an important vegetable crop in southern Ontario, Canada, where winter air and soil temperatures below 0°C are common. Consequently, cultivars growing in this area must possess winterhardiness and freezing tolerance for survival. Asparagus acquires freezing tolerance in the fall through cold acclimation and loses freezing tolerance in the spring through deacclimation. To understand the molecular bases of these processes, transcriptomic analysis (RNA-Seq) was conducted on two cultivars, one adapted, 'Guelph Millennium' (GM), and one unadapted, 'UC157' (UC), to the winter conditions of southern Ontario. RNA extracted from bud and rhizome tissues, sampled on three dates during early spring and late fall, was subjected to sequencing. In the fall, the numbers of differentially expressed (DE) genes at the second and third harvests increased, relative to the first harvest, in dormant buds and rhizomes as freezing tolerance of cultivars increased, and the majority of DE genes were downregulated. In spring, freezing tolerance decreased as plants deacclimated and most genes DE at second and third harvests were upregulated in both cultivars. GM had lower LT50 (lethal temperature at which 50% of plants die) values and hence higher freezing tolerance than UC on specific sampling dates during both spring and fall, and expression patterns of specific genes were correlated with LT50 differences. Functional analysis revealed that these genes were involved in carbohydrate metabolic process, plant hormone signal transduction (auxin and gibberellin), proline metabolism, biosynthesis of secondary metabolites, circadian rhythm, and late embryogenesis abundant proteins and could be associated with cold acclimation and deacclimation processes. These findings will help researchers understand the molecular mechanisms of freezing tolerance in asparagus, leading to breeding and genetic strategies to improve the trait.