- Research article
- Open Access
Transcriptional silencing of 35S driven-transgene is differentially determined depending on promoter methylation heterogeneity at specific cytosines in both plus- and minus-sense strands
BMC Plant Biologyvolume 19, Article number: 24 (2019)
De novo DNA methylation triggered by short interfering RNAs is called RNA-directed DNA methylation (RdDM). Transcriptional gene silencing (TGS) through RdDM can be induced using a viral vector. We have previously induced RdDM on the 35S promoter in the green fluorescent protein (GFP)-expressing Nicotiana benthamiana line 16c using the cucumber mosaic virus vector. The GFP fluorescence phenotype segregated into two types, “red” and “orange” in the first self-fertilized (S1) progeny plants by the difference in degree of recovery from TGS on GFP expression. In the second self-fertilized generation (S2 plants), the phenotypes again segregated. Explaining what generates the red and orange types could answer a very important question in epigenetics: How is the robustness of TGS maintained after RdDM induction?
In bisulfite sequencing analyses, we found a significant difference in the overall promoter hypermethylation pattern between the red and orange types in S1 plants but little difference in S2 plants. Therefore, we assumed that methylation at some specific cytosine residues might be important in determining the two phenotypes. To find the factor that discriminates stable, robust TGS from the unstable TGS with incomplete inheritance, we analyzed the direct effect of methylated cytosine residues on TGS. Because it has not yet been demonstrated that DNA methylation at a few specific cytosine residues on known sequence elements can indeed determine TGS robustness, we newly developed a method by which we can directly evaluate the effect of specific methylation on promoter activity. In this assay, we found that the effects of the specific cytosine methylation on TGS differed between the plus- and minus-strands.
We found two distinct phenotypes, the stable and unstable TGS in the progenies of virus-induced TGS plants. Our bisulfite sequencing analyses suggested that methylation at some specific cytosine residues in the 35S promoter played a role in determining whether stable or unstable TGSs are induced. Using the developed method, we inferred that DNA methylation heterogeneity in and between the plus- and minus-strands can differentially determine TGS.
DNA methylation is a highly conserved epigenetic hallmark in plant genomes that controls gene expression. Regulation of gene expression by transcriptional gene silencing (TGS) of transgenes and endogenous genes has been extensively studied [1,2,3,4]. In plants, DNA methylation occurs in three sequence contexts: CG, CHG and CHH (where H is A, C or T). The CG methylation is maintained by METHYLTRANSFERASE1 (MET1), whereas CHG and CHH methylation is maintained by plant-specific CHROMOMETHYLASE3 (CMT3) and CHROMOMETHYLASE 2 (CMT2), respectively [1,2,3,4,5]. In addition, the de novo DNA methyltrasferase, DOMAINS REARRANGED METHYLTRANSFERASE2 (DRM2) is involved in both maintenance and initiation of DNA methylation .
In plants, DRM2 is required for de novo DNA methylation, which is triggered by short-interfering RNAs (siRNAs) . siRNAs guide DRM2 to the target sequences on the genomic DNA and directs de novo DNA methylation through RNA-directed DNA methylation (RdDM) [7,8,9,10,11,12,13,14,15,16,17,18], in which two plant-specific RNA polymerases, Pol IV and Pol V play essential roles. The 24-nucleotide (nt) siRNAs are generated and amplified by Pol IV, RNA-DEPENDENT RNA POLYMERASE 2 (RDR2) and DICER-LIKE 3 (DCL3). The 24-nt siRNAs are then loaded onto Argonaute (AGO) for the subsequent DRM2-mediated DNA methylation. AGO4 eventually interacts with the Pol V subunit, NUCLEAR RNA POLYMERASE E1 (NRPE1) and makes a complex with DRM2 to initiate DNA methylation. Recently, RdDM has been classified into two pathways, canonical and non-canonical. In the canonical pathway, a cascade of enzyme reactions (Pol IV-RDR2-DCL3) produces 24-nt siRNAs in concert, whereas several additional small RNAs seem to be involved in RdDM through the non-canonical pathway .
Assuming that RdDM is triggered by siRNAs from a promoter sequence, there are two options to induce specific DNA methylation artificially: transgene-induced or virus-induced. The first involves the creation of transgenic plants where double-stranded RNAs (dsRNAs) against the target promoter sequence are generated, for example, by an inverted repeat construct [20,21,22], while the second involves the virus-induced transcriptional gene silencing (VITGS), which is quite convenient because we do not have to produce transgenic plants . However, there is a large difference between the two methods in that transgene-induced siRNAs are generated even in the progeny plants but virus-driven siRNAs are not detectably carried into the next generation without viral infection.
A substantial number of viruses has been developed for VITGS, e.g., apple latent spherical virus (ALSV) [24,25,26], tobacco rattle virus [27, 28], barley stripe mosaic virus  and potato virus X [30, 31]. We have previously developed the cucumber mosaic virus (CMV)-based vector, CMV-A1 for both PTGS and TGS [32,33,34,35,36,37]. Using CMV-A1, we could successfully induce VITGS against some endogenous genes as well as a transgene [33, 38]. When we inoculated the green fluorescent protein (GFP)-expressing transgenic Nicotiana benthamiana (16c) with the CMV-A1 constructs containing various sizes of the 35S promoter sequences, we found that the GFP expression levels were downregulated due to RdDM and that the sizes of the virus-integrated sequences were important for DNA methylation .
The CaMV 35S promoter is such a strong promoter that it has been widely used for gene expression in plants. Several sequence domains affecting promoter activity have already been identified. The activation sequence factors 1 and 2 (as-1 and as-2), which are located in subdomains A1 and B1, respectively, have been found as a transcription factor-binding site [40,41,42,43,44,45]. Subdomains B2 and B4 have been also reported to be important for the promoter activity . However, it has not yet been demonstrated that specific DNA methylation on those sequence elements can control the promoter activity.
In our observation of the VITGS against the 35S promoter in 16c plants, we noticed that the virus-free S1 progeny plants comprised two phenotypes when exposed to UV light: red coloration (RED) resulting from autofluorescence of chlorophylls due to loss of GFP fluorescence by TGS and orange coloration (ORN) from the combined chlorophyll autofluorescence and GFP fluorescence. Most of the S2 progeny plants from the S1 ORN plants were ORN plants, but occasionally RED plants appeared. Virus vector-derived siRNAs corresponding to the 35S promoter sequence would not be generated in the progeny plants because progenies are no longer infected by the recombinant virus. We raise the following questions regarding unknown mechanisms for maintenance and initiation of the DNA methylation: What discriminates the RED and ORN plants? How are stable and unstable TGS controlled? Our preliminary analyses of the overall DNA methylation status of the 35S promoter in the two phenotypes indicated that there was not much difference between the RED and ORN plants, suggesting that methylation at specific cytosine residues on the target sequence may be important for inheritance of TGS. To answer the above questions, we here analyzed the DNA methylation in our VITGS-induced 16c plants in detail.
Differential induction of VITGS depending on the sizes of target sequences
We previously induced TGS of the 35S promoter in GFP-expressing N. benthamiana plants (16c) by the CMV-A1 vector (Fig. 1a). In this study, we further characterized this VITGS from a mechanistic viewpoint. Depending on the sizes of the 35S promoter sequences integrated into the viral vector, different degrees of TGS induction were found. When the 345-nt fragment corresponding to almost the entire 35S promoter was cloned, the virus-infected plants developed efficient VITGS, which was well maintained in the first self-fertilized (S1) progeny plants (S1:345) that fluoresced red under UV light (Figs. 1b and 2). The red-color must be due to chlorophyll autofluorescence because GFP expression was greatly reduced. On the other hand, the 116-nt fragment could induce TGS in the inoculated plants, in which weak GFP fluorescence still remained on the leaf veins, but failed to maintain the TGS in the S1 progeny plants (S1:116); GFP fluorescence in S1:116 plant was recovered to the level in 16c plant, yielding a yellow phenotype in leaves (Figs. 1b and 2). For the 208-nt fragment, we observed VITGS induction in the A1–208-inoculated plants, but the plants showed somewhat GFP fluorescence in the petiole of leaves. The S1 progeny plants from A1–208-inoculated plants segregated into two phenotypes: red type, which lost GFP fluorescence (designated S1:208 RED) and orange type, which was intermediate between S1:345 and 16c (designated S1:208 ORN) (Figs. 1b and 2). ORN plants seem to retain a certain level of GFP fluorescence, which masks the red autofluorescence. None of the S1 progeny plants from the A1–208-inoculated plants had fully recovered GFP fluorescence.
Because VITGS is essentially operated through RdDM, we then investigated the methylation status on the 35S promoter in S1:345, S1:208 ORN and 16c plants. Our bisulfite sequencing analyses revealed that the methylation level of the entire 35S promoter was significantly higher in S1:345 than in S1:208 ORN, and that the 35S promoter in 16c had little methylation (Fig. 3).
Analysis of methylation of 35S promoter in second-generation self-fertilized progenies from VITGS plants
When we produced S2 progeny plants by a second self-fertilization of the S1:208 ORN plants, the GFP fluorescence phenotype again segregated into RED and ORN plants: S2:208 RED and S2:208 ORN (Figs. 1b and 4a). When the S1:208 RED plant was self-fertilized, we observed only RED phenotype plants in the S2 progenies. These results raised questions: Does the stable, inherited TGS depend on the methylation level? What is the threshold methylation level for the stable TGS? However, because it is often the case that a few methylated cytosine residues can inhibit promoter activity, the stable TGS may not necessarily be explained only by a difference in hypermethylation. To find further clues, we scrutinized the intermediate TGS in S2:208 progenies.
Among S2 progenies, we found that the leaf veins (Vein) fluoresced brighter yellow than the other parts (leaf without vein, LwV) (Fig. 4a and b). We separated Vein and LwV and compared the GFP expression levels in those tissues by qRT-PCR (Fig. 4c). As we expected from the GFP fluorescence intensity, GFP was expressed in S2:208 ORN at a level of about 1/3 of that in 16c, while it was barely detected in S2:208 RED. On the other hand, unlike our expectation, the GFP level differed little between Vein and LwV in the samples.
We next used bisulfite sequencing analyses to search for differences in the methylation pattern on the 35S promoter between S2:208 RED and ORN plants. As shown in Additional file 1: Figure S1, there were few methylations in the 35S promoter in 16c plants, and no large difference was observed in methylation status between Vein and LwV tissues. The overall DNA methylations on the plus-sense of the 35S promoter in S2:208 RED and S2:208 ORN were shown in Fig. 5 for Vein, and in Additional file 2: Figure S2 for LwV tissues. The total methylation in Vein tissues was higher in S2:208 ORN (Fig. 5b), although total methylation in LwV tissues was higher in S2:208 RED (Additional file 2: Figure S2b). On the other hand, we noticed that asymmetric methylation (CHH methylation) was significantly lower, but symmetric CG methylation was higher in the Vein of S2:208 RED than S2:208 ORN (Fig. 5b). Because asymmetric methylation on one strand does not occur on its complementary strand at the corresponding positions and is not inherited by progeny, we further analyzed the methylation status on the complementary, minus-strand. DNA methylations on the minus-strand of the 35S promoter were shown in Fig. 5c for Vein and in Additional file 2: Figure S2c for LwV. When the methylation status of minus-strand was statistically analyzed, symmetric (CG and CHG) methylations in S2:208 RED were significantly higher than S2:208 ORN in both Vein and LwV tissues (Fig. 5d and Additional file 2: Figure S2d). In addition, the overall CHH methylation was greatly reduced in S2:208 RED compared to S2:208 ORN. It is thus likely that symmetric methylation at specific sites and a decrease in CHH methylation may be correlated with the stable, inherited TGS in the VITGS progeny plants. To make it easier to find any differences in methylation frequency between RED and ORN, we subtracted the values in S2:208 ORN from those in S2:208 RED (Fig. 6 and Additional file 3: Figure S3). In the subdomain A1, the methylation levels of the cytosine residues involved in symmetric methylation (cytosines at − 45, − 65, − 77 and − 80) were significantly higher in S2:208 RED than in S2:208 ORN (Fig. 6b and Additional file 3: Figure S3b). Symmetric methylation in plus-sense tended to be higher at the specific cytosine (position − 82) in S2:208 RED than S2:208 ORN, although statistically not supported. These results suggest that methylation in the minus-strand would be more important than that in the plus-strand for stable TGS.
In parallel with the VITGS experiments, we were creating transgenic Arabidopsis plants to screen for plants that showed TGS of the GFP gene under the 35S promoter; we obtained a transgenic line that contained a direct repeat sequence of the 35S promoter fused to the GFP gene. In the progeny generations (T1 to T5) of the direct repeat line, we noted two types of GFP fluorescence: one showing GFP fluorescence at an early stage, then developed post-transcriptional gene silencing (PTGS) for GFP expression and eventually induced TGS at a late stage (PTGS-TGS plants); the other already had perfect TGS in small seedlings (TGS plants), and the TGS was stably inherited by the next generation. From the results of bisulfite sequencing to analyze the methylation status of the 35S promoter sequences in those plants, we found that CG or CHG methylations at positions − 46, − 66, − 78 and − 82 in the subdomain A1 were higher in the TGS plants than in the PTGS-TGS plants (Additional file 4: Figure S4). Furthermore, the overall CHH methylation was greatly reduced in the TGS plants compared to the PTGS-TGS plants. These observations therefore agree well with the results of VITGS.
CG methylation at a few cytosine residues in the minus-strand of the 105-bp sequence containing the subdomain A1 caused significant reduction in promoter activity
To verify that the specific CG methylation could indeed compromise the gene expression under the 35S promoter, we developed a new method, in which we used the 105-bp 35S promoter sequence made by hybridizing two chemically synthesized oligonucleotides with methylated cytosine residues in a Golden gate cloning strategy (Fig. 7a). To evaluate the ability of the synthetic promoter to induce the downstream gene expression, we ligated it to the reporter gene firefly luciferase (Fluc), and Fluc activity was measured in protoplasts after transfection of the ligation products (Fig. 7b). As shown in Fig. 8a, Fluc activity was greatly reduced when both plus- and minus-strands were methylated at the cytosine positions − 45, − 46, − 65, − 66, − 77 and − 78. However, when either plus- or minus-strand was methylated, we found that methylation on the plus-strand barely suppressed Fluc activity and sometimes even increased it. In contrast, cytosine methylations at the minus-strand greatly reduced Fluc activity to 1/3 of the control, 35S-Fluc (Fig. 8b). Intriguingly, when we introduced both strands-methylated DNA, Fluc was suppressed at least to a level similar to that of the minus-strand or seemed to be more effective in the reduction of promoter activity, although the difference was not statistically significant (Fig. 8b).
CHH methylation at a few cytosine residues in the subdomain A1 also contributed to reducing promoter activity
Because only three cytosine residues at CHH context (− 33, − 56, and − 84 in plus-sense) were found to be highly methylated within the domain A (− 90 ~ + 1) of the 35S promoter (Fig. 5a, Additional file 2: Figure S2a), we integrated methylated cytosine residues at those three CHH positions; two residues are located in the subdomain A1 and one residue is located in subdomain mp (region from − 45 to + 1). As shown in Fig. 8c, the introduced CHH methylation reduced the luciferase activity down to 1/3 of the control although the specific CG methylation tested was twice more effective (1/6 of the control).
Unstable TGS phenotype in the progeny from VITGS-16c plants
In our analyses of RdDM on the 35S promoter induced by VITGS, we found discrete phenotypes of the TGS patterns against the GFP gene downstream of the 35S promoter. The TGS seemed to depend on the sizes of the 35S promoter sequence inserted in the viral vector. The TGS in the 16c plants inoculated with A1–345 was maintained in the S1 progeny, whereas we observed two phenotypes (RED and ORN) in the S1 progeny plants after inoculation with A1–208. The RED, from chlorophyll autofluorescence, represents a strong GFP TGS, and the ORN seems to be an intermediate TGS between those of 16c and the RED plants. In the S1 progeny from the plants inoculated with A1–116, most of the S1 plants lost the TGS of the GFP gene that had been observed in the virus-inoculated plants. Based on these observations, we considered that by analyzing the 35S promoters in these phenotypes, we might be able to reveal some important factor(s), which can control the robustness of the RdDM-mediated TGS.
We then analyzed the S2:208 plants in detail to obtain a clue about fixing a stable TGS for next generation. As shown in Fig. 4, the S2:208 ORN plants clearly had stronger GFP fluorescence in the stems and leaf veins (Vein) than in the leaf areas without veins (LwV). The GFP mRNA levels in the S2:208 ORN plants were significantly higher than in the S2:208 RED plants. Although we initially expected that the GFP levels in Vein would be higher than those in LwV, GFP expression differed little between the two tissues, perhaps because chlorophyll content in Vein tissues is lower than in LwV and the chlorophyll autofluorescence interfered less with the GFP fluorescence.
Effect of cytosine methylation in the 35S promoter in a TSM assay on downstream gene expression
Here, we successfully developed an assay by which we can analyze the effect of specific cytosine methylation in a promoter sequence on the downstream gene expression. In addition, we can even discriminate DNA strand-sense when specific methylation is introduced into a target promoter. We here call this method strand-sense-specific methylation (triple S methylation, TSM) assay (Fig. 7). Previously, the promoter sequences that are important for TGS had been identified mainly using point-mutation and subdomain-deletion strategies [40,41,42,43,44,45,46]. However, without introducing actual methylation to specific cytosine residues, we can never conclude that methylation of the predicted sites is really responsible for the influence on the expression of downstream gene. Because we can directly manipulate cytosine methylation in a target promoter with a TSM assay, various promoter sequences can be analyzed to elucidate a link between specific cytosine methylation and promoter activities. The results of our TSM assays demonstrated that TGS levels of the 35S promoter might be differentially regulated depending on the heterogeneity of cytosine methylation in and between the plus- and minus-strands corresponding to the target domain sequences. Especially cytosine methylation at positions close to the as-1 element (− 83 to − 63) is very important for stable, inherited TGS, which agrees well with the previously predicted cytosine residues for the promoter activity [40, 43, 45]. To be more specific, we found that TGS of the 35S promoter occurred efficiently when the cytosine residues (positions − 78, − 66 and − 46) on the minus-strand were methylated. On the other hand, surprisingly, the methylation of the corresponding sites on the plus-strand had little effect on the TGS, but the maximum TGS was observed when both strands were methylated at the specific cytosine residues. Our TSM assay also indicated that CHH methylation in the subdomain A1 significantly reduced promoter activity. However, considering that CHH methylation is not maintained through the next generation, the stable inherited TGS in the RED phenotype will be driven mainly by some specific symmetric methylations. No report has previously discussed the effect of heterogeneity of DNA methylation in and between the plus- and minus-DNA strands in relation to TGS. Although we have always assumed that both strands would be simultaneously and uniformly methylated in a model for RdDM , we should rather consider that the two strands may not be necessarily uniformly methylated in nature. Whether our observation for the 35S promoter can be generally applied to other promoters or not should be independently tested using a similar assay. Although in this study, we focused only on methylation of the subdomain A1, subdomains B2 and B4 in the 35S promoter have also been proved to be essential for promoter activity . It would be worth evaluating the cytosine methylation on the other elements in the 35S promoter using the TSM assay.
Link between symmetric and asymmetric methylation in inherited DNA methylation
In our bisulfite sequencing analyses, we noticed an interesting phenomenon that was always associated with the stable, inherited TGS in many replicates. In the RED plants in S2:208, we found that the overall levels of asymmetric (CHH) methylation were significantly reduced in both plus- and minus-strands, while the levels of symmetric methylation (CG and CHG) were conversely increased (Fig. 5 and Additional file 2: Figure S2). Most of the S1 plants generated from A1–116-inoculated plants lost TGS of the GFP gene. This unstable TGS in next generation may be associated with a decrease in symmetric methylation and an increase of asymmetric methylation in the target sequence. For cytosine methylation of the 35S promoter in the TGS-induced progeny plants (S1 to S4 generations), similar observations have also been shown in the study of the ALSV vector . To understand the involvement of methylation status in TGS induction, we can test how the level of asymmetric methylation in the 35S promoter can affect TGS by a method like our TSM assay. Although we do not know the exact reason for a decrease in asymmetric methylation at this moment, it is conceivable that some factor during reproduction may have effects in the association between RdDM and inherited methylation. Considering the results together, we infer that the generation of the RED and ORN plants observed among the VITGS progeny plants may depend on the accumulation levels of symmetric methylation at some specific cytosine sites and a decrease in asymmetric methylation in the progeny plants.
We induced TGS against the GFP gene under the 35S promoter through RdDM by a viral vector containing various sizes of the promoter sequence. In the first self-fertilized generation, the stable TGS seemed to be correlated with high levels of methylation on the promoter. However, we observed an unstable TGS in the second generation; S2 progenies were segregated into RED and ORN. Our bisulfite sequencing analyses suggested that specific methylation at a few cytosine residues in the subdomain A1 and a decrease in frequency of asymmetric methylation might be important for efficient induction of TGS through RdDM. To probe the importance of the methylation at specific cytosine residues, we developed a method by which we can analyze a direct effect of methylated cytosine residues on TGS in a strand-specific manner. In this assay, we found that cytosine methylation on either or both of the plus- and minus-sense promoter sequences could differentially drive gene expression downstream of the promoter. We therefore infer that robust TGS may be determined by the heterogeneity of DNA methylation at specific sites in and between the plus- and minus-strands.
Nicotiana benthamiana line 16c was obtained from Dr. David Baulcombe. The plants were grown with 16-h light (8-h dark) at 24 °C.
Inoculation with virus and GFP observation
Fragments from the 35S promoter were amplified by PCR and inserted into the CMV-A1 vector. N. benthamiana line 16c plants were inoculated with in vitro transcripts of CMV as previously described by Otagaki et al. . Total RNA was isolated 15 days postinoculation from the infected tissues. Viral cDNAs were PCR-amplified using a primer pair of 2b-5up and R2–2814-R2 (Additional file 5: Table S1) and the TaKaRa One Step RNA PCR Kit AMV (TaKaRa) kit according to the supplier’s instruction. To confirm that the viruses did not have serious mutations during their replication, the amplified viral cDNAs were directly sequenced: the sequencing chromatograms were shown in Additional file 6: Figure S5. GFP fluorescence was examined with a Keyence Microscope VB-7000 with the GFP filter (OP-42313).
DNA methylation assay
DNA was extracted using the Illusta DNA extraction Kit PhytoPure (GE Healthcare) according to the supplier’s instructions. Bisulfite treatment of DNA was performed by using EZ DNA Methylation-Lightning Kit (Zymo Research) according to the supplier’s protocol. To amplify the target sequences, primer pair 35S-346F-bisuT/35S + 1A-bisuA were used for the plus-strand, and primer pair 35S (−)-5-BS/35S (−)-3-BS were used for the minus-strand (Additional file 5: Table S1). The primers were initially designed by the program, MethPrimer (http://www.urogene.org/cgi-bin/methprimer/methprimer.cgi). To amplify the entire (~ 345 bp) core sequence of the 35S promoter, we modified the designed primer so that they contained cytosine residues as less as possible, and synthesized degenerated primers. The PCR was conducted by Takara Epi-taq (TaKaRa) and the PCR conditions were as follows, 40 cycles of 94 °C for 30 s, 55 °C for 30 s and 72 °C for 30 s. The PCR products were then cloned into the pTAC1 vector using the Dyna Express TA PCR Cloning Kit (Bio Dynamics Laboratory). For clone-based bisulfite sequencing, we used 22–24 clones and the sequencing data were aligned using MEGA version 6 . The graphs were created in Microsoft Excel (Additional file 7: Table S2).
Quantitative RT-PCR analysis
Total RNA was isolated using TRIzol reagent (Invitrogen) according to the supplier’s instructions. For the DNA and RNA extraction in the S2 progenies, we separated class I and class II veins (Vein)  and leaf tissues without those veins (LwV). Extracted RNA was treated by DNase I recombinant RNase-free (Roche Applied Science), and then cDNA was synthesized from 0.2 μg of RNA using the PrimeScript RT reagent Kit (TaKaRa) according to the supplier’s instruction. qRT-PCR analysis was carried out by the Applied Biosystems StepOnePlus Real-time PCR system using the PowerUp SYBR Green Master Mix (Applied Biosystems). The PCR cycling conditions were 40 cycles of 95 °C for 15 s, 55 °C for 30 s and 72 °C for 30 s. The GFP expression levels were calculated using the comparative Ct method and normalized by the expression of the L23 gene. Target sequences were amplified using primer pair mGFP-5-160/mGFP3–160 for the GFP gene and primer pair Nb-L23–5-110/Nb-L23–3-110 for the internal control.
Golden gate cloning and TSM assay
All oligonucleotides used in this study are listed in Additional file 5: Table S1. The unmethylated 35S promoter fragment was amplified from pBI221 (Clontech) by PCR. The target sequences were amplified using primer pair 35S-GGC-5/35S-GGC-3. The Fluc PCR product was also amplified by PCR using primer pair LUC-GGC-5/LUC-GGC-3. The custom-synthesized oligonucleotides (oligo DNAs) with specific methylation at three CG sites (positions − 78, − 66 and − 46) and at three CHH sites (positions − 33, − 56 and − 84) were prepared by the Custom DNA Synthesis Service (Hokkaido System Science). The oligo DNAs (100 μM) were annealed in oligo annealing buffer (10 mM Tris-HCl, 1 mM EDTA, 100 mM NaCl) with gradually cooling after incubation at 95 °C for 4 min. To make a DNA fragment where either strand was methylated, 35S-CG-DNA (+), 35S-CG-DNA (−) and 35S-CHH-DNA (+) were annealed with either unmethylated oligonucleotide 35S-3-120 PM or 35S-5-120 PM and then amplified by PCR. The PCR cycling conditions were 40 cycles of 94 °C for 30 s, 55 °C for 30 s and 72 °C for 30 s. All four types of DNA constructs were first digested by BsaI and then ligated to the firefly luciferase reporter (Fluc) gene as essentially described as golden gate cloning by Engler et al. . We then analyzed the promoter activity in protoplasts. We call this method strand-sense-specific methylation (triple S methylation, TSM) assay.
Protoplast transfection and luciferase assay
Protoplasts were prepared from the leaves of N. benthamiana as described by Shimura et al. . Prepared Fluc-fused 35S promoter fragments (0.5 μg or 1 μg) were directly transfected to protoplasts together with pE-Rluc (0.45 μg) as an internal control in the presence of polyethylene glycol. After 18 h incubation, luciferase activities (Fluc and Rluc) were assayed using the Dual-Luciferase Reporter Assay System (Promega) and a fluorometer (Wallac 1420 ARVO MX) as previously described by Shimura et al. .
Apple latent spherical virus
Cucumber mosaic virus
DOMAINS REARRANGED METHYLTRANSFERASE2
double stranded RNA
Green fluorescent protein
NUCLEAR RNA POLYMERASE E1
Post-transcriptional gene silencing
RNA directed DNA methylation
RNA dependent RNA polymerase 2
short interfering RNA
Transcriptional gene silencing
Virus-induced transcriptional gene silencing
Baulcombe D. RNA silencing in plants. Nature. 2004;431:356–63.
Bender J. DNA methylation and epigenetics. Annu Rev Plant Biol. 2004;55:41–68.
Chan SW, Henderson IR, Jacobsen SE. Gardening the genome: DNA methylation in Arabidopsis thaliana. Nat Rev Genet. 2005;6:351–60.
Lucy AP, Guo HS, Li WX, Ding SW. Suppression of post-transcriptional gene silencing by a plant viral protein localized in the nucleus. EMBO J. 2000;19:1672–80.
Zhang Y, Harris CJ, Liu Q, Liu W, Ausin I, Long Y, Xiao L, Feng L, Chen X, Xie Y, Chen X, Zhan L, Feng S, Li JJ, Wang H, Zhai J, Jacobsen SE. Large-scale comparative epigenomics reveals hierarchical regulation of non-CG methylation in Arabidopsis. Proc Natl Acad Sci U S A. 2018;115:E1069–74.
Cao X, Jacobsen SE. Role of the Arabidopsis DRM methyltransferases in de novo DNA methylation and gene silencing. Curr Biol. 2002;12:1138–44.
Cao X, Aufsatz W, Zilberman D, Mette MF, Huang MS, Matzke M, Jacobsen SE. Role of the DRM and CMT3 methyltransferases in RNA-directed DNA methylation. Curr Biol. 2003;13:2212–7.
Haag JR, Pikaard CS. Multisubunit RNA polymerases IV and V: purveyors of non-coding RNA for plant gene silencing. Nat Rev Mol Cell Biol. 2011;12:483–92.
Law JA, Jacobsen SE. Establishing, maintaining and modifying DNA methylation patterns in plants and animals. Nat Rev Genet. 2010;11:204–20.
Matzke M, Aufsatz W, Kanno T, Daxinger L, Papp I, Mette MF. Genetic analysis of RNA-mediated transcriptional gene silencing. Biochem Biophys Acta. 2004;1677:129–41.
Matzke M, Kanno T, Daxinger L, Huettel B, Matzke AJ. RNA-mediated chromatin-based silencing in plants. Curr Opin Cell Biol. 2009;21:367–76.
Matzke MA, Kanno T, Matzke AJ. RNA-directed DNA methylation: the evolution of a complex epigenetic pathway in flowering plants. Annu Rev Plant Biol. 2015;66:243–67.
Mette MF, Aufsatz W, van der Winden J, Matzke M, Matzke AJM. Transcriptional silencing and promoter methylation triggered by double-stranded RNA. EMBO J. 2000;19:5194–201.
Pikaard CS, Haag JR, Pontes OM, Blevins T, Cocklin R. A transcription fork model for pol IV and pol V-dependent RNA-directed DNA methylation. Cold Spring Harb Symp Quant Biol. 2012;77:205–12.
Pontes O, Li CF, Costa Nunes P, Haag J, Ream T, Vitins A, Jacobsen SE, Pikaard CS. The Arabidopsis chromatin-modifying nuclear siRNA pathway involves a nucleolar RNA processing center. Cell. 2006;126:79–92.
Saze H, Tsugane K, Kanno T, Nishimura T. DNA methylation in plants: relationship to small RNAs and histone modifications, and functions in transposon inactivation. Plant Cell Physiol. 2012;53:766–84.
Sijen T, Vijn I, Rebocho A, van Blokland R, Roelofs D. Transcriptional and posttranscriptional gene silencing are mechanistically related. Curr Biol. 2001;11:436–40.
Zhang H, Zhu JK. RNA-directed DNA methylation. Curr Opin Plant Biol. 2011;14:142–7.
Cuerda-Gil D, Slotkin RK. Non-canonical RNA-directed DNA methylation. Nat Plants. 2016;2:16163.
Deng S, Chua NH. Inverted-repeat RNAs targeting FT intronic regions promote FT epression in arabidopsis. Plant Cell Physiol. 2015;56:1667–78.
Kinoshita Y, Saze H, Kinoshita T, Miura A, Soppe WJ, Koornneef M, Kakutani T. Control of FWA gene silencing in Arabidopsis thaliana by SINE-related direct repeats. Plant J. 2007;49:38–45.
Miki D, Shimamoto K. De novo DNA methylation induced by siRNA targeted to endogenous transcribed sequences is gene-specific and OsMet1-independent in rice. Plant J. 2008;56:539–49.
Burch-Smith TM, Anderson JC, Martin GB, Dinesh-Kumar SP. Applications and advantages of virus-induced gene silencing for gene function studies in plants. Plant J. 2004;39:734–46.
Kon T, Yoshikawa N. Induction and maintenance of DNA methylation in plant promoter sequences by apple latent spherical virus-induced transcriptional gene silencing. Front Microbiol. 2014;5:595.
Sasaki S, Yamagishi N, Yoshikawa N. Efficient virus-induced gene silencing in apple, pear and Japanese pear using Apple latent spherical virus vectors. Plant Methods. 2011;7:15.
Yamagishi N, Yoshikawa N. Virus-induced gene silencing in soybean seeds and the emergence stage of soybean plants with Apple latent spherical virus vectors. Plant Mol Biol. 2009;71:15–24.
Xu H, Xu L, Yang P, Cao Y, Tang Y, He G, Yuan S, Ming J. Tobacco rattle virus-induced PHYTOENE DESATURASE (PDS) and Mg-chelatase H subunit (ChlH) gene silencing in Solanum pseudocapsicum L. Peer J. 2018;6:e4424.
Jones L, Ratcliff F, Baulcombe D. RNA-directed transcriptional gene silencing in plants can be inherited independently of the RNA trigger and requires Met1 for maintenance. Curr Biol. 2001;11:747–57.
Holzberg S, Brosio P, Gross C, Pogue GP. Barley stripe mosaic virus-induced gene silencing in a monocot plant. Plant J. 2002;30:315–27.
Faivre-Rampant O, Gilroy EM, Hrubikova K, Hein I, Millam S, Loake GJ, Birch P, Taylor M, Lacomme C. Potato virus X-induced gene silencing in leaves and tubers of potato. Plant Physiol. 2004;134:1308–16.
Jones L, Hamilton A, Voinnet O, Thomas C, Maule A, Baulcombe D. RNA-DNA interactions and DNA methylation in post-transcriptional gene silencing. Plant Cell. 1999;11:2291–301.
Goto K, Kobori T, Kosaka Y, Natsuaki T, Masuta C. Characterization of silencing suppressor 2b of cucumber mosaic virus based on examination of its small RNA-binding abilities. Plant Cell Physiol. 2007;48:1050–60.
Kanazawa A, Inaba J, Shimura H, Otagaki S, Tsukahara S, Matsuzawa A. Virus-mediated efficient induction of epigenetic modifications of endogenous genes with phenotypic changes in plants. Plant J. 2011;65:156–68.
Nagamatsu A, Masuta C, Matsuura H, Kitamura K, Abe J, Kanazawa A. Downregulation of flavonoid 3′-hydroxylase gene expression by virus-induced gene silencing in soybean reveals the presence of a threshold mRNA level associated with pigmentation in pubescence. J Plant Physiol. 2009;166:32–9.
Nagamatsu A, Masuta C, Senda M, Matsuura H, Kasai A, Hong JS. Functional analysis of soybean genes involved in flavonoid biosynthesis by virus-induced gene silencing. Plant Biotechnol J. 2007;5:778–90.
Otagaki S, Arai M, Takahashi A, Goto K, Hong JS, Masuta C. Rapid induction of transcriptional and post-transcriptional gene silencing using a novel Cucumber mosaic virus vector. Plant Biotechnol. 2006;23:259–65.
Tasaki K, Terada H, Masuta C, Yamagishi M. Virus-induced gene silencing (VIGS) in Lilium leichtlinii using the Cucumber mosaic virus vector. Plant Biotechnol. 2016;33:373–81.
Kanazawa A, Inaba J, Kasai M, Shimura H, Masuta C. RNA-mediated epigenetic modifications of an endogenous gene targeted by a viral vector: a potent gene silencing system to produce a plant that does not carry a transgene but has altered traits. Plant Sig Behav. 2011;8:1090–3.
Otagaki S, Kawai M, Masuta C, Kanazawa A. Size and positional effects of promoter RNA segments on virus-induced RNA-directed DNA methylation and transcriptional gene silencing. Epigenetics. 2011;6:681–91.
Benfey PN, Chua NH. The cauliflower mosaic virus 35S promoter: combinatorial regulation of transcription in plants. Science. 1990;250:959–66.
Benfey PN, Ren L, Chua NH. Combinatorial and synergistic properties of CaMV 35S enhancer subdomains. EMBO J. 1990;9:1685–96.
Lam E, Benfey PN, Gilmartin PM, Fang RX, Chua NH. Site-specific mutations alter in vitro factor binding and change promoter expression pattern in transgenic plants. Proc Natl Acad Sci U S A. 1989;86:7890–4.
Lam E, Chua NH. ASF-2: a factor that binds to the cauliflower mosaic virus 35S promoter and a conserved GATA motif in cab promoters. Plant Cell. 1989;1:1147–56.
Fang RX, Nagy F, Sivasubramaniam S, Chua NH. Multiple cis regulatory elements for maximal expression of the cauliflower mosaic virus 35S promoter in transgenic plants. Plant Cell. 1989;1:141–50.
Kanazawa A, O’ Dell M, Hellens RP. The binding of nuclear factors to the as-1 element in the CaMV 35S promoter is affected by cytosine methylation in vitro. Plant Biol. 2007;9:435–41.
Bhullar S, Datta S, Advani S, Chakravarthy S, Gautam T, Pental D, Burma PK. Functional analysis of cauliflower mosaic virus 35S promoter: re-evaluation of the role of subdomains B5, B4 and B2 in promoter activity. Plant Biotechnol J. 2007;5:696–708.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.
Mekuria T, Bamunusinghe D, Payton M, Verchot-Lubicz J. Phloem unloading of potato virus X movement proteins is regulated by virus and host factors. Mol Plant-Microbe Interact. 2008;8:1106–17.
Engler C, Kandzia R, Marillonnet S. A one pot, one step, precision cloning method with high throughput capability. PLoS One. 2008;11:e3647.
Shimura H, Kogure Y, Goto K, Masuta C. Degree of RNA silencing and the ability of a viral suppressor vary depending on the cell species in a protoplast system. J Gen Plant Pathol. 2008;74:326–30.
We thank David Baulcombe (University of Cambridge, UK) for kindly providing Nicotiana benthamiana 16c. We are also grateful to Peter Palukaitis (Seoul Women’s University, Korea) for reviewing the manuscript.
This work was supported by the promotion services of the New Energy and Industrial Technology Development Organization (NEDO) (project code: P16009). The funding body was used for the design of the study; collection, analysis, and interpretation of data; or in writing the manuscript.
Availability of data and materials
All data generated and analyzed during this study are included in this published article (Additional file 1: Figure S1, Additional file 2: Figure S2, Additional file 3: Figure S3, Additional file 4: Figure S4, Additional file 5: Table S1, Additional file 6: Figure S5 and Additional file 7: Table S2).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Methylation status in the plus-strand of the 35S promoter of the 16c plants. Methylation frequency in the plus-strand of the 35S promoter in the Vein and LwV tissues was analyzed by bisulfite sequencing. N is the number of the clones used for sequencing. The x-axis is the position relative to the nucleotide distance from the transcription start sites. (TIFF 2598 kb)
Figure S2. Cytosine methylation frequency in the 35S promoter in the LwV tissues of S2:208 RED and S2:208 ORN. a Comparison of methylation status in the plus-sense of the 35S promoter between S2:208 RED and S2:208 ORN. b Summary of the results in a to show differences in CHH, CHG, CG and total methylation. c Comparison of methylation status in the minus-sense of the 35S promoter between S2:208 RED and ORN. Asterisk indicates cytosine residue that is significantly different in methylation frequency between RED and ORN as explained in Additional file 3: Figure S3b. d Summary of the results in c to show differences in CHH, CHG, CG and total methylation. N is the number of the clones used for the bisulfite sequencing. The x-axis shows the position relative to the transcription start site (+1). The asterisks in b and d indicate a statistical significance in methylation frequencies by two-tailed Fisher’s exact test (* P < 0.05, ** P < 0.01). (TIFF 6706 kb)
Figure S3. Difference in specific cytosine methylation between S2:208 RED and S2:208 ORN LwV tissues. a Difference in methylation frequencies in the overall plus- and minus-strand of the 35S promoter in S2:208 RED and S2:208 ORN. To make it easier to find any differences in methylation frequencies in each strand, values of S2:208 RED LwV (Additional file 2: Figure S2a and c, upper graphs) were subtracted from those of S2:208 ORN LwV (Additional file 2: Figure S2a and c, lower graphs). b Close-up of the subdomain A1 in a. The x-axis is the nucleotide sequence of the subdomain A1. Nucleotide sequences of the plus- and minus-strand of the subdomain A1 are indicated below the graph. Methylated CG and CHG sites are underlined. The asterisk in b indicates a statistical significance between S2:208 RED and S2:208 ORN by two-tailed Fisher’s exact test (* P < 0.05). (TIFF 7061 kb)
Figure S4. Methylation status in the plus-strand of the 35S promoter in transgenic plants. a, b Methylation frequencies of GFP-silenced lines. Lines 2–6–1-7-6 (a) and 2–6–1-7-2 (b) were the T5 progeny lines derived from the original line 2–6–1-7, which contained a direct repeat of the 35S promoter followed by the GFP gene sequence. GFP expression in 2–6–1-7 and 2–6–1-7-6 initially decreased as a result of post-transcriptional gene silencing (PTGS) and later by TGS, while GFP expression in 2–6–1-7-2 was stably suppressed by TGS. Twelve to fifteen clones were used for the bisulfite sequencing. c Values for a were subtracted from those for b to show differences in methylation frequencies between 2 and 6–1-7-2 and 2–6–1-7-6. The positions − 82, − 78, − 66 and − 46 indicated specific cytosine residues located in the subdomain A1. The x-axis shows the position relative to the transcription start site (+ 1). Positions of the subdomains A1, B2 and B4 are also indicated. (TIFF 4447 kb)
Table S1. List of primers used in this study. (XLSX 10 kb)
Figure S5. 35S promoter sequences in the viral genomes of the vector isolated from infected tissues. a-c Sequencing chromatograms of the inserts containing each of 345-, 208- and 116-bp portion (a to c, respectively). Total RNA was isolated 15 days postinoculation, and RT-PCR-amplified fragment were directly sequenced. We confirmed that the original sequences integrated into the viral vector did not change in the infected tissues. (TIFF 7535 kb)