Skip to main content

DNA methylation analysis of floral parts revealed dynamic changes during the development of homostylous Fagopyrum tataricum and heterostylous F. esculentum flowers

Abstract

Background

Proper flower development is essential for plant reproduction, a crucial aspect of the plant life cycle. This process involves precisely coordinating transcription factors, enzymes, and epigenetic modifications. DNA methylation, a ubiquitous and heritable epigenetic mechanism, is pivotal in regulating gene expression and shaping chromatin structure. Fagopyrum esculentum demonstrates anti-hypertensive, anti-diabetic, anti-inflammatory, cardio-protective, hepato-protective, and neuroprotective properties. However, the heteromorphic heterostyly observed in F. esculentum poses a significant challenge in breeding efforts. F. tataricum has better resistance to high altitudes and harsh weather conditions such as drought, frost, UV-B radiation damage, and pests. Moreover, F. tataricum contains significantly higher levels of rutin and other phenolics, more flavonoids, and a balanced amino acid profile compared to common buckwheat, being recognised as functional food, rendering it an excellent candidate for functional food applications.

Results

This study aimed to compare the DNA methylation profiles between the Pin and Thrum flower components of F. esculentum, with those of self-fertile species of F. tataricum, to understand the potential role of this epigenetic mechanism in Fagopyrum floral development. Notably, F. tataricum flowers are smaller than those of F. esculentum (Pin and Thrum morphs). The decline in DNA methylation levels in the developed open flower components, such as petals, stigmas and ovules, was consistent across both species, except for the ovule in the Thrum morph. Conversely, Pin and Tartary ovules exhibited a minor decrease in DNA methylation levels. The highest DNA methylation level was observed in Pin stigma from closed flowers, and the most significant decrease was in Pin stigma from open flowers. In opposition, the nectaries of open flowers exhibited higher levels of DNA methylation than those of closed flowers. The decrease in DNA methylation might correspond with the downregulation of genes encoding methyltransferases.

Conclusions

Reduced overall DNA methylation and the expression of genes associated with these epigenetic markers in fully opened flowers of both species may indicate that demethylation is necessary to activate the expression of genes involved in floral development.

Peer Review reports

Background

Buckwheat constitutes a small genus of 22 members [1]. Two of the most widely cultivated species of this genus are Fagopyrum esculentum Moench (common buckwheat) and Fagopyrum tataricum (Tartary buckwheat) [2, 3]. Common buckwheat is commonly cultivated in Russia, China, France, Poland, and the American regions [4, 5]. In terms of overall yield, buckwheat is a minor crop. However, its popularity among consumers steadily increases due to its remarkable properties. Common buckwheat is rich in beneficial phenolic compounds like rutin, quercetin and C-glycosyl flavones (orientin, isoorientin and vitexin). Additionally, it is gluten-free and provides amino acids, dietary fibre, resistant starch, and vitamins [6,7,8]. It was demonstrated that F. esculentum exhibits anti-hypertension, anti-diabetic, anti-inflammatory, cardio-protective, hepato-protective, and neuroprotective properties [7, 8]. Tartary buckwheat exhibits higher than common buckwheat resistance to high altitudes, harsh weather conditions such as drought, frost, UV-B radiation damage, and pests [9]. Additionally, compared to common buckwheat, F. tataricum contains more rutin and other phenolics and flavonoids, with well-balanced amino acids, making it a functional pseudocereal [10]. One significant advantage of Tartary buckwheat lies in the differences in flower structure and pollination methods, which directly impact yield. F. esculentum is an obligatory cross-pollinating, heterostylous species, while F. tataricum is self-pollinating, homostylous species eliminating reliance on external pollinators [11, 12]. F. esculentum has two floral morphs: Pin and Thrum, which differ in stamen and style length ratio, amount and size of the produced pollen grain, exine sculpturing, and nectar production [13, 14]. The distyly of F. esculentum is controlled by a cluster of genes, i.e. supergene S [13, 15]. Recently, the S-LOCUS EARLY FLOWERING 3 gene (S-ELF3), present within the S supergene locus, was shown to control style length and style incompatibility. The inactivation of S-ELF3 has led to the breakage of heterostyly and successful pollination [15]. While in the incompatible intra-morph pollination, the pollen tube growth within the style is inhibited [13, 16, 17].

Proper flower development ensures reproduction, one of the most vital parts of the plant’s life cycle. This process relies on fine-tuned machinery of transcription factors, enzymes and epigenetic modifications to function correctly. DNA methylation, a universal and heritable epigenetic mechanism, regulates gene expression and the chromatin structure [16, 18]. The extent of DNA methylation significantly shapes the plant genome’s structure and function. Various molecular mechanisms are impacted by DNA methylation, highlighting its pivotal role in regulating cellular processes within plants [16, 19,20,21,22]. DNA methylation takes place at CG, CHG and CHH sites (H = A, T or C); methylation at these three sites is achieved by DNA METHYLTRANSFERASE 1 (MET1), the plant-specific CHROMOMETHYLASE 3 (CMT3), and DOMAINS REARRANGED METHYLTRANSFERASEs (DRMs), respectively [23, 24]. Each type is crucial for development and the ability to respond to environmental stresses [16, 25]. An opposite phenomenon, demethylation, is dependent on four bifunctional 5-methylcytosine glycosylases: REPRESSOR OF SILENCING 1 (ROS1), DEMETER (DME), DME-LIKE 2 (DML2), and DML3, which are involved in the removal of methylated bases and cleavage of the DNA backbone at abasic sites [26]. ROS1 counteracts the DNA methylation pathway to prevent plant gene silencing [27]. Research demonstrated that DNA methylation plays an essential role in flower development, with demethylation dominating during early flower development [28]. It was also reported that decreased levels of DNA methylation induced early flowering and that demethylation activities are correlated with changes in the expression levels of DNA methylation genes [28,29,30]. On the other hand, experiments on Azalea japonica and Arabidopsis demonstrated an increase in global DNA methylation during the transition from the vegetative to the flowering stage [31, 32]. DNA methylation was determined to influence various aspects of flower development and morphology [33]. It affects floret closing in barley, and the development of bisexual flowers in Populus cathayana and Fraxinus mandshurica [34,35,36]. DNA methylation is also involved in forming double flowers and dichromatic petals [37, 38]. Moreover, it plays a crucial role in the expression of genes related to anthocyanin pigmentation in flower tissue [9, 39, 40]. Existing evidence suggests that methylation alterations accompany the flowering process in plants. Methylation regulates key flowering-related genes such as SOC1, AP1, and SPL and specific transcription factor genes, including WUS homeobox-containing (WOX) genes in apple [41].

While numerous processes related to DNA methylation in plant vegetative tissues have been extensively studied, there remains a scarcity of information regarding DNA methylation during flower development, particularly in flowers exhibiting distinct morphologies. Thus, this study aimed to compare the DNA methylation status of Pin and Thrum flower components of F. esculentum that participate in pollination and embryo production with self-sufficient species of F. tataricum.

Methods

Plant material

