Background: Cytosine methylation, the main type of DNA methylation, regulates gene expression in plant response to environmental stress. The winter rapeseed has high economic and ecological value in China's Northwest, but the DNA methylation pattern of winter rapeseed during freezing stress remains unclear. Result: This study integrated the methylome and transcriptome to explore the genome-scale DNA methylation pattern and its regulated pathway of winter rapeseed, using freezing-sensitive (NF) and freezing-resistant (NS) cultivars.The average methylation level decreased under freezing stress, and the decline in NF was stronger than NS after freezing stress. The CG methylation level was the highest among the three contexts of CG, CHG, and CHH. At the same time, the CHH proportion was high, and the methylation levels were highest 2 kb up/downstream, followed by the intron region. The C sub-genomes methylation level was higher than the A sub-genomes. The methylation levels of chloroplast and mitochondrial DNA were much lower than the B. napus nuclear DNA, the SINE methylation level was highest among four types of transposable elements (TEs), and the preferred sequence of DNA methylation did not change after freezing stress. A total of 1732 differentially expressed genes associated with differentially methylated genes (DMEGs) were identified in two cultivars under 12 h and 24 h in three contexts by combining whole-genome bisulfite sequencing( and RNA-Seq data. Function enrichment analysis showed that most DMEGs participated in linoleic acid metabolism, alpha-linolenic acid metabolism, carbon fixation in photosynthetic organisms, flavonoid biosynthesis, and plant hormone signal transduction pathways. Meanwhile, some DMEGs encode core transcription factors in plant response to stress. Conclusion: Based on the findings of DNA methylation, the freezing tolerance of winter rapeseed is achieved by enhanced signal transduction, lower lipid peroxidation, stronger cell stability, increased osmolytes, and greater reactive oxygen species (ROS) scavenging. These results provide novel insights into better knowledge of the methylation regulation of tolerance mechanism in winter rapeseed under freezing stress.

Two winter rapeseed cultivars, “NS” (cold tolerant) and “NF” (cold sensitive), were used to reveal the morphological, physiological, and proteomic characteristics in leaves of plants after treatment at -4°C for 12 h(T1) and 24 h(T2), and at room temperature(T0), to understand the molecular mechanisms of cold tolerance. Antioxidant activity and osmotic adjustment ability were higher, and plasma membrane injury was less obvious, in NS than in NF under cold stress. We detected different abundant proteins (DAPs) related to cold tolerance in winter rapeseed through data-independent acquisition (DIA). Compared with NF, A total of 1,235 and 1,543 DAPs were identified in the NSs under T1 and T2, respectively. Compared with NF, 911 proteins were more abundant in NS only after cold treatment. Some of these proteins were related to ROS scavenging through four metabolic pathways: lysine degradation; phenylalanine, tyrosine, and tryptophan; flavonoid biosynthesis; and ubiquinone and other terpenoid-quinone biosynthesis. Analysis of these proteins in the four candidate pathways revealed that they were rapidly accumulated to quickly enhance ROS scavenging and improve the cold tolerance of NS. These proteins were noticeably more abundant during the early stage of cold stress, which was critical for avoiding ROS damage.

Winter rapeseed is susceptible to low temperature during winter in Northwest China, which could lead to a severe reduction of crop production. The freezing temperature could stress the whole plant, especially the leaf, and ultimately harm the survival rate of winter rapeseed. However, the molecular mechanism underlying freezing tolerance is still unclear in winter rapeseed. In this study, a comprehensive investigation of winter rapeseed freezing tolerance was conducted at the levels of transcript, protein, and physiology and biochemistry, using a pair of freezing-sensitive and freezing-resistant cultivars NQF24 and 17NTS57. There were 4,319 unique differentially expressed genes (DEGs) and 137 unique differentially abundant proteins (DAPs) between two cultivars identified in leaf under freezing stress. Function enrichment analysis showed that most of the enriched DEGs and DAPs were involved in plant hormone signal transduction, alpha-linolenic/linoleic acid metabolism, peroxisome, glutathione metabolism, fatty acid degradation, and secondary metabolite biosynthesis pathways. Based on our findings, it was speculated that freezing tolerance formation is caused by increased signal transduction, enhanced biosynthesis of protein, secondary metabolites, and plant hormones, elevated energy supply, greater reactive oxygen species scavenging, and lower lipid peroxidation as well as stronger cell stability in leaf under freezing stress. These results provide a comprehensive profile of leaf response under freezing stress, which have potential to be used as selection indicators of breeding programs to improve freezing tolerance in rapeseed.

Flowering at the proper time is an important part of acclimation to the ambient environment and season and maximizes the plant yield. To reveal the genetic architecture and molecular regulation of flowering time in oilseed rape (Brassica napus), we performed an RNA-seq analysis of the two parents after vernalization at low temperature and combined this with quantitative trait loci (QTL) mapping in an F₂ population. A genetic linkage map that included 1,017 markers merged into 268 bins and covered 793.53 cM was constructed. Two QTLs associated with flowering time were detected in the F₂ population. qFTA06 was the major QTL in the 7.06 Mb interval on chromosome A06 and accounted for 19.3% of the phenotypic variation. qFTC08 was located on chromosome C06 and accounted for 8.6% of the phenotypic variation. RNA-seq analysis revealed 4,626 differentially expressed genes (DEGs) between two parents during vernalization. Integration between QTL mapping and RNA-seq analysis revealed six candidate genes involved in the regulation of flowering time through the circadian clock/photoperiod, auxin and ABA hormone signal, and cold signal transduction and vernalization pathways. These results provide insights into the molecular genetic architecture of flowering time in B. napus.