Skip to main content

Advertisement

HSI2/VAL1 PHD-like domain promotes H3K27 trimethylation to repress the expression of seed maturation genes and complex transgenes in Arabidopsis seedlings

Abstract

Background

The novel mutant allele hsi2-4 was isolated in a genetic screen to identify Arabidopsis mutants with constitutively elevated expression of a glutathione S-transferase F8::luciferase (GSTF8::LUC) reporter gene in Arabidopsis. The hsi2-4 mutant harbors a point mutation that affects the plant homeodomain (PHD)-like domain in HIGH-LEVEL EXPRESSION OF SUGAR-INDUCIBLE GENE2 (HSI2)/VIVIPAROUS1/ABI3-LIKE1 (VAL1). In hsi2-4 seedlings, expression of this LUC transgene and certain endogenous seed-maturation genes is constitutively enhanced. The parental reporter line (WTLUC) that was used for mutagenesis harbors two independent transgene loci, KanR and KanS. Both loci express luciferase whereas only the KanR locus confers resistance to kanamycin.

Results

Here we show that both transgene loci harbor multiple tandem insertions at single sites. Luciferase expression from these sites is regulated by the HSI2 PHD-like domain, which is required for the deposition of repressive histone methylation marks (H3K27me3) at both KanR and KanS loci. Expression of LUC and Neomycin Phosphotransferase II transgenes is associated with dynamic changes in H3K27me3 levels, and the activation marks H3K4me3 and H3K36me3 but does not appear to involve repressive H3K9me2 marks, DNA methylation or histone deacetylation. However, hsi2-2 and hsi2-4 mutants are partially resistant to growth inhibition associated with exposure to the DNA methylation inhibitor 5-aza-2′-deoxycytidine. HSI2 is also required for the repression of a subset of regulatory and structural seed maturation genes in vegetative tissues and H3K27me3 marks associated with most of these genes are also HSI2-dependent.

Conclusions

These data implicate HSI2 PHD-like domain in the regulation of gene expression involving histone modifications and DNA methylation-mediated epigenetic mechanisms.

Background

Transition from seed maturation to seed germination and seedling development involves a complex network of genetic and epigenetic mechanisms that down-regulate the expression of seed maturation genes in seedlings [1] [6]. Seed maturation is under the control of a group of transcriptional activators including LEAFY COTYLEDON1 (LEC1 [7]), LEC1-LIKE (L1L [8]), ABSCISIC ACID INSENSITIVE3 (ABI3 [9]), FUSCA3 (FUS3 [10]) and LEC2 [11], which are collectively called the “LAFL network” [3]. The B3-domain containing transcriptional repressors HIGH-LEVEL EXPRESSION OF SUGAR-INDUCIBLE GENE2 (HSI2) /VP1/ABI3-LIKE1 (VAL1) and its homolog HSI2-LIKE1 (HSL1)/VAL2 act redundantly to repress ectopic activation of embryonic traits during seed germination and seedling development by the “LAFL network” of transcriptional activators [12] [16]. HSI2 was also shown to negatively regulate the expression of β-glucuronidase (GUS) or luciferase (LUC) reporters under the control of seed-maturation specific gene promoters in transgenic Arabidopsis seedlings and vegetative organs [17],[18]. Since many of the genes repressed by HSI2 in vegetative tissues are involved in the maturation phase of seed development, including desiccation tolerance, knock-out hsi2 mutant seedlings show enhanced tolerance to water deficit whereas the overexpression of HSI2 resulted in hypersensitivity to desiccation stress [19]. Recently, it was shown that both fus3 and lec2 loss of function mutants can completely suppress the embryonic phenotype of hsi2/hsl1 double mutant seedlings, while it is partially suppressed in abi3, lec1 and l1l mutants [15]. These results indicate that HSI2 and HSL1 function redundantly to repress the expression of these regulatory genes in seedlings to prevent ectopic expression of embryonic traits during seed germination and vegetative development.

Developmental regulation of gene expression in plants is affected by chromatin mediated epigenetic mechanisms that include DNA methylation, chromatin remodeling, histone variants, and histone modifications [20],[21]. DNA methylation at the 5′position of cytosine plays important roles in transcriptional silencing of transposons, repeat sequences, transgenes and transcribed genes [22]. In addition to DNA methylation, histone modifications also play a vital role in the regulation of both transposons and transcribed genes in plants. Methylation of various lysine residues in the N-terminal tail of histone H3 is a well characterized epigenetic mechanism. In Arabidopsis, mono- (me1), di- (me2) or tri- (me3) methylation of histone H3 occurs mainly at lysine 4 (K4), lysine 9 (K9), lysine 27 (K27) and lysine 36 (K36) [23]. H3K4me3 and H3K36me3 are enriched on actively transcribed genes whereas H3K27me3 marks are associated with developmental repression of transcribed genes. H3K9me2/3 marks, which are associated with DNA methylation and small interfering RNAs (siRNAs), are enriched in heterochromatic regions known to be involved in transcriptional silencing of transposons, repeat sequences and transgenes [23],[24].

HSI2 and HSL1 proteins were predicted to contain a PHD-like domain, a B3-DNA binding domain, a conserved cysteine and tryptophan residue-containing (CW) domain and an ethylene-responsive element binding factor-associated amphiphilic repression (EAR) motif [3],[12],[14],[17],[25],[26]. Both CW and PHD protein domains are known to recognize methylated histone marks [23],[27] [30]. Hoppmann et al. [29] showed that the CW domain of HSI2 binds to H3K4me2 and H3K4me3 in vitro and, recently, it was reported that the HSL1 CW domain interacts with the histone deacetylase HDA19 to repress the “LAFL network” genes, including LEC1 and LEC2, by promoting histone deacetylation and the addition of H3K27me3 marks [31]. However, molecular and epigenetic mechanisms underlying the HSI2 PHD-like domain-mediated regulation of gene expression remain to be elucidated.

Previously, we reported a novel HSI2 allele, hsi2-4, in Arabidopsis that harbors a point mutation resulting in an amino acid substitution (C66Y) in the second zinc finger of the HSI2 PHD-like domain. The hsi2-4 mutant seedlings that carry a glutathione S-transferase F8::luciferase (GSTF8::LUC) reporter gene showed constitutively elevated transgene expression [14]. In addition to the LUC transgene, HSI2 PHD-like domain is required for the non-redundant repression of several seed-maturation genes in seedlings. These genes include those that encode both regulatory factors such as FUS3, and AGAMOUS-Like 15 (AGL15) and structural proteins that include cupin family storage protein, oleosins, late-embryogenesis-related proteins and seed storage albumins. Moreover, seed-specific genes that are de-repressed in hsi2-4 mutant seedlings are targets of H3K27me3 marks. Chromatin immunoprecipitation and quantitative PCR (ChIP-qPCR) analyses indicated that HSI2 PHD-like domain promotes H3K27me3 marks on transgene GSTF8 promoter and LUC coding sequences to repress transgene expression in parental GSTF8::LUC reporter (WTLUC) seedlings [14]. Both WTLUC and hsi2-4LUC mutant plants harbor two independent transgene loci [14]. One locus, located on chromosome IV, confers kanamycin resistance and luminescence, whereas the second locus, which is on chromosome V, confers only luminescence. Based on kanamycin sensitivity, the chromosome IV and chromosome V loci were named as KanR and KanS, respectively [14].

In this work, we show that HSI2 PHD-like domain represses LUC transgene expression from both KanR and KanS loci by promoting H3K27me3 marks but not DNA methylation and siRNA associated H3K9me2 marks. Expression of Neomycin Phosphotransferase II (NPTII) from the KanR locus is also partially suppressed in an HSI2-dependent mechanism. However, while our data indicate that DNA methylation and histone deacetylation are not involved in the transcriptional repression of transgene loci in WTLUC, the HSI2 PHD-like domain may play a role in the inhibition of seedling growth and development caused by DNA methylation inhibitor 5-aza-2′-deoxycytidine (5-azadC).

Results

Disruption of HSIPHD-like domain affects the expression of both KanRand KanStransgene loci

The GSTF8::LUC reporter construct contains a GSTF8 promoter sequence that controls the transcription of a luciferase expression cassette, along with an NPTII gene under control of the nopaline synthase promoter and terminator sequences, which confers kanamycin resistance in plants (Figure 1A). The parental WTLUC reporter line harbors two independent transgene insertion sites, KanR and KanS. The KanR locus was mapped to chromosome IV, while the KanS locus is located on chromosome V (Table 1) [14]. Active luciferase is expressed by both KanR and KanS loci, conferring a luminescent phenotype; however, only the KanR locus expresses NPTII; thus, plants that harbor only the KanR locus are resistant to kanamycin, while KanS plants are sensitive to this antibiotic.

Figure 1
figure1

Genomic structure of GSTF8::LUC transgene and luminescence imaging of WTLUCand hsi2-4 mutant seedlings harboring either KanRor KanStransgene locus or both. A. GSTF8::LUC transgene contain neomycin phosphotransferase (NPTII) coding sequences under the control of nopaline synthase (NOS) promoter and a modified luciferase (LUC+) coding sequences from firefly driven by glutathione S-transferase F8 (GSTF8) promoter conferring kanamycin resistance and luminescence expression respectively in plants. The 3′ ends of both NPTII and LUC+ coding sequences include NOS terminator sequences for transcriptional termination. B. Plants harboring either KanR or KanS transgene locus alone in the wild-type or in hsi2-4 mutant background were obtained by crossing of either WTLUC or hsi2-4LUC into Columbia-0 wild-type and homozygous lines were identified in F2 and F3 generations. Five days old seedlings of various genotypes grown on Murashige and Skoog media plates were imaged using cooled CCD camera after spraying with the substrate luciferin. Pseudocolor image indicates luminescence intensity from lowest (blue) to highest (white).

Table 1 Estimation of LUC copy number in Arabidopsis plants containing either Kan S or Kan R locus or both

