Genetics, Vol. 177, 749-760, October 2007, Copyright © 2007
doi:10.1534/genetics.107.072702

This Article Abstract Full Text (PDF) Data Supplement All Versions of this Article: genetics.107.072702v1 177/2/749    most recent Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via Google Scholar Google Scholar Articles by Makarevitch, I. Articles by Springer, N. M. Search for Related Content PubMed PubMed Citation Articles by Makarevitch, I. Articles by Springer, N. M.

Natural Variation for Alleles Under Epigenetic Control by the Maize Chromomethylase Zmet2

Irina Makarevitch^*, Robert M. Stupar^*, A. Leonardo Iniguez, William J. Haun^*, W. Brad Barbazuk, Shawn M. Kaeppler and Nathan M. Springer^*,1

^* Department of Plant Biology, Microbial and Plant Genomics Institute, University of Minnesota, Saint Paul, Minnesota 55108, Department of Agronomy, University of Wisconsin, Madison, Wisconsin 53706 and Donald Danforth Plant Sciences Center, St. Louis, Missouri 63132

¹ Corresponding author: Department of Plant Biology, 250 Biological Sciences Center, 1445 Gortner Ave., Saint Paul, MN 55108.
E-mail: springer{at}umn.edu

	ABSTRACT

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION ACKNOWLEDGEMENTS LITERATURE CITED

 
The contribution of epigenetic alterations to natural variation for gene transcription levels remains unclear. In this study, we investigated the functional targets of the maize chromomethylase ZMET2 in multiple inbred lines to determine whether epigenetic changes conditioned by this chromomethylase are conserved or variable within the species. Gene expression microarrays were hybridized with RNA samples from the inbred lines B73 and Mo17 and from near-isogenic derivatives containing the loss-of-function allele zmet2-m1. A set of 126 genes that displayed statistically significant differential expression in zmet2 mutants relative to wild-type plants in at least one of the two genetic backgrounds was identified. Analysis of the transcript levels in both wild-type and mutant individuals revealed that only 10% of these genes were affected in zmet2 mutants in both B73 and Mo17 genetic backgrounds. Over 80% of the genes with expression patterns affected by zmet2 mutations display variation for gene expression between wild-type B73 and Mo17 plants. Further analysis was performed for 7 genes that were transcriptionally silent in wild-type B73, but expressed in B73 zmet2-m1, wild-type Mo17, and Mo17 zmet2-m1 lines. Mapping experiments confirmed that the expression differences in wild-type B73 relative to Mo17 inbreds for these genes were caused by cis-acting regulatory variation. Methylation-sensitive PCR and bisulfite sequencing demonstrated that for 5 of these genes the CpNpG methylation in the wild-type B73 genetic background was substantially decreased in the B73 zmet2-m1 mutant and in wild-type Mo17. A survey of eight maize inbreds reveals that each of these 5 genes exhibit transcriptionally silent and methylated states in some inbred lines and unmethylated, expressed states in other inbreds, providing evidence for natural variation in epigenetic states for some maize genes.

Epigenetic states can be quite stable and faithfully transmitted through mitosis and meiosis (CUBAS et al. 1999; CHANDLER and STAM 2004). Epigenetic variation can result in formation of epigenetic alleles with altered gene expression patterns that can mimic traditional mutations. There are several examples of epigenetic variation that result in phenotypic variation (BENDER and FINK 1995; CUBAS et al. 1999; MORGAN et al. 1999; CHANDLER et al. 2000; STOKES et al. 2002; MANNING et al. 2006). Epigenetic variation has been classified as obligatory, facilitated, and pure according to the dependence on sequence changes (RICHARDS 2006). Obligatory epigenetic variation requires a genetic change, such as a nearby transposon insertion, that conditions the allele to regulation by epigenetic mechanisms. Facilitated epigenetic variation requires a genetic change that conditions an allele to stochastically fall under epigenetic control, such as the intracisternal A-particle (IAP) insertion at some agouti alleles (MORGAN et al. 1999). Pure epigenetic variation has no genetic changes associated with the epigenetic state. In many cases, the lack of detailed sequence comparisons makes it difficult to ascertain whether a particular instance is an example of an obligatory, facilitated, or pure epigenetic variation. Studies in Arabidopsis have provided evidence for natural variation for epigenetic states (RANGWALA et al. 2006) and factors that effect epigenetic states (RIDDLE and RICHARDS 2005), suggesting that epigenetic states may be a target for natural or artificial selection (KALISZ and PURUGGANAN 2004).

DNA methylation is frequently associated with epigenetic alterations in higher plants. In plants, DNA methylation is most commonly found in CpG dinucleotide and CpNpG trinucleotide contexts, which are maintained by different enzymes (CHAN et al. 2005). DNA methylation patterns are created and maintained by different families of DNA methyltransferase enzymes, including CpG maintenance methyltransferases (FINNEGAN and DENNIS 1993; KANKEL et al. 2003; XIAO et al. 2006), CpNpG maintenance/de novo methyltransferases (HENIKOFF and COMAI 1998; BARTEE et al. 2001), and de novo methyltransferases (CAO et al. 2000; CAO and JACOBSEN 2002b). Additional factors, including a SWI/SNF2 factor (DDM1), a histone deacetylase (HDA6), a histone methyltransferase (KYP), and several components of the RNA interference silencing pathway, also play roles in maintaining proper DNA methylation levels (VONGS et al. 1993; JEDDELOH et al. 1999; AUFSATZ et al. 2002; JACKSON et al. 2002; MALAGNAC et al. 2002; CHAN et al. 2006). Although unique targets of specific DNA methyltransferases have been identified, there is evidence for locus-specific, partial-to-complete redundancy for the function of DNA methyltransferases in the creation and maintenance of DNA methylation patterns (CAO and JACOBSEN 2002a,b).

