Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

A Role of Tomato UV-Damaged DNA Binding Protein 1 (DDB1) in Organ Size Control via an Epigenetic Manner

  • Jikai Liu ,

    Contributed equally to this work with: Jikai Liu, Xiaofeng Tang, Lanyang Gao

    Affiliations Ministry of Education Key Laboratory for Bio-resource and Eco-environment, State Key Laboratory of Hydraulics and Mountain River Engineering, College of Life Science, Sichuan University, Chengdu, China, School of Biotechnology and Food Engineering, Hefei University of Technology, Hefei, China

  • Xiaofeng Tang ,

    Contributed equally to this work with: Jikai Liu, Xiaofeng Tang, Lanyang Gao

    Affiliations Ministry of Education Key Laboratory for Bio-resource and Eco-environment, State Key Laboratory of Hydraulics and Mountain River Engineering, College of Life Science, Sichuan University, Chengdu, China, School of Biotechnology and Food Engineering, Hefei University of Technology, Hefei, China

  • Lanyang Gao ,

    Contributed equally to this work with: Jikai Liu, Xiaofeng Tang, Lanyang Gao

    Affiliation Ministry of Education Key Laboratory for Bio-resource and Eco-environment, State Key Laboratory of Hydraulics and Mountain River Engineering, College of Life Science, Sichuan University, Chengdu, China

  • Yongfeng Gao,

    Affiliation Ministry of Education Key Laboratory for Bio-resource and Eco-environment, State Key Laboratory of Hydraulics and Mountain River Engineering, College of Life Science, Sichuan University, Chengdu, China

  • Yuxiang Li,

    Affiliation Ministry of Education Key Laboratory for Bio-resource and Eco-environment, State Key Laboratory of Hydraulics and Mountain River Engineering, College of Life Science, Sichuan University, Chengdu, China

  • Shengxiong Huang,

    Affiliation Ministry of Education Key Laboratory for Bio-resource and Eco-environment, State Key Laboratory of Hydraulics and Mountain River Engineering, College of Life Science, Sichuan University, Chengdu, China

  • Xiaochun Sun,

    Affiliation Ministry of Education Key Laboratory for Bio-resource and Eco-environment, State Key Laboratory of Hydraulics and Mountain River Engineering, College of Life Science, Sichuan University, Chengdu, China

  • Min Miao,

    Affiliations Ministry of Education Key Laboratory for Bio-resource and Eco-environment, State Key Laboratory of Hydraulics and Mountain River Engineering, College of Life Science, Sichuan University, Chengdu, China, Department of Plant, Soil and Entomological Sciences, University of Idaho, Moscow, Idaho, United State of America

  • Hui Zeng,

    Affiliation Ministry of Education Key Laboratory for Bio-resource and Eco-environment, State Key Laboratory of Hydraulics and Mountain River Engineering, College of Life Science, Sichuan University, Chengdu, China

  • Xuefen Tian,

    Affiliation Ministry of Education Key Laboratory for Bio-resource and Eco-environment, State Key Laboratory of Hydraulics and Mountain River Engineering, College of Life Science, Sichuan University, Chengdu, China

  • Xiangli Niu,

    Affiliation School of Biotechnology and Food Engineering, Hefei University of Technology, Hefei, China

  • Lei Zheng,

    Affiliation School of Biotechnology and Food Engineering, Hefei University of Technology, Hefei, China

  • Jim Giovannoni,

    Affiliation United States Department of Agriculture-Agricultural Research Service, Robert Holly Center and Boyce Thompson Institute for Plant Research, Cornell University, Ithaca, New York, United State of America

  • Fangming Xiao ,

    fxiao@uidaho.edu (FX); liuyongsheng1122@yahoo.com.cn (Y. Liu)

    Affiliation Department of Plant, Soil and Entomological Sciences, University of Idaho, Moscow, Idaho, United State of America

  • Yongsheng Liu

    fxiao@uidaho.edu (FX); liuyongsheng1122@yahoo.com.cn (Y. Liu)

    Affiliations Ministry of Education Key Laboratory for Bio-resource and Eco-environment, State Key Laboratory of Hydraulics and Mountain River Engineering, College of Life Science, Sichuan University, Chengdu, China, School of Biotechnology and Food Engineering, Hefei University of Technology, Hefei, China

Abstract

Epigenetic modification generally refers to phenotypic changes by a mechanism other than changes in DNA sequence and plays a significant role in developmental processes. In this study, we found that overexpression of one alternatively spliced tomato DDB1 transcript, DDB1F that is prevalently present in all tested tissues, resulted in reduction of organ size. Transgenic plants constitutively expressing the DDB1F from a strong cauliflower mosaic virus (CaMV) 35S promoter displayed moderately reduced size in vegetative organs (leaves and stems) and radically decreased size in reproductive organs (flowers, seeds and fruits), in which several genes encoding negative regulators for cell division were upregulated. Significantly, reduction of organ size conferred by overexpression of DDB1F transgene appears not to segregate in the subsequent generations, suggesting the phenotypic alternations are manipulated in an epigenetic manner and can be transmitted over generations. This notion was further substantiated by analysis of DNA methylation level at the SlWEE1 gene (encoding a negative regulator of cell division), revealing a correlation between less methylation in the promoter region and elevated expression level of this gene. Thus, our results suggest DDB1 plays an important role in regulation of the epigenetic state of genes involved in organogenesis, despite the underlying mechanism remains to be elucidated.

Introduction

The UV-damaged DNA binding protein 1 (DDB1) was originally identified as a nuclear factor that binds to UV-damaged DNA and participates in versatile DNA repair pathways at the stage of binding and recognition [1]. A growing body of evidence suggests that, conserved from yeast to human, the DDB1 acts as an adapter linking the CUL4-ROC1 catalytic core to substrate receptors to form the DDB1-CUL4-ROC1 complex that is recently identified as a cullin-RING ubiquitin ligase [2]. This DDB1-CUL4-based ubiquitin E3 ligase involves in many physiological and developmental processes, such as transcription, cell cycle, cell death and embryonic development [3], [4]. Particularly in plants, the DDB1-CUL4-based E3 ligase complex (DDB1-CUL4-RBX1) has been found to function in regulation of photomorphogenesis [5], ABA signaling [6], flowering time control [7], parental imprinting [8], UV-B tolerance and genome integrity [9], [10].

Epigenetic modifications on chromatin structure without changes in DNA sequence, in many cases through DNA methylation or/and histone modification, usually result in repressing gene expression. It has been demonstrated that epigenetic modifications play a role in various physiological processes in plants ranging from plant growth and reproduction regulation to stress responses [11], [12], [13]. In plants, the epigenetic state of a particular gene can be inherited during cell propagation, through both mitosis and meiosis. In the case of latter, the epigenetic state can be transmitted over generations, despite the mechanisms underlying this transgenerational epigenetic inheritance are still largely unknown [14]. For example, the epigenetic state of enhanced homologous recombination in Arabidopsis induced by stress cues can be transmitted over 4 generations [15].