To estimate the number of LUC copies at both KanR and KanS loci, real-time quantitative PCR (qPCR) was performed using genomic DNA from WTLUC, KanR and KanS plants. Since both KanR and KanS loci confer luminescence expression, we used PCR primers that are specific to the LUC coding sequences to estimate the copy numbers. The results show that KanS plants contain 2 copies of LUC whereas the KanR locus harbors 5 LUC copies. Independent analysis of WTLUC plants, which contain both KanR and KanS loci, showed seven copies of the LUC transgene (Table 1). Therefore, both KanR and KanS loci are complex and contain multiple copies of the GSTF8::LUC transgene.

Previously, we showed that disruption of the HSI2 PHD-like domain affects the expression of the KanR transgene locus [14] but the effect of this mutation on the KanS locus was not evaluated. Therefore, to further investigate whether the KanS transgene locus is also regulated by the HSI2 PHD-like domain mutation and investigate potential interactions between KanR and KanS transgene loci in WTLUC and hsi2-4LUC mutant plants, these two loci were separated by crossing plants of the WTLUC reporter line and the hsi2-4LUC mutant line into Col-0 wild-type Arabidopsis and subsequent selection for homozygous WTLUC and hsi2-4 lines that carry either the KanR or KanS reporter gene locus.

Comparison of luciferase expression in seedlings homozygous for the isolated KanR and KanS transgene loci in the wild-type background showed that KanR seedlings had higher luminescence signals (Figure 1B) and steady state levels of LUC mRNA (Figure 2) than KanS seedlings. This is in agreement with the relative number of transgene copies at these loci. However, in spite of carrying more luciferase transgene copies than KanR seedlings, WTLUC seedlings, showed significantly lower luminescence signal and LUC transcript levels. On the other hand, analysis of the expression of these transgenes in the hsi2-4 background showed strongly enhanced luciferase expression in all of the lines and the relative levels of both luminescence signal and LUC transcripts corresponded with transgene copy number, with highest levels seen in hsi2-4LUC seedlings and lowest levels in hsi2-4-KanS samples (Figures 1B and 2). This could indicate that, in a wild-type background, the presence of both the KanR and Kans loci may lead to stronger suppression of transgene expression but disruption of the HSI2 PHD-like domain affects the expression of LUC transgenes at both KanR and KanS loci similarly. Thus, the more complete HSI2-mediated repression of the GSTF8::LUC transgenes in WTLUC plants results in stronger relative activation of their expression in the presence of the hsi2-4 mutation.

Figure 2
figure2

Transcript levels of endogenous GSTF8 and transgenes in WTLUCand hsi2 mutants carrying either KanRor KanStransgene locus or both. Real-time reverse transcription quantitative PCR was used to determine the relative transcript levels of endogenous GSTF8, LUC and NPTII genes in five day old seedlings of various genotypes. GSTF8 produces two different transcripts with different fragment lengths by alternative start sites namely GSTF8-Long and GSTF8-Short[34]. Expression of GSTF8-Total represents transcripts from both GSTF8-Long and GSTF8-Short versions whereas GSTF8-Long expression level corresponds to GSTF8-Long transcript. EF1α was used for normalization. Data represent means (’sD) of two biological replicates with three technical replicates each. Significant differences in LUC transcript levels between the three luciferase reporter lines in the wild-type background and the respective hsi2-4 mutant background, determined using two-tailed Student’s t-test assuming unequal variances, are indicated by letters (a = p < 0.001 and b = p < 0.0001).

Alternative transcriptional start sites of the endogenous GSTF8 gene result in two different transcripts with different sizes: GSTF8-Long (GSTF8-L) and GSTF8-Short (GSTF8-S) [34]. To determine whether endogenous GSTF8 expression is altered in hsi2 mutant alleles, we performed qRT-PCR using various wild type lines (WTLUC, KanR and KanS) and hsi2 mutant lines (hsi2-2LUC, hsi2-4-KanR, hsi2-4-KanS and hsi2-4LUC). hsi2-2 is a loss-of-function mutant allele that carries a T-DNA insertion in the seventh exon of HSI2 gene (SALK_088606) [12] [14],[17]. To obtain the hsi2-2LUC line, GSTF8::LUC transgenes were introgressed into the hsi2-2 mutant background by genetic crossing. Expression of GSTF8-Total (GSTF8-T) represents both GSTF8-L and GSTF8-S transcripts whereas GSTF8-Long expression represents GSTF8-L transcripts only. As shown in Figure 2, levels of endogenous GSTF8-L and GSTF8-T transcripts were not significantly affected in hsi2-2LUC and hsi2-4LUC plants and NPTII expression was not detected in the KanS reporter line, consistent with the kanamycin sensitivity of these plants. NPTII transcripts were expressed at similar levels in the WTLUC and hsi2-4LUC seedlings but expression of NPTII in KanR seedlings was responsive to the HSI2 PHD-like domain point mutation (Figure 2). Steady-state levels of NPTII transcripts from the KanR locus were about 3-fold higher in hsi2-4 seedlings than in the wild-type background. Taken together, these results indicate that both the GSTF8::LUC and NOS::NPTII transgenes of the T-DNA cassette are partially suppressed by HSI2 in KanR seedlings; however the NPTII genes at the KanS locus may be fully silenced and/or contain loss-of-function mutations. Furthermore, in WTLUC seedlings where the KanS locus is present, expression of both LUC and the NPTII genes at the KanR is more strongly suppressed.

Luciferase expression is not affected by DNA methylation or histone deacetylation inhibitors

Complex transgenes with tandem repeats in plants are often subjected to DNA methylation and histone deacetylation mediated transcriptional gene silencing [35]. If the GSTF8::LUC transgene loci in the WTLUC plants are targets of DNA methylation, treatment of these seedlings with an inhibitor of DNA methylation should derepress the luminescence expression similar to that seen in hsi2-4 seedlings. Previous reports showed that treatment with 5 μM/mL of the DNA methylation inhibitor 5-aza-2′-deoxycytidine (5-azadC) was effective in derepressing the transcriptional silencing of auxin-responsive ß-GUS reporter lines [35]. Treatment of Arabidopsis seedlings with 7 μM/mL 5-azadC also caused global changes in gene expression and derepression of silenced transgenes [36],[37]. To investigate whether DNA methylation is involved in repressing GSTF8::LUC transgene expression, WTLUC, hsi2-2LUC and hsi2-4LUC seedlings were grown on media containing 5-azadC at various concentrations. Luminescence expression in WTLUC seedlings was not affected at either 1 or 5 μM/mL 5-azadC concentrations (Figure 3A). Hence, DNA-methylation does not appear to be required for the repression of LUC expression in WTLUC seedlings.

Figure 3
figure3

Treatments of WTLUC, hsi2-2LUCand hsi2-4LUCseedlings with DNA methylation inhibitor 5-aza-2′-deoxycytidine (5-azadC) and histone deacteylase inhibitor Trichostatin A (TSA). Seeds were germinated and grown vertically on media plates containing the indicated concentrations of either 5-azadC (A) or TSA (B). Luminescence imaging of 10 days old seedlings was performed using cooled CCD camera after spraying with the substrate luciferin. Pseudocolor images show luminescence intensity from lowest (blue) to highest (white). The experiment was repeated with two technical replicates and representative images are shown.

Histone deacetylation is known to regulate gene expression by transcriptional repression in eukaryotes [38]. Recently, it was reported that treatment of Arabidopsis seedlings with the histone deacetylase inhibitor trichostatin A (TSA) or down-regulation of two histone deacetylase genes, HDA6 and HDA19 by RNA interference resulted in derepression of the embryonic program in germinating seeds and seedlings [39]. Arabidopsis seedlings treated with TSA and histone deacetylase mutants mimic the phenotypes of hsi2-2/hsl1double mutant seedlings [12],[13], indicating that HSI2- and HSL1-mediated repression of the embryonic program could involve histone deacetylation. HSL1 was shown to physically interact with HDA19 via its CW domain and disruption of HSL1 resulted in increased H3K4me3 and decreased H3K27me3 marks on genes that encode transcriptional activators involved in the embryonic program [31]. To test the effects of TSA on the luminescence expression of WTLUC seedlings, WTLUC, hsi2-2LUC and hsi2-4LUC seedlings were grown on media containing 0.1 and 1 μg/mL TSA. Since higher concentrations of TSA resulted in severe growth retardation and developmental delay in all seedlings tested (Figure 3B), only, 0.1 and 1 μg/mL of TSA was used in these assays. Luminescence imaging data showed that treatments of WTLUC seedlings with TSA did not affect their luminescence expression (Figure 3B), indicating that HSI2 PHD-like domain mediated repression of LUC transgene expression in WTLUC seedlings is not dependent on TSA-sensitive histone deacetylation.

hsi2-2LUCand hsi2-4LUCmutant seedlings are partially resistant to DNA methylation inhibitor 5-azadC induced growth inhibition

We noticed that the growth and development of WTLUC seedlings on plates that contained 5 μM/mL 5-azadC was more strongly inhibited than hsi2-2LUC and hsi2-4LUC mutant seedlings (Figure 4A). To further characterize the effects of 5-azadC on hypocotyl and root growth, WTLUC, hsi2-2LUC and hsi2-4LUC seeds were germinated on media containing 0, 1, 5, 10 and 20 μM 5-azadC. After 7 days of incubation on 5-azadC-containing media, all seedlings showed dose-dependent inhibition of growth and development. However, the most severe effects were seen with WTLUC seedlings whereas the growth of hsi2-2LUC and hsi2-4LUC mutant seedlings was less inhibited (Figure 4A). While WTLUC seeds germinated when incubated on media containing 20 μM 5-azadC, subsequent root growth and cotyledon development was almost completely abrogated while both and hsi2-4LUC and hsi2-2LUC mutant seedlings continued to grow and develop, albeit slowly, under these conditions (Figure 4A and B). Comparative measurements of hypocotyl and root growth indicated that hsi2-2LUC and hsi2-4LUC mutant seedlings were about one half as sensitive to 5-azadC-dependent inhibition as WTLUC seedlings at 5, 10 and 20 μM 5-azadC treatments (Figure 4B). These data indicate that, although 5-azadC does not affect the HSI2-dependent suppression of luciferase expression in WTLUC plants, HSI2 does affect sensitivity to 5-azadC-dependent inhibition of seedling development.

