Heritable changes of epialleles near genes in maize can be triggered in the absence of CHH methylation.

Trans-chromosomal interactions resulting in changes in DNA methylation during hybridization have been observed in several plant species. However, little is known about the causes or consequences of these interactions. Here, we compared DNA methylomes of F1 hybrids that are mutant for a small RNA biogenesis gene, Mop1 (Mediator of paramutation1), with that of their parents, wild-type siblings, and backcrossed progeny in maize (Zea mays). Our data show that hybridization triggers global changes in both trans-chromosomal methylation (TCM) and trans-chromosomal demethylation (TCdM), most of which involved changes in CHH methylation. In more than 60% of these TCM differentially methylated regions (DMRs) in which small RNAs are available, no significant changes in the quantity of small RNAs were observed. Methylation at the CHH TCM DMRs was largely lost in the mop1 mutant, although the effects of this mutant varied depending on the location of these DMRs. Interestingly, an increase in CHH at TCM DMRs was associated with enhanced expression of a subset of highly expressed genes and suppressed expression of a small number of lowly expressed genes. Examination of the methylation levels in backcrossed plants demonstrates that both TCM and TCdM can be maintained in the subsequent generation, but that TCdM is more stable than TCM. Surprisingly, although increased CHH methylation in most TCM DMRs in F1 plants required Mop1, initiation of a new epigenetic state of these DMRs did not require a functional copy of this gene, suggesting that initiation of these changes is independent of RNA-directed DNA methylation.

Drought-induced circular RNAs in maize roots: separating signal from noise.

Circular RNAs (CircRNAs) play an important role in diverse biological processes; however, their origin and functions, especially in plants, remain largely unclear. Here, we used two maize (Zea mays) inbred lines, as well as 14 of their derivative RILs with different drought sensitivity, to systematically characterize 8,790 circRNAs in maize roots under well-watered (WW) and water-stress (WS) conditions. We found that a diverse set of circRNAs expressed at significantly higher levels under WS. Enhanced expression of circRNAs was associated with longer flanking introns and an enrichment of long interspersed nuclear element (LINE) retrotransposable elements. The epigenetic marks found at the back-splicing junctions of circRNA-producing genes were markedly different from canonical splicing, characterized by increased levels of H3K36me3/H3K4me1, as well as decreased levels of H3K9Ac/H3K27Ac. We found that genes expressing circRNAs are subject to relaxed selection. The significant enrichment of trait-associated sites along their genic regions suggested that genes giving rise to circRNAs were associated with plant survival rate under drought stress, implying that circRNAs play roles in plant drought responses. Furthermore, we found that overexpression of circMED16, one of the drought-responsive circRNAs, enhances drought tolerance in Arabidopsis (Arabidopsis thaliana). Our results provide a framework for understanding the intricate interplay of epigenetic modifications and how they contribute to the fine-tuning of circRNA expression under drought stress.

Broad-spectrum fungal resistance in sorghum is conferred through the complex regulation of an immune receptor gene embedded in a natural antisense transcript.

Sorghum (Sorghum bicolor), the fifth most widely grown cereal crop globally, provides food security for millions of people. Anthracnose caused by the fungus Colletotrichum sublineola is a major disease of sorghum worldwide. We discovered a major fungal resistance locus in sorghum composed of the nucleotide-binding leucine-rich repeat receptor gene ANTHRACNOSE RESISTANCE GENE1 (ARG1) that is completely nested in an intron of a cis-natural antisense transcript (NAT) gene designated CARRIER OF ARG1 (CARG). Susceptible genotypes express CARG and two alternatively spliced ARG1 transcripts encoding truncated proteins lacking the leucine-rich repeat domains. In resistant genotypes, elevated expression of an intact allele of ARG1, attributed to the loss of CARG transcription and the presence of miniature inverted-repeat transposable element sequences, resulted in broad-spectrum resistance to fungal pathogens with distinct virulence strategies. Increased ARG1 expression in resistant genotypes is also associated with higher histone H3K4 and H3K36 methylation. In susceptible genotypes, lower ARG1 expression is associated with reduced H3K4 and H3K36 methylation and increased expression of NATs of CARG. The repressive chromatin state associated with H3K9me2 is low in CARG-expressing genotypes within the CARG exon and higher in genotypes with low CARG expression. Thus, ARG1 is regulated by multiple mechanisms and confers broad-spectrum, strong resistance to fungal pathogens.