It is generally thought there are two major mechanisms responsible for epigenetic modifications: DNA methylation and histone modification [13]. DNA methylation on the cyclic carbon-5 of cytosine can be asymmetric (mCpHpH) and symmetric (mCpG/mCpHpG) and often occurs in the promoter region, which results in repression of gene transcription [11]. One of the key questions regarding epigenetic regulations in developmental biology is how the epigenetic states are maintained faithfully through successive rounds of cell division. Recent investigations in metazoans, plants and microorganisms point to an important and conserved role of the DDB1-CUL4-containing ubiquitin E3 ligase in perpetuating epigenetic marks on chromatin, presumably via regulating the histone modification or/and DNA methylation. For example, in human, CUL4-DDB1 complex was shown to be essential for histone H3 methylation at K4 or K9 and K27, respectively, through interaction with multiple WD40-repeat proteins and Polycomb-group proteins [16], [17]. In mammalian cells, the histone H4 monomethylase PR-Set7 is shown to be a target of ubiquitin-conjugated degradation regulated by CUL4-DDB1Cdt2-mediated PCNA-dependent E3 ligase activity during DNA damage and replication [18], [19]. The Arabidopsis MSI1, a WD40 repeat protein, physically interacts with DDB1A and forms CUL4-DDB1A-MSI1 E3 ligase complex that is required to maintain MEDEA parental imprinting by interacting with the epigenetic regulatory Polycomb repressive complex 2 (PRC2) [8]. Another WD40 repeat protein, MSI4, is reported as a DDB1- and CUL4-associated factor that represses the expression of flowering locus C (FLC) and flowering locus T (FT) through its association with a CLF-Polycomb Repressive Complex 2 (PRC2) in Arabidopsis [7]. However, a very recent study suggests that MSI4 may not act as part of the PRC2-like complexes but act in concert to establish a repressive chromation environment at FLC for its transcriptional silencing [20]. Recently, several research groups have demonstrated that ubiquitin ligase components CUL4, DDB1 and DCAF (DDB1- and Cul4-associated factor) are essential for DNA methylation in the filamental fungus Neurospora crassa [21][23].

Tomato (Solanum lycopersicum) is an economically and experimentally important crop. Tomato fruits carrying the high pigment 1 mutations (hp-1 and hp-1w) are characterized by an enhanced biogenesis of plastids coupled with elevated levels of functional metabolites including carotenoids and flavonoids [24][27]. Previously, we have characterized the tomato HIGH PIGMENT 1 (HP1) gene encoding a protein homologous to human DDB1 and two alternatively spliced transcripts have been identified without further elucidation on their function [28]. Genetic analysis suggests this HP1/DDB1 plays a pivotal role in regulation of plastid division [28][30]. The HP1/DDB1 is also implicated in pathogenesis-related (PR) gene induction and basal defense in response to non-pathogenic Agrobacterium tumefaciens [31]. In addition, ectopic expression study infers that the DDB1 protein may be profoundly involved in regulating a variety of developmental aspects in tomato plants [27]. Nevertheless, how DDB1 regulates these various physiological processes is largely unknown.

In the current investigation, we found a significant role of DDB1 in the organogenesis in tomato, apparently via an epigenetic regulation. Transgenic tomato overexpressing the alternatively spliced transcript of tomato DDB1 exhibited reduction of organ size and enhanced expression of genes negatively regulating cell division probably due to less methylation. Significantly, these phenotypic alternations in both T1 and T2 generation plants were not necessarily associated with the DDB1F transgene, implicating an inheritable epigenetic state existing in these plants.

Results

Differential Expression Patterns of Two Alternatively Spliced DDB1 Transcripts (DDB1F and DDB1+15)

In eukaryotes, a large number of precursor mRNAs with introns can undergo alternative splicing (AS) to produce structurally and functionally different proteins from the same gene [32]. Previously, two different transcripts derived from HIGH PIGMENT 1 locus, designated as DDB1F and DDB1+15, respectively, (GenBank accession No. AY531660 and AY531661) were identified [28]. The ORF (open reading frame) of DDB1+15 contains additional 15 bp nucleotides (tgcatttgtctgcag) compared to that of DDB1F (Figure 1A). To determine the differential tissue-specific expression of the two alternatively spliced transcripts DDB1F and DDB1+15, we performed a detailed analysis of transcript composition and expression ratio in seedlings, roots, stems, leaves, flowers and fruits at different developmental stages by using a polymerase chain reaction (PCR)-based approach (described in detail in Materials and Methods). As shown in Figure 1B, in all the tissues tested, two transcripts were differentially accumulated: about 70–90% of the detected DDB1 transcripts were the isoform of DDB1F, suggesting DDB1F isoform is prevalently expressed in tomato plants.

thumbnail
Figure 1. Expression analysis of transcript variants of DDB1.

(A) Schematic diagram of DDB1 gene splicing pattern. Introns (I) and exons (E) are marked by lowercase and uppercase letters, respectively. Numbers above the letter and arrows indicate the position of the nucleotide in the genome sequence and the site of alternative splicing, respectively. The extra 15 bp nucleotides of DDB1+15 transcript are shown in grey. (B) The percentage of two DDB1 transcripts in different tissues. SD: seedlings; R: roots; SM: stems; L: leaves; F: flowers; IG1: immature green fruits at 7-DPA; IG2: immature green fruits at 14-DPA; MG: mature green fruits at 35-DPA. The light gray and black bars indicate the expression percentage of the DDB1+15 and DDB1F transcripts, respectively. Data presented are the means of three technical replicates. Similar results were obtained in at least two independent experiments. (C) Identification of the percentage of two DDB1 transcripts in seedlings under hormones or stress treatment. Total RNAs were extracted from 10-day-old seedlings with application of abscisic acid (ABA), auxin (indole acetic acid; IAA), cytokinin (6-benzylaminopurine, BA), gibberellins (GA) and mannitol. 10-day-old seedlings without any hormone and stress treatment were used as control. The light gray and black bars indicate the expression percentage of the DDB1+15 and DDB1F transcripts, respectively.

https://doi.org/10.1371/journal.pone.0042621.g001

To further investigate the possible effect of phytohormone and stress cues on the alternative splicing of the DDB1 transcripts, 10-day-old tomato seedlings were treated with various phytohormones or subjected to abiotic stress, including four phytohormones [abscisic acid (ABA), auxin (indole acetic acid; IAA), cytokinin (6-benzylaminopurine, BA), and gibberellins (GA)] and the stressing chemical mannitol. As shown in Figure 1C, no significant change was observed in untreated or treated seedlings and the DDB1F was still the dominant transcript.

Overexpression of DDB1F Isoform Affects Organ Size in an Epigenetic Manner

