Unraveling epigenetic and genetic variations in response to cold stress in two lotus ecotypes

Genetic variations accumulate over long evolutionary timescales, whereas epigenetic modifications can arise rapidly and be inherited across generations. However, the interplay between genetic and epigenetic variations in shaping ecotype-specific phenotypic plasticity remains elusive. Focusing on two lotus ecotypes that evolved under distinct winter temperature regions and display divergent annual growth cycles under cold stress, we generated DNA methylation landscapes across three sequence contexts (CG, CHG, and CHH, where H = A, T, or C) and identified single-cytosine methylation polymorphisms (SMPs) and single-nucleotide polymorphisms (SNPs). Interestingly, only CG methylation patterns mirror population-level genetic variations. Using epigenetic genome-wide association analysis, we identified differentially methylated CG sites that are either cis- or trans-regulated by SNP. Notably, we constructed a multifactorial regulatory network centered on the NnMKK4-NnCYCD5 module, linking cold response with cell cycle regulation. Temperature stress experiments conducted on lotus ecotypes and transgenic Arabidopsis (OE-NnMKK4 and OE-NnCYCD5) confirmed that NnMKK4 acts as a cold receptor and that NnCYCD5 promotes cell cycle progression and growth under cold conditions. Collectively, our findings provide novel insights into the co-evolutionary dynamics of epigenetic and genetic variations that are associated with different growth cycles of lotus ecotypes in response to cold stress.