ANTHRACNOSE RESISTANCE GENE2 confers fungal resistance in sorghum.

Sorghum is an important food and feed crop globally; its production is hampered by anthracnose disease caused by the fungal pathogen Colletotrichum sublineola (Cs). Here, we report identification and characterization of ANTHRACNOSE RESISTANCE GENE 2 (ARG2) encoding a nucleotide-binding leucine-rich repeat (NLR) protein that confers race-specific resistance to Cs strains. ARG2 is one of a cluster of several NLR genes initially identified in the sorghum differential line SC328C that is resistant to some Cs strains. This cluster shows structural and copy number variations in different sorghum genotypes. Different sorghum lines carrying independent ARG2 alleles provided the genetic validation for the identity of the ARG2 gene. ARG2 expression is induced by Cs, and chitin induces ARG2 expression in resistant but not in susceptible lines. ARG2-mediated resistance is accompanied by higher expression of defense and secondary metabolite genes at early stages of infection, and anthocyanin and zeatin metabolisms are upregulated in resistant plants. Interestingly, ARG2 localizes to the plasma membrane when transiently expressed in Nicotiana benthamiana. Importantly, ARG2 plants produced higher shoot dry matter than near-isogenic lines carrying the susceptible allele suggesting an absence of an ARG2 associated growth trade-off. Furthermore, ARG2-mediated resistance is stable at a wide range of temperatures. Our observations open avenues for resistance breeding and for dissecting mechanisms of resistance.

Conserved noncoding sequences and de novo Mutator insertion alleles are imprinted in maize.

Genomic imprinting is an epigenetic phenomenon in which differential allele expression occurs in a parent-of-origin-dependent manner. Imprinting in plants is tightly linked to transposable elements (TEs), and it has been hypothesized that genomic imprinting may be a consequence of demethylation of TEs. Here, we performed high-throughput sequencing of ribonucleic acids from four maize (Zea mays) endosperms that segregated newly silenced Mutator (Mu) transposons and identified 110 paternally expressed imprinted genes (PEGs) and 139 maternally expressed imprinted genes (MEGs). Additionally, two potentially novel paternally suppressed MEGs are associated with de novo Mu insertions. In addition, we find evidence for parent-of-origin effects on expression of 407 conserved noncoding sequences (CNSs) in maize endosperm. The imprinted CNSs are largely localized within genic regions and near genes, but the imprinting status of the CNSs are largely independent of their associated genes. Both imprinted CNSs and PEGs have been subject to relaxed selection. However, our data suggest that although MEGs were already subject to a higher mutation rate prior to their being imprinted, imprinting may be the cause of the relaxed selection of PEGs. In addition, although DNA methylation is lower in the maternal alleles of both the maternally and paternally expressed CNSs (mat and pat CNSs), the difference between the two alleles in H3K27me3 levels was only observed in pat CNSs. Together, our findings point to the importance of both transposons and CNSs in genomic imprinting in maize.

RNA-directed DNA methylation prevents rapid and heritable reversal of transposon silencing under heat stress in Zea mays.