The differential expression level of two alternative DDB1 transcript isoforms prompted us to investigate the biological significance of DDB1F or DDB1+15 transcript by using gain-of-function strategy. We generated transgenic tomato overexpressing DDB1F or DDB1+15 under the control of the constitutive CaMV 35S promoter. 9 and 10 primary transgenic T0 lines harboring 35S::DDB1F or 35S::DDB1+15 were created and verified by PCR using NPTII (Kanr)-specific primers. All 10 35S::DDB1+15 lines (named 35S::DDB1+15 a to j) were morphologically indistinguishable from the nontransformed wild type (WT) plants despite the overexpression of 35S::DDB1+15 transgene was verified by real-time RT-PCR analysis, whereas 6 of 9 35S::DDB1F lines displayed significantly small fruits, which weighed only 11–16% of the nontransformed WT fruits (Figure 2A, B, C). Further real-time RT-PCR assay indicated the 35S::DDB1F transcripts were dramatically overaccumulated in these 6 lines producing small fruits (named 35S::DDB1F a, b, c, f, g and i) but not in the other 3 lines (Figure 2A), suggesting the phenotype of reduced size fruits is caused by overexpression of 35S::DDB1F. Since the morphological appearance of T0 transgenic lines may be affected by the tissue culture process, we followed up the observed phenotype of 35S::DDB1F transgenic lines to T1 generation by self-pollination. To our surprise, the small fruit phenotype caused by 35S::DDB1F did not segregate with the 35S::DDB1F transgene. All T1 plants derived from the 6 lines of 35S::DDB1F transgenic plants produced small fruits (Figure 2D). Subsequently, we selected T1 plants carrying hemizygous 35S::DDB1F transgene derived from 3 T0 lines (lines a, b and c) to generate T2 generation. As expected, the 35S::DDB1F transgene exhibited a 3∶1 segregation ratio in T2 population, and all plants produced small fruits regardless harboring 35S::DDB1F transgene or not. In the mean time, we also selected nullizygous (without 35S::DDB1F transgene) T1 plants to generate the T2 and T3 nullizygous 35S::DDB1F tomato plants. As expected, no significant alteration in fruit size was observed among T1, T2 and T3 fruits (Figure 2E). Thus, we conclude: firstly, the small fruit phenotype elicited by overexpression of 35S::DDB1F transgene can be regulated in an epigenetic manner; secondly, this epigenetic state can be inherited over generations.

thumbnail
Figure 2. DDB1 mRNA level analysis and phenotypes in transgenic plants.

(A) Real-time PCR analysis of total DDB1 mRNA level in WT Ailsa Craig (WT), 35S::DDB1+15 and 35S::DDB1F T0 transgenic lines. a–j and a–i stand for 10 and 9 independent transgenic lines for 35S::DDB1+15 and 35S::DDB1F constructs, respectively. Each bar represents three repetitions from each RNA sample (derived from pools of at least three fruits per plant). Error bars representing standard errors are shown in each case. (B) Representative red ripe fruits from field-grown plants of WT Ailsa Craig (WT) and 3 independent T0 transgenic lines (a, b and c) containing 35S::DDB1+15 or 35S::DDB1F. Bar = 1 cm. (C) Comparison of red ripe fruit weight from WT Ailsa Craig (WT) and 3 independent T0 transgenic lines (a, b and c) containing 35S::DDB1+15 or 35S::DDB1F. Mean values from 10 fruits from each line and standard errors are shown. Statistical analysis was performed using Student's t-test. (*P<0.05, **P<0.01). (D) Mature green fruits from WT Ailsa Craig (WT), 35S::DDB1+15 and 35S::DDB1F T1 transgenic lines. (E) Comparison of red ripe fruit weight of T1, T2 and T3 nullizygous 35S::DDB1F plants segregated out from 3 T0 transgenic lines (lines a, b and c). Mean values from 10 fruits from each line and standard errors are shown. Error bars representing standard errors are shown in each case. p and n represent plants with or without transgene, respectively.

https://doi.org/10.1371/journal.pone.0042621.g002

The Altered Organogenesis Is Characterized by Substantial Reduction in Organ Size

The reduced size of fruits in the 35S::DDB1F transgenic plants was the morphological change that caught our attention at first. To investigate other possible altered morphological features in the progeny of the 35S::DDB1F transgenic plants, we examined the T2 generation plants displaying small fruits with or without 35S::DDB1F transgene more closely. As shown in Figure 3, besides reduced size in fruits (Figure 3A), a significant size reduction occurred in all the observed organs including seedling roots and stems (Figure 3B), stems (Figure 3C), mature leaves (Figure 3D), flowers (Figure 3E) and seeds (Figure 3F). In contrast, there's no significant alteration in organ/tissue size observed in T2 generation plants derived from transgenic overexpressing DDB1+15 driven by CaMV 35S promoter (Figure 3A, B, C, D, E, F).

thumbnail
Figure 3. Organ size reduction in DDB1F overexpression line.

(A) Representative red ripe fruits from field-grown plants of WT Ailsa Craig (WT), 35S::DDB1+15 and 35S::DDB1F T2 transgenic lines (a, b, c). Bar = 1 cm. (B) Representative 7-day-old seedlings of WT Ailsa Craig (WT), 35S::DDB1+15 and 35S::DDB1F T2 transgenic lines (a, b, c) germinated under white-light (16 h light/8 h dark photoperiod). Bar = 1 cm. (C) Comparison of stems from WT Ailsa Craig (WT), 35S::DDB1+15 and 35S::DDB1F T2 transgenic lines, showing the reduction in stem size. Bar = 1 cm. (D) Comparison of mature leaves from WT Ailsa Craig (WT), 35S::DDB1+15 and 35S::DDB1F T2 transgenic lines (a, b, c). Bar = 5 cm. (E) Comparison of flowers from WT Ailsa Craig (WT), 35S::DDB1+15 and 35S::DDB1F T2 transgenic lines. Bar = 1 cm. (F) Comparison of seeds from WT Ailsa Craig (WT), 35S::DDB1+15 and 35S::DDB1F transgenic lines. p and n represent plants with or without transgene, respectively. Bar = 1 cm.

https://doi.org/10.1371/journal.pone.0042621.g003

As both cell size and cell number could affect final organ size, we next sought to investigate the cytological basis of the drastically altered organ size, focusing on the fruit pericarp tissue due to a sharp reduction observed in fruit size. Paraffin sections were generated and then histochemically stained with fast green using pericarp tissues harvested from the T2 small-fruited plants (with or without 35S::DDB1F transgene) or non-transgenic WT plants at breaker stage. As shown in Figure 4, at the breaker stage wherein the cell division and cell expansion ceased completely, the pericarp tissue of small fruits was significantly thinner and displayed an incisive decrease of cell layers, while the cell size was slightly reduced as compared with that of the non-transgenic WT plants. These observations suggest that overexpression of the DDB1F isoform has influence on cell division.