F. tataricum seeds, sample k-17, were obtained from the N. I. Vavilov Institute of Plant Genetic Resources, Saint Petersburg, Russia. K-17 sample is a widely cultivated landrace of F. tataricum. Seeds are available upon request from the authors. F. esculentum seeds of the Panda cultivar are commercially available and were supplied from the Malopolska Plant Breeding, Poland. Both plant species were grown in pots with soil mixed with vermiculite (3:1, w/v) in a greenhouse at 20 ± 1 °C, under a 16/8 h light/dark photoperiod, provided by lamps emitting white light at the intensity of 90 µmol m − 2 s − 1. After approximately three to four weeks, the flowers started to appear and were gradually collected, photographed with the use of Keyence VHX-970 F digital microscope (Japan) equipped with an ultra-small high-performance zoom lens VH-Z20R/Z20T and wide-area illumination adapter OP-87,298 and designated for further procedures.

For analysis, the closed flowers from both species were fixed in the intact state. Open flowers were dissected into parts that involved excised petals, stigmas and remains that contained ovaries and nectaries.

Histological and immunostaining procedures

Closed and open flower components from F. tataricum, as well as closed and open Pin and Thrum flower components from F. esculentum were fixed in 4% paraformaldehyde (Sigma-Aldrich, USA) in 1x phosphate-buffered saline (PBS), pH 7.3 and placed in the vacuum desiccator for three hours (with 30 min intervals), after which the incubation at 4 °C overnight has followed. After approximately 24 h, the fixative was replaced with 1 × PBS (twice for 15 min) and followed by dehydration in a graded ethanol series diluted in 1xPBS solution in each concentration (10%, 30%, 50%, 70% and 90%) and 99,8% twice for 30 min each. Subsequently, the embedding procedure was performed according to Wolny et al., 2014 [42]. 5 μm thick sections were prepared using a HYRAX M40 rotary microtome (Zeiss, Oberkochen, Germany) and placed on polysine-coated microscope slides (Epredia, Netherlands). Next, the de-embedding procedure comprised of placing the slides in 99,8% ethanol three times for 10 min, followed by the rehydration in ethanol/ 1xPBS solutions: 90%, 50%v/v and, finally in 1xPBS for 10 min each. Such prepared slides were used for both, the immunostaining and histological analysis. For the histological analysis, slides were stained with 0.05% aqueous solution of Toluidine blue O (TBO, Sigma-Aldrich, USA) for approximately 10 min, rinsed and mounted with 50% glycerol/ distilled water solution (v/v). Observations and photographs of histological sections were performed with an Olympus BX43F microscope equipped with an Olympus XC50 digital camera.

The immunostaining method was previously established by Braszewska-Zalewska et al., 2013 [43]. After the wax de-embedding described above, samples underwent the 2 N HCl (Sigma-Aldrich, USA) digestion for 45 min to denature DNA. Following, the slides were incubated with 5% bovine serum albumin (BSA, Sigma-Aldrich, USA) in 1xPBS for 1 h in the humid chamber at room temperature. Next, the primary antibody diluted in 1% BSA in 1xPBS (1:100) was applied, and slides were incubated at 4o C overnight (Table 1). After the incubation, the slides were washed three times in 1xPBS. Subsequently, a secondary antibody diluted in 1% BSA in 1xPBS (1:100) was applied, and samples were incubated at 37 o C in the humid chamber in the dark for 1 h. After the incubation, the slides were again washed three times in 1xPBS and nuclei were counterstained with 4’,6-diamidyno-2-fenyloindol (DAPI, 2.5 g/ml in Vectashield).

Graphical depiction of the experimental design is available in Additional Files, File 1.

Table 1 List of antibodies used in the immunostaining

Fluorescence intensity measurements and statistical analysis

Images were captured with the Olympus FV1000 confocal system (Olympus, Poland) equipped with an Olympus IX81 inverted microscope. Fluorescence of Alexa488 (excitation 488 nm, emission 500–600 nm) was acquired from a 60x Plan Apo oil-immersion objective lens (NA 1.35), a 50 mW 405 nm diode laser and a 100 mW multi-line argon ion laser (Melles Griot BV, the Netherlands). The confocal laser scanning microscope offers a significant advantage in obtaining high-resolution images of optical sections through fluorescently labelled nuclei (z-stacks) [44, 45]. The pixel-by-pixel image generation produces exceptional quality images [45], enabling visualization of fluorescence distribution corresponding to DNA methylation within the three-dimensional chromatin architecture of a nucleus. The nucleus area was used to standardise the fluorescence intensity, preventing any artificial positive correlation between intensity and nucleus size. Consequently, the fluorescence intensity-to-nucleus area ratio indicated the level of DNA methylation [46]. An axial series of two-dimensional fluorescence images of the optical sections through the nuclei (z-stacks) was collected with the use of two separate photomultipliers (R6357, Hamamatsu Photonics, Hamamatsu, Japan) set to work in the integration mode at 4 µs pixel dwell time and 12-bit signal digitisation (4096 intensity levels). The fluorescence intensity levels of Alexa488 fluorochrome were subsequently measured in the ImageJ version 1.53s software (Wayne Rasband, National Institutes of Health, USA). ImageJ is a Java-based image processing program developed at the National Institutes of Health and the Laboratory for Optical and Computational Instrumentation (LOCI, University of Wisconsin) [47,48,49,50]. It performs various complex tasks such as editing, processing, and analysing 8-bit colour and grayscale images. ImageJ calculates area and pixel value statistics of selections and intensity-thresholded objects defined by the user and offers standard image processing [51,52,53]. Images were converted to eight bits and segmented with the threshold value parameter. Alexa488 fluorescence intensity was calculated as the mean values from the Integrated Density parameter per one nucleus, which depicted the sum of all pixels within the region of interest. Results are presented in relative units. Raw data obtained from the images were further analysed using statistical analysis software. This study utilised R Studio, an integrated development environment for R, a statistical programming language [54]. R is widely recognized as a statistical language, dominating other programming languages in developing statistical tools. It offers various packages suitable for scientific data analysis, including agricolae (Statistical Procedures for Agricultural Research), which can be used to determine significant differences between means [55], as well as other packages dedicated to data visualisation [56]. R facilitates the reproducible process of summarizing data after statistical analysis and visualizing it in graph form [57]. 500 nuclei were analysed from petals; 1000 nuclei from the nectary; 1500 nuclei from the ovary, and 1000 nuclei from the stigma in open and closed F. tataricum flowers and F. esculentum Pin open and closed and Thrum open and closed flower type. Numerical data is included in Additional File 2. Fluorescence data analysis and plotting was performed in R, a software environment for statistical computing and graphics in R Studio 2022.12.0 Build 353 (Script included in Additional File 3), an integrated development environment for R [58, 59]. One-way ANOVA test with the R Stats package [59] was followed by Tukey’s HSD test at the significance level p ≤ 0.05. Standard errors were likewise calculated with the stats package. The package agricolae (Statistical Procedures for Agricultural Research) calculated significant differences between means [55]. Subsequently, data was plotted with packages dedicated to visualising data: ggplot2 [60] and ggpubr [56]. Letters on the graphs indicate statistically significant differences between samples. Full immunostaining procedure and the fluorescence intensity analysis is described for Fagopyrum species [61, 62].

RNA isolation and real-time qPCR