In large complex plant genomes, RNA-directed DNA methylation (RdDM) ensures that epigenetic silencing is maintained at the boundary between genes and flanking transposable elements. In maize, RdDM is dependent on Mediator of Paramutation1 (Mop1), a gene encoding a putative RNA dependent RNA polymerase. Here we show that although RdDM is essential for the maintenance of DNA methylation of a silenced MuDR transposon in maize, a loss of that methylation does not result in a restoration of activity. Instead, heritable maintenance of silencing is maintained by histone modifications. At one terminal inverted repeat (TIR) of this element, heritable silencing is mediated via histone H3 lysine 9 dimethylation (H3K9me2), and histone H3 lysine 27 dimethylation (H3K27me2), even in the absence of DNA methylation. At the second TIR, heritable silencing is mediated by histone H3 lysine 27 trimethylation (H3K27me3), a mark normally associated with somatically inherited gene silencing. We find that a brief exposure of high temperature in a mop1 mutant rapidly reverses both of these modifications in conjunction with a loss of transcriptional silencing. These reversals are heritable, even in mop1 wild-type progeny in which methylation is restored at both TIRs. These observations suggest that DNA methylation is neither necessary to maintain silencing, nor is it sufficient to initiate silencing once has been reversed. However, given that heritable reactivation only occurs in a mop1 mutant background, these observations suggest that DNA methylation is required to buffer the effects of environmental stress on transposable elements.

The mop1 mutation affects the recombination landscape in maize.

Meiotic recombination is a fundamental process that generates genetic diversity and ensures the accurate segregation of homologous chromosomes. While a great deal is known about genetic factors that regulate recombination, relatively little is known about epigenetic factors, such as DNA methylation. In maize, we examined the effects on meiotic recombination of a mutation in a component of the RNA-directed DNA methylation pathway, Mop1 (Mediator of paramutation1), as well as a mutation in a component of the trans-acting small interference RNA biogenesis pathway, Lbl1 (Leafbladeless1). MOP1 is of particular interest with respect to recombination because it is responsible for methylation of transposable elements that are immediately adjacent to transcriptionally active genes. In the mop1 mutant, we found that meiotic recombination is uniformly decreased in pericentromeric regions but is generally increased in gene rich chromosomal arms. This observation was further confirmed by cytogenetic analysis showing that although overall crossover numbers are unchanged, they occur more frequently in chromosomal arms in mop1 mutants. Using whole genome bisulfite sequencing, our data show that crossover redistribution is driven by loss of CHH (where H = A, T, or C) methylation within regions near genes. In contrast to what we observed in mop1 mutants, no significant changes were observed in the frequency of meiotic recombination in lbl1 mutants. Our data demonstrate that CHH methylation has a significant impact on the overall recombination landscape in maize despite its low frequency relative to CG and CHG methylation.

Transposable elements employ distinct integration strategies with respect to transcriptional landscapes in eukaryotic genomes.

Transposable elements (TEs) are ubiquitous DNA segments capable of moving from one site to another within host genomes. The extant distributions of TEs in eukaryotic genomes have been shaped by both bona fide TE integration preferences in eukaryotic genomes and by selection following integration. Here, we compare TE target site distribution in host genomes using multiple de novo transposon insertion datasets in both plants and animals and compare them in the context of genome-wide transcriptional landscapes. We showcase two distinct types of transcription-associated TE targeting strategies that suggest a process of convergent evolution among eukaryotic TE families. The integration of two precision-targeting elements are specifically associated with initiation of RNA Polymerase II transcription of highly expressed genes, suggesting the existence of novel mechanisms of precision TE targeting in addition to passive targeting of open chromatin. We also highlight two features that can facilitate TE survival and rapid proliferation: tissue-specific transposition and minimization of negative impacts on nearby gene function due to precision targeting.

Patterns and Consequences of Subgenome Differentiation Provide Insights into the Nature of Paleopolyploidy in Plants.

