1 Department of Plant Systems Biology, Flanders Interuniversity Institute for Biotechnology (VIB), Ghent University, B-9052 Gent, Belgium
2 Laboratory of Plant Physiology, University of Liege, B-4000 Liège, Belgium
3 VIB MicroArray Facility, Gasthuisberg, Onderzoek en Navorsing, B-3000 Leuven, Belgium
* Author for correspondence (e-mail: dirk.inze{at}psb.ugent.be)
Accepted 19 June 2003
![]() |
Summary |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
Key words: Arabidopsis thaliana, Cell cycle, E2F, Microarray, Nitrogen assimilation
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
The mechanism of DNA replication seems to be conserved between mammals and plants, because E2F and DP genes have been isolated from different plant species, including wheat, tobacco, carrot, Arabidopsis and rice (Ramírez-Parra et al., 1999; Sekine et al., 1999
; Albani et al., 2000
; Magyar et al., 2000
; Ramirez-Parra and Gutierrez, 2000
; Kosugi and Ohashi, 2002a
). In the Arabidopsis genome there are three E2F (E2Fa, E2Fb and E2Fc) and two DP (DPa and DPb) genes (Vandepoele et al., 2002
). Recently, we have analyzed the phenotypes of plants co-overexpressing the E2Fa-DPa genes (De Veylder et al., 2002
). Transgenic plants were smaller than control plants, had curled leaves and cotyledons, and were arrested in growth at an early stage of development. Microscopic analysis revealed that E2Fa-DPa-overproducing cells underwent ectopic cell division or endoreduplication, depending on the cell type. Whereas extra cell divisions resulted in cells smaller than those seen in the same tissues of control plants, supplementary endoreduplication caused the formation of giant nuclei. By using reverse transcription (RT)PCR, we demonstrated that the expression levels of genes involved in DNA replication (CDC6, ORC1, MCM and DNA pol
) were strongly up-regulated (De Veylder et al., 2002
).
Physiologically important targets of the mammalian E2F-DP transcription factors have been identified by microarray hybridization experiments, chromatin immunoprecipitations and computer-assisted prediction (Ishida et al., 2001; Kel et al., 2001
; Müller et al., 2001
; Weinmann et al., 2001
; Ren et al., 2002
). E2F-DP-responsive genes can be found among genes involved in cell division, DNA repair and replication, mitotic progression, apoptosis and differentiation. Although little is known about the plant E2F-DP target genes, a database search has been published recently, in which the Arabidopsis genome was screened for genes harboring the TTTCCCGCC cis-acting element in their promoter (Ramirez-Parra et al., 2003
). However, it is still unclear whether this specific cis-acting element is the only one recognized by the plant E2F-DP complexes, or whether the presence of the TTTCCCGCC element is sufficient to mark a gene as a true E2F-DP target gene. In order to identify the functional classes of genes regulated by E2Fa-DPa and to understand the nature of the phenotype of the E2Fa-DPa-overexpressing plants, we designed a microarray experiment that compared the transcript levels of 4,571 genes of wild-type and transgenic lines. We found distinct classes of genes that were up- or down-regulated in the E2Fa-DPa plants. Promoter analysis allowed us to distinguish among the downstream expressed genes, the genes that were putatively under direct control of E2Fa-DPa. Furthermore, we found that the increased expression levels of E2Fa-DPa have a large impact on the expression levels of genes involved in nitrogen assimilation and metabolism.
![]() |
Materials and Methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
Construction of microarrays
The Arabidopsis thaliana (L.) Heynh. microarray consisted of 4,608 cDNA fragments spotted in duplicate, distant from each other, on Type V silane-coated slides (Amersham Biosciences, Little Chalfont, UK). The clone set included 4,571 Arabidopsis cDNAs from the unigene clone collection Arabidopsis Gem I (Incyte Genomics, Palo Alto, CA). The functional annotation of the genes related to the spotted cDNAs was retrieved by BLASTN against genomic sequences. To facilitate the analysis, a collection of genomic sequences was built each bearing only one gene. In each of these sequences, the upstream intergenic sequence was followed by the exon-intron structure of the gene and the downstream intergenic sequence, or, in other words, the whole genomic sequence between start and stop codons from neighboring protein-encoding genes. From the BLASTN output, the best hits were extracted and submitted to a BLASTX search against protein databases. From this analysis, the set of 4,571 cDNAs appeared to constitute 4,390 unique clones. To obtain more detailed information concerning the potential function of the genes, protein domains were searched using ProDom. The complete set can be found at http://www.psb.ugent.be/E2F/. The cDNA inserts were amplified by PCR with M13 primers, purified with MultiScreen- PCR plate (Millipore, Bedford, MA), and arrayed on slides using a Generation III printer (Amersham Biosciences). Slides were blocked in 3.5% SSC (1xSSC, 150 mM NaCl, 15 mM sodium citrate, pH 7.0), 0.2% sodium dodecyl sulfate (SDS), 1% bovine serum albumin for 10 minutes at 60°C.
RNA amplification and labeling
Antisense RNA was amplified with a modified protocol of in vitro transcription (Puskás et al., 2002). For the first-strand cDNA synthesis, 5 µg of total RNA was mixed with 2 µg of a HPLC-purified anchored oligo(dT) + T7 promoter (5'-GGCCAGTGAATTGTAATACGACTCACTATAGGGAGGCGG-T24(ACG)-3') (Eurogentec, Seraing, Belgium), 40 units of RNAseOUT (Invitrogen, Gaithersburg, MD) and 0.9 M D(+)trehalose (Sigma-Aldrich, St. Louis, MO) in a total volume of 11 µl and heated to 75°C for 5 minutes. To this mixture, 4 µl 5xfirst-strand buffer (Invitrogen), 2 µl 0.1 M dithiothreitol, 1 µl 10 mM dNTP mix, 1 µl 1.7 M D(+)trehalose (Sigma-Aldrich), and 1 µl SuperScript II (Invitrogen) were added to 20 µl final volume. The sample was incubated in a UnoII thermocycler (Whatman Biometra, Göttingen, Germany) at 37°C for 5 minutes, at 45°C for 10 minutes, 10 cycles at 60°C for 2 minutes and at 55°C for 2 minutes. To the first-strand reaction mix, 103.8 µl water, 33.4 µl 5xsecond-strand synthesis buffer (Invitrogen), 3.4 µl 10 mM dNTP mix, 1 µl of 10 U/µl DNA ligase (Invitrogen), 4 µl 10 U/µl DNA Polymerase I (Invitrogen), and 1 µl 2 U/µl RNAse H (Invitrogen) were added and incubated at 16°C for 2 hours. The synthesized double-stranded cDNA was purified with Qiaquick (Qiagen, Hilden, Germany). Antisense RNA was synthesized by AmpliScribe T7 high-yield transcription kit (Epicentre Technologies, Madison, WI) in a total volume of 20 µl according to the manufacturer's instructions. The RNA was purified with the RNeasy purification kit (Qiagen). From this RNA, 5 µg was labeled by reverse transcription using random nonamer primers (Genset, Paris, France), 0.1 mM d(G/T/A)TPs, 0.05 mM dCTP (Amersham Biosciences), 0.05 mM Cy3-dCTP or Cy5-dCTP (Amersham Biosciences), 1xfirst-strand buffer, 10 mM dithiothreitol, and 200 U of SuperScript II (Invitrogen) in 20 µl total volume. The RNA and primers were denatured at 75°C for 5 minutes and cooled on ice before the remaining reaction components were added. After 2 hours incubation at 42°C, mRNA was hydrolyzed in 250 mM NaOH for 15 minutes at 37°C. The sample was neutralized with 10 µl of 2 M 3-(N-morpholino)propanesulfonic acid and purified with Qiaquick (Qiagen).
Array hybridization and post-hybridization processes
The probes were resuspended in 30 µl hybridization solution (50% formamide, 5xSSC, 0.1% SDS, 100 µg/ml salmon sperm DNA) and prehybridized with 1 µl poly(dT) (1 mg/ml) at 42°C for 30 minutes to block hybridization on the polyA/T tails of the cDNA on the arrays. Mouse COT DNA (1 mg/ml) (Invitrogen) was added to the mixture and placed on the array under a glass coverslip. Slides were incubated for 18 hours at 42°C in a humid hybridization cabinet (Amersham Biosciences). Post-hybridization washing was performed for 10 minutes at 56°C in 1xSSC, 0.1% SDS, twice for 10 minutes at 56°C in 0.1xSSC, 0.1% SDS, and for 2 minutes at 37°C in 0.1xSSC.
Scanning and data analysis
Arrays were scanned at 532 nm and 635 nm using a Generation III scanner (Amersham Biosciences). Image analysis was performed with ArrayVision (Imaging Research Inc, St. Catharines, Ontario, Canada). Spot intensities were measured as artifact-removed total intensities (ARVol) without correction for background. We first addressed withinslide normalization by plotting for each single slide a `MA-plot' (Yang et al., 2002), where M=log2 (R/G) and
. Dye intensity differences were corrected with the `LOWESS' normalization. Subsequently, between-slide normalization and differentially expressed gene identification between the two genotypes were performed by sequential analysis of variances (ANOVAs), according to the method of Wolfinger et al. (Wolfinger et al., 2001
). (i) The base- 2 logarithm of the `LOWESS'-transformed measurements for all 73,136 spots (yiklm) were subjected to a normalization model yiklm=µ+Ak +AkDlRm+
iklm, where µ is the sample mean, Ak the effect of the kth array (k=1-4), AkDlRm the channel effect (AD) for the mth replication (m=2; left and right) of the total collection of i (i=1,..., 4571) cDNA fragments, and
iklm the stochastic error. (ii) The residuals from this model were subjected to 4,571 gene-specific models rijkl=µ+GiAk+GiDl+GiCj+
ijkl, where GiAk is the spot effect, GiDl the gene-specific dye effect, GiCj the signal intensity for genes that can specifically be attributed to the genotypes (effect of interest), and
ijkl the stochastic error. All effects were assumed to be fixed, except for
iklm and
ijkl. T-tests for differences between the GiCj effects were performed, all based on n1+n2-6 degrees of freedom, where n1 and n2 correspond to the number of wild-type and E2Fa-DPa hybridizations, respectively. Bonferroni adjustment for the 4,571 tests to assure an experiment-wise false positive rate of 0.05 results in a P-value cut-off of 1e-5.0, which is certainly too conservative. Thus, no further adjustments for multiple testing were done. Therefore, we chose to set the P-value cut-off arbitrarily at the 0.05 level. We used Genstat for both the normalization and gene model fits.
RT-mediated PCR analysis
RNA was isolated from plants 8 days after sowing with the Trizol reagent (Amersham Biosciences). First-strand cDNA was prepared from 3 µg of total RNA with the Superscript RT II kit (Invitrogen) and oligo(dT)18 according to the manufacturer's instructions. A 0.25 µl aliquot of the total RT reaction volume (20 µl) was used as a template in a semi-quantitative RT-mediated PCR amplification, ensuring that the amount of amplified product remained in linear proportion to the initial template present in the reaction. From the PCR reaction, 10 µl was separated on a 0.8% agarose gel and transferred onto Hybond N+ membranes (Amersham Biosciences) that were hybridized at 65°C with fluorescein-labeled probes (Gene Images random prime module; Amersham Biosciences). The hybridized bands were detected with the CDP Star detection module (Amersham Biosciences). Primers used were 5'-AAAAAGCAGGCTGTGTCGTACGATCTTCTCCCGG-3' and 5'-AGAAAGCTGGGTCATGTGATAGGAGAACCAGCG-3' for E2Fa, 5'-ATAGAATTCGCTTACATTTTGAAACTGATG-3' and 5'-ATAGTCGACTCAGCGAGTATCAATGGATCC-3' for DPa, 5'-CAGATCTTGTTAACCTTGACATCTCAG-3' and 5'-GGGTCAAAAGATACAACCACACCAG-3' for glutamine synthetase (GS), 5'-GGTTTACGAGCTACATGGCCC-3' and 5'-GAGCAATCCGTTCAGCCTCC-3' for glutamate synthase (GOGAT), 5'-GCGTTTGACCACTCTTGGAGAC-3' and 5'-GAACGCCATTGAGAAAGTCCGC-3' for histone acetylase HAT B, 5'-GTTACCGGCTCGACTTGAAGATC-3' and 5'-GAATCGGAGGGAAAGTCTGACG-3' for LOB domain protein 41, 5'-GTGTGGTTTCCAAGCTTTCCTACG-3' and 5'-GGTGAAGGGACTAGCCTTGTGG-3' for isocitrate lyase, 5'-GGGATCAATCCTCAGGAGAAGG-3' and 5'-CCGTCCATCTTTATTAGCGGCATG-3' for nitrite reductase (NiR), and 5'-TTACCGAGGCTCCTCTTAACCC-3' and 5'-ACCACCGATCCAGACACTGTAC-3' for actin 2 (ACT2).
Promoter analysis
The intergenic sequence corresponding to the promoter area of each gene spotted on the microarray was deduced from genomic sequences. From these intergenic sequences, up to 500 bp upstream of the ATG start codon were extracted and subjected to motif searches to retrieve potential E2F elements. Of the 4,571 expressed sequence tags (ESTs) spotted on the microarray, we could retrieve the genomic sequence of 4,390. This difference is due to the presence of duplicate genes and mitochondrial or chloroplast DNA on the microarray. Both the position and frequency of occurrence were determined with the publicly available MatInspector (version 2.2) by using matrices extracted from PlantCARE and matrices made especially for this particular analysis (Lescot et al., 2002). The relevance of each motif was evaluated against a background consisting of all the sequences from the dataset by using the Fisher exact test.
![]() |
Results and Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
Fluorescence levels were analyzed to establish whether the expression level of each gene varied according to the overexpression of the E2Fa-DPa transcription factor. Two sequential ANOVA models were used, as proposed by Wolfinger et al. (Wolfinger et al., 2001). First, the model called `normalization model' accounts for experiment-wise systematic effects, such as array and channel effects, which could bias inferences made on the data from the individual genes. The residuals from this model represent normalized values and are the input data for the second model, called the `gene' model. The gene models are fitted separately to the normalized data from each gene (see Materials and Methods). In this procedure, normalized expression levels rather than ratios are used as units.
For each of the 4,571 genes on the arrays the genotype-specific signal intensity was determined and t-tested for significant differences (P<0.05). Fig. 1 presents the P values obtained (as the negative log10 of the P value) against the magnitude of the effect (log2 of estimated fold change). This so-called volcano plot illustrates the substantial difference of significance testing as opposed to cut-offs strictly based on the fold change. The two vertical reference lines indicate a twofold cut-off for either repression or induction, whereas the horizontal reference line refers to the P-value cut-off at 0.05. These reference lines divide the plot into six meaningful sectors. The 3,126 genes in the lower middle sector have low significance and low fold change, and both methods are in agreement that the corresponding changes are not significant. The 188 genes in the upper left and right sectors have high significance (P<0.05) and high fold change (2); 84 of these genes show a significant two-or-more-fold induction of expression, whereas the remaining 104 genes show a significant two-or-more-fold repression of expression in the E2Fa-DPa plants. The identity of these genes was confirmed by sequencing, and the induction of a random set of selected genes was confirmed by RT-PCR analysis (Fig. 2). Finally, the 1,257 genes in the upper middle sector represent significant (P<0.05) up- or down-regulated genes, but with a low (
2) fold change. The full dataset of genes can be viewed at http://www.psb.ugent.be/E2F/.
|
|
DNA replication and cell cycle genes
Genes up- or down-regulated in the E2Fa-DPa transgenic plants can be classified into clear groups according to their function (Tables 1 and 2). Among the genes that are twofold or more up-regulated, 14 belong to the class of DNA replication and modification, correlating with the observation that E2Fa-DPa-overexpressing plants undergo extensive endoreduplication. Most of these genes have previously been shown to be up-regulated by E2F-DP overexpression in mammalian cells, including a putative thymidine kinase, replication factor c, adenosylhomocysteinase, DNA (cytosine-5)-methyltransferase, and histone genes (Ishida et al., 2001; Müller et al., 2001
; Ren et al., 2002
). Other E2Fa-DPa-induced S phase genes include a linker histone protein, the topoisomerase 6 subunit A, and two subunits of the histone acetyltransferase HAT B complex, namely HAT B and Msi3. The HAT B complex is responsible for the specific diacetylation of newly synthesized histone H4 during nucleosome assembly on newly synthesized DNA (Lusser et al., 1999
).
|
|
In addition to the overexpressed E2Fa gene (90-fold more abundant in transgenic than in control plants), only one cell cycle gene (CDKB1;1) has a twofold or more change in expression level upon E2Fa-DPa overexpression. CDKB1;1 had already been predicted to be a candidate E2F-DP target by the presence of a consensus E2F-DP-binding site in its promoter (de Jager et al., 2001). Whereas CDKB1;1 activity is highest at the G2 to M transition, its transcript levels start to increase during S phase (Porceddu et al., 1999
; Menges and Murray, 2002
). Therefore, up-regulation of CDKB1;1 might be a mechanism linking DNA replication with the following mitosis. That other cell cycle genes modulated in the E2Fa-DPa plants are not detected can be explained by the lack of many important E2F-DP target genes on the microarray and the putative difficulty in identifying changes in expression levels of lowly expressed genes in microarray hybridizations.
Cell wall biogenesis genes
Four members of the xyloglucan endotransglucosylase (XET) gene family are found to be twofold or more up-regulated in the E2Fa-DPa plants, one of them identical to the previously described Meri-5 gene (Medford et al., 1991). XETs are enzymes that modify cell wall components and are presumed to play a role in altering size, shape and physical properties of plant cells. Reversal breakage of the xyloglucan tethers by XETs has been proposed as a mechanism for allowing cell wall loosening in turgor-driven cell expansion (Campbell and Braam, 1999
). However, there are several reasons for believing that E2Fa-DPa-induced XETs are not required for cell expansion. First, cells divide more frequently in the E2Fa-DPa plants, but the overall cell size is smaller in transgenic than in control plants; so, no overall increase in expansion rates is needed. Second, no induction is seen of genes with a known role in cell expansion, such as expansins. Therefore, the hydrolytic activity of the XETs might rather be required to incorporate the newly synthesized cell walls formed during cytokinesis into the existing cell wall structure. Alternatively, because XET activity has been shown to be involved in the postgerminative mobilization of xyloglucan storage reserves in Nasturtium cotyledons (Farkas et al., 1992
; Fanutti et al., 1993
), induction of XETs in E2Fa-DPa plants might be related to polysaccharide breakdown to serve the metabolic and energy needs that are required to synthesize new nucleotides (see below).
Interestingly, two XETs can be identified in the set of twofold-or-more down-regulated genes. These XETs are more related to each other than to the induced XET genes. This differential response of XETs toward the E2Fa-DPa-induced phenotypes suggests that plant XETs can be classified into at least two different functional classes.
Genes involved in metabolism and biogenesis
A relatively large number of genes involved in metabolism and biogenesis were found in both the up-regulated and down-regulated gene groups. Most remarkable is the induction of genes involved in nitrogen assimilation, such as nitrate reductase (NIA2), glutamine synthetase (GS), and glutamate synthase (GOGAT) (Fig. 3). Although not present on the microarray, the nitrite reductase (NiR) gene was found to be induced as well in the transgenic lines, as demonstrated by RT-mediated PCR analysis (Fig. 2). Nitrogen and nitrite reductase catalyze the first two steps in the nitrogen assimilation pathway, whereas GS and GOGAT are involved both in the primary assimilation of nitrogen and the reassimilation of free ammonium. This mechanism supplies the plant with all nitrogen needed for the biosynthesis of amino acids and other nitrogen-containing compounds.
|
There are other indications that the nitrogen metabolism is altered in the E2Fa-DPa plants; these include the modification of genes homologous to genes expressed during the formation of nitrogen-fixing nodules in Medicago sativa (MTN3 and a nodulin-like gene), and the down-regulation of genes involved in sulfur assimilation (two different genes encoding adenylylsulfate reductase [APR] and a putative adenine phosphosulfate kinase). Genes involved in sulfur assimilation have been shown before to be transcriptionally down-regulated during nitrogen deficiency (Koprikova et al., 2000).
The altered expression of genes involved in nitrogen assimilation and metabolism in the E2Fa-DPa transgenic plants might reflect the need for nitrogen for the nucleotide biosynthesis, because purine and pyrimidine bases are rich in nitrogen. If nitrogen assimilation were indeed stimulated by E2Fa-DPa overexpression, two requirement should be fulfilled. Firstly, there should be enough -ketoglutarate to act as an acceptor molecule for ammonium (Lancien et al., 2000
) and secondly, because assimilation of nitrogen is energy consuming, the rate of reductant production should be higher in the E2Fa-DPa transgenic than in the wild-type plants.
Our microarray data suggest that in the accumulation of -ketoglutarate in E2Fa-DPa-overexpressing plants is stimulated in different ways. First,
-ketoglutarate production is improved by increased photosynthetic activity, as indicated by the 4.7-fold up-regulation of the large subunit of ribulose-1,5-bisphosphate carboxylase/oxygenase (Fig. 3), with accumulation of glyceraldehyde-3-phosphate as a result. Glyceraldehyde-3-phosphate can be converted into fructose-1,6-bisphosphate by fructose bisphosphate aldolase. However, a sixfold down-regulation of the fructose bisphosphate aldolase gene rather suggests the conversion of glyceraldehyde-3-phosphate into pyruvate, which can be converted into
-ketoglutarate in the citrate cycle. The preferential conversion of glyceraldehyde-3-phosphate to pyruvate fits the increased need for amino acids rather than for sugars to drive nucleotide biosynthesis (Fig. 3).
A second source of -ketoglutarate can be provided by the glyoxylate cycle. In E2Fa-DPa-overproducing plants we observed a 3.1-fold increase in expression of isocitrate lyase, suggesting an increased lipid turnover. Isocitrate lyase activity cleaves isocitrate into glyoxylate and succinate (Fig. 3). Whereas the produced glyoxylate can be converted into glycine, which is also required for de novo nucleotide biosynthesis, succinate can be converted into
-ketoglutarate in the citrate cycle. A 2.3-fold decrease in the expression of the fumarase gene presumably stimulates the subsequent conversion of
-ketoglutarate to glutamate by triggering an accumulation of succinate and fumarate, which are also side products formed during de novo nucleotide biosynthesis (Fig. 3). Reductant in plants mainly originates from photosynthetic electron transport in leaves. Corresponding with the increased need for reductant, several components of the chloroplast electron transport chain and associated ATP-synthesizing apparatus, such as cytochrome B6, a photosystem II subunit, and the ATPase
subunit, are up-regulated in the E2Fa-DPa transgenic plants. Increased expression of the protochlorophyllide reductase precursor even indicates an increase in chlorophyll biosynthesis.
E2Fa-DPa plants may suffer from nitrogen starvation that has an impact on amino acid biosynthesis. Three different amino acid aminotransferases are down-regulated in the E2Fa-DPa plants. Shortage of nitrogen-rich amino acids is also evident from the reduced expression of genes encoding vegetative storage proteins (VSP1 and VSP2) and ERD10, a protein with a compositional bias toward glutamate (Kiyosue et al., 1994). Additional evidence for amino acid shortage comes from the down-regulation of a myrosinase-binding protein and the cytochrome P450 monooxygenase CYP83A1. Both proteins are involved in the biosynthesis of glucosinolates, nitrogen- and sulfur-containing products derived from amino acids (Wittstock and Halkier, 2002
).
Promoter analysis of E2Fa-DPa-regulated genes
The DNA-binding domains of the E2F and DP proteins are highly conserved between plants and mammals and, correspondingly, plant E2F-DP proteins have been shown by the technique of electrophoresis mobility shift assay to bind to the same canonical DNA-binding site as their mammalian counterparts (Albani et al., 2000; Ramirez-Parra and Gutierrez, 2000
; de Jager et al., 2001
). Furthermore, these E2F-binding sites regulate the expression of several plant genes involved in DNA synthesis (Kosugi and Ohashi, 2002a
; Chabouté et al., 2000
; Castellano et al., 2001
; Egelkrout et al., 2001
; Stevens et al., 2002
).
To distinguish between the putatively direct target genes of E2Fa-DPa and the secondarily induced genes, the first 500 bp upstream of the ATG start codon of the genes with 2-fold or higher change in expression were scanned for the presence of an E2F-like-binding site matching the (A/T)TT(G/C)(G/C)C(G/C)(G/C) sequence, which corresponds to all the different E2F-DP-binding motifs that have been described in plants. Of all the different permutations only the TTTCCCGC and TTTGGCGG elements were enriched significantly (P<0.01) in the set of E2Fa-DPaupregulated genes, suggesting these are the preferred binding site of the E2Fa-DPa complex (Table 3). Moreover, six out of eight target genes containing one or more of these elements belong to the group of genes involved in DNA replication and modification. The observation that not all genes that enclase this DNA sequence in their promoter are induced upon E2Fa-DPa overexpression suggests that the presence of the TTTCCCGC or TTTGGCGG motif is not the only element to make a gene responsive toward E2Fa-DPa, and that E2Fa-DPa may cooperate with other factors to activate transcription. Alternatively, the promoters of non-responsive genes might be shielded with other transcription factor complexes. A putative candidate is the E2Fc protein which, in analogy with the mammalian E2F6 protein, lacks a strong transactivation domain (del Pozo et al., 2002). Alternative candidates are the recently discovered DEL proteins, proven to bind as monomers to the canonical E2F-binding site (Kosugi and Ohashi, 2002b
; Mariconti et al., 2002
). Because of a lack of transcriptional activation domain, the DEL proteins are postulated to act as repressors of E2F-DP-regulated genes by competing for the same binding site.
|
It is not excluded that genes without an E2F-like-binding site are not directly activated by E2Fa-DPa. Chromatin immunoprecipitation experiments have shown that mammalian E2F factors can bind to promoters without a clear E2F recognition motif (Kiyosue et al., 1994), suggesting that E2FDP might recognize non-canonical binding sites, or might be recruited by promoters through the association of other factors. In this respect, the Chlorella vulgaris nitrate reductase gene, of which the Arabidopsis homologue was shown here to be induced by E2F-DPa, binds an E2F-DP complex, although a clear consensus binding site is lacking (Cannons and Shiflett, 2001
).
E2Fs can activate as well as repress promoter activity (Trimarchi and Lees, 2002). In the PCNA, MCM3 and RNR2 promoters, E2F sequences have been identified that act as a negative regulatory element during development (Chabouté et al., 2000
; Egelkrout et al., 2001
; Stevens et al., 2002
). In the set of down-regulated genes, no particular enrichment of a specific E2F sequence could be seen (Table 3). Therefore, the data suggest that the E2Fa-DPa complex works as a transcriptional activator and that other E2F-DP complexes are involved in E2F-mediated transcriptional repression.
Conclusions
Microarray analysis of E2Fa-DPa-overexpressing lines identified a cross-talking genetic network between DNA replication, nitrogen assimilation and photosynthesis. The putatively direct E2Fa-DPa target genes as identified by the presence of an E2F-DP-binding site, belong to the group of genes involved in DNA synthesis, whereas the secondarily induced genes are mainly linked to nitrogen assimilation. In a recently published microarray experiment in which the periodic expression of genes during the cell cycle was monitored, genes with a role in nitrogen assimilation (aspartate aminotransferase and a nitrate transporter) were found to be specifically expressed during the S phase (Menges et al., 2002). Because purine and pyrimidine bases are nitrogen rich, we postulate that induction of nitrogen assimilation genes during DNA synthesis in wild-type and E2Fa-DPa transgenic plants is required to supply enough nitrogen for nucleotide biosynthesis. However, in the EFa-DPa transgenic plants, increased nitrogen assimilation most probably does not meet all the nucleotide biosynthesis needs, as seen by the expression modulation of many genes involved in nitrogen and carbohydrate metabolism. The drain of nitrogen from essential biosynthetic pathways to the nucleotide biosynthesis pathway is expected to affect other aspects of plant metabolism, as can be seen from the reduced expression of vegetative storage protein genes and genes involved in amino acid biosynthesis. This altered metabolism might, at least in part, contribute to the growth arrest observed in E2Fa-DPa transgenic plants.
The exact regulatory pathways and factors controlling the nitrogen assimilation pathway in plants are still unknown. In addition to the genes involved in DNA replication and metabolism, our data contain a relatively large number of genes with unspecified function (Tables 1 and 2). For instance, a GATA zinc-finger-encoded gene with a still unknown function is found between the up-regulated regulatory genes. This gene might encode the ortholog of the Neurospora crassa nit-2 protein that has been shown to positively regulate expression of the nitrate reductase gene (Fu and Marzluf, 1990). Other regulatory genes modified in the E2Fa-DPa plants encode protein kinases and several putative receptor kinases. These genes might include some novel key regulatory components in the process of nitrogen assimilation or regulation of efficient nitrogen usage. It will be of great interest to analyze their role in nitrogen assimilation, metabolism and plant growth.
![]() |
Acknowledgments |
---|
![]() |
References |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
Albani, D., Mariconti, L., Ricagno, S., Pitto, L., Moroni, C., Helin, K. and Cella, R. (2000). DcE2F, a functional plant E2F-like transcriptional activator from Daucus carota. J. Biol. Chem. 275, 19258-19267.
Campbell, P. and Braam, J. (1999). Xyloglucan endotransglycosylases: diversity of genes, enzymes and potential wall-modifying functions. Trends Plant Sci. 4, 361-366.[CrossRef][Medline]
Cannons, A. C. and Shiflett, S. D. (2001). Transcriptional regulation of the nitrate reductase gene in Chlorella vulgaris: identification of regulatory elements controlling expression. Curr. Genet. 40, 128-135.[CrossRef][Medline]
Castellano, M. M., del Pozo, J. C., Ramirez-Parra, E., Brown, S. and Gutierrez, C. (2001). Expression and stability of Arabidopsis CDC6 are associated with endoreplication. Plant Cell 13, 2671-2686.
Chabouté, M.-E., Clément, B., Sekine, M., Philipps, G. and Chaubet-Gigot, N. (2000). Cell cycle regulation of the tobacco ribonucleotide reductase small subunit gene is mediated by E2F-like elements. Plant Cell 12, 1987-1999.
de Jager, S. M., Menges, M., Bauer, U.-M. and Murray, J. A. H. (2001). Arabidopsis E2F1 binds a sequence present in the promoter of S-phase-regulated gene AtCDC6 and is a member of a multigene family with differential activities. Plant Mol. Biol. 47, 555-568.[CrossRef][Medline]
De Veylder, L., Beeckman, T., Beemster, G. T. S., de Almeida Engler, J., Ormenese, S., Maes, S., Naudts, M., Van Der Schueren, E., Jacqmard, A., Engler, G. and Inzé, D. (2002). Control of proliferation, endoreduplication and differentiation by the Arabidopsis E2Fa/DPa transcription factor. EMBO J. 21, 1360-1368.
del Pozo, J. C., Boniotti, M. B. and Gutierrez, C. (2002). Arabidopsis E2Fc functions in cell division and is degraded by the ubiquitin-SCFAtSKP2 pathway in response to light. Plant Cell 14, 3057-3071.
Egelkrout, E. M., Robertson, D. and Hanley-Bowdoin, L. (2001). Proliferating cell nuclear antigen transcription is repressed through an E2F consensus element and activated by geminivirus infection in mature leaves. Plant Cell 13, 1437-1452.
Fanutti, C., Gidley, M. J. and Reid, J. S. G. (1993). Action of a pure xyloglucan endo-transglycosylase (formerly called xyloglucan-specific endo-(14)-ß-D-glucanase) from the cotyledons of germinated nasturtium seeds. Plant J. 3, 691-700.[CrossRef][Medline]
Farkas, V., Sulova, Z., Stratilova, E., Hanna, R. and Maclachlan, G. (1992). Cleavage of xyloglucan by nasturtium seed xyloglucanase and transglycosylation to xyloglucan subunit oligosaccharides. Arch. Biochem. Biophys. 298, 365-370.[Medline]
Fu, Y.-H. and Marzluf, G. A. (1990). nit-2, the major positive-acting nitrogen regulatory gene of Neurospora crassa, encodes a sequence-specific DNA-binding protein. Proc. Natl. Acad. Sci. USA 87, 5331-5335.[Abstract]
Ishida, S., Huang, E., Zuzan, H., Spang, R., Leone, G., West, M. and Nevins, J. R. (2001). Role for E2F in control of both DNA replication and mitotic functions as revealed from DNA microarray analysis. Mol. Cell. Biol. 21, 4684-4699.
Kel, A. E., Kel-Margoulis, O. V., Farnham, P. J., Bartley, S. M., Wingender, E. and Zhang, M. Q. (2001). Computer-assisted identification of cell cycle-related genes: new targets for E2F transcription factors. J. Mol. Biol. 309, 99-120.[CrossRef][Medline]
Kiyosue, T., Yamaguchi-Shinozaki, K. and Shinozaki, K. (1994). Cloning of cDNAs for genes that are early-responsive to dehydration stress (ERDs) in Arabidopsis thaliana L.: identification of three ERDs as HSP cognate genes. Plant Mol. Biol. 25, 791-798.[Medline]
Koprivova, A., Suter, M., Op den Camp, R., Brunold, C. and Kopriva, S. (2000). Regulation of sulfate assimilation by nitrogen in Arabidopsis. Plant Physiol. 122, 737-746.
Kosugi, S. and Ohashi, Y. (2002a). E2F sites that can interact with E2F proteins cloned from rice are required for meristematic tissue-specific expression of rice and tobacco proliferating cell nuclear antigen promoters. Plant J. 29, 45-59.[CrossRef][Medline]
Kosugi, S. and Ohashi, Y. (2002b). E2Ls, E2F-like repressors of Arabidopsis that bind to E2F sites in a monomeric form. J. Biol. Chem. 277, 16553-16558.
Lancien, M., Gadal, P. and Hodges, M. (2000). Enzyme redundancy and the importance of 2-oxoglutarate in higher plant ammonium assimilation. Plant Physiol. 123, 817-824.
Lescot, M., Déhais, P., Thijs, G., Marchal, K., Moreau, Y., Van de Peer, Y., Rouzé, P. and Rombauts, S. (2002). PlantCare: a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 30, 325-327.
Lusser, A., Eberharter, A., Loidl, A., Goralik-Schramel, M., Horngacher, M., Haas, H. and Loidl, P. (1999). Analysis of the histone acetyltransferase B complex of maize embryos. Nucleic Acids Res. 27, 4427-4435.
Magyar, Z., Atanassova, A., De Veylder, L., Rombauts, S. and Inzé, D. (2000). Characterization of two distinct DP-related genes from Arabidopsis thaliana. FEBS Lett. 486, 79-87.[CrossRef][Medline]
Mariconti, L., Pellegrini, B., Cantoni, R., Stevens, R., Bergounioux, C., Cella, R. and Albani, D. (2002). The E2F family of transcription factors from Arabidopsis thaliana. Novel and conserved components of the retinoblastoma/E2F pathway in plants. J. Biol. Chem. 277, 9911-9919.
Medford, J. I., Elmer, J. S. and Klee, H. J. (1991). Molecular cloning and characterization of genes expressed in shoot apical meristems. Plant Cell 3, 359-370.
Menges, M. and Murray, J. A. H. (2002). Synchronous Arabidopsis suspension cultures for analysis of cell-cycle gene activity. Plant J. 20, 203-212.
Menges, M., Hennig, L., Gruissem, W. and Murray, J. A. H. (2002). Cell cycle-regulated gene expression in Arabidopsis. J. Biol. Chem. 277, 41987-42002.
Müller, H., Bracken, A. P., Vernell, R., Moroni, M. C., Christians, F., Grassilli, E., Prosperini, E., Vigo, E., Oliner, J. D. and Helin, K. (2001). E2Fs regulate the expression of genes involved in differentiation, development, proliferation, and apoptosis. Genes Dev. 15, 267-285.
Porceddu, A., De Veylder, L., Hayles, J., Van Montagu, M., Inzé, D. and Mironov, V. (1999). Mutational analysis of two Arabidopsis thaliana cyclin-dependent kinases in fission yeast. FEBS Lett. 446, 182-188.[CrossRef][Medline]
Puskás, L. G., Zvara, A., Hackler, L. J. and Van Hummelen, P. (2002). RNA amplification results in reproducible microarray data with slight ratio bias. Biotechniques 32, 1330-1341.[Medline]
Ramirez-Parra, E. and Gutierrez, C. (2000). Characterization of wheat DP, a heterodimerization partner of the plant E2F transcription factor which stimulates E2F-DNA binding. FEBS Lett. 486, 73-78.[CrossRef][Medline]
Ramirez-Parra, E., Fründt, C. and Gutierrez, C. (2003). A genome-wide identification of E2F-regulated genes in Arabidopsis. Plant J. 33, 801-811.[CrossRef][Medline]
Ramírez-Parra, E., Xie, Q., Boniotti, M. B. and Gutierrez, C. (1999). The cloning of plant E2F, a retinoblastoma-binding protein, reveals unique and conserved features with animal G1/S regulators. Nucleic Acids Res. 27, 3527-3533.
Ren, B., Cam, H., Takahashi, Y., Volkert, T., Terragni, J., Young, R. A. and Dynlacht, B. D. (2002). E2F integrates cell cycle progression with DNA repair, replication, and G2/M checkpoints. Genes Dev. 16, 245-256.
Sekine, M., Ito, M., Uemukai, K., Maeda, Y., Nakagami, H. and Shinmyo, A. (1999). Isolation and characterization of the E2F-like gene in plants. FEBS Lett. 460, 117-122.[CrossRef][Medline]
Stevens, R., Mariconti, L., Rossignol, P., Perennes, C., Cella, R. and Bergounioux, C. (2002). Two E2F sites in the Arabidopsis MCM3 promoter have different roles in cell cycle activation and meristematic expression. J. Biol. Chem. 277, 32978-32984.
Trimarchi, J. M. and Lees, J. A. (2002). Sibling rivalry in the E2F family. Nature Rev. Mol. Cell. Biol. 3, 11-20.[CrossRef][Medline]
Valvekens, D., Van Montagu, M. and Van Lijsebettens, M. (1988). Agrobacterium tumefaciens-mediated transformation of Arabidopsis thaliana root explants by using kanamycin selection. Proc. Natl. Acad. Sci. USA 85, 5536-5540.[Abstract]
Vandepoele, K., Raes, J., De Veylder, L., Rouzé, P., Rombauts, S. and Inzé, D. (2002). Genome-wide analysis of core cell cycle genes in Arabidopsis. Plant Cell 14, 903-916.
Weinmann, A., Bartley, S. M., Zhang, T., Zhang, M. Q. and Farnham, P. J. (2001). Use of chromatin immunoprecipitation to clone novel E2F target promoters. Mol. Cell. Biol. 21, 6820-6832.
Wittstock, U. and Halkier, B. A. (2002). Glucosinolate research in the Arabidopsis era. Trends Plant Sci. 7, 263-270.[CrossRef][Medline]
Wolfinger, R. D., Gibson, G., Wolfinger, E. D., Bennett, L., Hamadeh, H., Bushel, P., Afshari, C. and Paules, R. S. (2001). Assessing gene significance from cDNA microarray expression data via mixed models. J. Comput. Biol. 8, 625-637.[CrossRef][Medline]
Yang, Y. H., Dudoit, S., Luu, P., Lin, D. M., Peng, V., Ngai, J. and Speed, T. P. (2002). Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucleic Acids Res. 30, e15.