RNA-directed DNA methylation enforces boundaries between heterochromatin and euchromatin in the maize genome.

The maize genome is relatively large (∼ 2.3 Gb) and has a complex organization of interspersed genes and transposable elements, which necessitates frequent boundaries between different types of chromatin. The examination of maize genes and conserved noncoding sequences revealed that many of these are flanked by regions of elevated asymmetric CHH (where H is A, C, or T) methylation (termed mCHH islands). These mCHH islands are quite short (∼ 100 bp), are enriched near active genes, and often occur at the edge of the transposon that is located nearest to genes. The analysis of DNA methylation in other sequence contexts and several chromatin modifications revealed that mCHH islands mark the transition from heterochromatin-associated modifications to euchromatin-associated modifications. The presence of an mCHH island is fairly consistent in several distinct tissues that were surveyed but shows some variation among different haplotypes. The presence of insertion/deletions in promoters often influences the presence and position of an mCHH island. The mCHH islands are dependent upon RNA-directed DNA methylation activities and are lost in mop1 and mop3 mutants, but the nearby genes rarely exhibit altered expression levels. Instead, loss of an mCHH island is often accompanied by additional loss of DNA methylation in CG and CHG contexts associated with heterochromatin in nearby transposons. This suggests that mCHH islands and RNA-directed DNA methylation near maize genes may act to preserve the silencing of transposons from activity of nearby genes.

Accessible DNA and relative depletion of H3K9me2 at maize loci undergoing RNA-directed DNA methylation.

RNA-directed DNA methylation (RdDM) in plants is a well-characterized example of RNA interference-related transcriptional gene silencing. To determine the relationships between RdDM and heterochromatin in the repeat-rich maize (Zea mays) genome, we performed whole-genome analyses of several heterochromatic features: dimethylation of lysine 9 and lysine 27 (H3K9me2 and H3K27me2), chromatin accessibility, DNA methylation, and small RNAs; we also analyzed two mutants that affect these processes, mediator of paramutation1 and zea methyltransferase2. The data revealed that the majority of the genome exists in a heterochromatic state defined by inaccessible chromatin that is marked by H3K9me2 and H3K27me2 but that lacks RdDM. The minority of the genome marked by RdDM was predominantly near genes, and its overall chromatin structure appeared more similar to euchromatin than to heterochromatin. These and other data indicate that the densely staining chromatin defined as heterochromatin differs fundamentally from RdDM-targeted chromatin. We propose that small interfering RNAs perform a specialized role in repressing transposons in accessible chromatin environments and that the bulk of heterochromatin is incompatible with small RNA production.

Differential nuclease sensitivity profiling of chromatin reveals biochemical footprints coupled to gene expression and functional DNA elements in maize.

The eukaryotic genome is organized into nucleosomes, the fundamental units of chromatin. The positions of nucleosomes on DNA regulate protein-DNA interactions and in turn influence DNA-templated events. Despite the increasing number of genome-wide maps of nucleosome position, how global changes in gene expression relate to changes in nucleosome position is poorly understood. We show that in nucleosome occupancy mapping experiments in maize (Zea mays), particular genomic regions are highly susceptible to variation introduced by differences in the extent to which chromatin is digested with micrococcal nuclease (MNase). We exploited this digestion-linked variation to identify protein footprints that are hypersensitive to MNase digestion, an approach we term differential nuclease sensitivity profiling (DNS-chip). Hypersensitive footprints were enriched at the 5' and 3' ends of genes, associated with gene expression levels, and significantly overlapped with conserved noncoding sequences and the binding sites of the transcription factor KNOTTED1. We also found that the tissue-specific regulation of gene expression was linked to tissue-specific hypersensitive footprints. These results reveal biochemical features of nucleosome organization that correlate with gene expression levels and colocalize with functional DNA elements. This approach to chromatin profiling should be broadly applicable to other species and should shed light on the relationships among chromatin organization, protein-DNA interactions, and genome regulation.

Genetic perturbation of the maize methylome.

DNA methylation can play important roles in the regulation of transposable elements and genes. A collection of mutant alleles for 11 maize (Zea mays) genes predicted to play roles in controlling DNA methylation were isolated through forward- or reverse-genetic approaches. Low-coverage whole-genome bisulfite sequencing and high-coverage sequence-capture bisulfite sequencing were applied to mutant lines to determine context- and locus-specific effects of these mutations on DNA methylation profiles. Plants containing mutant alleles for components of the RNA-directed DNA methylation pathway exhibit loss of CHH methylation at many loci as well as CG and CHG methylation at a small number of loci. Plants containing loss-of-function alleles for chromomethylase (CMT) genes exhibit strong genome-wide reductions in CHG methylation and some locus-specific loss of CHH methylation. In an attempt to identify stocks with stronger reductions in DNA methylation levels than provided by single gene mutations, we performed crosses to create double mutants for the maize CMT3 orthologs, Zmet2 and Zmet5, and for the maize DDM1 orthologs, Chr101 and Chr106. While loss-of-function alleles are viable as single gene mutants, the double mutants were not recovered, suggesting that severe perturbations of the maize methylome may have stronger deleterious phenotypic effects than in Arabidopsis thaliana.