Polyploidy is an important feature of plant genomes, but the nature of many polyploidization events remains to be elucidated. Here, we demonstrate that the evolutionary fates of the subgenomes in maize (Zea mays) and soybean (Glycine max) have followed different trajectories. One subgenome has been subject to relaxed selection, lower levels of gene expression, higher rates of transposable element accumulation, more small interfering RNAs and DNA methylation around genes, and higher rates of gene loss in maize, whereas none of these features were observed in soybean. Nevertheless, individual gene pairs exhibit differentiation with respect to these features in both species. In addition, we observed a higher number of chromosomal rearrangements and higher frequency of retention of duplicated genes in soybean than in maize. Furthermore, soybean "singletons" were found to be more frequently tandemly duplicated than "duplicates" in soybean, which may, to some extent, counteract the genome imbalance caused by gene loss. We propose that unlike in maize, in which two subgenomes were distinct prior to the allotetraploidization event and thus experienced global differences in selective constraints, in soybean, the two subgenomes were far less distinct prior to polyploidization, such that individual gene pairs, rather than subgenomes, experienced stochastic differences over longer periods of time, resulting in retention of the majority of duplicates.

Creating Order from Chaos: Epigenome Dynamics in Plants with Complex Genomes.

Flowering plants have strikingly distinct genomes, although they contain a similar suite of expressed genes. The diversity of genome structures and organization is largely due to variation in transposable elements (TEs) and whole-genome duplication (WGD) events. We review evidence that chromatin modifications and epigenetic regulation are intimately associated with TEs and likely play a role in mediating the effects of WGDs. We hypothesize that the current structure of a genome is the result of various TE bursts and WGDs and it is likely that the silencing mechanisms and the chromatin structure of a genome have been shaped by these events. This suggests that the specific mechanisms targeting chromatin modifications and epigenomic patterns may vary among different species. Many crop species have likely evolved chromatin-based mechanisms to tolerate silenced TEs near actively expressed genes. These interactions of heterochromatin and euchromatin are likely to have important roles in modulating gene expression and variability within species.

How important are transposons for plant evolution?

For decades, transposable elements have been known to produce a wide variety of changes in plant gene expression and function. This has led to the idea that transposable element activity has played a key part in adaptive plant evolution. This Review describes the kinds of changes that transposable elements can cause, discusses evidence that those changes have contributed to plant evolution and suggests future strategies for determining the extent to which these changes have in fact contributed to plant adaptation and evolution. Recent advances in genomics and phenomics for a range of plant species, particularly crops, have begun to allow the systematic assessment of these questions.

Natural antisense transcripts are significantly involved in regulation of drought stress in maize.

Natural antisense transcripts (NATs) are a prominent and complex class of regulatory RNAs. Using strand-specific RNA sequencing, we identified 1769 sense and antisense transcript pairs (NAT pairs) in two maize inbreds with different sensitivity to drought, as well as in two derivative recombination inbred lines (RILs). A significantly higher proportion of NATs relative to non-NATs are specifically expressed under water stress (WS). Surprisingly, expression of sense and antisense transcripts produced by NAT pairs is significantly correlated, particularly under WS. We found an unexpected large proportion of NATs with protein coding potential, as estimated by ribosome release scores. Small RNAs significantly accumulate within NAT pairs, with 21 nt smRNA particularly enriched in overlapping regions of these pairs of genes. The abundance of these smRNAs is significantly altered in the leafbladeless1 mutant, suggesting that these genes may be regulated by the tasiRNA pathway. Further, NATs are significantly hypomethylated and include fewer transposable element sequences relative to non-NAT genes. NAT gene regions also exhibit higher levels of H3K36me3, H3K9ac, and H3K4me3, but lower levels of H3K27me3, indicating that NAT gene pairs generally exhibit an open chromatin configuration. Finally, NAT pairs in 368 diverse maize inbreds and 19 segregating populations were specifically enriched for polymorphisms associated with drought tolerance. Taken together, the data highlight the potential impact of that small RNAs and histone modifications have in regulation of NAT expression, and the significance of NATs in response to WS.

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.

Origin, inheritance, and gene regulatory consequences of genome dominance in polyploids.