Figure 4
figure4

Effects of DNA methylation inhibitor 5-aza-2′-deoxycytidine (5-azadC) on the growth and development of WTLUC, hsi2-2LUCand hsi2-4LUCseedlings. Seeds were germinated vertically on media plates containing various concentrations of 5-azadC and pictures were taken 7 days after germination. A. Morphology of seedlings. B. Measurements of hypocotyl and root growths. Hypocotyl and root lengths were measured using ImageJ software. Data represent mean values (’sD) from 10 seedlings. The experiments were repeated with two technical replicates. Letters indicate significant differences between WTLUC and hsi2-2LUC or hsi2-4LUC at each time point (a = p< 0.005, b = p< 0.0001).

LUC and NPTIItransgene expression is associated with changes in histone methylation marks

To examine the histone methylation properties along the transgene cassette and the role of HSI2 PHD-like domain in regulating those marks and transgene expression, ChIP-qPCR analyses were performed using 5 day old seedlings of various genotypes. Antibodies specific to H3K4me3, H3K9me2, H3K27me3 and H3K36me3 marks were used, along with PCR primers that specifically amplify sequences from the endogenous (native) and transgene GSTF8 promoters and LUC and NPTII coding regions. For the specific amplifications of E1 and E2 PCR fragments only from the endogenous GSTF8 promoter sequence during ChIP-qPCR, at least one PCR primer that binds outside of the −495 bp endogenous GSTF8 region that is not part of the GSTF8::LUC transgene cassette was used. Also, to make sure the PCR products of T1 and T2 fragments are only amplified from the GSTF8 transgene promoter sequence, at least one primer that binds outside of the −495 bp region in the GSTF8::LUC transgene cassette was used ([14], Figure 5A). PCR amplification specificities of E1, E2, T1 and T2 fragments were confirmed using Col-0 wild-type and WTLUC. Among the histone methylation marks, H3K4me3 and H3K36me3 are associated with actively transcribed genes, while H3K9me2 is a repressive mark commonly enriched on transposable elements and repetitive sequences [24]. H3K27me3 is a repressive mark associated with transcribed genes that are under tissue-specific or developmental regulation [40] [42]. Preimmune immunoglobulin G (IgG) was used as a negative control for non-specific binding and all genomic DNA fragments tested show very low background levels of enrichment when chromatin samples were immunoprecipitated with IgG (Figure 5B). FUS3 was used as a positive control for H3K27me3, while actin2/7 (ACT2/7) was used as a negative control for H3K27me3 and as a positive control for H3K4me3 and H3K36me3. TA2 was used as a positive control for H3K9me2 and as a negative control for H3K4me3 and H3K36me3. In agreement with our previous report [14], chromatin from the transgene GSTF8 promoter region, and both 5′and 3′ ends of the LUC coding sequences was highly enriched in H3K27me3 marks in WTLUC seedlings (Figure 5B) while endogenous GSTF8 promoter sequences showed consistently low levels of H3K27me3 (Figure 5B). Chromatin from the 5′region of the NPTII coding sequence was also highly H3K27me3 enriched in WTLUC seedlings (Figure 5B). A substantial decrease in H3K27me3 levels was detected on chromatin from the transgene GSTF8 promoter sequences and LUC and NPTII coding sequences in hsi2-4LUC seedlings that carry a point mutation in HSI2 PHD-like domain (Figure 5B). Though the transgene sequences tested showed considerable H3K9 dimethylation, unlike H3K27me3, no significant differences in H3K9me2 enrichment were seen between chromatin from WTLUC and hsi2-4LUC seedlings at any of the sites tested (Figure 5B). Therefore, among the histone methylation marks associated with transcriptional suppression, only H3K27me3 was dependent on the HSI2 PHD-like domain.

Figure 5
figure5

Chromatin immunoprecipitation (ChIP) and quantitative PCR (qPCR) analyses of H3K4me3, H3K9me2, H3K27me3 and H3K36me3 levels on endogenous GSTF8 promoter and transgene chromatin in WTLUCand hsi2-4LUCmutant. A. Genomic structures of endogenous GSTF8 gene and GSTF8::LUC transgene showing the locations of amplified regions by ChIP-qPCR. B. qPCR analyses of chromatin samples from 5-day old seedlings of WTLUC and hsi2-4LUC that were immunoprecipitated using either specific antibodies recognizing indicated histone methylation marks or IgG (non-specific binding control). Data is expressed as percentage of immunoprecipitated DNA relative to input DNA. ACT2/7 (H3K4me3 and H3K36me3), FUS3 (H3K27me3) and TA2 (H3K9me2) serve as positive controls whereas TA2 (H3K4me3 and H3K36me3) and ACT2/7 (H3K9me2 and H3K27me3) were used as negative controls. Data represent means (’sD) from two biological replicates with three qPCR replicates each. Significant differences in enrichment between WTLUC and hsi2-4LUC for each genomic region tested were determined using two-tailed Student’s t-test assuming unequal variances and P values are indicated by letters (a = p < 0.05, b = p < 0.005, c = p < 0.0005, d = p < 0.0001).

Histone methylation marks H3K4me3 and H3K36me3, which are associated with chromatin from actively transcribed genes, were enriched at all of the transgene sequences assayed in hsi2-4LUC seedlings, relative to WTLUC (Figure 5B). These marks were particularly abundant at the proximal transgene GSTF8 promoter and 5′LUC coding sequences but significant enrichment was also seen at the endogenous GSTF8 promoter and NPTII coding sequence.

To examine whether the hsi2-4-dependent changes in histone methylation marks are associated with both the KanR and KanS transgene loci, ChIP-qPCR analyses were performed on various regions of the endogenous GSTF8 gene and the GSTF8::LUC transgene in KanR, KanS, hsi2-4-KanR and hsi2-4-KanS seedlings (Figure 6). As in chromatin from WTLUC seedlings, higher levels of H3K27me3 marks at transgene GSTF8 promoter sequences and at LUC and NPTII coding sequences were detected in wild-type seedlings carrying either the KanR or KanS transgene locus than in corresponding hsi2-4-KanR or hsi2-4-KanS mutant seedlings. Thus, the significant decrease in H3K27me3 levels at the GSTF8::LUC transgene associated with homozygosity for the hsi2-4 allele was seen at both insertion sites. While chromatin from transgene sequences generally had higher levels of H3K9me2 marks than the endogenous GSTF8 gene, no significant change was seen between these genotypes.

Figure 6
figure6

ChIP-qPCR analyses of H3K4me3, H3K9me2, H3K27me3 and H3K36me3 enrichments on endogenous GSTF8 promoter and GSTF8::LUC transgene regions in WTLUCand hsi2-4 carrying either KanSor KanRlocus. Chromatin samples were prepared using five day old seedlings from various genotypes. ChIP-qPCR data represents mean values (’sD) of three PCR reactions obtained from each of two independent immunoprecipitations. Significant differences between wild-type lines and corresponding mutant lines (KanS versus hsi2-4-KanS, KanR versus hsi2-4-KanR) were determined using two-tailed Student’s t-test assuming unequal variances and P values are indicated by letters (a = p < 0.05, b = p < 0.005, c = p < 0.0005, d = p < 0.0001).

Enrichment of H3K4me3 and H3K36me3 was seen in chromatin at both KanR and KanS loci in hsi2-4 seedlings (Figure 6). This enrichment was most pronounced at LUC coding sequences rather than in promoter regions and significant enrichment was also seen in chromatin of the NPTII gene at the KanR locus. Therefore, disruption of HSI2 PHD-like domain resulted in increased activation marks on 5′and 3′ end of LUC coding sequences in both KanR and KanS backgrounds. However, increased H3K36me3 marks on the NPTII coding sequences were observed only in hsi2-4-KanR seedlings (Figure 6).

H3K27me3 levels are significantly decreased on a subset of seed-maturation genes in hsi2LUCmutant seedlings

Some members of the LAFL clade of regulatory genes that control the expression of seed maturation genes [3] are misregulated in hsi2 mutant seedlings [14]. LEC1 and ABI3 are ectopically expressed in hsi2-2 but not in hsi2-4 seedlings, while FUS3 is upregulated in both hsi2-2 and hsi2-4 lines [14]. These results suggested to us that the HSI2-dependent negative regulation of LEC1 and ABI3 in seedlings does not require the PHD-like domain, while suppression of FUS3 could be dependent on the PHD-like domain of HSI2. To determine if correlations exist between these expression patterns and histone modifications, ChIP-qPCR analysis of these genes was carried out on chromatin samples from WTLUC, hsi2-2LUC and hsi2-4LUC seedlings using antiH3K27me3 (Figure 7). Consistent with our hypothesis, significant reductions in H3K27me3 chromatin marks were detected in association with ABI3 and LEC1 genomic sequences only in chromatin from hsi2-2LUC but not hsi2-4LUC mutant seedlings and genes such as LEC2 and L1L, which are not misregulated in either hsi2 mutant allele alone, also showed no change in H3K27me3 marks in these mutant backgrounds. However, H3K27me3 marks associated with FUS3 sequences were not altered in either mutant background. Thus, the effects of HSI2 on FUS3 expression do not appear to depend on alterations in H3K27me3. On the other hand, the enrichment of H3K27me3 marks detected on both 5′and 3′ coding sequences of AGL15 (At5g13790) in WTLUC seedlings was significantly decreased in both hsi2-2LUC and hsi2-4LUC backgrounds, which is consistent with the increased expression of AGL15 in these mutants [14]. Therefore, the HSI2 PHD-like domain does appear to be required for both the repressed expression and H3K27 hypermethylation of the AGL15 locus.

Figure 7
figure7