The plant-specific chromomethylase (CMT) enzymes contain a chromo-domain embedded within the methyltransferase motifs (HENIKOFF and COMAI 1998; GENGER et al. 1999; BARTEE et al. 2001). The chromomethylases are primarily responsible for methylation found in CpNpG contexts, although they may also play a role in de novo methylation (LINDROTH et al. 2001; PAPA et al. 2001; CAO and JACOBSEN 2002a,b; CAO et al. 2003). Although chromomethylases are important for silencing of endogenous sequences such as SUP and PAI loci (BARTEE et al. 2001; CAO and JACOBSEN 2002a), knockout mutations in CMT genes do not result in notable morphological phenotypes in Arabidopsis or maize (BARTEE et al. 2001; PAPA et al. 2001; CAO and JACOBSEN 2002a). In addition to controlling the methylation state and expression level for SUP and PAI genes, CMT3 has been shown to methylate transposable elements (LINDROTH et al. 2001). A series of genetic and biochemical studies have provided evidence for the targeting of CMT methylation through histone methylation (LINDROTH et al. 2001; JACKSON et al. 2002; MALAGNAC et al. 2002; EBBS and BENDER 2006) and short interfering RNA (siRNA) pathways (ZILBERMAN et al. 2003; CHAN et al. 2004, 2006).