thumbnail
Figure 4. Cellular parameters of pericarp of fruits at breaker stage from WT Ailsa Craig (WT) and DDB1 overexpression lines.

(A) Number of cell layers across pericarp. (B) Mean cross-sectional area of one pericarp cell. (C) Pericarp thickness. (D) Cross section of fruit (left) and pericarp (right) are shown for WT Ailsa Craig (WT), 35S::DDB1+15 and 35S::DDB1F T2 transgenic lines. p and n represent plants with or without transgene, respectively. Bar: 1 cm (F, left sections), 500 µm (F, right sections).

https://doi.org/10.1371/journal.pone.0042621.g004

Reduction of Organ Size Is Correlated With the Enhanced Expression of Genes Encoding Negative Regulators of Mitotic Division

It has been demonstrated in recent publications that several negative cell-division regulators are involved in controlling meristem size and final organ size in plants, including WEE1, CCS52A and DA1 [33][37]. It is possible that the reduction of organ size originally elicited by overexpression of DDB1F is implemented through regulating these genes. To test this hypothesis, fruit tissues from the small-fruited T2 plants with or without 35S::DDB1F transgene or from non-transgenic WT plants were harvested at cell division stage (7-DAP) for total RNA extraction. Gene-specific primers were designed to determine the expression of tomato orthologs of these 3 genes, SlWEE1 (GenBank accession No. AM180939.1), SlCCS52A (SGN-U570015) and SlDA1 (SGN-U584698). Real-time PCR analysis using the gene-specific primers revealed the transcription of all the 3 tested genes were significantly enhanced in fruit tissues from the small-fruited T2 plants, regardless of containing 35S::DDB1F transgene or not, using SlUBI3 (GenBank accession No. X58253) as internal control (Figure 5). We also performed real-time PCR using SlGAPDH (GenBank accession No. U97257.1) as the second reference gene, and obtained the similar results (Figure S1). Thus, our results suggest the reduced organ size phenotype is probably due to upregulation of genes involved in cell-division and, significantly, this is not dependent on the original cue – the 35S::DDB1F transgene.

thumbnail
Figure 5. Real-time RT-PCR analysis of mRNA levels of genes regulating cell division in 7-DPA fruits from WT Ailsa Craig (WT) and 35S::DDB1F T2 transgenic lines (a, b, c).

p and n represent plants with or without transgene, respectively. Each bar represents three repetitions from each RNA sample (derived from pools of at least three fruits per plant). Error bars representing standard errors are shown in each case. Statistical analysis was performed using Student's t-test (*P<0.05, **P<0.01). Similar results were obtained in at least two independent experiments.

https://doi.org/10.1371/journal.pone.0042621.g005

Correlation Between Elevated Expression of the SlWEE1 Gene and Less Methylation in Its Promoter Region

It has been demonstrated that DDB1 acts as a CUL4-based ubiquitin E3 ligase component and plays a pivotal role in regulation of chromatin compaction during cell division by recruiting enzymes involved in histone methylation for ubiquitin/proteasome-mediated degradation [38], [39]. In addition, there is a growing body of evidence indicating histone modification is a prerequisite for DNA methylation and consequently represses gene expression [40]. We thus predicted that overexpression of DDB1F would affect above genes via regulation of their methylation states and this epigenetic regulation itself could be implemented via an epigenetic manner, in another word, may indirectly impacting methylation of genes of the transgenic progenies not containing the 35S::DDB1F transgene. To test this notion, we examined the correlation between the enhanced SlWEE1 expression and its possibly reduced DNA methylation during the cell division stage of fruit development. Using the approach of bisulfite genomic sequencing, DNA methylation levels of the promoter region of the SlWEE1 gene (genomic sequences of −457 to −147 bp upstream the start codon) were determined between the immature fruits (7-DPA) of nontransgenic WT plants and both T1 and T2 plants with reduction of organ size derived from three individual primary 35S::DDB1F transgenic lines. As shown in Figure 6, in fruits from the T1 progenies with or without the 35S::DDB1F transgene, both CpG dinucleotide motifs and non-CG motifs in the tested region of the SlWEE1 genomic sequences displayed less methylation level than in fruits from nontransgenic WT plants (Figure 6B). Similar result was also found in T2 progenies (Figure 6C). These results not only indicate a correlation between enhanced expression of the SlWEE1 and decreased methylation in its promoter region, but also suggest that, at least in part, the reduced organ size phenotype is possibly due to upregulation of genes controlling cell division via epigenetic demethylation, which is originally triggered by overexpression of DDB1F transgene.

thumbnail
Figure 6. Analysis of SlWEE1 methylation through bisulfite genomic sequencing.

(A) A schematic diagram of SlWEE1 gene structure, with +1 indicating the transcription start site and rectangle representing the methylation region −457 to −147. (B) and (C) Histograms of the percentage of CpG, CpHpG, and CpHpH methylation at region −457 to −147 of the SlWEE1 gene from WT Ailsa Craig (WT) and 35S::DDB1F T1 (B) and T2 (C) transgenic lines (a, b, c). p and n represent plants with or without transgene, respectively. H represents A, T or C.

https://doi.org/10.1371/journal.pone.0042621.g006

Discussion

One of the key questions in developmental biology is how epigenetic states of gene activity are maintained faithfully through successive rounds of cell division. Recent investigations in metazoans, plants and microorganisms point to an important and conserved role of the CUL4-DDB1-containing ubiquitin ligase in perpetuating epigenetic marks on chromatin [7], [16], [18], [19], [22]. Thus, it is reasonable to hypothesize that genetic manipulation of DDB1 gene, either by knock-down or up-regulation, may result in impact on the epigenetic state of certain genes. In this study, we have demonstrated that overexpression of DDB1F isoform, but not DDB1+15 isoform, can significantly affect the organogenesis in tomato plants, apparently via an epigenetic manner.

DDB1F Is the Prevalent Form and Plays a Role in Organogenesis

We have previously showed that tomato DDB1 gene possesses two alternatively spliced transcripts [28], designated as DDB1F and DDB1+15, respectively. In the current investigation, we demonstrate the DDB1F isoform accounts for the overwhelming majority of the expressed DDB1 mRNAs in all the tissues tested (Figure 1B). Unlike many other examples where alternative splicing is regulated by hormones and other abiotic stress [41], the ratio of these two transcripts remains largely unchanged in response to 4 different hormones and mannitol, again suggesting the DDB1F isoform plays a dominant role in the physiological processes in tomato plants. This hypothesis is further supported by our data from experiments in which overexpression of DDB1F, but not DDB1+15, exhibited reduced-size organs, including fruits, seedling roots and stems, stems, leaves, flowers and seeds (Figure 3A, B, C, D, E, F). These results are consistent with a recent observation that overexpression of the DDB1 gene from the CaMV 35S promoter in transgenic tomato does not cause morphological changes [42]. Nevertheless, although our results do not provide much information on the physiological significance of DDB1+15 transcript, we cannot rule out the possible role of DDB1+15 transcript in other physiological processes. Further experiments such as specific knockdown of DDB1+15 transcript will be helpful to determine its function.