Whole-genome duplications happen repeatedly in a typical flowering plant lineage. Following most ancient tetraploidies, the two subgenomes are distinguishable because one subgenome, the dominant subgenome, tends to have more genes than the other subgenome. Additionally, among retained pairs, the gene on the dominant subgenome tends to be expressed more than its recessive homeolog. Using comparative genomics, we show that genome dominance is heritable. The dominant subgenome of one postpolyploidy event remains dominant through a subsequent polyploidy event. We show that transposon-derived 24-nt RNAs target and cover the upstream region of retained genes preferentially when located on the recessive subgenome, and with little regard for a gene's level of expression. We hypothesize that small RNA (smRNA)-mediated silencing of transposons near genes causes position-effect down-regulation. Unlike 24-nt smRNA coverage, transposon coverage tracks gene expression, so not all transposons behave identically. We propose that successful ancient tetraploids begin as wide crosses between two lines, each evolved for different tradeoffs between transposon silencing and negative position effects on gene expression. We hypothesize that following a chaotic wide-cross/new tetraploid period, genes acquire their new expression balances based on differences in transposon coverage in the parents. We envision patches of silenceable transposon as quantitative cis-regulators of baseline transcription rate. Attractive solutions to heterosis and the C-value paradox are mentioned.

Pack-Mutator-like transposable elements (Pack-MULEs) induce directional modification of genes through biased insertion and DNA acquisition.

In monocots, many genes demonstrate a significant negative GC gradient, meaning that the GC content declines along the orientation of transcription. Such a gradient is not observed in the genes of the dicot plant Arabidopsis. In addition, a lack of homology is often observed when comparing the 5' end of the coding region of orthologous genes in rice and Arabidopsis. The reasons for these differences have been enigmatic. The presence of GC-rich sequences at the 5' end of genes may influence the conformation of chromatin, the expression level of genes, as well as the recombination rate. Here we show that Pack-Mutator-like transposable elements (Pack-MULEs) that carry gene fragments specifically acquire GC-rich fragments and preferentially insert into the 5' end of genes. The resulting Pack-MULEs form independent, GC-rich transcripts with a negative GC gradient. Alternatively, the Pack-MULEs evolve into additional exons at the 5' end of existing genes, thus altering the GC content in those regions. We demonstrate that Pack-MULEs modify the 5' end of genes and are at least partially responsible for the negative GC gradient of genes in grasses. Such a unique and global impact on gene composition and gene structure has not been observed for any other transposable elements.

Heritable transposon silencing initiated by a naturally occurring transposon inverted duplication.

It has been suggested that gene silencing evolved as a defense against genomic parasites such as transposons. This idea is based on analysis of mutations that reactivate transposons that are stably silenced: they affect maintenance rather than initiation of silencing. Here we describe the cloning and characterization of a naturally occurring locus able to heritably silence the otherwise highly active MuDR transposon in maize. This locus, Mu killer (Muk), results from the inverted duplication of a partially deleted autonomous MuDR element located at the breakpoint of a genomic deletion. Muk produces a hybrid hairpin transcript that is processed into small RNAs, which are amplified when the target MuDR transcript is present. Muk provides the first example of a naturally occurring transposon derivative capable of initiating the heritable silencing of an active transposon family. Further, transposon-generated inverted duplications may be important for the generation of double-stranded RNAs used in gene silencing.

Following tetraploidy in maize, a short deletion mechanism removed genes preferentially from one of the two homologs.

Previous work in Arabidopsis showed that after an ancient tetraploidy event, genes were preferentially removed from one of the two homologs, a process known as fractionation. The mechanism of fractionation is unknown. We sought to determine whether such preferential, or biased, fractionation exists in maize and, if so, whether a specific mechanism could be implicated in this process. We studied the process of fractionation using two recently sequenced grass species: sorghum and maize. The maize lineage has experienced a tetraploidy since its divergence from sorghum approximately 12 million years ago, and fragments of many knocked-out genes retain enough sequence similarity to be easily identifiable. Using sorghum exons as the query sequence, we studied the fate of both orthologous genes in maize following the maize tetraploidy. We show that genes are predominantly lost, not relocated, and that single-gene loss by deletion is the rule. Based on comparisons with orthologous sorghum and rice genes, we also infer that the sequences present before the deletion events were flanked by short direct repeats, a signature of intra-chromosomal recombination. Evidence of this deletion mechanism is found 2.3 times more frequently on one of the maize homologs, consistent with earlier observations of biased fractionation. The over-fractionated homolog is also a greater than 3-fold better target for transposon removal, but does not have an observably higher synonymous base substitution rate, nor could we find differentially placed methylation domains. We conclude that fractionation is indeed biased in maize and that intra-chromosomal or possibly a similar illegitimate recombination is the primary mechanism by which fractionation occurs. The mechanism of intra-chromosomal recombination explains the observed bias in both gene and transposon loss in the maize lineage. The existence of fractionation bias demonstrates that the frequency of deletion is modulated. Among the evolutionary benefits of this deletion/fractionation mechanism is bulk DNA removal and the generation of novel combinations of regulatory sequences and coding regions.