ChIP-qPCR analyses of H3K27me3 levels on the master transcriptional regulators of seed-maturation program and seed maturation related target genes of HSI2 PHD-like domain in WTLUCand hsi2LUCmutant seedlings. Chromatin samples were prepared using five day old WTLUC, hsi2-2LUC and hsi2-4LUC seedlings. Immunoprecipitations were performed using either IgG (non-specific binding control) or anti-H3K27me3 antibody. Data represent means (’sD) obtained from three qPCR reactions each from three independent immunoprecipitations from three different biological replicates. Significant differences in H3K27me3 enrichments between WTLUC and hsi2LUC mutants (WTLUC versus hsi2-2LUC, WTLUC versus hsi2-4LUC) are indicated by letters (a = p < 0.05, b = p < 0.005, c = p < 0.0005, d = p < 0.0001). ABSCISIC ACID-INSENSITIVE3 (ABI3, At3g24650), FUSCA3 (FUS3, At3g26790), LEAFY COTYLEDON1 (LEC1, At1g21970), LEAFY COTYLEDON2 (LEC2, At1g28300), LEC1-LIKE (L1L, At5g47670), AGAMOUS-Like 15 (AGL15, At5g13790), At2g41260 (Late-embryogenesis-abundant protein), At3g22640 (cupin family seed storage protein), At1g04660 (glycine-rich protein), At2g29370 (NAD(P)-binding Rossmann-fold superfamily protein), At5g45830 (DOG1, DELAY OF GERMINATION1), At2g34700 (putative proline-rich glycoprotein).

As previously reported, a number of seed maturation-related structural genes are derepressed in hsi2-4 mutant seedlings [14]. To investigate whether the increased expression of these genes in hsi2 mutant seedlings correlates with changes in associated H3K27me3 marks, ChIP-qPCR analyses were performed on chromatin from WTLUC, hsi2-2LUC and hsi2-4LUC mutant seedlings using primers specific to At2g41260 (Late-embryogenesis-abundant (LEA) protein), At3g22640 (cupin family seed storage protein), At1g04660 (glycine-rich protein), At2g29370 (NAD(P)-binding Rossmann-fold superfamily protein), At5g45830 (Delay of Germination1, DOG1) and At2g34700 (proline-rich glycoprotein). These seed maturation genes were previously reported to be targets of H3K27 hypermethylation [41] and our results confirmed that H3K27me3 marks were enriched, relative to IgG, on the 5′coding sequences of these genes (Figure 7). In comparison to WTLUC, significant decreases in H3K27me3 enrichments were detected in both hsi2-2LUC and hsi2-4LUC mutant seedlings on chromatin associated with 5′genomic sequences of DOG1 and At2g34700, and at both 5′and 3′ sequences of At3g22640 (Figure 7). On the other hand, H3K27me3 levels on At1g04660 and At2g29370 were significantly decreased only in the hsi2-2LUC background and not in chromatin from hsi2-4LUC seedlings (Figure 7). These data indicate that, as with the regulatory genes described above, derepression of seed maturation-specific gene expression in hsi2 mutant seedlings often corresponds with decreased accumulation of H3K27me3 marks that may or may not depend on the presence of an intact HSI2 PHD-like domain. However, as with FUS3, the LEA-like protein gene At2g41260, which is expressed at elevated levels in both hsi2-2 and hsi2-4 mutant plants [14], is strongly enriched for H3K27me3 but these marks are not significantly reduced in either hsi2 mutant.

Discussion

Despite detailed genetic and functional characterization, the molecular mechanisms that underlie HSI2- and HSL1-mediated repression of seed maturation program in seedlings are still not fully understood [12] [16]. HSI2 contains a PHD-like domain [3],[12],[14],[26] and PHD domains can act as “readers” of the histone methylation status of target genes to regulate their expression [30]. Through characterization of a novel mutant allele, hsi2-4, which affects the expression of GSTF8::LUC transgenes and certain seed maturation genes [14], we provide evidence that the HSI2 PHD-like domain is involved in regulating the expression of some genes by altering histone modifications.

Quantitative PCR data indicates that both KanR and KanS transgene loci in the WTLUC reporter gene line contain multiple copies of the GSTF8::LUC transgene. The KanR locus is more complex than KanS, harboring five copies of the transgene, while the KanS locus includes two (Table 1). Transgene loci in plants that harbor multiple and complex transgene repeats at a single locus were frequently targeted by DNA methylation-associated H3K9me2 marks and also histone deacetylation mediated transcriptional gene silencing [35],[43] [45]. However, treatment of WTLUC seedlings with either the DNA methylation inhibitor 5-azadC or the histone deacetylase inhibitor TSA failed to derepress the LUC expression in WTLUC seedlings (Figure 3A and B). Also, ChIP-qPCR analyses showed no differences in DNA methylation-associated H3K9me2 histone methylation marks on transgene sequences between WTLUC and hsi2-4 mutants that harbor either individual KanR and KanS loci or both (Figures 5B and 6). Based on these data, it appears that DNA methylation and histone deacetylation mechanisms are not involved in HSI2 PHD-like domain-mediated repression of transgene expression in WTLUC seedlings.

Although 5-azadC did not affect LUC expression in WTLUC seedlings (Figure 3A), hsi2 mutant seedlings treated with various concentrations of 5-azadC maintained root and hypocotyl growth better than WTLUC seedlings under the same conditions (Figure 4A and B). These results could indicate that HSI2 is somehow involved in the inhibition of seedling growth and development caused by 5-azadC and, since this effect was observed with both the PHD-like domain mutant allele hsi2-4 and the hsi2-2 T-DNA knock-out allele, it appears that the PHD-like domain may be required for this 5-azadC-dependent inhibition of growth. LEC1, an HSI2 and HSL1 target gene and member of the “LAFL network”, was shown to be regulated by DNA methylation [46],[47] and the embryonic phenotypes of gain-of-function lec1 mutants were enhanced by treatment with a DNA methylation inhibitor 5-azacytidine [46]. Hence, the partial rescue of seedling growth in hsi2 mutants in the presence of 5-azadC could be an indirect effect of changes in the DNA methylation status of HSI2-targeted regulatory genes, including LEC1.

Data presented here clearly indicate that the HSI2 PHD-like domain is involved in suppressing the expression of GSTF8::LUC transgenes in both KanR and KanS transgene loci (Figure 2). Furthermore, the levels of LUC expression seen in hsi2-4 plants that carry these reporter complexes correlate with transgene copy number, with relatively low levels of LUC expression seen in hsi2-4, KanS plants that contain two GSTF8::LUC copies and correspondingly higher levels expression in hsi2-4, KanR or hsi2-4LUC lines with five and seven total reporter gene copies, respectively. The direct correlation between derepressed LUC expression and the number of reporter gene copies means that the luminescence of mutant plants with compromised HSI2-dependent repression will be amplified in a high copy number reporter gene line, making these mutants far more apparent in a luminescence-based mutant screen. On the other hand, LUC expression in wild-type plants that carry the KanR locus alone showed higher levels of LUC transcripts than did WTLUC seedlings (Figure 2). These results appear to indicate that the presence of both KanR and KanS loci in the same genome may lead to stronger transcriptional suppression of the GSTF8::LUC reporter genes than when the KanR locus is present alone. No NPTII transcripts were detected in either KanS or hsi2-4-KanS seedlings (Figure 2), which is consistent with the kanamycin sensitivity of these plants. However, hsi2-4-KanR seedlings showed 3-fold higher expression of NPTII transcripts relative to wild-type KanR seedlings but this was not seen in corresponding hsi2-4LUC and WTLUC seedlings (Figure 2). Thus, as with the LUC reporter gene, the co-existence of KanR and KanS loci is associated with stronger transcriptional repression of NPTII gene expression. To better understand the function of the HSI2 PHD-like domain, interactions between activation-associated and repressive histone methylation marks at the KanR and KanS transgene loci were evaluated by ChIP-qPCR assays. Transgene GSTF8 promoter sequences, along with LUC and NPTII coding sequences were highly enriched in H3K27me3 marks in WTLUC seedlings (Figures 5B and 6) and significantly lower amounts of H3K27me3 were observed on these transgene sequences in the hsi2-4 mutant background (Figure 5B). Similar histone modification patterns were observed in seedlings harboring individual KanR or KanS loci (Figure 6). Thus, the PHD-like domain of HSI2, which is required to repress the expression of these transgene complexes, is also necessary for the appearance of H3K27me3 marks on these loci. In contrast, H3K4me3 and H3K36me3 histone methylation marks, which are associated with active gene expression and have been shown to inhibit H3K27me3 marks on transcribed genes in both animals and plants [48] [50], were enriched on these transgene sequences in hsi2-4 seedlings, relative to those with the wild-type HSI2 allele. Thus, the decrease in H3K27me3 marks on transgene sequences in both KanR and KanS loci in the hsi2-4 mutant background is associated with both increased expression (Figure 2) and increased accumulation of H3K4me3 and H3K36me3 marks (Figures 5B and 6). Developmental repression of transcribed genes is often associated with H3K27me3 marks [41],[42] but emerging evidence also suggests that H3K27me3 may act as an alternative to DNA methylation-associated H3K9me2 in transposable elements and repetitive sequence silencing [40],[51] [54]. Turck et al. [40] showed that the chromodomain-containing H3K27me3 “reader” protein LHP1 (LIKE HETEROCHROMATIN PROTEIN 1) is enriched on tandemly duplicated genes, such as the nine closely linked chitinase/glucosylase-18 genes (At4g19720-At4g19820) on chromosome 4 of Arabidopsis, but not on segmentally duplicated genes. Expressed genes that flank tandemly duplicated gene loci are also not associated with LHP1. Both KanR and KanS loci contain multiple GSTF8::LUC transgenes at individual loci (Table 1) and expression of genes that flank the KanR transgene locus does not differ between WTLUC and hsi2-4LUC seedlings [14]. However, similarities in HSI2-dependent H3K27me3 accumulation at the KanR and KanS loci and the corresponding relative changes in LUC expression in hsi2-4 plants lead us to speculate that transcriptional repression is mediated by the GSTF8::LUC transgene itself and is not dependent on tandem T-DNA insertions. The presence of multiple transgene copies results in high levels of expression that accentuates the apparent repressive effect of HSI2 when its activity is compromised by mutation. However, since LUC expression in WTLUC seedlings is lower than in KanR seedlings, the presence of these two unlinked loci appears to have synergistic effects on reporter gene silencing.