Correlation Between Reduction of Organ Size and Enhanced Expression of Genes Encoding Negative Regulators of Cell-Division

Overexpression of the DDB1F resulted in reduced size of many organs, which is supported by the histological analysis showing that the decreased thickness of fruit pericarp is mainly due to the attenuation of cell division as evident from the reduced cell layers and cell numbers (Figure 4). This is consistent with observation of slightly faster cell division rate in hp1 mutant plants in response to cytokinin [43]. Moreover, in T2 generation plants derived from the primary 35S::DDB1F transgenic lines, correlation was found between reduction of organ size and enhanced expression of several genes encoding negative regulators of cell-division (SlWEE1, SlCCS52A and SlDA1) (Figure 5). It is possible that the upregulation of these negative regulator genes leads to premature termination of the mitotic division due to more production of the encoded proteins, resulting in proliferative growth arrest in organs. In fact, both WEE1, a kinase partially responsible for the functional abrogation of CDK/CYC complex whose activity is required for mitotic division, and CCS52A, acting as a substrate-specific activator of the anaphase-promoting complex (APC) ubiquitin ligase, are involved in the transition of mitotic cycles to endoreduplication cycles [37], [44][46]. Increases in either WEE1 or CCS52A activity result in inhibition of mitosis and reduced cell number, but still allow DNA replication to occur with cell enlargement as a consequence of endoreduplication [36], [37]. Thus, further investigation of the possible altered enzyme activity of WEE1 in the transgenic tomato plants will be helpful to understand the complexity of these cell division-related factors, which will be one of the topics of our future research. Nevertheless, since the final organ size, determined by cell number and size, must be compromised by these two counteracting cellular processes, it is contradictory but understandable that altering tomato CCS52A (SlCCS52A) expression in either a negative or positive manner would result in reduction of fruit size [37]. The expression levels of several mitotic-associated genes, including CDKB (RefSeq: NM_001246970.1), CYCB2 (RefSeq: NM_001246857.1) and CYCA (RefSeq: NM_001246833.1) were also determined. But there was no significant alterations observed between WT plants and transgenic plants (Figure S2), suggesting the specificity of DDB1 effect on genes involved in regulation of cell division.

Inheritable Epigenetic State Initiated by Overexpression of DDB1F Is Not Dependent on the 35S::DDB1F Transgene

Based on the observation that reduction of organ size phenotype did not segregate in T1 population (containing homozygous, hemizygous or nullizygous 35S::DDB1F allele), in which the 35S::DDB1F transgene must show segregation, we conclude that this unique phenotype of reduction of organ size is controlled by an epigenetic manner. By the same token, T2 plants (without harboring 35S::DDB1F transgene) exhibiting organ size reduction were found not necessarily associated with the 35S::DDB1F transgene, also suggesting this epigenetically phenotypic change can be transmitted over generations.

We speculate that overexpression of 35S::DDB1F could cause epigenetic change on certain genes and this epigenetic alternation is inheritable. The DNA methylation analysis on the SlWEE1 promoter region suggests the overexpression of 35S::DDB1F probably leads to decrease of DNA methylation level. However, how DDB1F involves in manipulation of DNA methylation is largely unknown. Recently, several research groups have demonstrated that ubiquitin ligase components CUL4, DDB1 and DCAF (DDB1- and Cul4-associated factor), are essential for DNA methylation in the filamentous fungus Neurospora crassa [21][23]. Genetic analysis by impairing components of the E3 complex suggests that the CUL4-DDB1 E3 ubiquitin ligase indirectly regulates DNA methylation via histone H3K9 trimethylation [22]. Coimmunoprecipitation-based biochemical analysis also showed CUL4-DDB1 ligase recruits histone methyltransferase DIM-5 for histone H3K9 trimethylation [21]. Moreover, interaction between CUL4-DDB1 ligase components and the gene-specific chromatin was shown to donate epigenetic repression or parental imprinting of target genes by associations with the epigenetic regulatory Polycomb Repressive Complex 2 in Arabidopsis [7], [8]. All these observations point to a conserved positive role for DDB1 protein in maintenance of DNA methylation and/or repressed gene activities. However, in both human cells and Xenopus egg extract, biochemical and molecular evidence demonstrated that Set8 protein, a methyltransferase that monomethylates histone H4 on lysine 20 (H4K20me1) [47], [48], is targeted for proteolysis in S phase during cell division by the CUL4-DDB1 E3 ubiquitin ligase [38], [39]. Similarly, in our tomato system, it is possible that DDB1-CUL4 E3 ligase manipulates cellular factors involved in methylation of genes (such as SlWEE1) via ubiquitination-mediated proteasome-dependent degradation. These factors likely are positive regulators of DNA methylation, such as enzymes involved in methylation of histones or factors negatively regulating the DNA demethylation. In the tomato system, overexpression of DDB1F may result in expediting the degradation of these positive regulators, consequently ablating the methylation reaction of target genes for their up-regulation. Nevertheless, the decreased methylation of the SlWEE1 gene is apparently initiated by, but not dependent on, the 35S::DDB1F transgene, despite the underlying mechanism remains to be elucidated, which will be the subject of our future research.

Materials and Methods

Plant Growth Conditions

Tomato plants were germinated and grown in the greenhouse under standard conditions (26°C day, 18°C night; 16 h light, 8 h dark). Primary transformants (T0), transgenic generation 1 and 2 (T1 and T2) plants were planted in the greenhouse and transplanted into the field 35 days later.

Analysis of DDB1 Transcripts

Total RNAs were extracted from seedlings, roots, stems, leaves, flowers and fruits pericarps at various developmental stages (7, 14, 35-DPA). The ratio of alternatively spliced forms of DDB1 mRNA was evaluated by PCR-based analysis. Briefly, the first-strand cDNA prepared from these organs/tissues was used in RT-PCR. Gene-specific primer set (forward primer: 5′-CATTCCTGTGTTGGCAATTTC-3′; reverse primer: 5′-AGATTATTTTGAGCCCATGGC-3′) was used to amplify both DDB1+15 and DDB1F. The RT-PCR products containing mixed fragments derived from either DDB1F or DDB1+15 were cloned into the pGM-T vector and subjected to a second PCR reaction to verify the identity of DDB1 isoform and quantify the original amount of mRNA of DDB1F or DDB1+15. This second PCR was conducted using an alternative splicing-specific primer set (forward, 5′- GGTTTACCAGTGCATTTGTCTG-3′; reverse, 5′-AGATTATTTTGAGCCCATGGC-3′) that specifically target DDB1+15,but not DDB1F without containing the extra 15 bp nucleotides, and to unambiguously distinguish DDB1+15 from DDB1F. The accuracy of this second PCR method was confirmed by the subsequent sequencing analysis for the randomly selected clones.