Epigenetic reprogramming during vegetative phase change in maize.

An important step during plant development is the transition from juvenile to adult growth. It is only after this transition that plants are reproductively competent. Given the great danger that transposon activity represents to the germ line, this may also be an important period during development with respect to transposon regulation and silencing. We demonstrate that a change in expression of a key component of the RNA silencing pathway is associated with both vegetative phase change and shifts in epigenetic regulation of a maize transposon.

Heritable changes of epialleles in maize can be triggered in the absence of DNA methylation.

Trans-chromosomal interactions resulting in changes in DNA methylation during hybridization have been observed in several plant species. However, very little is known about the causes or consequences of these interactions. Here, we compared DNA methylomes of F1 hybrids that are mutant for a small RNA biogenesis gene, Mop1 (mediator of paramutation1) with that of their parents, wild type siblings, and backcrossed progeny in maize. Our data show that hybridization triggers global changes in both trans-chromosomal methylation (TCM) and trans-chromosomal demethylation (TCdM), most of which involved changes in CHH methylation. In more than 60% of these TCM differentially methylated regions (DMRs) in which small RNAs are available, no significant changes in the quantity of small RNAs were observed. Methylation at the CHH TCM DMRs was largely lost in the mop1 mutant, although the effects of this mutant varied depending on the location of the CHH DMRs. Interestingly, an increase in CHH at TCM DMRs was associated with enhanced expression of a subset of highly expressed genes and suppressed expression of a small number of lowly expressed genes. Examination of the methylation levels in backcrossed plants demonstrates that TCM and TCdM can be maintained in the subsequent generation, but that TCdM is more stable than TCM. Surprisingly, although increased CHH methylation in F1 plants did require Mop1, initiation of the changes in the epigenetic state of TCM DMRs did not require a functional copy of this gene, suggesting that initiation of these changes is not dependent on RNA-directed DNA methylation.

Production and processing of siRNA precursor transcripts from the highly repetitive maize genome.

Mutations affecting the maintenance of heritable epigenetic states in maize identify multiple RNA-directed DNA methylation (RdDM) factors including RMR1, a novel member of a plant-specific clade of Snf2-related proteins. Here we show that RMR1 is necessary for the accumulation of a majority of 24 nt small RNAs, including those derived from Long-Terminal Repeat (LTR) retrotransposons, the most common repetitive feature in the maize genome. A genetic analysis of DNA transposon repression indicates that RMR1 acts upstream of the RNA-dependent RNA polymerase, RDR2 (MOP1). Surprisingly, we show that non-polyadenylated transcripts from a sampling of LTR retrotransposons are lost in both rmr1 and rdr2 mutants. In contrast, plants deficient for RNA Polymerase IV (Pol IV) function show an increase in polyadenylated LTR RNA transcripts. These findings support a model in which Pol IV functions independently of the small RNA accumulation facilitated by RMR1 and RDR2 and support that a loss of Pol IV leads to RNA Polymerase II-based transcription. Additionally, the lack of changes in general genome homeostasis in rmr1 mutants, despite the global loss of 24 nt small RNAs, challenges the perceived roles of siRNAs in maintaining functional heterochromatin in the genomes of outcrossing grass species.