Direct and Indirect Transcriptional Effects of Abiotic Stress in Zea mays Plants Defective in RNA-Directed DNA Methylation.

Plants respond to abiotic stress stimuli, such as water deprivation, through a hierarchical cascade that includes detection and signaling to mediate transcriptional and physiological changes. The phytohormone abscisic acid (ABA) is well-characterized for its regulatory role in these processes in response to specific environmental cues. ABA-mediated changes in gene expression have been demonstrated to be temporally-dependent, however, the genome-wide timing of these responses are not well-characterized in the agronomically important crop plant Zea mays (maize). ABA-mediated responses are synergistic with other regulatory mechanisms, including the plant-specific RNA-directed DNA methylation (RdDM) epigenetic pathway. Our prior work demonstrated that after relatively long-term ABA induction (8 h), maize plants homozygous for the mop1-1 mutation, defective in a component of the RdDM pathway, exhibit enhanced transcriptional sensitivity to the phytohormone. At this time-point, many hierarchically positioned transcription factors are differentially expressed resulting in primary (direct) and secondary (indirect) transcriptional outcomes. To identify more immediate and direct MOP1-dependent responses to ABA, we conducted a transcriptomic analysis using mop1-1 mutant and wild type plants treated with ABA for 1 h. One h of ABA treatment was sufficient to induce unique categories of differentially expressed genes (DEGs) in mop1-1. A comparative analysis between the two time-points revealed that distinct epigenetically-regulated changes in gene expression occur within the early stages of ABA induction, and that these changes are predicted to influence less immediate, indirect transcriptional responses. Homology with MOP1-dependent siRNAs and a gene regulatory network (GRN) were used to identify putative immediate and indirect targets, respectively. By manipulating two key regulatory networks in a temporal dependent manner, we identified genes and biological processes regulated by RdDM and ABA-mediated stress responses. Consistent with mis-regulation of gene expression, mop1-1 homozygous plants are compromised in their ability to recover from water deprivation. Collectively, these results indicate transcriptionally and physiologically relevant roles for MOP1-mediated regulation of gene expression of plant responses to environmental stress.

Epigenetic Regulation of ABA-Induced Transcriptional Responses in Maize.

Plants are subjected to extreme environmental conditions and must adapt rapidly. The phytohormone abscisic acid (ABA) accumulates during abiotic stress, signaling transcriptional changes that trigger physiological responses. Epigenetic modifications often facilitate transcription, particularly at genes exhibiting temporal, tissue-specific and environmentally-induced expression. In maize (Zea mays), MEDIATOR OF PARAMUTATION 1 (MOP1) is required for progression of an RNA-dependent epigenetic pathway that regulates transcriptional silencing of loci genomewide. MOP1 function has been previously correlated with genomic regions adjoining particular types of transposable elements and genic regions, suggesting that this regulatory pathway functions to maintain distinct transcriptional activities within genomic spaces, and that loss of MOP1 may modify the responsiveness of some loci to other regulatory pathways. As critical regulators of gene expression, MOP1 and ABA pathways each regulate specific genes. To determine whether loss of MOP1 impacts ABA-responsive gene expression in maize, mop1-1 and Mop1 homozygous seedlings were subjected to exogenous ABA and RNA-sequencing. A total of 3,242 differentially expressed genes (DEGs) were identified in four pairwise comparisons. Overall, ABA-induced changes in gene expression were enhanced in mop1-1 homozygous plants. The highest number of DEGs were identified in ABA-induced mop1-1 mutants, including many transcription factors; this suggests combinatorial regulatory scenarios including direct and indirect transcriptional responses to genetic disruption (mop1-1) and/or stimulus-induction of a hierarchical, cascading network of responsive genes. Additionally, a modest increase in CHH methylation at putative MOP1-RdDM loci in response to ABA was observed in some genotypes, suggesting that epigenetic variation might influence environmentally-induced transcriptional responses in maize.

Function of Succinoglycan Polysaccharide in Sinorhizobium meliloti Host Plant Invasion Depends on Succinylation, Not Molecular Weight.