Plasmid Construction and Tomato Transformation

DNA manipulations were carried out by using standard procedures [49]. Sequences from DDB1F and DDB1+15 cDNA (accession No. AY531660 and AY531661) for construction of overexpression vectors were amplified from proper plasmids by PCR with primers DDB1OE-F (5′-GGATCCATGAGTGTATGGAACTACGTGG-3′), DDB1OE-R (5′-GAATTCCTAATGCAACCTTGTCAACTC-3′). The fragments were inserted into vector PBI121 (driven by the 35S promoter) at the BamHI and EcoRI restriction enzyme sites. Then, NOS terminator fragment from PBI121 vector was inserted into EcoRI restriction site by PCR using primers NOS-F (5′-GAATTCGCTCGAATTTCCCCGATC-3′), NOS-R (5′-AGTGAATTCCCGATCTAGTAAC-3′). The direction of insertion was confirmed by sequencing.

Transgenic plants were generated by Agrobacterium tumefaciens-mediated transformation according to the method described by Fillatti et al. (1987) [50], and transformed lines were first selected for kanamycin (70 mg l−1) resistance and then analyzed by PCR to determine the presence of T-DNA. The primers designed to the NPTII (Kanr) marker of PBI121 for confirmation of integration were 5′-ATTGAACAAGATGGATTGCACG-3′ and 5′-CTCGTCAAGAAGGCGATAGAAG -3′.

Molecular Analyses