In contrast to the GSTF8::LUC reporter genes, expression of native GSTF8 transcripts derived from either the “long” or “short” transcriptional start sites show no significant increase in the hsi2-4LUC mutant background (Figure 2). The most parsimonious explanation for the discrepancy between native GSTF8 expression and the expression of the GSTF8::LUC reporter gene is that the isolated GSTF8 promoter sequence used in the GSTF8::LUC gene construct, which corresponds with the short promoter as defined by Thatcher et al. [34], could contain cis-acting suppressor elements that are masked in the context of the native gene. Support for this explanation can be seen in Figure 4B of Veerappan et al. [14]. In this experiment, LUC expression (measured as luminescence) was assayed in WT Col-0 and mutant hsi2-4 Arabidopsis plants newly transformed with either short-GSTF8::LUC or long-GSTF8::LUC gene constructs (Col-S, Col-L and hsi2-4-S and hsi2-4-L, respectively). No significant differences in luminescence were apparent between Col-L and hsi2-4-L plants but LUC expression driven by the short GSTF8 promoter in hsi2-4-S plants was substantially elevated relative to that detected in Col-S plants.

Polycomb group (PcG) proteins are evolutionarily conserved multi-protein complexes required for developmental repression of gene expression by chromatin based mechanisms. PcG proteins in plants comprise of two major complexes: Polycomb Repressive Complex 1 (PRC1) and PRC2 [4] [6],[55]. Arabidopsis PRC1 proteins BMI1 and RING1 were shown to have histone H2A mono ubiquitination (H2Aub) activity in vitro[16],[56],[57], whereas PRC2 complex proteins catalyze the deposition of H3K27me3 marks to promote developmental repression in animals and plants [58] [60].

The GSTF8 promoter sequence used in the GSTF8::LUC reporter contains an octopine synthase (OCS) sequence element at −460 (Additional file 1) that is known to be required for transcriptional activation in response to a variety of biotic and abiotic stress signals [61]. This OCS element is flanked by OCS element binding factor 5 (OBF5)and OCS element binding proteins 1 (OBP1) elements that were shown to bind proteins of the DNA binding with One Finger (DOF) family of transcription factors and are reported to act either as positive or negative regulatory factors in various plant genes [62]. A putative myeloblastosis2 (MYB2) binding element is also located at −311 but its potential function is unknown. Sequence elements with a potential role in PRC2-based silencing can also be identified in the GSTF8 promoter. By analyzing the co-distribution of the ubiquitous Arabidopsis PRC2 protein FERTILIZATION INDEPENDENT ENDOSPERM (FIE) and H3K27me3 marks, Deng et al. [63] identified four sequence motifs that could act as PRC2 binding sites. Two of these motifs are found in the GSTF8 promoter, a putative GAGA box located between −123 and −131 relative to the transcription start site and a TTC repeat element located in the 5′untranslated region between +36 and +56. GAGA elements were found to be specifically associated with FIE and H3K27me3 enriched sites and these sequences are also found in polycomb response elements (PREs) of Drosophila [63]. However, while GAA (reverse complement of TTC) sequence motifs were found to be associated with genomic regions that bind both FIE and H3K27me3, they were also enriched in random promoter sequences and were, therefore, not considered to be specific PRC2 binding sites [63].

The role of the GAGA element in the transcriptional regulation of the PRC2-repressed LEC2 gene was confirmed by Berger et al. [64]. However, mutational analysis showed that, in this context, the GAGA element acted as a required cis-activating element, which was associated with a distinct cis-repressing element, termed repressive LEC2 element (RLE) that apparently consists of two component sequences. Comparison of the LEC2 RLE sequence with the GSTF8 promoter identified a duplicated element identical to the 5′component of the LEC2 RLE. As in the LEC2 gene, the putative GSTF8 RLE-like sequence is located immediately downstream of the GAGA element (Additional file 1). Whether this putative GAGA-RLE motif plays a role in the transcriptional regulation of the GSTF8 promoter in a transgene context is not known.

ChIP-qPCR analyses were carried out to investigate whether LAFL network genes and other seed maturation genes that are up-regulated in hsi2 mutant seedlings are also associated with HSI2-dependent changes in H3K27me3 marks. With the exception of L1L, all the tested seed maturation genes are enriched with H3K27me3 marks in WTLUC seedlings (Figures 7) in agreement with previous reports [41],[58]. Significant reductions in H3K27me3 levels, relative to WTLUC seedlings, were observed at some gene loci in both hsi2-2LUC and hsi2-4LUC seedlings, while other genes showed reductions only in hsi2-2LUC seedlings and a few showed no changes in either hsi2 mutant background. In general, these HSI2-dependent differences in H3K27me3 marks correlate well with HSI2-dependent changes in gene expression. For example, among the regulatory genes tested, ABI3 and LEC1, which are expressed at elevated levels in hsi2-2 seedlings but not in hsi2-4 seedlings [14], showed correspondingly decreased levels of H3K32me3 enrichment in chromatin from hsi2-2LUC plants but not from hsi2-4LUC plants (Figure 7). On the other hand, H3K27me3 marks on AGL15 were strongly decreased in both hsi2-2LUC and hsi2-4LUC mutant plants (Figure 7) and expression of this gene is also upregulated in both hsi2 mutant backgrounds. These results can be interpreted to indicate that HSI2-dependent transcriptional repression and H3K27 hypermethylation of AGL15 is mediated by the PHD-like domain while that of ABI3 and LEC1 may be mediated by other HSI2 domains, such as the CW domain, which was shown to interact with HDA19 to promote histone deacetylation and H3K27me3 marks to repress seed maturation genes [31].

Similar patterns can be seen in the structural (non-regulatory) seed maturation gene sample. The putative glycine-rich protein gene (At1g04660) is more strongly expressed in hsi2-2LUC plants than in hsi2-4LUC plants [14] and H3K27me3 marks at this site are correspondingly reduced in hsi2-2LUC but not hsi2-4LUC mutants. On the other hand, genes for a cupin-like protein (At3g22640), DOG1 (At5g45830) and a proline-rich glycoprotein (At2g34700) are similarly up-regulated in hsi2-2LUC and hsi2-4LUC plants and decreased H3K27me3 marks are also apparent in both mutant lines.

Two of the genes in our sample group do not show correlations between HSI2-dependent transcriptional repression and H3K27me3 marks. These genes, FUS3 and the LEA protein gene At2g41260, are expressed at elevated levels in both hsi2-2LUC and hsi2-4LUC backgrounds and accumulate substantial H3K27me3 in WTLUC plants but these marks are not diminished in either hsi2 mutant. Therefore, derepression of these genes in hsi2 mutant plants does not appear to require the depletion of H3K27me3. It seems likely that, in addition to the PRC2-mediated accumulation of H3K27me3 marks, other mechanisms are involved in the repression of these genes.

AGL15 is a member of the MIKC subfamily of MADS domain transcription factors that is preferentially expressed in developing embryos. Ectopic overexpression of AGL15 results in enhanced somatic embryogenesis [65],[66]. AGL15 acts upstream of LAFL network genes and several LAFL genes are direct regulatory targets of AGL15 [67],[68]. DOG1 is a seed-specific gene and plays a critical role in promoting seed dormancy by integrating environmental signals [69],[70]. Recently, AGL15 and DOG1 were shown to be targets of H3K27me3 marks, which could be mediated by PRC1 proteins [58],[71]. Our data shows that HSI2 regulates the expression of AGL15 and DOG1 in seedlings by promoting H3K27me3 marks possibly via PRC1-PRC2 complex which requires HSI2 PHD-like domain.

HSI2 was shown to directly interact with PRC1 complex proteins AtBMI1A/B/C and is required for the deposition of H2Aubi and H3K27me3 marks on “LAFL network” seed-maturation genes including LEC1, FUS3 and ABI3[16]. Disruption of PRC2 complex genes in Arabidopsis led to decreased H3K27me3 levels, activation of “LAFL network” transcription factor genes and ectopic expression of embryonic traits during seed germination and vegetative development [58],[72],[73]. Thus, it is possible that HSI2 interacts with PRC1 proteins like AtBMI1 to recruit PRC2 proteins such as the histone methyltransferase CURLY LEAF (CLF), to deposit H3K27me3 marks on the GSTF8::LUC transgene loci. AtRING1a was also shown to physically interact with the PRC2 core component CLF [74] and several reports have demonstrated the involvement of PHD-PRC2 and PHD-PRC1-PRC2 complexes in deposition of H3K27me3 marks to promote transcriptional repression of gene expression in plants [16],[71],[75]. Similarly, HSI2 could be part of a repression complex that involves the HSI2 PHD-like domain, PRC1 and PRC2 complex proteins to promote high levels of H3K27me3 marks on native seed maturation genes and GSTF8::LUC transgene loci to repress their expression during the seed to seedling developmental phase transition.

Conclusions

HSI2 contains a putative PHD domain, which could act as a “reader” of histone methylation marks. In this work, we show that HSI2 PHD-like domain regulates both LUC and NPTII transgenes from two independent transgene loci. Transcriptional repression of both of these transgene loci by HSI2 PHD-like domain is associated with repressive histone methylation marks H3K27me3 but not siRNA and DNA methylation associated H3K9me2 marks. In addition to the transgenes, HSI2 is also required for the repression of a subset of seed maturation genes in seedlings by promoting H3K27me3 marks in a PHD-like domain dependent and independent manner.

Methods

Plant materials, growth conditions and chemical treatments