Total RNA was isolated from the closed and open flowers and Thrum and Pin types of F. esculentum and F. tataricum. Total RNA was isolated using a FastPure Plant Total RNA Isolation Kit (Polysaccharides and polyphenolics-rich) (Vazyme Biotech, Red Maple Hi-tech Industry Park, Nanjing, PRC). RNA concentrations were measured using a Nano-Drop ND-1000 (NanoDrop Technologies, Wilmington, DE, USA). The DNA was removed from the RNA samples by digesting them with an RNase-free DNase Set (Qiagen, Hilden, Germany). The oligo-dT primers and a Maxima H Minus First Strand cDNA Synthesis Kit (Thermo Fisher Scientific, Waltham, MA, USA) were used to produce the cDNA. The obtained cDNA was diluted four-fold with water and used at a volume of 2 µl in a qPCR reaction. Analyses were performed in a 10 µl volume using a LightCycler® 480 SYBR Green I Master (Roche, Basel, Switzerland). The primers were designed based on Fagopyrum esculentum “Pintian4” and Fagopyrum tataricum “Pinku1” references genomes with Primer3Plus (Additional File 4) [63]. The control genes (SAND, ACTIN) had a constant expression level in all tissue samples. Analyses were performed using a LightCycler 480 (Roche, Basel, Switzerland) under the following reaction conditions: initial denaturation of 5 min at 95 °C, followed by 10 s at 95 °C, 20 s at a temperature specific for the primers, 10 s at 72 °C, repeated in 40 cycles. Denaturation for the melt curve analysis was conducted for 5 s at 95 °C, followed by 1 min at 65 °C and heating to 98 °C (0.1 °C/s for the fluorescence measurement). The Ct values were calculated using LinRegPCR software (version 11, Academic Medical Centre, Amsterdam, The Netherlands). The plant tissues for the Real-Time qPCR analysis were produced in three biological repetitions, and two technical replicates of each repetition were analysed. The relative expression level was calculated using 2–∆∆CT, where ∆∆CT represents ∆CTreference condition − ∆CTcompared condition.

Statistical analysis

The Student t-test and one-way ANOVA (p < 0.05) followed by Tukey’s honestly-significant-difference test (Tukey HSD-test) (p < 0.05) were used to calculate any significant differences between the experimental combinations. The graphs show the average values with the standard error (SE) in Fig. 4 and the standard deviation (SD) in Fig. 5.

Graphics