Total RNAs were extracted using Trizol reagent according to the protocol provided by the manufacturer (Invitrogen, http://www.Invitrogen.com/) and treated with DNaseI (TaKaRa; http://www.takara-bio.com). Primers for real-time RT-PCR were designed for DDB1 (DDB1 rt-F, 5′-CCGAGAATTGCAGACAGAATGT-3′; DDB1 rt-R, 5′-CCCTCTTCATGCTTGAAAATCA-3′), SlWEE1 (WEE1 rt-F, 5′-GAGCAAATCGGTAGTGGGAAC-3′; WEE1 rt-R, 5′-CCATCAAAGCCTGTCTCCTATC-3′), SlCCS52A (CCS52A rt-F, 5′-TTCAGGAAGCCGAGACAAGAG-3′; CCS52A rt-R, 5′-AAGCCGATTATCATTTCCACCT-3′), SlDA1 (DA1 rt-F, 5′-TTTCAATGTCAGATAACCGTCCA-3′; DA1 rt-R, 5′-AATTAGACCAGCCGGATTTGTT-3′), and the control SlUBI3 (UBI3 rt-F, 5′-AGGTTGATGACACTGGAAAGGTT-3′; UBI3 rt-R, 5′-AATCGCCTCCAGCCTTGTTGTA-3′), SlGAPDH (GAPDH rt-F, 5′-AGGACTGGAGAGGTGGAAGAGC-3′; GAPDH rt-R, 5′-CGACAACGGAGACATCAGCAGT-3′). The real-time PCR was performed using an SsoFast EvaGreen Supermix (Bio-rad catalog #172-5203). Each sample was amplified in triplicate and all PCR reactions were performed on the Applied Biosystems StepOne Real-Time PCR System (Applied biosystems, http://www.appliedbiosystems.com.cn/). Dissociation curve analysis was performed at the end of each run to ensure that unique products were amplified. The tomato SlUBI3 gene was used as an internal reference. The RT-PCR conditions were as follows: 95°C for 30 s, followed by 40 cycles of 95°C for 5 s, 60°C for 20 s. The expression level was normalized to the SlUBI3 control, and relative expression values were determined against the buffer-treated sample or the WT Ailsa Craig (WT) sample using the 2−ΔΔCt method. To confirm the specificity of the PCR reaction, PCR products were verified on a 1% agarose gel for the accurate amplification product size.

Phenotypic and Cytological Analysis

After fruit ripening, fruits from individual plants were collected and weighed. At least 10 fruits from individual plants were measured. After weighing, fruits were transversally cut in half and imaged.

Cytological analysis was performed according to the methods described previously [51]. Fruits at breaker stage of WT Ailsa Craig and transgenic lines were prepared for cytological analysis by a paraffin-embedding method. An equatorial slice was excised and cut into fragments less than 4 mm wide before immersion in the fixative. During fixation, a partial vacuum was applied to extract intercellular gas. Samples were rinsed, dehydrated through an ethanol series, and embedded in paraffin. Sections (10 µm thick) were made with LEICA RM2235 microtome, stained with fast green, and photographed on OLYMPUS BX51 microscope with an OLYMPUS DP71 camera.

Bisulfite Sequencing

Bisulfite sequencing was performed according to the previously described methods [52]. 4 µg of genomic DNA was heated at 94°C for 5 min, and quenched on ice. 5.5 µL NaOH (3 M; freshly prepared) was added and incubated at 42°C for 30 min, followed by the addition of 416 µL of bisulfite solution to the denatured DNA. Bisulfite solution was prepared as follows: 3.76 g of sodium bisulfite (Fisher S654-500) was dissolved in 4 mL of water, adjusted to pH 5.0, 330 µL of 20 mM hydroquinone was added (Sigma-Aldrich H-9003), and the volume was adjusted to 10 mL water. Samples were incubated in a PCR machine for five cycles of 55°C for 3 h and 95°C for 5 min. After bisulfite conversion, the universal DNA Purification Kit was u sed to remove extra salt (Tiangen). Then, NaOH was added to a final concentration of 0.3 M and incubated at 37°C for 20 min. 28 µL 5 M ammonium acetate was added to the bisulfite-treated DNA. Finally, the bisulfite-treated DNA was precipitated with 3 volumes of ethanol. DNA was dissolved in 50 µL of water. Semi-nested PCR was performed, and the PCR products were cloned into pGM-T vector (Tiangen) and sequenced. The primers used for SlWEE1 were WEE1-577F (5′-TGTtATTTGGTTGTGAATAAGTTAT-3′), WEE1-457F (5′-TTATTATGATATGAGTGGTGGAA-3′) and WEE1-147R (5′-AAAATTTCATATATAACAAACTTATTAC-3′), with lowercase letters representing C-to-T or G-to-A substitution in SlWEE1, respectively.

Supporting Information

Figure S1.

Real-time RT-PCR analysis of mRNA levels of genes regulating cell division using SlGAPDH as the reference gene. Real-time RT-PCR analysis of mRNA levels of genes regulating cell division in 7-DPA fruits from WT Ailsa Craig (WT) and 35S::DDB1F T2 transgenic lines (a, b, c). p and n represent plants with or without transgene, respectively. Each bar represents three repetitions from each RNA sample (derived from pools of at least three fruits per plant). Error bars representing standard errors are shown in each case. Statistical analysis was performed using Student's t-test (*P<0.05, **P<0.01).

https://doi.org/10.1371/journal.pone.0042621.s001

(TIF)

Figure S2.

Semi-quantitative RT-PCR analysis of mRNA levels of SlWEE1, SlCDKB, SlCYCB2 and SlCYCA genes at different fruit development stages. Semi-quantitative RT-PCR analysis of mRNA levels of SlWEE1 gene and cell cycle controlling genes in 7-DPA, 21-DPA and 35-DPA fruits from WT Ailsa Craig (WT) and 35S::DDB1F T2 transgenic lines (a, b, c). p and n represent plants with or without transgene, respectively. Since the expression level of each gene was reduced with the fruit development, we used different PCR cycles at each stage. 7-DPA: 25 cycles; 21-DPA: 28 cycles; 35-DPA: 32 cycles.

https://doi.org/10.1371/journal.pone.0042621.s002

(TIF)

Author Contributions

Conceived and designed the experiments: Y. Liu FX JG. Performed the experiments: JL XT LG YG Y. Li SH XS MM HZ XT LZ XN. Analyzed the data: JL XT LG YG. Contributed reagents/materials/analysis tools: Y. Liu FX JG. Wrote the paper: Y. Liu FX.

References

  1. 1. Chu G, Chang E (1988) Xeroderma pigmentosum group E cells lack a nuclear factor that binds to damaged DNA. Science 242: 564–567.
  2. 2. Jackson S, Xiong Y (2009) CRL4s: the CUL4-RING E3 ubiquitin ligases. Trends Biochem Sci 34: 562–570.
  3. 3. Petroski MD, Deshaies RJ (2005) Function and regulation of cullin-RING ubiquitin ligases. Nat Rev Mol Cell Biol 6: 9–20.
  4. 4. O'Connell BC, Harper JW (2007) Ubiquitin proteasome system (UPS): what can chromatin do for you? Curr Opin Cell Biol 19: 206–214.
  5. 5. Chen H, Huang X, Gusmaroli G, Terzaghi W, Lau OS, et al. (2010) Arabidopsis CULLIN4-damaged DNA binding protein 1 interacts with CONSTITUTIVELY PHOTOMORPHOGENIC1-SUPPRESSOR OF PHYA complexes to regulate photomorphogenesis and flowering time. Plant Cell 22: 108–123.
  6. 6. Lee JH, Yoon HJ, Terzaghi W, Martinez C, Dai M, et al. (2010) DWA1 and DWA2, two Arabidopsis DWD protein components of CUL4-based E3 ligases, act together as negative regulators in ABA signal transduction. Plant Cell 22: 1716–1732.
  7. 7. Pazhouhandeh M, Molinier J, Berr A, Genschik P (2011) MSI4/FVE interacts with CUL4-DDB1 and a PRC2-like complex to control epigenetic regulation of flowering time in Arabidopsis. Proc Natl Acad Sci USA 108: 3430–3435.
  8. 8. Dumbliauskas E, Lechner E, Jaciubek M, Berr A, Pazhouhandeh M, et al. (2011) The Arabidopsis CUL4–DDB1 complex interacts with MSI1 and is required to maintain MEDEA parental imprinting. EMBO J 30: 731–743.
  9. 9. Biedermann S, Hellmann H (2010) The DDB1a interacting proteins ATCSA-1 and DDB2 are critical factors for UV-B tolerance and genomic integrity in Arabidopsis thaliana. Plant J 62: 404–415.
  10. 10. Castells E, Molinier J, Benvenuto G, Bourbousse C, Zabulon G, et al. (2011) The conserved factor DE-ETIOLATED 1 cooperates with CUL4-DDB1DDB2 to maintain genome integrity upon UV stress. EMBO J 30: 1162–1172.
  11. 11. Chinnusamy V, Zhu JK (2009) Epigenetic regulation of stress responses in plant. Curr Opin Plant Biol 12: 133–139.
  12. 12. Liu C, Lu F, Cui X, Cao X (2010) Histone methylation in higher plants. Annu Rev Plant Biol 61: 395–420.
  13. 13. Feng S, Jacobsen SE, Reik W (2010) Epigenetic reprogramming in plant and animal development. Science 330: 622–627.
  14. 14. Hauser MT, Aufsatz W, Jonak C, Luschnig C (2011) Transgenerational epigenetic inheritance in plants. Biochim Biophys Acta 1809: 459–468.
  15. 15. Molinier J, Ries G, Zipfel C, Hohn B (2006) Transgeneration memory of stress in plants. Nature 442: 1046–1049.
  16. 16. Higa LA, Wu M, Ye T, Kobayashi R, Sun H, et al. (2006) CUL4-DDB1 ubiquitin ligase interacts with multiple WD40-repeat proteins and regulates histone methylation. Nat Cell Biol 8: 1277–1283.
  17. 17. Kotake Y, Zeng Y, Xiong Y (2009) DDB1-CUL4 and MLL1 mediate oncogene-induced p16INK4a activation. Cancer Res 69: 1809–1814.
  18. 18. Oda H, Hübner MR, Beck DB, Vermeulen M, Hurwitz J, et al. (2010) Regulation of the histone H4 monomethylase PR-Set7 by CRL4Cdt2-mediated PCNA-dependent degradation during DNA damage. Mol Cell 40: 364–376.
  19. 19. Tardat M, Brustel J, Kirsh O, Lefevbre C, Callanan M, et al. (2010) The histone H4 Lys 20 methyltransferase PR-Set7 regulates replication origins in mammalian cells. Nat Cell Biol 12: 1086–1093.
  20. 20. Gu X, Jiang D, Yang W, Jacob Y, Michaels SD, et al. (2011) Arabidopsis homologs of Retinoblastoma-Associated Protein 46/48 associate with a histone deacetylase to act redundantly in chromatin silencing. PLoS Genet 7: e1002366.
  21. 21. Zhao Y, Shen Y, Yang S, Wang J, Hu Q, et al. (2010) Ubiquitin ligase components Cullin4 and DDB1 are essential for DNA methylation in Neurospora crassa. J Biol Chem 285: 4355–4365.
  22. 22. Xu H, Wang J, Hu Q, Quan Y, Chen H, et al. (2010) DCAF26, an adaptor protein of Cul4-based E3, is essential for DNA methylation in Neurospora crassa. PLoS Genet 6: e1001132.
  23. 23. Lewis ZA, Adhvaryu KK, Honda S, Shiver AL, Knip M, et al. (2010) DNA methylation and normal chromosome behavior in Neurospora depend on five components of a histone methyltransferase complex, DCDC. PLoS Genet 6: e1001196.
  24. 24. Yen HC, Shelton BA, Howard LR, Lee S, Vrebalov J, et al. (1997) The tomato high-pigment (hp) locus maps to chromosome 2 and influences plastome copy number and fruit quality. Theor Appl Genet 95: 1069–1079.
  25. 25. Bino RJ, Ric de Vos CH, Lieberman M, Hall RD, Bovy A, et al. (2005) The light-hyperresponsive high pigment-2dg mutation of tomato: alterations in the fruit metabolome. New Phytol 166: 427–438.
  26. 26. Kolotilin I, Koltai H, Tadmor Y, Bar-Or C, Reuveni M, et al. (2007) Transcriptional profiling of high pigment-2dg tomato mutant links early fruit plastid biogenesis with its overproduction of phytonutrients. Plant Physiol 145: 389–401.
  27. 27. Azari R, Reuveni M, Evenor D, Nahon S, Shlomo H, et al. (2010) Overexpression of UV-DAMAGED DNA BINDING PROTEIN 1 links plant development and phytonutrient accumulation in high pigment-1 tomato. J Exp Bot 61: 3627–3637.
  28. 28. Liu Y, Roof S, Ye Z, Barry C, van Tuinen A, et al. (2004) Manipulation of light signal transduction as a means of modifying fruit nutritional quality in tomato. Proc Natl Acad Sci USA 101: 9897–9902.
  29. 29. Lieberman M, Segev O, Gilboa N, Lalazar A, Levin I (2004) The tomato homolog of the gene encoding UV DAMAGED DNA BINDING protein 1 (DDB1) underlined as the gene that causes the high pigment-1 mutant phenotype. Theor Appl Genet 108: 1574–1581.
  30. 30. Wang S, Liu J, Feng Y, Niu X, Giovannoni J, et al. (2008) Altered plastid levels and potential for improved fruit nutrient content by downregulation of the tomato DDB1-interacting protein CUL4. Plant J 55: 89–103.
  31. 31. Liu J, Li H, Miao M, Tang X, Giovannoni J, et al. (2011) The tomato UV-damaged DNA-binding protein-1 (DDB1) is implicated in pathogenesis-related (PR) gene expression and resistance to Agrobacterium tumefaciens. Mol plant pathol DOI: 10.1111/J.1364-3703.2011.00735.X.
  32. 32. Blencowe BJ (2006) Alternative splicing: new insights from global analyses. Cell 126: 37–47.
  33. 33. Li Y, Zheng L, Corke F, Smith C, Bevan MW (2008) Control of final seed and organ size by the DA1 gene family in Arabidopsis thaliana. Genes Dev 22: 1331–1336.
  34. 34. Vanstraelen M, Baloban M, Da Ines O, Cultrone A, Lammens T, et al. (2009) APC/CCCS52A complexes control meristem maintenance in the Arabidopsis root. Proc Natl Acad Sci USA 106: 11806–11811.
  35. 35. Inzé D, De Veylder L (2006) Cell cycle regulation in plant development. Annu Rev Genet 40: 77–105.
  36. 36. Gonzalez N, Gévaudant F, Hernould M, Chevalier C, Mouras A (2007) The cell cycle-associated protein kinase WEE1 regulates cell size in relation to endoreduplication in developing tomato fruit. Plant J 51: 642–655.
  37. 37. Mathieu-Rivet E, Gévaudant F, Sicard A, Salar S, Do PT, et al. (2010) Functional analysis of the anaphase promoting complex activator CCS52A highlights the crucial role of endo-reduplication for fruit growth in tomato. Plant J 62: 727–741.
  38. 38. Abbas T, Shibata E, Park J, Jha S, Karnani N, et al. (2010) CRL4(Cdt2) regulates cell proliferation and histone gene expression by targeting PR-Set7/Set8 for degradation. Mol Cell 40: 9–21.
  39. 39. Centore RC, Havens CG, Manning AL, Li JM, Flynn RL, et al. (2010) CRL4(Cdt2)-mediated destruction of the histone methyltransferase Set8 prevents premature chromatin compaction in S phase. Mol Cell 40: 22–23.
  40. 40. Fuks F (2005) DNA methylation and histone modifications: teaming up to silence genes. Curr Opin Genet Dev 15: 490–495.
  41. 41. König H, Ponta H, Herrlich P (1998) Coupling of signal transduction to alternative pre-mRNA splicing by a composite splice regulator. EMBO J 17: 2904–2913.
  42. 42. Azari R, Tadmor Y, Meir A, Reuveni M, Evenor D, et al. (2010) Light signaling genes and their manipulation towards modulation of phytonutrient content in tomato fruits. Biotechnol Adv 28: 108–118.
  43. 43. Caspi N, Levin I, Chamovitz DA, Reuveni M (2008) A mutation in the tomato DDB1 gene affects cell and chloroplast compartment size and CDT1 transcript. Plant Signal Behav 3: 641–649.
  44. 44. Sun Y, Dilkes BP, Zhang C, Dante RA, Carneiro NP, et al. (1999) Characterization of maize (Zea mays L.) Wee1 and its activity in developig endosperm. Proc Natl Acad Sci USA 96: 4180–4185.
  45. 45. Gonzalez N, Hernould M, Delmas F, Gévaudant F, Duffe P, et al. (2004) Molecular characterization of a WEE1 gene homologue in tomato (Lycopersicon esculentum Mill.). Plant Mol Biol 56: 849–861.
  46. 46. De Schutter K, Joubès J, Cools T, Verkest A, Corellou F, et al. (2007) Arabidopsis WEE1 kinase controls cell cycle arrest in response to activation of the DNA integrity checkpoint. Plant Cell 19: 211–225.
  47. 47. Fang J, Feng Q, Ketel CS, Wang H, Cao R, et al. (2002) Purification and functional characterization of SET8, a nucleosomal histone H4-lysine 20-specific methyltransferase. Curr Biol 12: 1086–1099.
  48. 48. Nishioka K, Rice JC, Sarma K, Erdjument-Bromage H, Werner J, et al. (2002) PR-Set7 is a nucleosome-specific methyltransferase that modifies lysine 20 of histone H4 and is associated with silent chromatin. Mol Cell 9: 1201–1213.
  49. 49. Sambrook J, Russell WR (2001) Molecular Cloning: A Laboratory Manual, 3rd edn. Cold Spring Harbor, NY: Cold Spring Harbor Laboratory Press.
  50. 50. Fillatti JJ, Kiser J, Rose R, Comai L (1987) Efficient transfer of a glyphosate tolerance gene into tomato using a binary Agrobacterium tumefaciens vector. BioTechnology 5: 726–730.
  51. 51. Cheniclet C, Rong WY, Causse M, Frangne N, Bolling L, et al. (2005) Cell expansion and endoreduplication show a large genetic variability in pericarp and contribute strongly to tomato fruit frowth. Plant Physiol 139: 1984–1994.
  52. 52. Ding Y, Wang X, Su L, Zhai J, Cao S, et al. (2007) SDG714, a histone H3K9 methyltransferase, is involved in Tos17 DNA methylation and transposition in rice. Plant Cell 19: 9–22.