Arabidopsis thaliana Columbia-0 (Col-0; CS60000) wild-type was obtained from Arabidopsis Biological Resources Center. WTLUC and hsi2-4LUC, which contain both KanR and KanS transgene loci in Col-0 background were described before [14]. The other genotypes that were used in this study including KanR, KanS, hsi2-4-KanR and hsi2-4-KanS were obtained by crossing the WTLUC and hsi2-4LUC into the Col-0 wild-type. The hsi2-2LUC line harbors GSTF8::LUC transgenes in the HSI2 T-DNA knock-out allele hsi2-2 (SALK_088606) background. For all the experiments described here, plants were grown under continuous illumination at 24°C on 0.3% Phytagel plates containing 0.5X Murashige and Skoog (MS) salt, 0.5 g/L MES (2-(N-morpholino) ethanesulfonic acid, 1X Gamborg vitamin mix and 1 % sucrose (pH adjusted to 5.7). 5-Aza-2′-deoxycytidine (5-azadC, A3656; Sigma) and Trichostatin A (T8552; Sigma) stocks were prepared using dimethyl sulfoxide and methanol respectively, and added directly to the MS media plates. Hypocotyl and root length measurements were made using Image J software (http://rsbweb.nih.gov/ij/). Digital photos of seedlings grown on 5-azadC plates were taken along with a ruler of known length. Hypocotyl length was measured from the tip of the apical meristem to the junction between hypocotyl and root, while root length was measured from the hypocotyl/root junction to the tip of the primary root.

Luminescence imaging, genetic crosses and genotyping

Luminescence imaging was performed using Andor iKON-M DU934N-BV CCD camera (Andor Technology). After spraying with 1 mM D-luciferin potassium salt (Gold Biotechnology) containing 0.01% Triton X-100 solution, seedlings were kept in the dark for 5 minutes and imaging was performed with a 5 minute exposure. Andor SOLIS (I) imaging software (Andor Technology) was used for the acquisition of luminescence images and processing. To separate KanR and KanS loci, WTLUC and hsi2-4LUC mutant were crossed into Col-0 wild-type plants. Successful crosses were identified based on luciferase imaging in the F1 generation and plants were allowed to self-pollinate. Progeny lines homozygous for the KanS locus were identified based on kanamycin sensitivity whereas plants homozygous for the KanR locus were identified by PCR genotyping using T-DNA and genomic primers. To genotype hsi2-4 mutation, a previously described CAPS marker [14] was used.

Preparation of total RNA, cDNA synthesis and real-time reverse transcription quantitative PCR

Total RNA extraction and real-time reverse transcription quantitative-PCR (RT-qPCR) analysis was performed as described in Veerappan et al. [14]. Primers used in RT-qPCR are listed in Additional file 2.

Estimation of LUCtransgene copy numbers by real-time quantitative PCR

To determine the copy numbers of LUC transgenes in WTLUC, KanR and KanS lines, real-time quantitative PCR (qPCR) was performed as described before [32],[33]. All qPCR reactions were performed using AB StepOnePlus Real-Time PCR System (Applied Biosystems) in 10 μl volume containing different amounts of DNA, 0.2 μM of each primers and 5 μl iTaqTM SYBR green supermix (Bio-Rad). Several sets of primers were tested for optimal performance. Temperature cycling conditions were 95°C for 10 minutes, 40 cycles for 15 seconds at 95°C and 1 minute at 60°C. Each DNA sample was tested in triplicates with three different DNA concentrations. Calibration curves were also performed in triplicates with five different DNA concentrations. Ct values were calculated using StepOne Software v2.1 (Applied Biosysytems). Concentrations of DNA samples were measured using Nanodrop 2000 (Thermoscientific) and the exact copy numbers of the template genome in the reactions were calculated using the following website: http://cels.uri.edu/gsc/cndna.html), applying the formula: number of copies = (amount * 6.022×1023)/(length * 1×109 * 650). The calibration curves for LUC were created using the plasmid DNA pBI121-GSTF8::LUC as a template. At5g47480, a single copy gene from Arabidopsis, was used as an internal control for normalization of the data. PCR primers used in the estimation of transgene copy numbers can be found in Additional file 2.

Chromatin immunoprecipitation and quantitative PCR analyses

Chromatin immunoprecipitation (ChIP) and quantitative PCR (qPCR) analyses were performed as described by Veerappan et al. [14]. Percentage of immunoprecipitated DNA relative to the total chromatin input was calculated for various samples using qPCR. Antibodies used for ChIP: normal rabbit IgG (Millipore, 12–370), anti-H3K4me3 (Millipore, 07–473), anti-H3K9me2 (Abcam, ab1220), anti-H3K27me3 (Millipore, 07–449) and anti-H3K36me3 (Abcam, ab9050). Primers used for ChIP PCR analyses are listed in Additional file 2.

Authors’ contributions

VV designed and performed genetic crosses, luminescence imaging, qRT-PCR, DNA methylation and histone deacetylase inhibitor experiments, chromatin immunopreciptation and qPCR analyses. NC and AR designed and performed the estimation of LUC copy numbers in various GSTF8::LUC transgene reporter lines. VV and RDA coordinated all the experiments, wrote and edit the manuscript. All authors read and approved the final manuscript.

Additional files

Abbreviations

HSI2:

High-level expression of sugar inducible gene2

HSL1:

HSI2-Like1

VAL1:

Viviparous1/ABI3-Like1

PHD:

Plant homeodomain

LUC:

Luciferase

5-azadC:

5-aza-2′-deoxycytidine

GSTF8:

:LUC: Glutathione S-transferase F8::luciferase

NPTII:

Neomycin phosphotransferase II

LEC1:

Leafy cotyledon1

LEC2:

Leafy cotyledon2

L1L:

LEC1-Like

ABI3:

Abscisic acid insensitive3

FUS3:

FUSCA3

LAFL:

AFL (ABI3/ FUS3/ LEC2)/LEC (LEC1/L1L)

GUS:

β-glucuronidase

CW:

Cysteine and tryptophan residue-containing

EAR:

Ethylene-responsive element binding factor-associated amphiphilic repression

siRNAs:

Small interfering RNAs

ChIP:

Chromatin immunoprecipitation

qPCR:

Quantitative PCR

WTLUC:

Wild-type luciferase

KanR:

Kanamycin resistant

KanS:

Kanamycin sensitive

AGL15:

AGAMOUS-Like15

GSTF8-L:

GSTF8-long

GSTF8-S:

GSTF8-short

TSA:

Trichostatin A

ACT2/7:

Actin2/7

IgG:

Immunoglobulin G

H3K27me3:

Trimethylation of histone H3 at lysine 27

H3K4me3:

Trimethylation of histone H3 at lysine 4

H3K36me3:

Trimethylation of histone H3 at lysine 36

H3K9me2:

Dimethylation of histone H3 at lysine 9

LEA:

Late embryogenesis abundant

DOG1:

Delay of germination1

LHP1:

Like heterochromatin protein1

PcG:

Polycomb group

PRC1:

Polycomb repressive complex 2

PRC2:

Polycomb repressive complex 1

H2Aub:

H2A mono ubiquitination

OCS:

Octopine synthase

OBF5:

OCS element binding factor 5

OBP1:

OCS element binding proteins 1

DOF:

DNA binding with one finger

FIE:

Fertilization independent endosperm

MYB2:

Myeloblastosis2

PREs:

Polycomb response elements

RLE:

Repressive LEC2 element

CLF:

Curly leaf

References

  1. 1.

    Zhang H, Ogas J: An epigenetic perspective on developmental regulation of seed genes. Mol Plant. 2009, 2: 470-627. 10.1093/mp/ssp027.

  2. 2.

    Muller K, Bouyer D, Schnittger A, Kermode AR: Evolutionarily conserved histone methylation dynamics during seed life-cycle transitions. PLoS One. 2012, 7: e51532-10.1371/journal.pone.0051532.

  3. 3.

    Jia H, Suzuki M, McCarty DR: Regulation of the seed to seedling developmental phase transition by the LAFL and VAL transcription factor networks. Wires Dev Biol. 2014, 3: 135-145. 10.1002/wdev.126.

  4. 4.

    Calonje M: PRC1 marks the difference in plant PcG repression. Mol Plant. 2014, 7: 459-471. 10.1093/mp/sst150.

  5. 5.

    Derkacheva M, Hennig L: Variations on a theme: polycomb group proteins in plants. J Exp Bot. 2014, 65: 2769-2784. 10.1093/jxb/ert410.

  6. 6.

    Molitor A, Shen W-H: The polycomb complex PRC1: composition and function in plants. J Genet Genomics. 2013, 40: 231-238. 10.1016/j.jgg.2012.12.005.

  7. 7.

    Lotan T, Ohto M, Yee KM, West MA, Lo R, Kwong RW, Yamagishi K, Fischer RL, Goldberg RB, Harada JJ: Arabidopsis LEAFY COTYLEDON1 is sufficient to induce embryo development in vegetative cells. Cell. 1998, 93: 1195-1205. 10.1016/S0092-8674(00)81463-4.

  8. 8.

    Kwong RW, Bui AQ, Lee H, Kwong LW, Fischer RL, Goldberg RB, Harada JJ: LEAFY COTYLEDON1-LIKE defines a class of regulators essential for embryo development. Plant Cell. 2003, 15: 5-18. 10.1105/tpc.006973.

  9. 9.

    Giraudat J, Hauge BM, Valon C, Smalle J, Parcy F, Goodman HM: Isolation of the Arabidopsis ABI3 gene by positional cloning. Plant Cell. 1992, 4: 1251-1261. 10.1105/tpc.4.10.1251.

  10. 10.

    Luerssen H, Kirik V, Herrmann P, Misera S: FUSCA3 encodes a protein with a conserved VP1/AB13-like B3 domain which is of functional importance for the regulation of seed maturation in Arabidopsis thaliana. Plant J. 1998, 15: 755-764. 10.1046/j.1365-313X.1998.00259.x.

  11. 11.

    Stone SL, Kwong LW, Yee KM, Pelletier J, Lepiniec L, Fischer RL, Goldberg RB, Harada JJ: LEAFY COTYLEDON2 encodes a B3 domain transcription factor that induces embryo development. Proc Natl Acad Sci U S A. 2001, 98: 11806-11811. 10.1073/pnas.201413498.

  12. 12.

    Suzuki M, Wang HHY, McCarty DR: Repression of the LEAFY COTYLEDON 1/B3 regulatory network in plant embryo development by VP1/ABSCISIC ACID INSENSITIVE 3-LIKE B3 genes. Plant Physiol. 2007, 143: 902-911. 10.1104/pp.106.092320.

  13. 13.

    Tsukagoshi H, Morikami A, Nakamura K: Two B3 domain transcriptional repressors prevent sugar-inducible expression of seed maturation genes in Arabidopsis seedlings. Proc Natl Acad Sci U S A. 2007, 104: 2543-2547. 10.1073/pnas.0607940104.

  14. 14.

    Veerappan V, Wang J, Kang M, Lee J, Tang Y, Jha AK, Shi H, Palanivelu R, Allen RD: A novel HSI2 mutation in Arabidopsis affects the PHD-like domain and leads to derepression of seed-specific gene expression. Planta. 2012, 236: 1-17. 10.1007/s00425-012-1630-1.

  15. 15.

    Jia H, McCarty DR, Suzuki M: Distinct roles of LAFL network genes in promoting the embryonic seedling fate in the absence of VAL repression. Plant Physiol. 2013, 163: 1293-1305. 10.1104/pp.113.220988.

  16. 16.

    Yang C, Bratzel F, Hohmann N, Koch M, Turck F, Calonje M: VAL- and AtBMI1-mediated H2Aub initiate the switch from embryonic to postgerminative growth in Arabidopsis. Curr Biol. 2013, 23: 1324-1329. 10.1016/j.cub.2013.05.050.

  17. 17.

    Tsukagoshi H, Saijo T, Shibata D, Morikami A, Nakamura K: Analysis of a sugar response mutant of Arabidopsis identified a novel B3 domain protein that functions as an active transcriptional repressor. Plant Physiol. 2005, 138: 675-685. 10.1104/pp.104.057752.

  18. 18.

    Tang X, Hou A, Babu M, Nguyen V, Hurtado L, Lu Q, Reyes JC, Wang A, Keller WA, Harada JJ, Tsang EW, Cui Y: The Arabidopsis BRAHMA chromatin-remodeling ATPase is involved in repression of seed maturation genes in leaves. Plant Physiol. 2008, 147: 1143-1157. 10.1104/pp.108.121996.

  19. 19.

    Sharma N, Bender Y, Boyle K, Fobert PR: High-level expression of sugar inducible gene2 (HSI2) is a negative regulator of drought stress tolerance in Arabidopsis. BMC Plant Biol. 2013, 13: 170-10.1186/1471-2229-13-170.

  20. 20.

    Ahmad A, Zhang Y, Cao XF: Decoding the epigenetic language of plant development. Mol Plant. 2010, 3: 719-728. 10.1093/mp/ssq026.

  21. 21.

    Feng S, Jacobsen SE, Reik W: Epigenetic reprogramming in plant and animal development. Science. 2010, 330: 622-627. 10.1126/science.1190614.

  22. 22.

    Law JA, Jacobsen SE: Establishing, maintaining and modifying DNA methylation patterns in plants and animals. Nat Rev Genet. 2010, 11: 204-220. 10.1038/nrg2719.

  23. 23.

    Liu C, Lu F, Cui X, Cao X: Histone methylation in higher plants. Annu Rev Plant Biol. 2010, 47: 395-420. 10.1146/annurev.arplant.043008.091939.

  24. 24.

    He G, Elling AA, Deng XW: The epigenome and plant development. Annu Rev Plant Biol. 2011, 62: 411-435. 10.1146/annurev-arplant-042110-103806.

  25. 25.

    Perry J, Zhao YD: The CW domain, a structural module shared amongst vertebrates, vertebrate-infecting parasites and higher plants. Trends Biochem Sci. 2003, 28: 576-580. 10.1016/j.tibs.2003.09.007.

  26. 26.

    Suzuki M, McCarty DR: Functional symmetry of the B3 network controlling seed development. Curr Opin Plant Biol. 2008, 11: 548-553. 10.1016/j.pbi.2008.06.015.

  27. 27.

    Lee WY, Lee D, Chung WI, Kwon CS: Arabidopsis ING and Alfin1-like protein families localize to the nucleus and bind to H3K4me3/2 via plant homeodomain fingers. Plant J. 2009, 58: 511-524. 10.1111/j.1365-313X.2009.03795.x.

  28. 28.

    He F, Umehara T, Saito K, Harada T, Watanabe S, Yabuki T, Kigawa T, Takahashi M, Kuwasako K, Tsuda K: Structural insight into the zinc finger CW domain as a histone modification reader. Structure. 2010, 18: 1127-1139. 10.1016/j.str.2010.06.012.

  29. 29.

    Hoppmann V, Thorstensen T, Kristiansen PE, Veiseth SV, Rahman MA, Finne K, Aalen RB, Aasland R: The CW domain, a new histone recognition module in chromatin proteins. EMBO J. 2011, 30: 1939-1952. 10.1038/emboj.2011.108.

  30. 30.

    Sanchez R, Zhou MM: The PHD finger: a versatile epigenome reader. Trends Biochem Sci. 2011, 36: 364-372.

  31. 31.

    Zhou Y, Tan B, Luo M, Li Y, Liu C, Chen C, Yu CW, Yang S, Dong S, Ruan J, Yuan L, Zhang Z, Zhao L, Li C, Chen H, Cui Y, Wu K, Huang S: HISTONE DEACETYLASE19 interacts with HSL1 and participates in the repression of seed maturation genes in Arabidopsis seedlings. Plant Cell. 2013, 25: 134-148. 10.1105/tpc.112.096313.

  32. 32.

    Mason G, Provero P, Varia AM, Acotto GP: Estimating the number of integrations in transformed plants by quantitative real-time PCR. BMC Biotechnol. 2003, 2: 20-10.1186/1472-6750-2-20.

  33. 33.

    Gadaleta A, Giancaspro A, Cardone MF, Blanco A: Real-time PCR for the detection of precise transgene copy number in durum wheat. Cell Mol Biol Lett. 2011, 16: 652-668. 10.2478/s11658-011-0029-5.

  34. 34.

    Thatcher LF, Carrie C, Andersson CR, Sivasithamparam K, Whelan J, Singh KB: Differential gene expression and subcellular targeting of Arabidopsis glutathione S-transferase F8 is achieved through alternative transcription start sites. J Biol Chem. 2007, 282: 28915-28928. 10.1074/jbc.M702207200.

  35. 35.

    Murfett J, Wang XJ, Hagen G, Guilfoyle TJ: Identification of Arabidopsis histone deacetylase HDA6 mutants that affect transgene expression. Plant Cell. 2001, 13: 1047-1047-10.1105/tpc.13.5.1047.

  36. 36.

    Chang S, Pikaard CS: Transcript profiling in Arabidopsis reveals complex responses to global inhibition of DNA methylation and histone deacetylation. J Biol Chem. 2005, 280: 796-804. 10.1074/jbc.M409053200.

  37. 37.

    Dinh TT, O’Leary M, Won SY, Li S, Arroyo L, Liu X, Defries A, Zheng B, Cutler SR, Chen X: Generation of a luciferase-based reporter for CHH and CG DNA methylation in Arabidopsis thaliana. Silence. 2013, 4: 1-10.1186/1758-907X-4-1.

  38. 38.

    Hollender C, Liu Z: Histone deacetylase genes in Arabidopsis development. J Integrative Plant Biol. 2008, 50: 875-885. 10.1111/j.1744-7909.2008.00704.x.

  39. 39.

    Tanaka M, Kikuchi A, Kamada H: The Arabidopsis histone deacetylases HDA6 and HDA19 contribute to the repression of embryonic properties after germination. Plant Physiol. 2008, 146: 149-161. 10.1104/pp.107.111674.

  40. 40.

    Turck F, Roudier F, Farrona S, Martin-Magniette ML, Guillaume E, Buisine N, Gagnot S, Martienssen RA, Coupland G, Colot V: Arabidopsis TFL2/LHP1 specifically associates with genes marked by trimethylation of histone H3 lysine 27. PLoS Genet. 2007, 3: e86-10.1371/journal.pgen.0030086.

  41. 41.

    Zhang X, Clarenz O, Cokus S, Bernatavichute YV, Pellegrini M, Goodrich J, Jacobsen SE: Whole-genome analysis of histone H3 lysine 27 trimethylation in Arabidopsis. PLoS Biol. 2007, 5: e129-10.1371/journal.pbio.0050129.

  42. 42.

    Lafos M, Kroll P, Hohenstatt ML, Thorpe FL, Clarenz O, Schubert D: Dynamic regulation of H3K27 trimethylation during Arabidopsis differentiation. PLoS Genet. 2011, 7: e1002040-10.1371/journal.pgen.1002040.

  43. 43.

    Gong ZH, Morales-Ruiz T, Ariza RR, Roldan-Arjona T, David L, Zhu J-K: ROS1, a repressor of transcriptional gene silencing in Arabidopsis, encodes a DNA glycosylase/lyase. Cell. 2002, 111: 803-814. 10.1016/S0092-8674(02)01133-9.

  44. 44.

    To TK, Kim JM, Matsui A, Kurihara Y, Morosawa T, Ishida J, Tanaka M, Endo T, Kakutani T, Toyoda T, Kimura H, Yokoyama S, Shinozaki K, Seki M: Arabidopsis HDA6 regulates locus-directed heterochromatin silencing in cooperation with MET1. PLoS Genet. 2011, 7: e1002055-10.1371/journal.pgen.1002055.

  45. 45.

    Weinhold A, Kallenbach M, Baldwin IT: Progressive 35S promoter methylation increases rapidly during vegetative development in transgenic Nicotiana attenuate plants. BMC Plant Biol. 2013, 13: 99-10.1186/1471-2229-13-99.

  46. 46.

    Casson SA, Lindsey K: The turnip mutant of Arabidopsis reveals that LEAFY COTYLEDON1 expression mediates the effects of auxin and sugars to promote embryonic cell identity. Plant Physiol. 2006, 142: 526-541. 10.1104/pp.106.080895.

  47. 47.

    Shibukawa T, Yazawa K, Kikuchi A, Kamada H: Possible involvement of DNA methylation on expression regulation of carrot LEC1 gene in its 5′-upstream region. Gene. 2009, 437: 22-31. 10.1016/j.gene.2009.02.011.

  48. 48.

    Buzas DM, Robertson M, Finnegan EJ, Helliwell CA: Transcription-dependence of histone H3 lysine 27 trimethylation at the Arabidopsis polycomb target gene FLC. Plant J. 2011, 65: 872-881. 10.1111/j.1365-313X.2010.04471.x.

  49. 49.

    Schmitges FW, Prusty AB, Faty M, Stützer A, Lingaraju GM, Aiwazian J, Sack R, Hess D, Li L, Zhou S, Bunker RD, Wirth U, Bouweester T, Bauer A, Ly-Hartig N, Zhao K, Chan H, Gu J, Gut H, Fischle W, Müller J, Thomä NH: Histone methylation by PRC2 is inhibited by active chromatin marks. Mol Cell. 2011, 42: 330-341. 10.1016/j.molcel.2011.03.025.

  50. 50.

    Yun JY, Tamada Y, Kang YE, Amasino RM: Arabidopsis trithorax-related3/SET domain GROUP2 is required for the winter-annual habit of Arabidopsis thaliana. Plant Cell Physiol. 2012, 53: 834-846. 10.1093/pcp/pcs021.

  51. 51.

    Mathieu O, Probst AV, Paszkowski J: Distinct regulation of histone H3 methylation at lysines 27 and by CpG methylation in Arabidopsis. EMBO J. 2005, 24: 2783-2791. 10.1038/sj.emboj.7600743.

  52. 52.

    Vaillant I, Paszkowski J: Role of histone and DNA methylation in gene regulation. Curr Opin Plant Biol. 2007, 10: 528-533. 10.1016/j.pbi.2007.06.008.

  53. 53.

    Weinhofer I, Hehenberger E, Roszak P, Hennig L, Köhler C: H3K27me3 profiling of the endosperm implies exclusion of polycomb group protein targeting by DNA methylation. PLoS Genet. 2010, 6: e1001152-10.1371/journal.pgen.1001152.

  54. 54.

    Deleris A, Stroud H, Bernatavichute Y, Johnson E, Klein G, Shubert D, Jacobsen SE: Loss of the DNA methyltransferase MET1 Induces H3K9 hypermethylation at PcG target genes and redistribution of H3K27 trimethylation to transposons in Arabidopsis thaliana. PLoS Genet. 2012, 8: e1003062-10.1371/journal.pgen.1003062.

  55. 55.

    Simon JA, Kingston RE: Polycomb mechanisms for getting to genomic targets, stopping transcriptional traffic, and staying put. Mol Cell. 2013, 49: 808-824. 10.1016/j.molcel.2013.02.013.

  56. 56.

    Bratzel F, López-Torrejón G, Koch M, Del Pozo JC, Calonje M: Keeping cell identity in Arabidopsis requires PRC1 RING-finger homologs that catalyze H2A monoubiquitination. Curr Biol. 2010, 20: 1853-1859. 10.1016/j.cub.2010.09.046.

  57. 57.

    Li W, Wang Z, Li J, Yang H, Cui S, Wang X, Ma L: Overexpression of AtBMI1C, a polycomb group protein gene, accelerates flowering in Arabidopsis. PLoS One. 2011, 6: e21364-10.1371/journal.pone.0021364.

  58. 58.

    Bouyer D, Roudier F, Heese M, Andersen ED, Gey D, Nowack MK, Goodrich J, Renou J-P, Grini PE, Colot V, Schnittger A: Polycomb repressive complex 2 controls the embryo-to-seedling phase transition. PLoS Genet. 2011, 7: e1002014-10.1371/journal.pgen.1002014.

  59. 59.

    Schubert D, Primavesi L, Bishopp A, Roberts G, Doonan J, Jenuwein T, Goodrich J: Silencing by plant polycomb-group genes requires dispersed trimethylation of histone H3 at lysine 27. EMBO J. 2006, 25: 4638-4649. 10.1038/sj.emboj.7601311.

  60. 60.

    Margueron R, Reinberg D: The polycomb complex PRC2 and its mark in life. Nature. 2011, 469: 343-349. 10.1038/nature09784.

  61. 61.

    Ellis JG, Tokuhisa JG, Llewellyn DJ, Bouchez D, Singh K, Dennis ES, Peacock WJ: Does the ocs-element occur as a functional component of the promoters of plant genes. Plant J. 1993, 4: 433-443. 10.1046/j.1365-313X.1993.04030433.x.

  62. 62.

    Chen WQ, Chao G, Singh KB: The promoter of a H 2 O 2 -inducible, Arabidopsis glutathione S -transferase gene contains closely linked OBF- and OBP1-binding sites. Plant J 1996, 10:955-966.,

  63. 63.

    Deng W, Buzas DM, Ying H, Robertson M, Taylor J, Peacock WJ, Dennis ES, Helliwell C: Arabidopsis polycomb repressive complex 2 binding sites contain putative GAGA factor binding motifs within coding regions of genes. BMC Genomics. 2013, 14: 593-10.1186/1471-2164-14-593.

  64. 64.

    Berger N, Dubreucq B, Roudier F, Dubos C, Lepiniec L: Transcriptional regulation of Arabidopsis LEAFY COTYLEDON2 involves RLE, a cis-element that regulates trimethylation of histone H3 at lysine 27. Plant Cell. 2011, 23: 4065-4078. 10.1105/tpc.111.087866.

  65. 65.

    Heck GR, Perry SE, Nichols KW, Fernandez DE: AGL15, a MADS domain protein expressed in developing embryos. Plant Cell. 1995, 7: 1271-1282. 10.1105/tpc.7.8.1271.

  66. 66.

    Harding EW, Tang W, Nichols KW, Fernandez DE, Perry SE: Expression and maintenance of embryogenic potential is enhanced through constitutive expression of AGAMOUS-Like 15. Plant Physiol. 2003, 133: 653-663. 10.1104/pp.103.023499.

  67. 67.

    Zheng Y, Ren N, Wang H, Stromberg AJ, Perry SE: Global identification of targets of the Arabidopsis MADS domain protein AGAMOUSLike15. Plant Cell. 2009, 21: 2563-2577. 10.1105/tpc.109.068890.

  68. 68.

    Wang F, Perry SE: Identification of direct targets of FUSCA3, a key regulator of Arabidopsis seed development. Plant Physiol. 2013, 161: 1251-1264. 10.1104/pp.112.212282.

  69. 69.

    Bentsink L, Hanson J, Hanhart CJ, Blankestijn-deVries H, Coltrane C, Keizer P: Natural variation for seed dormancy in Arabidopsis is regulated by additive genetic and molecular pathways. Proc Natl Acad Sci U S A. 2010, 107: 4264-4269. 10.1073/pnas.1000410107.

  70. 70.

    Nonogaki H: Seed dormancy and germination-emerging mechanisms and new hypotheses. Frontiers Plant Sci. 2014, 5: 233-10.3389/fpls.2014.00233.

  71. 71.

    Molitor AM, Bu Z, Yu Y, Shen W-H: Arabidopsis AL PHD-PRC1 complexes promote seed germination through H3K4me3-to-H3K27me3 chromatin state switch in repression of seed developmental genes. PLoS Genet. 2014, 10: e1004091-10.1371/journal.pgen.1004091.

  72. 72.

    Chanvivattana Y, Bishopp A, Schubert D, Stock C, Moon YH, Sung ZR, Goodrich J: Interaction of polycomb-group proteins controlling flowering in Arabidopsis. Development. 2004, 131: 5263-5276. 10.1242/dev.01400.

  73. 73.

    Aichinger E, Villar CB, Farrona S, Reyes JC, Hennig L, Kohler C: CHD3 proteins and polycomb group proteins antagonistically determine cell identity in Arabidopsis. PLoS Genet. 2009, 5: e1000605-10.1371/journal.pgen.1000605.

  74. 74.

    Xu L, Shen WH: Polycomb silencing of KNOX genes confines shoot stem cell niches in Arabidopsis. Curr Biol. 2008, 18: 1966-1971. 10.1016/j.cub.2008.11.019.

  75. 75.

    De Lucia F, Crevillen P, Jones AM, Greb T, Dean C: A PHD-polycomb repressive complex 2 triggers the epigenetic silencing of FLC during vernalization. Proc Natl Acad Sci U S A. 2008, 105: 16831-16836. 10.1073/pnas.0808687105.

Download references

Acknowledgements

The authors would like to thank Drs. Mohamed Fokar, Miyoung Kang and Million Tadege, for critically reading and providing helpful comments on the manuscript. We also thank Ms. Katie Pranger for her careful editing. This work was supported by the Oklahoma Agricultural Experiment Station, a grant from the Samuel Roberts Noble Foundation and an endowment from the Walter Sitlington Foundation to RDA.

Author information

Correspondence to Randy D Allen.

Additional information

Competing interests

The authors declare that they have no competing interests.

Electronic supplementary material

Additional file 1: Analysis of GSTF8 promoter sequence for the identification of putative cis -elements.(PDF 138 KB)

Additional file 2: List of PCR primers used in this study.(PDF 66 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Authors’ original file for figure 4

Authors’ original file for figure 5

Authors’ original file for figure 6

Authors’ original file for figure 7

Authors’ original file for figure 8

Rights and permissions

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Keywords

  • HSI2
  • VAL1
  • AGL15
  • DOG1
  • Transgene silencing
  • Seed-maturation
  • DNA methylation
  • Histone methylation
  • H3K27me3
  • 5-aza-2′-deoxycytidine