UNLABELLED: The acidic polysaccharide succinoglycan produced by the rhizobial symbiont Sinorhizobium meliloti 1021 is required for this bacterium to invade the host plant Medicago truncatula and establish a nitrogen-fixing symbiosis. S. meliloti mutants that cannot make succinoglycan cannot initiate invasion structures called infection threads in plant root hairs. S. meliloti exoH mutants that cannot succinylate succinoglycan are also unable to form infection threads, despite the fact that they make large quantities of succinoglycan. Succinoglycan produced by exoH mutants is refractory to cleavage by the glycanases encoded by exoK and exsH, and thus succinoglycan produced by exoH mutants is made only in the high-molecular-weight (HMW) form. One interpretation of the symbiotic defect of exoH mutants is that the low-molecular-weight (LMW) form of succinoglycan is required for infection thread formation. However, our data demonstrate that production of the HMW form of succinoglycan by S. meliloti 1021 is sufficient for invasion of the host M. truncatula and that the LMW form is not required. Here, we show that S. meliloti strains deficient in the exoK- and exsH-encoded glycanases invade M. truncatula and form a productive symbiosis, although they do this with somewhat less efficiency than the wild type. We have also characterized the polysaccharides produced by these double glycanase mutants and determined that they consist of only HMW succinoglycan and no detectable LMW succinoglycan. This demonstrates that LMW succinoglycan is not required for host invasion. These results suggest succinoglycan function is not dependent upon the presence of a small, readily diffusible form. IMPORTANCE: Sinorhizobium meliloti is a bacterium that forms a beneficial symbiosis with legume host plants. S. meliloti and other rhizobia convert atmospheric nitrogen to ammonia, a nutrient source for the host plant. To establish the symbiosis, rhizobia must invade plant roots, supplying the proper signals to prevent a plant immune response during invasion. A polysaccharide, succinoglycan, produced by S. meliloti is required for successful invasion. Here, we show that the critical feature of succinoglycan that allows infection to proceed is the attachment of a "succinyl" chemical group and that the chain length of succinoglycan is much less important for its function. We also show that none of the short-chain versions of succinoglycan is produced in the absence of two chain-cleaving enzymes.

Chromatin structure and gene expression changes associated with loss of MOP1 activity in Zea mays.

Though the mechanisms governing nuclear organization are not well understood, it is apparent that epigenetic modifications coordinately modulate chromatin organization as well as transcription. In maize, MEDIATOR OF PARAMUTATION1 (MOP1) is required for 24 nt siRNA-mediated epigenetic regulation and transcriptional gene silencing via a putative Pol IV- RdDM pathway. To elucidate the mechanisms of nuclear chromatin organization, we investigated the relationship between chromatin structure and transcription in response to loss of MOP1 function. We used a microarray based micrococcal nuclease sensitivity assay to identify genome-wide changes in chromatin structure in mop1-1 immature ears and observed an increase in chromatin accessibility at chromosome arms associated with loss of MOP1 function. Within the many genes misregulated in mop1 mutants, we identified one subset likely to be direct targets of epigenetic transcriptional silencing via Pol-IV RdDM. We found that target specificity for MOP1-mediated RdDM activity is governed by multiple signals that include accumulation of 24 nt siRNAs and the presence of specific classes of gene-proximal transposons, but neither of these attributes alone is sufficient to predict transcriptional misregulation in mop1-1 homozygous mutants. Our results suggest a role for MOP1 in regulation of higher-order chromatin organization where loss of MOP1 activity at a subset of loci triggers a broader cascade of transcriptional consequences and genome-wide changes in chromatin structure.

Identification of epigenetic regulators of a transcriptionally silenced transgene in maize.

Transcriptional gene silencing is a gene regulatory mechanism essential to all organisms. Many transcriptional regulatory mechanisms are associated with epigenetic modifications such as changes in chromatin structure, acetylation and methylation of core histone proteins, and DNA methylation within regulatory regions of endogenous genes and transgenes. Although several maize mutants have been identified from prior forward genetic screens for epigenetic transcriptional silencing, these screens have been far from saturated. Herein, the transcriptionally silent b1 genomic transgene (BTG-silent), a stable, epigenetically silenced transgene in Zea mays (maize), is demonstrated to be an effective phenotype for a forward genetic screen. When the transgene is reactivated, a dark purple plant phenotype is evident because the B1 transcription factor activates anthocyanin biosynthesis, making loss of silencing mutants easy to identify. Using BTG-silent, ten new putative mutants were identified and named transgene reactivated1 through 11 (tgr1-6 and tgr8-11). Three of these mutants have been examined in more detail, and molecular and genetic assays demonstrated that these mutants have both distinct and overlapping phenotypes with previously identified maize mutants that relieve epigenetic transcriptional silencing. Linkage analysis suggests that tgr2 and tgr3 do not correspond to a mutation at previously identified maize loci resulting from other forward genetic screens, while tgr1 shows linkage to a characterized gene. These results suggest that the mutants are a valuable resource for future studies because some of the mutants are likely to reveal genes that encode products required for epigenetic gene regulation in maize but are not currently represented by sequenced mutations.