In this study, we investigated the targets of the maize chromomethylase ZMET2 (PAPA et al. 2001). Mutations in the maize Zmet2 chromomethylase gene result in a 10–15% reduction in genomewide methylation levels, mainly affecting CpNpG methylation (PAPA et al. 2001). However, there is little evidence for any morphological phenotypes in homozygous mutant plants following multiple generations of introgression into several inbred lines (PAPA et al. 2001). Surprisingly, in this study we found that ZMET2 was involved in the maintenance of epigenetic states that displayed natural variation in different maize inbreds. The expression and methylation status of epigenetically regulated alleles conditioned by CpNpG methylation was found to vary between different maize inbred lines, suggesting the potential for this epigenetic variation to contribute to transcriptional variation and phenotypic variation among inbred lines.

	MATERIALS AND METHODS

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION ACKNOWLEDGEMENTS LITERATURE CITED

 
Plant materials and tissue collection:
The zmet2-m1 allele, a loss-of-function allele resulting from an insertion of a Mu transposon into the 18th exon of the gene (PAPA et al. 2001), was backcrossed into the Mo17, B73, and W22 r-g inbred lines for six generations and then self-pollinated for two generations to produce the homozygous near-isogenic lines of Mo17, B73, and W22 r-g that are homozygous for the zmet2-m1 allele. Biological replicates in the B73 background (six replicates) and Mo17 background (three replicates) were grown using standard greenhouse conditions; RNAs from these plants were used for microarray analyses. Each biological replicate was planted at a different time and represents the pooled aboveground tissue from eight 11-day-old seedlings (for full details on plant growth and tissue collection, see supporting online text at http://genetics.org/supplemental/). W22 r-g wild-type and zmet2-m1 plants, as well as plants of five other inbred lines (B14, B37, B84, Oh43, and Wf9) used for analysis of epigenetic variation in other inbreds, were grown under the same conditions. Plants segregating for the zmet2-m1 allele were genotyped with primers Zmet2F23, ZmetR23, and 9242 (see supplemental Table 1 at http://www.genetics.org/supplemental/ for all primer sequences) using the PCR conditions described below.

RNA isolation and microarray hybridization:
RNA isolation and Affymetrix (Santa Clara, CA) microarray hybridizations were performed as described (STUPAR and SPRINGER 2006) for six biological replicates of wild-type and mutant plants in the B73 genetic background and for three biological replicates of wild-type and mutant plants in the Mo17 genetic background. Spotted long oligonucleotide array hybridizations were performed for three biological replicates of the B73 and B73 zmet2-m1 genotypes using the same RNA stock used for Affymetrix hybridizations. For this experiment, dual channel (Cy3 and Cy5 per slide) studies were performed with complete dye swapping similar to previously described methods (GARDINER et al. 2005). Total RNA (1.5 µg) from each biological replicate was indirectly labeled using Ambion's (Austin, TX) Aminoallyl Message Amp II kit according to the manufacturer's recommendations with several modifications (see the supporting online text at http://www.genetics.org/supplemental/ for details of microarray hybridizations). The microarray data are available under GEO series accession no. GSE8188.

Microarray data analysis:
Affymetrix microarray data analysis was performed as described (STUPAR and SPRINGER 2006). Briefly, the GCOS software package v1.2 (Affymetrix) was used for signal acquisition and initial analysis. GeneSpring (Agilent Technologies, Palo Alto, CA) software was used for GC content robust multi-array analysis (GC-RMA) processing of the .cel files that involved normalization between the arrays and a subsequent per gene normalization of the resulting values (STUPAR and SPRINGER 2006). Genes differentially expressed in B73 zmet2-m1 relative to B73 or in Mo17 zmet2-m1 relative to Mo17 genotypes were identified by performing separate one-way ANOVA tests on the GC-RMA values using a parametric test with no assumption of equal variance for each genetic background. A Benjamin and Hochberg multiple testing correction was applied using a false-discovery-rate significance threshold of 0.25, such that 25% of the genes identified in a test are likely to be falsely identified. The genes that were identified in at least one of the comparisons were pooled together into a list of differentially expressed genes. The resulting list of the genes was further filtered on the basis of the criteria of expression level (i.e., the expression level must be >50 units for GC-RMA values in at least one of the genotypes) and significant t-test (P < 0.05 in B73 vs. B73 zmet2-m1, Mo17 vs. Mo17 zmet2-m1, or mutant vs. wild-type comparisons). The GeneSpring software was used to perform a hierarchical clustering analysis using a Pearson correlation method to create gene or condition trees based on specified gene lists, conditions, and genotypes.

For the long oligonucleotide spotted arrays, the data were extracted from each slide using GenePix Pro software (http://www.moleculardevices.com/pages/software/gn_genepix_pro.html). The analysis for differentially expressed genes was conducted using the limma package (SMYTH 2004). The complete dye swapping was accounted for by determining the correlation values between the same samples coupled to both dyes. The data were normalized using loess and quantile methods to account for the variation within arrays and between arrays, respectively. Once the data had been background corrected and normalized and the dye bias had been accounted for, the data were fitted into a linear model and differentially expressed genes were estimated with a false discovery rate of 0.05 (SMYTH 2004).

RT–PCR validation:
Primers were designed using Primer 3.0 software (ROZEN and SKALETSKY 2000) for a subset of the differentially expressed genes on the basis of corresponding tentative consensus contigs (http://www.tigr.org) and maize assembled genomic islands (FU et al. 2005) and are listed in supplemental Table SOM1 at http://www.genetics.org/. RNA samples used for cDNA synthesis were DNAse treated (Promega, Madison, WI) and reverse transcribed using Superscript III reverse transcriptase (Invitrogen, Carlsbad, CA) according to the manufacturer's instructions. PCR reactions were performed on genomic DNA or cDNA templates in a 15-µl total volume containing 50 ng of DNA, 2 pmol of each primer, 0.4 units of HotStar Taq polymerase (Eppendorf, Westbury, NY), 1.58 µl of 10x reaction buffer, and 0.1 µl of 50 mM dNTPs. Cycling conditions of the PCR reactions were as follows: 94° for 15 min, 35 cycles of 94° for 30 sec, 60° for 30 sec, 72° for 90 sec, followed by 72° for 5 min. Amplified products were separated in a 1% agarose TBE gel and visualized by ethidium bromide staining. The concentrations of cDNAs were calibrated and made approximately equal by amplifying a mez1 gene expressed at approximately the same level in all four genotypes analyzed.

Mapping of expression differences:
The genomic location responsible for expression differences for seven genes was mapped using a set of 51 recombinant inbred lines from the Syn4F7:9 IBM population (LEE et al. 2002). RNA isolation, cDNA synthesis, and RT–PCR reactions were performed as described above and each line was scored as "present" or "absent" for expression of the gene being amplified. These scores were mapped using "U Map It!" software (http://www.plantgenomics.iastate.edu/maize/).

Methylation assays:
Genomic DNA was extracted from 11-day-old seedlings for B73, B73 zmet2-m1, Mo17, and Mo17 zmet2-m1 genotypes using the CTAB method (SAGHAI-MAROOF et al. 1984). Two plants were sampled for each of the four genotypes. For methylation-sensitive PCR assays, 500 ng of genomic DNA was digested with the methylation-sensitive enzymes PstI, DdeI, MspI, HpaII, or EcoRII according to the manufacturer's recommendations. The glycerol mock treatment was performed as a control for each of the DNA samples. The digested DNA was extracted by phenol–chloroform and precipitated with 0.1x vol sodium acetate (pH 5.5) and 2.5x vol ethanol. PCR reactions were performed on digested DNA templates as described above (see supplemental Table 1 at http://www.genetics.org/supplemental/ for all primer sequences). Bisulfite sequencing analysis was performed as described (HAUN et al. 2007).

Bioinformatics analysis:
Annotations for differentially expressed genes were based on information available at the The Institute for Genomic Research Maize Gene Index (http://www.tigr.org/tigr-scripts/tgi/T_index.cgi?species=maize). The gene ontology (GO) annotations were obtained on the basis of assignment of the best Arabidopsis hits from the The Arabidopsis Information Resource website (http://www.arabidopsis.org/tools/bulk/go/index.jsp). The genetic map positions for Affymetrix array probe sets were predicted on the basis of identity with genetically mapped sequences or inferred on the basis of identity with BAC contig sequences that contained genetically mapped markers.

	RESULTS

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION ACKNOWLEDGEMENTS LITERATURE CITED

 
The expression of a limited number of genes is affected in zmet2 mutants:
To study the effects of disruption of CpNpG methylation in maize, we performed global profiling of gene expression levels in maize zmet2-m1 mutants using the Affymetrix maize GeneChip, which contains 17,622 probe sets corresponding to 13,500 genes. Although post-transcriptional regulation and protein modifications could affect the function of ZMET2 targets, in this study we investigated the gene transcript levels and used the word "expression" to refer to the abundance of RNA transcripts, not protein. Transcript levels of zmet2-m1 mutants and wild-type plants were compared in 11-day-old seedlings from two genetic backgrounds, B73 and Mo17, using separate statistical analyses of the GC-RMA processed signals from B73 and Mo17 genetic backgrounds. A series of statistical tests and expression criteria (see MATERIALS AND METHODS) were employed to identify 126 differentially expressed genes (0.72% of all GeneChip probe sets) with 65 genes upregulated and 61 genes downregulated in zmet2-m1 mutants (Table 1; supplemental Table 2 at http://www.genetics.org/supplemental/). Many of the differentially expressed genes showed only a transcriptional response to the zmet2-m1 mutation in one of the two inbred backgrounds. Notably, one of the genes identified as downregulated in the zmet2-m1 mutant relative to wild type in both B73 and Mo17 was the zmet2 gene itself. The microarray results for 33 genes were validated by semiquantitative RT–PCR with 30 (91%) genes showing a strong correlation between microarray and RT–PCR results (Table 1; supplemental Table 2). The genetic map positions were predicted for 50% of the Affymetrix array probe sets based on identity with genetically mapped sequences or with BAC contig sequences containing genetically mapped markers. Analysis of the map positions for 52 of the 126 differentially expressed genes did not display any chromosomal bias (supplemental Table 2). Additionally, we did not detect evidence for allelic variants linked to the Zmet2 locus on chromosome 10 that may result in differential expression not dependent upon effects of the zmet2-m1 mutation.

On the basis of studies from Arabidopsis that demonstrated a role for chromomethylases in maintaining DNA methylation of transposable elements, we might expect to see reactivation of silenced transposable elements in zmet2 mutant plants (LINDROTH et al. 2001; TOMPA et al. 2002; KATO et al. 2003). However, the Affymetrix maize microarray platform contains relatively few probe sets designed for repetitive features or transposable elements. Therefore, we hybridized the subset of RNA samples from B73 and B73 zmet2-m1 plants to a long oligonucleotide microarray platform with >50,000 features, including at least 355 probe sets derived from transposable elements (see supplemental online text, supplemental Table 3, and supplemental Figures 2 and 3 at http://www.genetics.org/supplemental/ for details on this experiment). Only 1 of the 179 differentially expressed genes identified using this platform was annotated as a transposable element (supplemental Table 3 at http://www.genetics.org/supplemental/), providing no evidence for extensive transcriptional activation of repetitive elements in the zmet2-m1 background. It is possible that the ZMET2 protein is involved in the methylation of repetitive elements in the maize genome (PAPA et al. 2001) but that the transcriptional silencing of these elements is not solely dependent upon CpNpG methylation.

Complete loss of ZMET2 activity is required to cause changes in gene expression:
The altered DNA methylation phenotype of the zmet2-m1 mutant displays codominant inheritance on the basis of changes in CpNpG methylation at different genomic loci (PAPA et al. 2001). To further validate the examples of differential expression and to determine if the altered expression patterns would exhibit codominant inheritance, we tested the expression of eight genes differentially expressed in B73 zmet2-m1 mutants in the sibling plants segregating for the zmet2-m1 allele (Figure 1). In all eight cases, regardless of whether the gene was up- or downregulated in the zmet2-m1 mutant, the heterozygous plants had the same expression levels as their wild-type siblings, suggesting that a complete loss of ZMET2 function is required for expression changes in these genes. The expression level of these eight genes was also tested in seedlings that were segregating for the zmet2-m2 allele, a second mutant allele of Zmet2 caused by a Mu insertion into exon 1 of the gene (data not shown). We found that the expression levels of these eight genes were similarly affected in plants homozygous for zmet2-m1 or zmet2-m2 alleles.

View larger version (46K): [in this window] [in a new window] [Download PPT slide]   FIGURE 1.— Complete loss of ZMET2 is required to cause changes in gene expression. Gene expression was analyzed by RT–PCR in sibling plants segregating for the zmet2-m1 allele in the B73 genotype. The Mez1 (AF443596) gene was used as a control to equalize cDNA concentrations and as a control for the absence of DNA in cDNA samples. The plant genotypes are shown as +/+ for wild-type, +/– for a heterozygous plant, and –/– for a homozygous zmet2-m1 plant. The no template control (NTC) lane did not include any cDNA. Gene CK827767 was silent in the wild-type or heterozygous backgrounds, but was expressed in the zmet2-m1 homozygous mutant background. Gene CA827096 was expressed in the wild-type or heterozygous backgrounds, but was silenced in the zmet2-m1 homozygous mutant background. We also confirmed gene expression changes specific to the homozygous mutant individuals for genes BM336416, CF629745, CO522679, CA402353, AW787466, and AI649924 (supplemental Table SOM2 at http://www.genetics.org/supplemental/). The "no-RT" controls were negative in all of the cases and are not shown.

View larger version (32K): [in this window] [in a new window] [Download PPT slide]   FIGURE 2.— Different sets of genes are affected in zmet2-m1 mutants in different genetic backgrounds. (A) The 126 genes affected in zmet2-m1 mutants were divided into five groups according to their relative expression in the four genotypes. Group 1 (yellow) genes were differentially expressed in both B73 and Mo17. Group 2 (red) genes were differentially expressed in B73 and were expressed in Mo17 and Mo17 zmet2-m1. Group 3 (orange) genes were differentially expressed in B73 but did not have detectable expression in Mo17 and Mo17 zmet2-m1. The remaining two groups display differential expression in Mo17 genetic background with detectable expression in B73 and B73 zmet2-m1 (group 4, green) or no detectable expression in B73 and B73 zmet2-m1 (group 5, blue). (B) The expression ratios of the wild-type vs. zmet2-m1 mutant plants were compared for the B73 and Mo17 genetic backgrounds to assess whether genes were differentially expressed in both genotypes or in only one of the genotypes. Colors indicate the group that each gene was assigned to, as described in A.

View larger version (55K): [in this window] [in a new window] [Download PPT slide]   FIGURE 3.— Genes affected in zmet2-m1 mutant plants revealed different patterns of gene expression changes. The semiquantitative RT–PCR results for the four genotypes are shown for five of the genes from different groups. The first two lanes on the left show the results obtained with B73 (B) and Mo17 (M) genomic DNA. The next four lanes show the products amplified using cDNA from wild-type B73 (B), B73 zmet2-m1 (Bz), wild-type Mo17 (M), or Mo17 zmet2-m1 (Mz). The "no-RT" controls were negative in all of the cases and are not shown. The Mez1 (AF443596) gene was used as a control to equalize cDNA concentrations and to check for the absence of DNA in cDNA samples.

To further assess the prevalence of inbred-variable differential expression conditioned by ZMET2 among maize inbred lines, we compared the expression of 27 genes in wild-type and mutant plants in a third genetic background, W22 (Table 2; Figure 4). All four group 1 genes that were tested in W22 zmet2-m1 plants displayed similar expression change patterns in all three genetic backgrounds. Of the 23 genes from groups 2–5, only 4 showed evidence of differential expression in the zmet2-m1 mutant relative to wild-type plants in W22. This finding suggests that the majority of the genes from these groups (19/23) were affected in only one of the three genotypes tested. The widespread prevalence of differential effects in different genetic backgrounds suggested that there are few conserved targets of ZMET2 (likely genes from group 1). Instead, the genes regulated by ZMET2 may represent genes that independently fell under control of ZMET2 in different genetic backgrounds.

View larger version (42K): [in this window] [in a new window] [Download PPT slide]   FIGURE 4.— Comparison of expression changes in zmet2-m1 mutants in B73, Mo17, and W22 genetic backgrounds. RT–PCR products are shown for the following genotypes: lane B, B73; lane Bz, B73 zmet2-m1; lane M, Mo17; lane Mz, Mo17 zmet2-m1; lane W, W22; lane Wz, W22 zmet2-m1; lane NTC, no template control (no cDNA). The "no-RT" controls were negative in all of the cases and are not shown. The Mez1 (AF443596) gene was used as a control to equalize cDNA concentrations and to check for the absence of DNA in cDNA samples (data not shown).

View larger version (45K): [in this window] [in a new window] [Download PPT slide]   FIGURE 5.— Expression, sequence, and methylation analysis for five CDE genes. In A–E, the expression level (based on normalized Affymetrix expression values) is displayed as a bar graph for each of the five genes. There is also a schematic for each of the five genes that indicates the position of exons (open boxes), transcription start sites (dashed arrow), primers used for sequencing and methylation analyses (arrows), restriction sites (vertical lines), B73-Mo17 SNPs (diamonds), insertion/deletions (triangles), and a LINE element (shaded box). The GenBank accession numbers for the Mo17 sequences are EU000403–EU000409. The line above the schematic for CDE4 indicates the region used for bisulfite sequencing (see Figure 6). Note that the scale for each of the schematics is different as indicated by the bars. The agarose gel electrophoresis images on the right show the results of methyl-sensitive PCR analyses for each gene. Control reactions with no restriction enzyme (No RE) were performed on two biological replicates of DNA from each of the four genotypes (B, B73; Bz, B73 zmet2-m1; M, Mo17; Mz, Mo17 zmet2-m1). The lack of amplified products following digestion with the methylation-insensitive enzyme FokI demonstrates complete digestion using these conditions. The presence of amplified products following digestion with the methylation-sensitive enzymes indicates the presence of CpNpG methylation for the methylation-sensitive enzymes EcoRII, MspI, and DdeI. The presence of amplified products following digestion with HpaII indicates the presence of CpG and/or CpCpG methylation. All five genes display a correlation between the low levels of gene expression and the high level of cytosine methylation in the B73 wild type compared to the other three genotypes. (F) For the gene CDE1 (CA402353), the presence of transcripts in seedling tissue from eight inbred lines was tested. In addition, the same methyl-sensitive PCR assay as in A was used to test the methylation level in genomic DNA isolated from eight inbred lines. The inbred lines are indicated above the lanes and the final lane is a no template control (NTC).

View larger version (25K): [in this window] [in a new window] [Download PPT slide]   FIGURE 6.— Analysis of DNA methylation for the upstream region of CDE4 (CK827767). The methylation status of the cytosines in a 413-bp region of the promoter of CDE4 (Figure 5D indicates the region of the promoter that was analyzed) was determined using bisulfite sequencing. Bisulfite sequencing was used to compare the methylation profile of B73 (20 cloned PCR products), B73 zmet2-m1 (10 cloned PCR products), and Mo17 (10 cloned PCR products) genotypes. The x-axis indicates the position of the cytosine while the y-axis indicates the percentage of methylation. The sequence context for each cytosine residue is indicated by the type of marker used (circle for CpG sites, triangle for CpNpG sites and "x" for asymmetric sites). The region sequenced included 10 CpG sites, 21 CpNpG sites, and 48 asymmetric sites. The CpNpG methylation is higher in B73 wild type than in the other two genotypes while the level of CpG and asymmetric methylation is similar in all three genotypes.

To determine whether the gene silencing conditioned by CpNpG methylation for ZMET2 targets is widely distributed among inbred lines or restricted to a single line, we tested the expression levels and promoter methylation states for CDE1–5 genes in seedling tissue from eight maize inbreds (Figure 5F; Table 3). For each of the genes we found that when the gene was expressed it was unmethylated and when the gene was not expressed there was CpNpG methylation at the restriction site(s) tested. Each of the five genes that we tested existed in a methylated and silenced state in at least two different inbred lines. One of the genes, CDE3, is methylated and silent in five of the eight inbreds that were surveyed. This finding, in combination with the ability to map the expression differences in recombinant inbred lines, suggests that the methylated and silenced allelic state of CDE genes is relatively stable.

	DISCUSSION

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION ACKNOWLEDGEMENTS LITERATURE CITED

 
DNA methylation patterns in plants are the result of a complex process. In addition to maintenance (MET) and domains rearranged methyltransferase (DRM) DNA methyltransferases, plants possess chromomethylases, responsible for maintenance of CpNpG methylation, which is specific to the plant kingdom. Chromomethylases are highly conserved across a variety of plant species and, therefore, should be functionally important. However, the exact role of CpNpG methylation has not been determined. Previous studies in maize (PAPA et al. 2001) and Arabidopsis (BARTEE et al. 2001; LINDROTH et al. 2001; CAO and JACOBSEN 2002a) have found little evidence for morphological phenotypes in CMT mutants. While mutations in either cmt3 or drm2 have very limited morphological phenotypes in Arabidopsis, the double mutant displays severe abnormalities, including developmental retardation and partial sterility, suggesting that, at least in Arabidopsis, chromomethylases may have partially redundant functions with the DRM methylases (LINDROTH et al. 2001; CAO and JACOBSEN 2002a,b). It is likely that partial redundancy also exists in maize. However, despite such functional redundancy, CMT3 has specific targets; for example, mutations in the cmt3 gene can have profound effects on the expression of the epigenetically controlled alleles of PAI or SUP (BARTEE et al. 2001; LINDROTH et al. 2001; CAO and JACOBSEN 2002a).

Functional targets of the maize chromomethylase ZMET2:
To determine the functional targets of CpNpG methylation in maize, we used two different genomewide gene expression profiling platforms. We found that relatively few genes displayed altered expression. These genes are distributed throughout the genome and do not display marked trends in their putative function. Moreover, many of the genes affected by the zmet2 mutation display large differences in expression level in different maize inbred lines with some of the genes being completely silent in some inbreds, suggesting that the functional targets of ZMET2 are primarily nonessential genes that can exhibit altered expression levels without severe phenotypes. We found very little evidence for upregulation of transposable elements in zmet2 mutant plants. In Arabidopsis, the silencing of some transposable elements was shown to be a redundant function of CMT3 and MET1 (KATO et al. 2003) Therefore, it is possible that ZMET2 does play a role in the methylation of transposable elements but that the silencing of these elements is redundant with methylation catalyzed by other enzymes, preventing reactivation of mobile elements in maize zmet2 mutants.

One caveat in identifying the functional targets of ZMET2 is that a gene with altered expression may represent a primary target of ZMET2 (a gene that is directly methylated by ZMET2) or a secondary target (a gene that is regulated by a primary target). In the five genes that we tested, we found evidence for direct CpNpG methylation changes in promoter or intron sequences correlating to gene expression levels, but we expect that the set of 126 genes identified represents a combination of primary and secondary targets of ZMET2. The genes upregulated in zmet2 mutants when CpNpG methylation is released are more likely to represent primary targets of ZMET2. The finding that expression changes occur only in homozygous mutant individuals suggests that ZMET2 affects individual gene targets in a qualitative, rather than quantitative, manner.

Natural variation for the functional targets of ZMET2:
It is notable that the majority of functional targets of ZMET2 identified in our study represent genotype-specific targets. It is possible that there are additional conserved ZMET2 targets that were not identified in our study. Our inability to detect such conserved ZMET2 targets may be explained by the redundancy among CMT, MET, and DRM methyltransferases in regulation of the conserved functional targets of ZMET2. Therefore, the analysis of plants containing mutations in several methyltransferases of different types is required to identify such targets. Alternatively, our study of whole seedlings might have missed important cases of conserved regulation by ZMET2 that occur in specific tissues or at specific developmental points.

Initially, we were quite surprised by the intraspecific variation for the targets of ZMET2. However, after extensive validation of this variation and reconsideration of previous studies, this finding is not as surprising. The previously documented targets of CMT in Arabidopsis are PAI and SUP. Both of these cases represent examples of epigenetic states that are present in some ecotypes but absent in others. Similarly, many of the maize genes affected by ZMET2 identified in this study also displayed epigenetic states that appear to be variable among inbred lines. We found that multiple instances of differential expression in the zmet2-m1 plants in B73 or Mo17 genetic backgrounds represent genes that are present in the genome but not expressed at detectable levels in one of the inbreds despite the absence of major rearrangements within coding or regulatory regions.

The seven CDE genes represent examples of putative epigenetic variation due to the dependence of their expression state upon ZMET2, the cis-regulatory variation, and methylation changes associated with expression states. However, we cannot determine whether these represent potential obligatory, facilitated, or pure epigenetic variation. We do not have enough sequence from the B73 and Mo17 alleles to assess contributions of linked genetic changes. There is a high level of intraspecific variation for the content and identity of repetitive sequences and transposons surrounding maize genes (FU and DOONER 2002; BRUNNER et al. 2005). We did not find evidence for variation in the content of transposable elements in the vicinity of the epigenetically controlled alleles identified in this study. It is likely that, at least for some of the alleles, alterations for the presence of specific repetitive sequences and/or transposable elements cause genes to become susceptible to regulation by ZMET2 in certain genetic backgrounds.

The ability to map the basis of the B73–Mo17 regulatory variation in recombinant inbred lines suggests that the expression differences exhibited by the B73 and Mo17 alleles are quite stable. In addition, it suggests that the expression differences cannot be communicated in a paramutation-like manner. If the silenced CDE allele were able to heritably alter the expressed CDE allele, we would not have been able to map the expression variation. Interestingly, if the allelic variation was caused by siRNAs, we might expect that the expressed CDE allele would be silenced when it is heterozygous with a silenced CDE allele. However, we did not find evidence supporting this mechanism. Moreover, we have evidence that one of the CDE genes exhibits mono-allelic expression in hybrids (R. M. STUPAR and N. M. SPRINGER, unpublished results). The presence of silenced alleles in multiple inbred lines suggests that such silent allelic states are stable and naturally exist within the population. The existence of naturally occurring epigenetically regulated alleles within maize suggests that these epigenetic changes may be selected upon and contribute to phenotypic variation. Naturally occurring epigenetic variation has been shown to produce different plant phenotypes and can have important implications for phenotypic variation and microevolution in natural plant populations (KALISZ and PURUGGANAN 2004). Thus, ZMET2 participates in creating novel epigenetic variation through modulation of gene expression. Further experimentation is necessary to understand the mechanisms that produce epigenetically regulated alleles and their stability within maize lines.

	ACKNOWLEDGEMENTS

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION ACKNOWLEDGEMENTS LITERATURE CITED

 
The authors gratefully acknowledge the assistance of Peter Hermanson in assisting with plant growth, nucleic acid isolations, and validation PCRs. We are also grateful to several anonymous reviews for their helpful comments on this manuscript. This work was supported by grants from the National Science Foundation DBI-0227310 to N.M.S. and DBI-0321663 to S.M.K.

	FOOTNOTES

 
Sequence data from this article have been deposited with the GenBank Data Libraries under accession nos. EU000403–EU000409. Microarray data from this article have been deposited at GEO database under the series accession no. GSE8188.

	LITERATURE CITED

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION ACKNOWLEDGEMENTS LITERATURE CITED

 

AUFSATZ, W., M. F. METTE, J. VAN DER WINDEN, M. MATZKE and A. J. MATZKE, 2002 HDA6, a putative histone deacetylase needed to enhance DNA methylation induced by double-stranded RNA. EMBO J. 21: 6832–6841.[CrossRef][Medline]

BARTEE, L., F. MALAGNAC and J. BENDER, 2001 Arabidopsis cmt3 chromomethylase mutations block non-CG methylation and silencing of an endogenous gene. Genes Dev. 15: 1753–1758.[Abstract/Free Full Text]

BENDER, J., and G. R. FINK, 1995 Epigenetic control of an endogenous gene family is revealed by a novel blue fluorescent mutant of Arabidopsis. Cell 83: 725–734.[CrossRef][Medline]

BRUNNER, S., K. FENGLER, M. MORGANTE, S. TINGEY and A. RAFALSKI, 2005 Evolution of DNA sequence nonhomologies among maize inbreds. Plant Cell 17: 343–360.[Abstract/Free Full Text]

BUCKLER, E. S., B. S. GAUT and M. D. MCMULLEN, 2006 Molecular and functional diversity of maize. Curr. Opin. Plant Biol. 9: 172–176.[CrossRef][Medline]

CAO, X., and S. E. JACOBSEN, 2002a Locus-specific control of asymmetric and CpNpG methylation by the DRM and CMT3 methyltransferase genes. Proc. Natl. Acad. Sci. USA 99(Suppl. 4): 16491–16498.[Abstract/Free Full Text]

CAO, X., and S. E. JACOBSEN, 2002b Role of the arabidopsis DRM methyltransferases in de novo DNA methylation and gene silencing. Curr. Biol. 12: 1138–1144.[CrossRef][Medline]

CAO, X., N. M. SPRINGER, M. G. MUSZYNSKI, R. L. PHILLIPS, S. KAEPPLER et al., 2000 Conserved plant genes with similarity to mammalian de novo DNA methyltransferases. Proc. Natl. Acad. Sci. USA 97: 4979–4984.[Abstract/Free Full Text]

CAO, X., W. AUFSATZ, D. ZILBERMAN, M. F. METTE, M. S. HUANG et al., 2003 Role of the DRM and CMT3 methyltransferases in RNA-directed DNA methylation. Curr. Biol. 13: 2212–2217.[CrossRef][Medline]

CHAN, S. W., D. ZILBERMAN, Z. XIE, L. K. JOHANSEN, J. C. CARRINGTON et al., 2004 RNA silencing genes control de novo DNA methylation. Science 303: 1336.[Free Full Text]

CHAN, S. W., I. R. HENDERSON and S. E. JACOBSEN, 2005 Gardening the genome: DNA methylation in Arabidopsis thaliana. Nat. Rev. Genet. 6: 351–360.[CrossRef][Medline]

CHAN, S. W., I. R. HENDERSON, X. ZHANG, G. SHAH, J. S. CHIEN et al., 2006 RNAi, DRD1, and histone methylation actively target developmentally important non-CG DNA methylation in Arabidopsis. PLoS Genet. 2: e83.[CrossRef][Medline]

CHANDLER, V. L., and M. STAM, 2004 Chromatin conversations: mechanisms and implications of paramutation. Nat. Rev. Genet. 5: 532–544.[CrossRef][Medline]

CHANDLER, V. L., W. B. EGGLESTON and J. E. DORWEILER, 2000 Paramutation in maize. Plant Mol. Biol. 43: 121–145.[CrossRef][Medline]

CUBAS, P., C. VINCENT and E. COEN, 1999 An epigenetic mutation responsible for natural variation in floral symmetry. Nature 401: 157–161.[CrossRef][Medline]

EBBS, M. L., and J. BENDER, 2006 Locus-specific control of DNA methylation by the Arabidopsis SUVH5 histone methyltransferase. Plant Cell 18: 1166–1176.[Abstract/Free Full Text]

FINNEGAN, E. J., and E. S. DENNIS, 1993 Isolation and identification by sequence homology of a putative cytosine methyltransferase from Arabidopsis thaliana. Nucleic Acids Res. 21: 2383–2388.[Abstract/Free Full Text]

FU, H., and H. K. DOONER, 2002 Intraspecific violation of genetic colinearity and its implications in maize. Proc. Natl. Acad. Sci. USA 99: 9573–9578.[Abstract/Free Full Text]

FU, Y., S. J. EMRICH, L. GUO, T. J. WEN, D. A. ASHLOCK et al., 2005 Quality assessment of maize assembled genomic islands (MAGIs) and large-scale experimental verification of predicted genes. Proc. Natl. Acad. Sci. USA 102: 12282–12287.[Abstract/Free Full Text]

GARDINER, J. M., C. R. BUELL, R. ELUMALAI, D. W. GALBRAITH, D. A. HENDERSON et al., 2005 Design, production, and utilization of long oligonucleotide microarrays for expression analysis in maize. Maydica 50: 425–435.

GENGER, R. K., K. A. KOVAC, E. S. DENNIS, W. J. PEACOCK and E. J. FINNEGAN, 1999 Multiple DNA methyltransferase genes in Arabidopsis thaliana. Plant Mol. Biol. 41: 269–278.[CrossRef][Medline]

HAUN, W. J., S. LAOUEILLE-DUPRAT, J. O'CONNELL, M. C. SPILLANE, U. GROSSNIKLAUS et al., 2007 Genomic imprinting, methylation and molecular evolution of maize Enhancer of zeste (Mez) homologs. Plant J. 49: 325–337.[CrossRef][Medline]

HENIKOFF, S., and L. COMAI, 1998 A DNA methyltransferase homolog with a chromodomain exists in multiple polymorphic forms in Arabidopsis. Genetics 149: 307–318.[Abstract/Free Full Text]

JACKSON, J. P., A. M. LINDROTH, X. CAO and S. E. JACOBSEN, 2002 Control of CpNpG DNA methylation by the KRYPTONITE histone H3 methyltransferase. Nature 416: 556–560.[CrossRef][Medline]

JEDDELOH, J. A., T. L. STOKES and E. J. RICHARDS, 1999 Maintenance of genomic methylation requires a SWI2/SNF2-like protein. Nat. Genet. 22: 94–97.[CrossRef][Medline]

KALISZ, S., and M. D. PURUGGANAN, 2004 Epialleles via DNA methylation: consequences for plant evolution. Trends Ecol. Evol. 19: 309–314.[CrossRef][Medline]

KANKEL, M. W., D. E. RAMSEY, T. L. STOKES, S. K. FLOWERS, J. R. HAAG et al., 2003 Arabidopsis MET1 cytosine methyltransferase mutants. Genetics 163: 1109–1122.[Abstract/Free Full Text]

KATO, M., A. MIURA, J. BENDER, S. E. JACOBSEN and T. KAKUTANI, 2003 Role of CG and non-CG methylation in immobilization of transposons in Arabidopsis. Curr. Biol. 13: 421–426.[CrossRef][Medline]

LEE, M., N. SHAROPOVA, W. D. BEAVIS, D. GRANT, M. KATT et al., 2002 Expanding the genetic map of maize with the intermated B73 x Mo17 (IBM) population. Plant Mol. Biol. 48: 453–461.[CrossRef][Medline]

LINDROTH, A. M., X. CAO, J. P. JACKSON, D. ZILBERMAN, C. M. MCCALLUM et al., 2001 Requirement of CHROMOMETHYLASE3 for maintenance of CpXpG methylation. Science 292: 2077–2080.[Abstract/Free Full Text]

MALAGNAC, F., L. BARTEE and J. BENDER, 2002 An Arabidopsis SET domain protein required for maintenance but not establishment of DNA methylation. EMBO J. 21: 6842–6852.[CrossRef][Medline]

MANNING, K., M. TOR, M. POOLE, Y. HONG, A. J. THOMPSON et al., 2006 A naturally occurring epigenetic mutation in a gene encoding an SBP-box transcription factor inhibits tomato fruit ripening. Nat. Genet. 38: 948–952.[CrossRef][Medline]

MORGAN, H. D., H. G. SUTHERLAND, D. I. MARTIN and E. WHITELAW, 1999 Epigenetic inheritance at the agouti locus in the mouse. Nat. Genet. 23: 314–318.[CrossRef][Medline]

PAPA, C. M., N. M. SPRINGER, M. G. MUSZYNSKI, R. MEELEY and S. M. KAEPPLER, 2001 Maize chromomethylase Zea methyltransferase2 is required for CpNpG methylation. Plant Cell 13: 1919–1928.[Abstract/Free Full Text]

RANGWALA, S. H., R. ELUMALAI, C. VANIER, H. OZKAN, D. W. GALBRAITH et al., 2006 Meiotically stable natural epialleles of Sadhu, a novel Arabidopsis retroposon. PLoS Genet. 2: e36.[CrossRef][Medline]

RICHARDS, E. J., 2006 Inherited epigenetic variation: revisiting soft inheritance. Nat. Rev. Genet. 7: 395–401.[CrossRef][Medline]

RIDDLE, N. C., and E. J. RICHARDS, 2005 Genetic variation in epigenetic inheritance of ribosomal RNA gene methylation in Arabidopsis. Plant J. 41: 524–532.[CrossRef][Medline]

ROZEN, S., and H. SKALETSKY, 2000 Primer3 on the WWW for general users and for biologist programmers. Methods Mol. Biol. 132: 365–386.[Medline]

SAGHAI-MAROOF, M. A., K. M. SOLIMAN, R. A. JORGENSEN and R. W. ALLARD, 1984 Ribosomal DNA spacer-length polymorphisms in barley: Mendelian inheritance, chromosomal location, and population dynamics. Proc. Natl. Acad. Sci. USA 81: 8014–8018.[Abstract/Free Full Text]

SMYTH, G. K., 2004 Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Stat. Appl. Genet. Mol. Biol. 3: Article3.[Medline]

STOKES, T. L., B. N. KUNKEL and E. J. RICHARDS, 2002 Epigenetic variation in Arabidopsis disease resistance. Genes Dev. 16: 171–182.[Abstract/Free Full Text]

STUPAR, R. M., and N. M. SPRINGER, 2006 Cis-transcriptional variation in maize inbred lines B73 and Mo17 leads to additive expression patterns in the F₁ hybrid. Genetics 173: 2199–2210.[Abstract/Free Full Text]

TOMPA, R., C. M. MCCALLUM, J. DELROW, J. G. HENIKOFF, B. VAN STEENSEL et al., 2002 Genome-wide profiling of DNA methylation reveals transposon targets of CHROMOMETHYLASE3. Curr. Biol. 12: 65–68.[CrossRef][Medline]

VONGS, A., T. KAKUTANI, R. A. MARTIENSSEN and E. J. RICHARDS, 1993 Arabidopsis thaliana DNA methylation mutants. Science 260: 1926–1928.[Abstract/Free Full Text]

XIAO, W., K. D. CUSTARD, R. C. BROWN, B. E. LEMMON, J. J. HARADA et al., 2006 DNA methylation is critical for Arabidopsis embryogenesis and seed viability. Plant Cell 18: 805–814.[Abstract/Free Full Text]

ZILBERMAN, D., X. CAO and S. E. JACOBSEN, 2003 ARGONAUTE4 control of locus-specific siRNA accumulation and DNA and histone methylation. Science 299: 716–719.[Abstract/Free Full Text]

This Article Abstract Full Text (PDF) Data Supplement All Versions of this Article: genetics.107.072702v1 177/2/749    most recent Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via Google Scholar Google Scholar Articles by Makarevitch, I. Articles by Springer, N. M. Search for Related Content PubMed