4.1. ABA catabolism
ABA plays vital roles in seed development and maturation, encompassing the accumulation of storage compounds, acquisition of desiccation tolerance, induction of dormancy, and suppression of precocious germination [
12,
17,
18,
19,
57,
58,
59]. However, to break dormancy and initiate germination, ABA needs to be catabolized, primarily through hydroxylation and conjugation. The primary ABA hydroxylation route is the ABA catabolic pathway (
Figure 5), which relies on the activities of CYP707A cytochromes P450, notably ABA 8'-hydroxylases [
60].
Initially, ABA is catalyzed by 8'-hydroxylase, converting it to 8'-hydroxy ABA (8'-OH ABA), an unstable intermediate [
61,
62]. This intermediate is then spontaneously rearranged into PA and subsequently reduced by PA reductase (PAR) to DPA [
24,
63]. The 9'-hydroxylation pathway, similar to 8'-hydroxylation, involves CYP707As enzymes and converts 9'-hydroxy ABA (9'-OH ABA) to neoPA with both 8'-C and 9'-C hydroxylations catalyzed by the same enzyme [
16]. Recently, Bai et al. (2022) [
24] identified a downstream catabolite of neoPA in the 9'-hydroxylation pathway as epi-neodihydrophaseic acid (
epi-neoDPA) and discovered the responsible enzyme, neoPA reductase 1 (NeoPAR1) (
Figure 5).
Our study examined ABA and ABA-related catabolites in embryonic axes of
P. sativum seeds before and after radicle protrusion. We found a decline in ABA content with a concurrent rise in its catabolites (PA, DPA, and neoPA) (
Figure 2). Intriguingly, PA, similar to ABA, can regulate stomatal closure and suppress seed germination [
64,
65]. Weng et al. (2016) Weng
et al. (2016) demonstrated that PA functions as a signaling molecule through ABA receptors. Similar ABA-like hormonal activity was observed for neoPA, but not for
epi-neoDPA [
24]. Additionally, altered seed germination patterns were noted in neo-PAR1 mutant and overexpression lines, implicating the ABA catabolic pathway as a critical regulatory mechanism during the seed-to-seedling transition [
24]. Despite reduced ABA levels, the accumulation of its catabolic products (PA, neoDPA) in embryonic axes suggests a continued regulatory influence
via ABA receptors.
4.2. Annotation of ABA-associated DEGs
In our prior RNA sequencing-based transcriptomic analysis of pea embryonic axes isolated from seeds before and after radicle protrusion [
4], we identified 24,184 DEGs, with 2,101 showing notably higher expression. This work extends that analysis by focusing on ABA-associated DEGs (ABA-DEGs). Of the 70 ABA-DEGs annotated, 46 genes up-regulated and 24 genes down-regulated more than 4-fold after radicle protrusion (
Figure 3).
The up-regulated ABA-DEGs predominantly pertained to cellular signaling, stress resistance, membrane transporters, and transcription factors that regulate seedling development. For instance,
Psat6g199400, encoding RD29B/LTI65 which responds to water deprivation, was up-regulated 4.5-fold. This gene's promoter region contains two ABA-responsive elements (ABREs) that required as
cis-acting elements for the dehydration-responsive expression of RD29B/LTI65 [
66,
67] Similarly,
Psat4g146960, encoding ANNEXIN4, a calcium-binding protein involved in drought and other stress responses [
68,
69], showed a 9-fold increase in expression (
Table S1).
Among the up-regulated genes were those coding for membrane transporters like
Psat4g117800 (P-ATPase) and
Psat4g184760 (potassium channel AKT2/3). P-type ATPases play a role in ion transport across membranes, utilizing ATP for transmembrane conformational changes [
70,
71]. Additionally,
Psat2g121520, encoding TCP15, a transcription factor implicated in cell expansion and proliferation [
72,
73] was up-regulated 8-fold. The TCP proteins, known as TEOSINTE BRANCHED 1 (TB1) in maize, CYCLOIDEA (CYC) in
Anthirrinum majus, and PCF in rice [
74] have been linked to various developmental processes, including light-induced cotyledon opening in
Arabidopsis [
75].
Conversely,
down-regulated ABA-DEGs included genes central to ABA signaling (
ABI3, ABI4, ABI5) and those involved in water deprivation response (
LEA14, RD22, HVA22, PER1, and
LTI65) (
Table S1). Seed germination is governed by the antagonistic balance of ABA/GA, with ABA catabolism preceding GA synthesis and activation [
5,
7,
17]. Key ABA signaling genes
ABI3, ABI4, and
ABI5 encode the TFs featuring B3, AP2, and bZIP domains which control the expression of ABA-responsive genes crucial for seed maturation, dormancy, longevity, germination, and post-germination growth ([
12,
16,
76,
77,
78]. .
ABI5 encodes a member of the basic leucine zipper TF family and involved in ABA signaling in seeds by acting as a signal integrator between ABA and other hormones [
41,
79,
80]. The Arabidopsis
abi5 mutants have pleiotropic defects in ABA response, including reduced sensitivity to ABA, inhibition of germination, and altered expression of some ABA-regulated genes [
81,
82]. Notably,
Psat3g033680, encoding ABI5, exhibited a 22-fold downregulation after radicle protrusion.
ABI4 was shown to be a key integration node for multiple signals, participating in critical transition steps during plant ontogenesis [
83,
84,
85]. In dormant seeds, ABI4 acts as a repressor of ABA catabolism by binding to the promoter of CYP707A, being the main enzyme of ABA catabolism [
86]. Thereby, ABA and GAs can antagonistically modify the expression and stability of
ABI4, suggesting the existence of regulatory loops [
83]. In germinating seeds, ABI4 can regulate both ABA synthesis and catabolism. Some authors suggest that ABI4 plays as a key regulator of the balance between ABA and GAs in seeds at post-germination stages [
83,
85]. In our study, the level of the
Psat2g031240 gene encoding ABI4 was decreased 21-fold.
ABI3 encoded AP2/B3-like transcriptional factor family protein [
87]. ABI3 belongs to the LAFL regulatory network where interacts with LEAFY COTYLEDON1 (LEC1), ABSCISIC ACID INSENSITIVE3 (ABI3), FUSCA3 (FUS3), and LEC2 [
31,
32]. The LAFL network is a positive regulator of seed dormancy and need to be suppressed for seed germination. Together, ABI3, FUS3, and LEC1 are involved in sensitivity of seeds to ABA, and regulate the expression of the 12S storage protein gene family [
88]. In addition, both FUS3 and LEC1 positively regulate the ABI3 protein abundance in the seeds [
89]. The expression of
Psat3g142040 encoded ABI3 was decreased 21-fold.
We also found the downregulation of ABA-dependent genes involved in response to water deprivation (
LEA14,
RD22,
HVA22,
PER1 and
LTI65) (
Table S1). In accordance with our findings,
Psat7g085840 encoding peroxiredoxin1 (PER1),
Psat0s2227g0040 encoding protein LTI65/78, and
Psat0s2780g0040 encoding late embryogenesis abundant (LEA) protein were down-regulated 20-30-fold. Peroxiredoxins are thiol-dependent antioxidants containing one (1-Cys) or two (2-Cys) conserved Cys residues [
90].
PER1 encodes a 1-Cys peroxiredoxin (PER1) protein that accumulates during seed development but rapidly disappears upon germination [
91]. PER1 is involved in the quenching reactive oxygen species (ROS) during late maturation, dormancy, and early germination, thereby maintaining seed viability [
91,
92,
93]. The low-temperature-induced (LTI) protein family is associated with responses to abiotic stresses. In Arabidopsis, the homologous genes
RD29A (LTI78) and
RD29B (LTI65) are induced by cold, drought, salt, and abscisic acid [
66]. Most LEA genes have ABA response elements in promoters and their expression can be induced not only by ABA, but also by cold or drought. The desiccation-related protein LEA14 belongs to the group II LEA proteins, also known as dehydrins [
94]. LEA14 was induced in response to salt and low temperature [
95].
4.3. Epigenetic regulation of PsABI3, PsABI4, and PsABI5 genes based on DNA promoter methylation.
Major transitions in the plant life cycle require fine-tuned regulation at the molecular and cellular levels. Epigenetic regulation, particularly DNA methylation, is crucial for maintaining genome stability in plants by inhibiting transposable element movement and modulating gene expression during development and stress responses [
8,
96]. DNA methylation patterns in seeds undergo significant changes during development and germination [
8,
46,
47,
50,
52,
97,
98]. DNA methylation (mC) occurs in three sequence contexts (CG, CHG, and CHH) and refers as addition a methyl group to the C5 position of cytosine to form 5-methylcytosine [
96]. Methylation of CHH sites notably increases from early to late stages of seed development, then decreases during germination [
8,
47,
48]. Two DNA methylases, RdDM (RNA directed DNA methylation) and CMT2 (DOMAINS REARRANGED METHYLTRANSFERASE 2), responsible for methylating CHH sites in developing seeds, are inactivated during germination [
51,
52]. In contrast, CG and CHG methylation patterns are relatively stable throughout seed development [
46,
98,
99]. Therefore, monitoring the level of 5-methylcytosine (m5C) is considered as a universal marker for seeds at the different stages of their ontogenesis [
97].
Our study reveals that during the transition from germination to post-germination, expression of key ABA signaling pathway genes (
ABI3, ABI4, and
ABI5) is markedly suppressed. We analyzed the DNA promoter methylation profiles of
PsABI3, PsABI4, and
PsABI5 to understand their epigenetic regulation. Contrary to our expectations of low promoter methylation levels based on their expression before radicle protrusion [
4], we observed high methylation levels both before and after this developmental stage (
Figure 4). Notably, approximately one-third of the
PsABI3 gene promoter region showed reduced methylation. However, this region might belong to the 5'-UTR as per the Pea Genome International Consortium version 1a (
Figure S3, pink).
We further investigated the coincidence of epigenetic marks with transcription factor binding sites in the promoters of these genes, using PlantPAN 3.0 and PCBase, focusing on stress and hormone response motifs.
PsABI5 showed numerous potential binding sites for LAFL network proteins, along with motifs associated with cold and water deprivation responses (
Figure S4). This finding aligns with the role of ABI5 as a major regulator of seed maturation and longevity in legumes [
41]. Our results suggest that epigenetic modifications, impacting the binding ability of
ABI3,
ABI4, and
ABI5 to DNA promoters, occur prior to the initiation of the seed transition from germination to post-germination.
Thus, our study provided an insight into the involvement of ABA in the transition in
P. sativum from germination to the post-germination stage when seeds turn into seedlings. The initiation of embryonic axis growth corresponds with changes in the abscisate profile: a decrease in ABA levels and an accumulation of its catabolites (PA, DPA, and neoPA), which possess hormonal activity similar to ABA [
24,
100]. Our in-depth analysis of ABA-DEGs revealed 46 up-regulated and 24 down-regulated genes with more than 4-fold changes. Most up-regulated ABA-DEGs were related to the regulation of seedling development. Most notably, the expression of
ABI3, ABI4, and
ABI5 was significantly downregulated, and their promoters exhibited a high level of methylation both before and after radicle protrusion. While ABA continues to be important, other regulators appear to be involved in the seed-to-seedling transition.