Photographs were cropped, adjusted (brightness, contrast) and arranged into figures in the Corel Draw 2020 program. A graphical depiction of the experimental design (Additional File 1) was prepared with BioRender (available online at: https://app.biorender.com/). Publication licences are included as Additional File 5.

Results

Anatomy and histology of the flowers

F. tataricum flowers, closed and open, were smaller than those of F. esculentum (Pin, Thrum) (Fig. 1). Closed flowers were coiled in petals of a greenish colour (Fig. 1a–c). In open flowers, the differences between two morphs were disclosed: Pin flowers featured a long style and shorter filaments (Fig. 1d), while Thrum had a short style and longer filaments (Fig. 1e). Open flowers of F. tataricum had green petals, with style and filaments of similar length (Fig. 1f).

Fig. 1
figure 1

Morphology of analysed flowers. F. esculentum (a, d) Pin, (b, e) Thrum, and (c, f) F. tataricum. (ac) closed flowers were tightly coiled in petals. (df) In open flowers, the anatomy is easily visible. Arrows – stigmas, double arrows – nectaries, open arrows – anthers, o – ovaries, p – petals. Scale bars: (a–f) = 300 μm

In anthers of closed flowers, the developmental stage of pollen was advanced, with microspores enwrapped in egzine sheath, tapetum cells were already undergoing programmed cell death (Additional File 6a-c). The anthers of open flowers were fully developed and open, with mature pollen and no tapetum cell remains (Additional File 6d-f). Thus, the anthers were not considered in DNA methylation analysis due to the lack of reference (closed vs. open).

Petals of closed flowers were rich in polyphenols, which were detected in the epidermis and mesophyll cells (Fig. 2a–c). Although phenolic compounds were abundantly present in papillae and subpapillae cells of the stigma, the nuclei were visible (Fig. 2d–f). Within ovaries, the ovules comprised a double-layered integument with tightly adhered nucellus parenchyma cells and embryo sac (Fig. 2g–l). The nectaries were composed of the epidermis, nectary parenchyma, and secretory trichomes in which polyphenolic compounds, were mainly detected (Fig. 2m–o). Similarly to closed flowers, petals of open flowers contained large amounts of phenolic compounds, and the same for stigma papillae and subpapillae cells, where nuclei were “masked” with abundant polyphenol presence (Fig. 3a-f). Part of the nucellus parenchyma cells degenerated, leaving space within the ovule (Fig. 3g–i). The histology of nectaries remained unchanged compared to closed flowers (Fig. 3j–l).

Fig. 2
figure 2

Histology of closed flower components (F. esculentum Pin, Thrum; F. tataricum). Schematic diagrams represent analysed flower components: I petals, II stigma, III ovary with ovule, IV nectary, marked in pink on the diagram. (a - c) petals, abundant polyphenols occurrence in epidermis (arrowheads) and mesophyll cells (arrows). (d - f) stigmas, papillae and subpapillae cells rich in polyphenols (arrowheads) but with visible nuclei (arrows). (g - l) ovaries and ovules, polyphenols present in basal part of ovule and in outer integument cells (arrowheads), arrows – nucellus parenchyma cells, embryo sacs (asterisks); double arrow stands for the extent of cells that were taken into account in the measurement of Alexa488 fluorescence intensity. (m - o) nectaries, epidermis rich in polyphenols (arrowheads), secretory trichomes (arrows), np nectary parenchyma. Scale bars: a, d, e, j, k = 10 μm; b, f, g, h, l, m – o = 20 μm; c, i = 50 μm

Fig. 3
figure 3

Histology of open flower components (F. esculentum Pin, Thrum; F. tataricum). Schematic diagrams represent analysed flower components: I petals, II stigma, III ovary with ovule, IV nectary, marked in pink on the diagram. (ac) petals, abundant occurrence of polyphenols in the epidermis (arrowheads) and mesophyll cells (arrows). (df) stigmas, papillae and subpapillae cells containing polyphenols (arrowheads) and hardly visible nuclei (arrows). (gi) ovaries and ovules, high polyphenolic content in the basal part of the ovule and in outer integument cells (arrowheads), arrows – nucellus parenchyma cells, embryo sacs (asterisks); double arrow stands for the extent of cells that were taken into account in the measurement of Alexa488 fluorescence intensity. (jl) nectaries, epidermis rich in polyphenols (arrowheads), secretory trichomes (arrows), np nectary parenchyma. Scale bars: b = 10 μm; a, d – f, i, j = 20 μm; c, k, l = 50 μm; g, h = 100 μm

DNA methylation analysis

The highest DNA methylation level in petals from closed flowers was in the Thrum morph; the Pin morph exhibited slightly lower values, whereas the F. tataricum open flowers presented the lowest level of this modification. In open flowers, the methylation values decreased in each of the analysed floral parts, but the decrease was the most prominent in Thrum petals (Fig. 4a). Summarising, the F. esculentum Thrum morph petals were characterised by the highest starting DNA methylation level and the most eminent decrease; F. tataricum, exhibited the lowest DNA methylation value in closed flowers and the lowest decrease during petal development.

Fig. 4
figure 4

Fluorescence intensity measurements of DNA methylation levels. (a) petal; n = 500, (b) stigma, n = 1000, (c) ovule, n = 1500, and (d) nectary, n = 1000 parts of the closed and the open flowers of F. esculentum and F. tataricum. Bars represent standard error; results presented in relative units, letters indicate statistically significant differences between the groups

Stigma is a pistil component crucial for the first pollination stages. In this floral part, we observed the highest values of the DNA methylation levels in closed flowers among analysed specimens and each flower part and the most considerable reduction in the open flower stage. The highest output DNA methylation level in closed flowers, as well as the decrease (approximately five-fold drop), was noticed in Pin morph (Fig. 4b). Thrum morph was characterised by a lower DNA methylation level concurrently with the lower decrease of the methylation values in comparison to Pin morph (approximately two-fold drop; Fig. 4b). F. tataricum stigmas exhibited the lowest starting DNA methylation level among each variant at the closed flower stage, and around two-fold drop in an open flower, which is similar to value decrease observed in Thrum morph (Fig. 4b).

Among the analysed ovules, the DNA methylation level was highest in closed flowers of F. tataricum, then in the Pin morph and the lowest in the Thrum morph (Fig. 4c). During flower development, a slight reduction of DNA methylation level was observed for both Pin and Tartary (Fig. 4c) on the contrary to Thrum morph ovules, in which the values of methylation level increased in open flowers (Fig. 4c).

The highest DNA methylation level in nectaries from closed flowers was observed in Tartary buckwheat, lower in Pin, and the lowest in the Thrum morph (Fig. 4d). In contrast to other flower parts (except for the Thrum morph ovule), the nectaries from open flowers had higher values of DNA methylation (compared to values of closed flowers) in all three specimens (Fig. 4d). However, the increase in DNA methylation level was minor, not as prominent as in Pin or Thrum morph. In the opened flowers stage, the F. esculentum Pin morph exhibited the highest DNA methylation level, following Tartary, and the lowest values were noted in the Thrum morph.

In conclusion, F. esculentum Pin and Thrum morphs differ significantly in DNA methylation levels in the analysed parts of the flower. In the ovule, stigma, and petals of open flowers, the Thrum morph exhibited higher values than the Pin. Results from Tartary buckwheat show some similarities to either Pin or Thrum, depending on the analysed part and stadium of the flower development.

The graphs show global DNA methylation analysis, but within parts such as ovules and nectaries, the DNA methylation level was not uniform in all nuclei. For example, in ovules, the nuclei from the embryo sac showed very low or no DNA methylation in contrast to other ovule nuclei (Additional File 7).

Gene expression analysis

We analysed the expression level of genes related to DNA methylation (MET1, MET2, CMET3) and demethylation (DME1, DME3, ROS1) in closed and open Pin and Thrum flowers of F. esculentum and flowers of F. tataricum. The results showed the down-regulation of all genes in Thrum flowers (Fig. 5a). The most intense decrease in expression was observed for MET1, which was over 45 times lower in the open than in closed Thrum flowers. In contrast to Thrum flowers, the expression of MET1 in open Pin flowers was the same as in closed ones (Fig. 5b). The decreased transcript level of the rest of the analysed genes was observed in open Pin flowers. The highest reduction was characteristic for the DME3 gene (14 times lower in open than in closed Pin flowers). The decreased transcript level of almost all analysed genes was observed in open flowers of F. tataricum (Fig. 5c). In general, during the development of flowers, the expression of genes encoding DNA methylases and demethylases was decreased.

Fig. 5
figure 5

Expression level of MET1, MET2, CMET3, DME1, DME3, and ROS1 in closed and open flowers, (a, d) Thrum and (b, d) Pin of F. esculentum and (c, d) F. tataricum. The expression level of genes in open flowers was calibrated to expression in closed flowers of the same type and species (a, b, c). * - values significantly different from closed flowers of the same type and species (a, b,c) (p < 0.05; n = 3; means ± SD are given). The expression level of genes in all types and species was calibrated to expression in closed flowers of the F. tataricum (d). Different letters indicate a significant difference between flower type and species according to Tukey’s HSD test (p < 0.05; n = 3; means ± SD are given)

In addition, we analysed the expression of the genes between flower types and species (Fig. 5d). The expression of all genes is higher (from 2.5 to over 9 times higher) in closed flowers of F. esculentum than in F. tataricum. Moreover, we observed the difference in expression of genes between Pin and Thrum flowers, both closed and open. MET2, DME1, and ROS1 expression were lower in Pin than in Thrum flowers. The other three analysed genes (MET1, CMET3, DME3) are characterised by higher transcript levels in Thrum than in Pin flowers.

Discussion

The objective of this study was to compare the DNA methylation status between the Pin and Thrum flower components of F. esculentum, with that of self-pollinating species like F. tataricum in order to shed light on how this epigenetic mark may have a role in the floral development of the Fagopyrum species. F. tataricum flowers are smaller than those of F. esculentum (Pin, Thrum). In F. esculentum open flowers, Pin flowers feature a long style and shorter anther stamens, while Thrum has a short style and longer anther stamens. Open flowers of F. tataricum were characterised by green petals, style, and anthers of similar length. The heteromorphic heterostyly of F. esculentum is a significant limitation in breeding [15, 64]. Around a few thousand phenolic structures in the plant kingdom were reported – ranging from single, aromatic-ringed compounds to complex structures [65]. In F. esculentum, a total of 60 different phenolic substances were identified, with the highest number found in the flower part of the plant [66]. In tea, another species rich in this type of secondary metabolites, the accumulation of phenolics in leaves was developmentally regulated during bud or first leaf expansion [67]. Moreover, in tobacco, the subsequent stages of flower development were correlated with tissue contents of polyphenols and activities of L-phenylalanine ammonia-lyase, polyphenoloxidase and peroxidase [68]. Also, a large amount of polyphenols was accumulated in tobacco pistil [68]– both reports are consistent with our observations. Although the occurrence of phenolic compounds was not the main objective of our study, it cannot be excluded that DNA methylation changes also affect polyphenolic accumulation, as was reported for Salvia sp [69].

DNA methylation, one of the most researched epigenetic modifications, has been linked with regulating chromatin structure and gene expression [16, 25]. It is involved in various developmental processes in flowers, including developing the floral organs, regulating the flowering time, and floral patterning [19]. DNA methylation impacts floral development by modulating the expression of flowering genes via epigenetic alterations. Promoter methylation typically correlates with the suppression of gene expression, while the effects of gene body methylation remain ambiguous, with studies reporting both positive and negative associations [20, 70, 71]. A conducted study reported that DNA methylation in non-promoter, intergenic regions and gene bodies facilitates gene expression. Non-promoter DNA methylation appears crucial for preserving the active chromatin states of genes [26, 72, 73]. It has been previously demonstrated that this epigenetic mark is essential in flower development and that its enrichment fluctuates in various developmental stages [28, 35, 41]. Here, we demonstrated the differences in the level of DNA methylation in flowers at various developmental stages in two buckwheat species, F. tataricum with homostylous flower type and F. esculentum with heterostylous flowers. Only the nectaries from open flowers had higher values of DNA methylation in comparison to results from closed flowers. DNA methylation values decreased in petals, stigmas and ovules (except for Thrum morph) of each analysed variant in open, developed flowers. In studies conducted on lotus (Nelumbo nucifera), the stamen petaloid exhibited a global decrease in DNA methylation levels [74]. Studies comparing DNA methylation patterns in ovules of the female-sterile rice and the wild-type, displayed a slightly lower whole-genome methylation level [75]. Research on hazel ovaries reported the reduction in DNA methylation levels of the ovule after pollination, indicating that the epigenetic mark is a crucial player in post-pollination stages [76]. It is known that stigma and pollen molecular cross-talk are crucial for successful pollination [77, 78]. During flower maturation, the papillae and sub-papillae cells of the stigma undergo histo- and biochemical changes to prepare for pollen reception, which involves cell loosening or exudate synthesis [79]. In our research, we observed a significant decrease in DNA methylation levels between stigmas of closed and open F. esculentum and F. tataricum flower; however, the divergence between DNA methylation levels in Pin and Thrum was the most prominent. This trait may point to differences between homo- and distylous species.

Studies conducted on Lilium longiflorum cv. Gelria focused on the analysis of the DNA methylation in pollen and reported that the vegetative nucleus in mature pollen grains was heavily methylated and that dramatic nonreplicative demethylation occurred during the pollen tube’s development. In conclusion, it was established that DNA hypomethylation ensures the survival of pollen grains without external sources of nutrients until they reach the stigma [80]. Since open flowers in Fagopyrum already have produced pollen grains, this might have contributed to the global reduction in DNA methylation levels. The reduction in DNA methylation is correlated with the down-regulation of genes coding DNA methyltransferases. In the Thrum open flower, the expression of MET1 was over 45 times lower in the open than in closed Thrum flowers.

In contrast to Thrum flowers, the expression of MET1 in open Pin flowers was the same as in closed ones. It seems that in Pin open flowers, the decrease in DNA methylation level might be related to the down-regulation of genes other than MET1 methyltransferase, such as MET2 or CMET3. The differential expression of methyltransferases is observed between various tissues [81]; for example, in rice, the transcript of MET1b accumulates more abundantly than those of MET1a in many different tissues [82]. In our research, the decreased transcript level of the rest of the analysed genes was observed in open Pin flowers. ROS1 encodes a nuclear protein containing an endonuclease III domain, exhibiting bifunctional DNA glycosylase/lyase activity specifically targeting methylated DNA while leaving unmethylated DNA unaffected [27]. It might explain why the levels of ROS1 expression are reduced in open Fagopyrum flowers, which exhibit reduced DNA methylation levels. Lower transcript accumulation of genes engaged in DNA demethylation was observed in open flowers of all analysed flower types and species. The DNA methylation level of analysed tissue is an interplay between DNA methylation and active demethylation. Methylation strongly affects flowering-related genes’ expression and mobility of transposons [83]. Effective demethylation occurs on the maternal central cell before fertilisation, and in the endosperm, maternal alleles are less methylated than paternal ones [84]. Among others, imprinted genes such as FLOWERING WAGENINGEN (FWA), MEDEA (MEA), and FERTILIZATION-INDEPENDENT SEED 2 (FIS2) are methylated and silenced in the spermatic cell, so after fertilisation, only maternal alleles are expressed in the endosperm [85]. Thus, the methylation and demethylation of specific DNA sequences are essential in regulating genes for the proper development of flowers and embryos.

Acquisition of the flowering abilities and subsequent development is precisely governed by the synchronised and specific expression patterns of microRNAs associated with flowering in plants [86, 87]. MicroRNAs represent a new category of intrinsic molecules that modulate gene expression, particularly during the pathways of flower development in plants [83, 88, 89]. Research has shown that the relative expression levels of several potential miRNAs remain notably consistent while governing the flowering time phenotype in a spatial and temporal context [90].

Conclusions

The current study illustrates the distinct DNA methylation patterns observed between the Pin and Thrum flower components of F. esculentum, as compared to self-pollinating species like F. tataricum, offering insights into the potential role of this epigenetic modification in the floral development of Fagopyrum species. Decreased overall DNA methylation and expression of genes connected with that epigenetic mark in open, developed flowers of both species might suggest that the demethylation is required to activate the expression of genes involved in floral development. Understanding the role of DNA methylation in flowers is crucial for unravelling the molecular mechanisms underlying floral development, adaptation to environmental changes, and the regulation of floral traits. It provides insights into how epigenetic modifications contribute to the diversity and complexity of flower morphology and function. Comparing the global methylation of DNA and expression of genes in both morphs of F. esculentum, as well as with a closely related, homostylous species, F. tataricum can lead to deciphering and understanding the differences between them and influence the self-incompatibility phenomenon to create self-compatible varieties. Our analyses offer insights into the potential roles of DNA methylation in gene expression and serve as valuable resources for further exploration of the genetic pathways governing flower development.

Data availability

All data generated or analysed during this study are included in this published article (and its additional information files).

Abbreviations

BSA:

Bovine albumin serum

CMT3:

CHROMOMETHYLASE 3

DAPI:

4’,6-diamidyno-2-fenyloindol

DME :

Demeter, 5-methylcytosine glycosylase

DML2 :

DME-like 2, 5-methylcytosine glycosylase

DML3 :

DME-like 3, 5-methylcytosine glycosylase

DRMs :

DOMAINS REARRANGED METHYLTRANSFERASEs

MET1 :

DNA METHYLTRANSFERASE 1

PBS:

1x phosphate-buffered saline

ROS1 :

Repressor of silencing 1, 5-methylcytosine glycosylase

S-ELF3 :

S-LOCUS EARLY FLOWERING 3

TBO:

Toluidine blue O

References

  1. Jha R, Zhang K, He Y, Mendler-Drienyovszki N, Magyar-Tábori K, Quinet M, Germ M, Kreft I, Meglič V, Ikeda K et al. Global nutritional challenges and opportunities: Buckwheat, a potential bridge between nutrient deficiency and food security. Trends Food Sci Technol 2024, 145.

  2. Tomasiak A, Zhou M, Betekhtin A. Buckwheat in tissue culture research: current status and future perspectives. Int J Mol Sci 2022, 23(4).

  3. Kwon SJ, Roy SK, Choi JY, Park JH, Cho SW, Sarker K, Woo SH. Recent research updates on functional components in buckwheat. J Agr Sci. 2018;34(1):1–8.

    Google Scholar 

  4. Sinkovič L, Kokalj D, Vidrih R, Meglič V. Milling fractions fatty acid composition of common (Fagopyrum esculentum Moench) and tartary (Fagopyrum tataricum (L.) Gaertn) buckwheat. J Stored Prod Res 2020, 85.

  5. Ji X, Han L, Liu F, Yin S, Peng Q, Wang M. A mini-review of isolation, chemical properties and bioactivities of polysaccharides from buckwheat (Fagopyrum Mill). Int J Biol Macromol. 2019;127:204–9.

    Article  CAS  PubMed  Google Scholar 

  6. Bonafaccia MM, Kreft I. Composition and technological properties of the flour and bran from common and tartary buckwheat. Food Chem. 2003;80:9–15.

    Article  CAS  Google Scholar 

  7. Gonçalves FMF, Debiage RR, Gonçalves da Silva RM, Porto PP, Yoshihara E, Cosendey E. Toledo De Mello Peixoto EC: Fagopyrum esculentum Moench: a crop with many purposes in agriculture and human nutrition. Afr J Agric Res. 2016;11(12):983–9.

    Article  Google Scholar 

  8. Huda MN, Lu S, Jahan T, Ding M, Jha R, Zhang K, Zhang W, Georgiev MI, Park SU, Zhou M. Treasure from garden: bioactive compounds of buckwheat. Food Chem. 2021;335:127653.

    Article  CAS  PubMed  Google Scholar 

  9. Liu M, Sun W, Ma Z, Guo C, Chen J, Wu Q, Wang X, Chen H. Integrated network analyses identify MYB4R1 neofunctionalization in the UV-B adaptation of Tartary buckwheat. Plant Commun. 2022;3(6):100414.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Sofi SA, Ahmed N, Farooq A, Rafiq S, Zargar SM, Kamran F, Dar TA, Mir SA, Dar BN, Mousavi Khaneghah A. Nutritional and bioactive characteristics of buckwheat, and its potential for developing gluten-free products: an updated overview. Food Sci Nutr. 2023;11(5):2256–76.

    Article  CAS  PubMed  Google Scholar 

  11. Matsui K, Yasui Y. Genetic and genomic research for the development of an efficient breeding system in heterostylous self-incompatible common buckwheat (Fagopyrum esculentum). Theor Appl Genet. 2020;133(5):1641–53.

    Article  CAS  PubMed  Google Scholar 

  12. Yasui Y, Hirakawa H, Ueno M, Matsui K, Katsube-Tanaka T, Yang SJ, Aii J, Sato S, Mori M. Assembly of the draft genome of buckwheat and its applications in identifying agronomically useful genes. DNA Res. 2016;23(3):215–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Cawoy V, Deblauwe V, Halbrecq B, Ledent JF, Kinet JM, Jacquemart AL. Morph differences and honeybee morph preference in the distylous species Fagopyrum esculentum Moench. IJPS. 2006;167(4):853–61.

    Google Scholar 

  14. Quinet M, Cawoy V, Lefevre I, Van Miegroet F, Jacquemart AL, Kinet JM. Inflorescence structure and control of flowering time and duration by light in buckwheat (Fagopyrum esculentum Moench). J Exp Bot. 2004;55(402):1509–17.

    Article  CAS  PubMed  Google Scholar 

  15. Fawcett JA, Takeshima R, Kikuchi S, Yazaki E, Katsube-Tanaka T, Dong Y, Li M, Hunt HV, Jones MK, Lister DL, et al. Genome sequencing reveals the genetic architecture of heterostyly and domestication history of common buckwheat. Nat Plants. 2023;9(8):1236–51.

    Article  CAS  PubMed  Google Scholar 

  16. Zang Y, Xie L, Su J, Luo Z, Jia X, Ma X. Advances in DNA methylation and demethylation in medicinal plants: a review. Mol Biol Rep. 2023;50(9):7783–96.

    Article  CAS  PubMed  Google Scholar 

  17. Zhang D, Li YY, Zhao X, Zhang C, Liu DK, Lan S, Yin W, Liu ZJ. Molecular insights into self-incompatibility systems: from evolution to breeding. Plant Commun. 2024;5(2):100719.

    Article  CAS  PubMed  Google Scholar 

  18. Lin PY, Chang YT, Huang YC, Chen PY. Estimating genome-wide DNA methylation heterogeneity with methylation patterns. Epigenet Chromatin. 2023;16(1):44.

    Article  CAS  Google Scholar 

  19. Lucibelli F, Valoroso MC, Aceto S. Plant DNA methylation: an epigenetic mark in development, environmental interactions, and evolution. Int J Mol Sci 2022, 23(15).

  20. Zemach A, McDaniel I, Silva P, Zilberman D. Genome-wide evolutionary analysis of eukaryotic DNA methylation. Science. 2010;328(5980):916–9.

    Article  CAS  PubMed  Google Scholar 

  21. Omony J, Nussbaumer T, Gutzat R. DNA methylation analysis in plants: review of computational tools and future perspectives. Brief Bioinform. 2020;21(3):906–18.

    Article  CAS  PubMed  Google Scholar 

  22. Bednarek PT, Orłowska R. Plant tissue culture environment as a switch-key of (epi)genetic changes. PCTOC. 2019;140(2):245–57.

    Article  Google Scholar 

  23. Chan S, Henderson I, Jacobsen S. Gardening the genome: DNA methylation in Arabidopsis thaliana. Nat Rev Genet. 2005;6:351–60.

    Article  CAS  PubMed  Google Scholar 

  24. Bennett M, Cleaves K, Hewezi T. Expression patterns of DNA methylation and demethylation genes during plant development and in response to phytohormones. Int J Mol Sci 2021, 22(18).

  25. Bird A. DNA methylation patterns and epigenetic memory. Genes Dev. 2002;16(1):6–21.

    Article  CAS  PubMed  Google Scholar 

  26. He XJ, Chen T, Zhu JK. Regulation and function of DNA methylation in plants and animals. Cell Res. 2011;21(3):442–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Gong Z, Morales-Ruiz T, Ariza RR, Roldan-Arjona T, David L, Zhu JK. ROS1, a repressor of transcriptional gene silencing in Arabidopsis, encodes a DNA glycosylase/lyase. Cell. 2002;111:803–14.

    Article  CAS  PubMed  Google Scholar 

  28. Yang H, Chang F, You C, Cui J, Zhu G, Wang L, Zheng Y, Qi J, Ma H. Whole-genome DNA methylation patterns and complex associations with gene structure and expression during flower development in Arabidopsis. Plant J. 2015;81(2):268–81.

    Article  CAS  PubMed  Google Scholar 

  29. Finnegan EJ, Genger RK, Kovac K, Peacock WJ, Dennis ES. DNA methylation and the promotion of flowering by vernalization. Plant Biol. 1998;95:5824–9.

    CAS  Google Scholar 

  30. Kumari P, Khan S, Wani IA, Gupta R, Verma S, Alam P, Alaklabi A. Unravelling the role of epigenetic modifications in development and reproduction of angiosperms: a critical appraisal. Front Genet. 2022;13:819941.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Meijón M, Feito I, Valledor L, Rodríguez R, Cañal MJ. Dynamics of DNA methylation and histone H4 acetylation during floral bud differentiation in azalea. BMC Plant Biol 2010, 10(10).

  32. Ruiz-Garcia L, Cervera MT, Martinez-Zapater JM. DNA methylation increases throughout Arabidopsis development. Planta. 2005;222(2):301–6.

    Article  CAS  PubMed  Google Scholar 

  33. Fulnecek J, Matyasek R, Votruba I, Holy A, Krizova K, Kovarik A. Inhibition of SAH-hydrolase activity during seed germination leads to deregulation of flowering genes and altered flower morphology in tobacco. Mol Genet Genomics. 2011;285(3):225–36.

    Article  CAS  PubMed  Google Scholar 

  34. Wang N, Ning S, Wu J, Tagiri A, Komatsuda T. An epiallele at cly1 affects the expression of floret closing (cleistogamy) in barley. Genetics. 2015;199(1):95–104.

    Article  CAS  PubMed  Google Scholar 

  35. Zou JJ, Cai X, Yang J, Zeng X, Liu DX, Huang S, Chen X, Yang QY, Wang C, Chen H. DNA hypomethylation mediates flower opening and senescence in sweet osmanthus through auxin and ethylene responsive pathways. Postharvest Biol Technol 2023, 198.

  36. Zhang S, Jiang H, Zhao H, Korpelainen H, Li C. Sexually different physiological responses of Populus cathayana to nitrogen and phosphorus deficiencies. Tree Physiol. 2014;34(4):343–54.

    Article  CAS  PubMed  Google Scholar 

  37. Sharifi A, Oizumi K, Kubota S, Bagheri A, Shafaroudi SM, Nakano M, Kanno A. Double flower formation in Tricyrtis macranthopsis is related to low expression of AGAMOUS ortholog gene. Sci Hortic. 2015;193:337–45.

    Article  CAS  Google Scholar 

  38. Wang Y, Zheng PC, Liu PP, Song XW, Guo F, Li YY, Ni DJ, Jiang CJ. Novel insight into the role of withering process in characteristic flavor formation of teas using transcriptome analysis and metabolite profiling. Food Chem. 2019;272:313–22.

    Article  CAS  PubMed  Google Scholar 

  39. Liu XJ, Chuang YN, Chiou CY, Chin DC, Shen FQ, Yeh KW. Methylation effect on chalcone synthase gene expression determines anthocyanin pigmentation in floral tissues of two Oncidium orchid cultivars. Planta. 2012;236(2):401–9.

    Article  CAS  PubMed  Google Scholar 

  40. Wu X, Zhou Y, Yao D, Iqbal S, Gao Z, Zhang Z. DNA methylation of LDOX gene contributes to the floral colour variegation in peach. J Plant Physiol 2020, 246–7:153116.

  41. Xing L, Li Y, Qi S, Zhang C, Ma W, Zuo X, Liang J, Gao C, Jia P, Shah K, et al. Comparative RNA-sequencing and DNA methylation analyses of apple (Malus domestica Borkh.) Buds with diverse flowering capabilities reveal novel insights into the regulatory mechanisms of flower bud formation. Plant Cell Physiol. 2019;60(8):1702–21.

    Article  CAS  PubMed  Google Scholar 

  42. Wolny E, Braszewska-Zalewska A, Hasterok R. Spatial distribution of epigenetic modifications in Brachypodium distachyon embryos during seed maturation and germination. PLoS ONE. 2014;9(7):e101246.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Braszewska-Zalewska AJ, Wolny EA, Smialek L, Hasterok R. Tissue-specific epigenetic modifications in root apical meristem cells of Hordeum vulgare. PLoS ONE. 2013;8(7):e69204.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Pawley JB. Handbook of biological confocal microscopy: Third edition. 2006.

  45. Wilson T. Resolution and optical sectioning in the confocal microscope. J Microsc 2011, 244.

  46. Inacio V, Martins MT, Graca J, Morais-Cecilio L. Cork oak young and traumatic periderms show PCD typical chromatin patterns but different chromatin-modifying genes expression. Front Plant Sci. 2018;9:1194.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Bankhead P. Analyzing fluorescence microscopy images with ImageJ. ImageJ; 2014.

  48. Abràmoff MD, Magalhães PJ, Ram SJ. Image processing with imageJ. Biophotonics Int 2004, 11.

  49. Bourne R. ImageJ. Fundamentals of Digital Imaging in Medicine 2010.

  50. Collins TJ. ImageJ for microscopy. Biotechniques 2007, 43.

  51. Hartig SM. Basic image analysis and manipulation in imageJ. Curr Protoc Mol Biol 2013.

  52. Schroeder AB, Dobson ETA, Rueden CT, Tomancak P, Jug F, Eliceiri KW. The ImageJ ecosystem: open-source software for image visualization, processing, and analysis. Prot Sci 2021, 30.

  53. Ferreira T, Rasband W. ImageJ user guide user guide ImageJ. Image J user Guide 2012, 146r.

  54. Martin G. R Studio. In: An Introduction to Programming with R 2021.

  55. Mendiburu F. Statistical procedures for agricultural research. In: Package ‘agricolae’ 2022.

  56. Kassambara A. ‘ggplot2’ Based Publication Ready Plots. In., 0.5.0 edn; 2022.

  57. Toro ID, Dickson K, Hakes AS, Newman SL. Undergraduate Biostatistics & Data Science Introduction Using R, R Studio & the Tidyverse. Am Biol Teach 2022.

  58. Team RC. R: A Language and environment for statistical computing. In. Vienna, Austria: R foundation for statistical computing; 2022.

  59. Team RS. RStudio: Integrated development environment for R. In: RStudio Boston, Massachusetts: R Studio; 2022.

  60. Wickham H. Easily Install and Load the ‘Tidyverse’. In: Package ‘tidyverse’ 2022.

  61. Tomasiak A, Sala K, Braszewska A. Immunostaining for epigenetic modifications in Fagopyrum calli. In: Buckwheat Methods in Molecular Biology Edited by Betekhtin A, Pinski, A., vol. 2791. Humana, New York: Springer 2024: 15–22.

  62. Tomasiak A, Berg LS, Sala K, Braszewska A. Quantitative analysis of epigenetic modifications in Fagopyrum nuclei with confocal microscope, ImageJ, and R Studio. In: Buckwheat Methods in Molecular Biology Edited by Betekhtin A, Pinski, A., vol. 2791. Humana, New York: Springer; 2024: 23–33.

  63. He Q, Ma D, Li W, Xing L, Zhang H, Wang Y, Du C, Li X, Jia Z, Li X, et al. High-quality Fagopyrum esculentum genome provides insights into the flavonoid accumulation among different tissues and self-incompatibility. J Integr Plant Biol. 2023;65(6):1423–41.

    Article  CAS  PubMed  Google Scholar 

  64. Ueno M, Yasui Y, Aii J, Matsui K, Sato S, Ota T. Genetic analyses of the heteromorphic self-incompatibility (S) lcus in buckwheat. In: Mol Breed Nutritional Aspects Buckwheat 2016: 411–21.

  65. Crozier A, Jaganath I, Clifford MN. Phenols, polyphenol and tannins: an overview. In: Plant secondary metabolites: occurence, structure and role in human diet 2007: 1–24.

  66. Nesovic M, Gasic U, Tosti T, Horvacki N, Nedic N, Sredojevic M, Blagojevic S, Ignjatovic L, Tesic Z. Distribution of polyphenolic and sugar compounds in different buckwheat plant parts. RSC Adv. 2021;11(42):25816–29.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  67. Jiang X, Liu Y, Li W, Zhao L, Meng F, Wang Y, Tan H, Yang H, Wei C, Wan X, et al. Tissue-specific, development-dependent phenolic compounds accumulation profile and gene expression pattern in tea plant [Camellia sinensis]. PLoS ONE. 2013;8(4):e62315.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Sheen SJ. Changes in amount of polyphenols and activity of related enzymes during growth of tobacco flower and capsule. Plant Physiol. 1973;51:839–44.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. He X, Chen Y, Xia Y, Hong X, You H, Zhang R, Liang Z, Cui Q, Zhang S, Zhou M, et al. DNA methylation regulates biosynthesis of tanshinones and phenolic acids during growth of Salvia miltiorrhiza. Plant Physiol. 2024;194(4):2086–100.

    Article  CAS  PubMed  Google Scholar 

  70. Zhang X, Yazaki J, Sundaresan A, Cokus S, Chan SW, Chen H, Henderson IR, Shinn P, Pellegrini M, Jacobsen SE, et al. Genome-wide high-resolution mapping and functional analysis of DNA methylation in Arabidopsis. Cell. 2006;126(6):1189–201.

    Article  CAS  PubMed  Google Scholar 

  71. Zilberman D, Gehring M, Tran R, Ballinger T, Henikoff S. Genome-wide analysis of Arabidopsis thaliana DNA methylation uncovers an interdependence between methylation and transcription. Nat Genet. 2007;39:61–9.

    Article  CAS  PubMed  Google Scholar 

  72. Wu H, Coskun V, Tao J, Xie W, Ge W, Yoshikawa K, Li E, Zhang Y, Sun YE. Dnmt3a-dependent nonpromoter DNA methylation facilitates transcription of neurogenic genes. Science. 2010;329(5990):444–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Feng S, Cokus SJ, Zhang X, Chen PY, Bostick M, Goll MG, Hetzel J, Jain J, Strauss SH, Halpern ME, et al. Conservation and divergence of methylation patterning in plants and animals. Proc Natl Acad Sci U S A. 2010;107(19):8689–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Lin Z, Liu M, Damaris RN, Nyong’a TM, Cao D, Ou K, Yang P. Genome-wide DNA methylation profiling in the lotus (Nelumbo nucifera) flower showing its contribution to the stamen petaloid. Plants 2019, 8(5).

  75. Liu H, Wu Y, Cao A, Mao B, Zhao B, Wang J. Genome-wide analysis of DNA methylation during ovule development of female-sterile rice fsv1. G3 2017, 7(11):3621–3635.

  76. Cheng YQ, Wang J, Liu JF, Zhao YX, Geng WT. Z. H. D.: Analysis of ovary DNA methylation during delayed fertilization in hazel using the methylation-sensitive amplification technique. Acta Physiol Plant 2015, 37(231).

  77. Cawoy V, Ledent JF, Kinet JM, Jacquemart AL. Floral biology of common buckwheat (Fagopyrum esculentum Moench). EJPSB. 2009;3(1):1–9.

    Google Scholar 

  78. Gao X-Q, Zhu D, Zhang X. Stigma factors regulating self-compatible pollination. Front Biol. 2010;5(2):156–63.

    Article  Google Scholar 

  79. Çetinbaş-Genç A, Yanık F, Vardar F. Histochemical and biochemical alterations in the stigma of Hibiscus syriacus (Malvaceae) during flower development. Caryologia. 2019;72(4):3–13.

    Google Scholar 

  80. Janousek B, Zluvova J, Vyskot B. Histone H4 acetylation and DNA methylation dynamics during pollen development. Protoplasma. 2000;211:116–22.

    Article  CAS  Google Scholar 

  81. Zhang H, Lang Z, Zhu JK. Dynamics and function of DNA methylation in plants. Nat Rev Mol Cell Biol. 2018;19(8):489–506.

    Article  CAS  PubMed  Google Scholar 

  82. Yamauchi T, Moritoh S, Johzuka-Hisatomi Y, Ono A, Terada R, Nakamura I, Iida S. Alternative splicing of the rice OsMET1 genes encoding maintenance DNA methyltransferase. J Plant Physiol. 2008;165(17):1774–82.

    Article  CAS  PubMed  Google Scholar 

  83. Kim WJ, Ahn JH. MicroRNA-target interactions: important signaling modules regulating flowering time in diverse plant species. Crit Rev Plant Sci. 2014;33(6):470–85.

    Article  CAS  Google Scholar 

  84. Klosinska M, Picard CL, Gehring M. Conserved imprinting associated with unique epigenetic signatures in the Arabidopsis genus. Nat Plants. 2016;2:16145.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  85. Luo MB, Dennis P, Peacock ES, Chaudhury WJ. Expression and parent-of-origin effects for FIS2, MEA, and FIE in the endosperm and embryo of developing Arabidopsis seeds. Proc Natl Acad Sci. 2000;97:10637–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  86. Arora S, Chaudhary B. Global expression dynamics and miRNA evolution profile govern floral/fiber architecture in the modern cotton (Gossypium). Planta. 2021;254(3):62.

    Article  CAS  PubMed  Google Scholar 

  87. Arora S, Pandey DK, Chaudhary B. Target-mimicry based diminution of miRNA167 reinforced flowering-time phenotypes in tobacco via spatial-transcriptional biases of flowering-associated miRNAs. Gene. 2019;682:67–80.

    Article  CAS  PubMed  Google Scholar 

  88. Hong Y, Jackson S. Floral induction and flower formation–the role and potential applications of miRNAs. Plant Biotechnol J. 2015;13(3):282–92.

    Article  CAS  PubMed  Google Scholar 

  89. Spanudakis E, Jackson S. The role of microRNAs in the control of flowering time. J Exp Bot. 2014;65(2):365–80.

    Article  CAS  PubMed  Google Scholar 

  90. Jin D, Wang Y, Zhao Y, Chen M. MicroRNAs and their cross-talks in plant development. J Genet Genomics. 2013;40(4):161–70.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The authors thank Lea Berg for preparing the R Studio script used in our analysis.

Funding

This research was funded by the National Science Centre, Poland. Research project OPUS-19 (no. reg. 2020/37/B/NZ9/01499 awarded to AB) and the Research Excellence Initiative of the University of Silesia in Katowice.

Author information

Authors and Affiliations

Authors

Contributions

Conceptualisation: AB; Methodology: KS-C, AT, KN, AP; Formal analysis: KS-C, AT, AB, KN, AP; Investigation: KS-C, AT, AB, KN, AP; Resources: AB; Writing—original draft: KS-C, AT, KN, AP; Writing—review & editing: KS-C, AT, KN, AP, AB; Visualisation: KS-C, AT, KN; Supervision: AB; Project administration: AB; Founding acquisition: AB. All authors have read and approved the final manuscript.

Corresponding authors

Correspondence to Katarzyna Sala-Cholewa or Alexander Betekhtin.

Ethics declarations

Ethics approval and consent to participate

The use of all plant materials in this study complies with relevant institutional, national, and international guidelines and legislation. Seeds of F. tataricum (sample k-17) are from the N. I. Vavilov Institute of Plant Genetic Resources collections, Saint Petersburg, Russia. The Plant Cytogenetic and Molecular Biology Group Institute of Biology, Biotechnology and Environmental Protection, Faculty of Natural Sciences, University of Silesia in Katowice, Poland, multiplied the obtained seeds. F. tataricum sample k-17 is a common cultivated landrace of F. tataricum, and seeds are available upon request from the publication’s authors. F. esculentum seeds of the Panda cultivar are commercially available and were supplied from the Malopolska Plant Breeding, Poland.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Additional File 1

: Graphical depiction of the experimental design

Additional File 2

: Excel table containing raw numerical data for DNA methylation

Additional File 3

: R-Studio script for statistical analysis and graphical depiction of the results

Additional File 4

: Excel table containing forward and reverse primers for the RT-qPCR analysis

Additional File 5

: BioRender certificate confirming the publication rights for Additional File

Additional File 6

: Histology of anthers from closed (a – c) and open (d – f) flowers (not included in DNA methylation analysis)

Additional File 7

: Selected images captured with the Olympus FV1000 confocal system

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Sala-Cholewa, K., Tomasiak, A., Nowak, K. et al. DNA methylation analysis of floral parts revealed dynamic changes during the development of homostylous Fagopyrum tataricum and heterostylous F. esculentum flowers. BMC Plant Biol 24, 448 (2024). https://doi.org/10.1186/s12870-024-05162-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-024-05162-w

Keywords