Extending the Genotype in Brachypodium by Including DNA Methylation Reveals a Joint Contribution with Genetics on Adaptive Traits

Epigenomic changes have been considered a potential missing link underlying phenotypic variation in quantitative traits but is potentially confounded with the underlying DNA sequence variation. Although the concept of epigenetic inheritance has been discussed in depth, there have been few studies attempting to directly dissect the amount of epigenomic variation within inbred natural populations while also accounting for genetic diversity. By using known genetic relationships between Brachypodium lines, multiple sets of nearly identical accession families were selected for phenotypic studies and DNA methylome profiling to investigate the dual role of (epi)genetics under simulated natural seasonal climate conditions. Despite reduced genetic diversity, appreciable phenotypic variation was still observable in the measured traits (height, leaf width and length, tiller count, flowering time, ear count) between as well as within the inbred accessions. However, with reduced genetic diversity there was diminished variation in DNA methylation within families. Mixed-effects linear modeling revealed large genetic differences between families and a minor contribution of DNA methylation variation on phenotypic variation in select traits. Taken together, this analysis suggests a limited but significant contribution of DNA methylation toward heritable phenotypic variation relative to genetic differences.

Excess Light Priming in Arabidopsis thaliana Genotypes with Altered DNA Methylomes

Plants must continuously react to the ever-fluctuating nature of their environment. Repeated exposure to stressful conditions can lead to priming, whereby prior encounters heighten a plant's ability to respond to future events. A clear example of priming is provided by the model plant Arabidopsis thaliana (Arabidopsis), in which photosynthetic and photoprotective responses are enhanced following recurring light stress. While there are various post-translational mechanisms underpinning photoprotection, an unresolved question is the relative importance of transcriptional changes toward stress priming and, consequently, the potential contribution from DNA methylation - a heritable chemical modification of DNA capable of influencing gene expression. Here, we systematically investigate the potential molecular underpinnings of physiological priming against recurring excess-light (EL), specifically DNA methylation and transcriptional regulation: the latter having not been examined with respect to EL priming. The capacity for physiological priming of photosynthetic and photoprotective parameters following a recurring EL treatment was not impaired in Arabidopsis mutants with perturbed establishment, maintenance, or removal of DNA methylation. Importantly, no differences in development or basal photoprotective capacity were identified in the mutants that may confound the above result. Little evidence for a causal transcriptional component of physiological priming was identified; in fact, most alterations in primed plants presented as a transcriptional 'dampening' in response to an additional EL exposure, likely a consequence of physiological priming. However, a set of transcripts uniquely regulated in primed plants provide preliminary evidence for a novel transcriptional component of recurring EL priming, independent of physiological changes. Thus, we propose that physiological priming of recurring EL in Arabidopsis occurs independently of DNA methylation; and that the majority of the associated transcriptional alterations are a consequence, not cause, of this physiological priming.

HOME: a histogram based machine learning approach for effective identification of differentially methylated regions

Background

The development of whole genome bisulfite sequencing has made it possible to identify methylation differences at single base resolution throughout an entire genome. However, a persistent challenge in DNA methylome analysis is the accurate identification of differentially methylated regions (DMRs) between samples. Sensitive and specific identification of DMRs among different conditions requires accurate and efficient algorithms, and while various tools have been developed to tackle this problem, they frequently suffer from inaccurate DMR boundary identification and high false positive rate.

Results

We present a novel Histogram Of MEthylation (HOME) based method that takes into account the inherent difference in the distribution of methylation levels between DMRs and non-DMRs to discriminate between the two using a Support Vector Machine. We show that generated features used by HOME are dataset-independent such that a classifier trained on, for example, a mouse methylome training set of regions of differentially accessible chromatin, can be applied to any other organism’s dataset and identify accurate DMRs. We demonstrate that DMRs identified by HOME exhibit higher association with biologically relevant genes, processes, and regulatory events compared to the existing methods. Moreover, HOME provides additional functionalities lacking in most of the current DMR finders such as DMR identification in non-CG context and time series analysis. HOME is freely available at https://github.com/ListerLab/HOME.

Conclusion

HOME produces more accurate DMRs than the current state-of-the-art methods on both simulated and biological datasets. The broad applicability of HOME to identify accurate DMRs in genomic data from any organism will have a significant impact upon expanding our knowledge of how DNA methylation dynamics affect cell development and differentiation.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2845-y) contains supplementary material, which is available to authorized users.

RNA Polymerase II Read-Through Promotes Expression of Neighboring Genes in SAL1-PAP-XRN Retrograde Signaling

In plants, the molecular function(s) of the nucleus-localized 5'-3' EXORIBONUCLEASES (XRNs) are unclear; however, their activity is reported to have a significant effect on gene expression and SAL1-mediated retrograde signaling. Using parallel analysis of RNA ends, we documented a dramatic increase in uncapped RNA substrates of the XRNs in both sal1 and xrn2xrn3 mutants. We found that a major consequence of reducing SAL1 or XRN activity was RNA Polymerase II 3' read-through. This occurred at 72% of expressed genes, demonstrating a major genome-wide role for the XRN-torpedo model of transcription termination in Arabidopsis (Arabidopsis thaliana). Read-through is speculated to have a negative effect on transcript abundance; however, we did not observe this. Rather, we identified a strong association between read-through and increased transcript abundance of tandemly orientated downstream genes, strongly correlated with the proximity (less than 1,000 bp) and expression of the upstream gene. We observed read-through in the proximity of 903 genes up-regulated in the sal1-8 retrograde signaling mutant; thus, this phenomenon may account directly for up to 23% of genes up-regulated in sal1-8 Using APX2 and AT5G43770 as exemplars, we genetically uncoupled read-through loci from downstream genes to validate the principle of read-through-mediated mRNA regulation, providing one mechanism by which an ostensibly posttranscriptional exoribonuclease that targets uncapped RNAs could modulate gene expression.

Maintenance of pre-existing DNA methylation states through recurring excess-light stress

The capacity for plant stress priming and memory and the notion of this being underpinned by DNA methylation-mediated memory is an appealing hypothesis for which there is mixed evidence. We previously established a lack of drought-induced methylome variation in Arabidopsis thaliana (Arabidopsis); however, this was tied to only minor observations of physiological memory. There are numerous independent observations demonstrating that photoprotective mechanisms, induced by excess-light stress, can lead to robust programmable changes in newly developing leaf tissues. Although key signalling molecules and transcription factors are known to promote this priming signal, an untested question is the potential involvement of chromatin marks towards the maintenance of light stress acclimation, or memory. Thus, we systematically tested our previous hypothesis of a stress-resistant methylome using a recurring excess-light stress, then analysing new, emerging, and existing tissues. The DNA methylome showed negligible stress-associated variation, with the vast majority attributable to stochastic differences. Yet, photoacclimation was evident through enhanced photosystem II performance in exposed tissues, and nonphotochemical quenching and fluorescence decline ratio showed evidence of mitotic transmission. Thus, we have observed physiological acclimation in new and emerging tissues in the absence of substantive DNA methylome changes.

The Arabidopsis DNA Methylome Is Stable under Transgenerational Drought Stress

Improving the responsiveness, acclimation, and memory of plants to abiotic stress holds substantive potential for improving agriculture. An unresolved question is the involvement of chromatin marks in the memory of agriculturally relevant stresses. Such potential has spurred numerous investigations yielding both promising and conflicting results. Consequently, it remains unclear to what extent robust stress-induced DNA methylation variation can underpin stress memory. Using a slow-onset water deprivation treatment in Arabidopsis (Arabidopsis thaliana), we investigated the malleability of the DNA methylome to drought stress within a generation and under repeated drought stress over five successive generations. While drought-associated epi-alleles in the methylome were detected within a generation, they did not correlate with drought-responsive gene expression. Six traits were analyzed for transgenerational stress memory, and the descendants of drought-stressed lineages showed one case of memory in the form of increased seed dormancy, and that persisted one generation removed from stress. With respect to transgenerational drought stress, there were negligible conserved differentially methylated regions in drought-exposed lineages compared with unstressed lineages. Instead, the majority of observed variation was tied to stochastic or preexisting differences in the epigenome occurring at repetitive regions of the Arabidopsis genome. Furthermore, the experience of repeated drought stress was not observed to influence transgenerational epi-allele accumulation. Our findings demonstrate that, while transgenerational memory is observed in one of six traits examined, they are not associated with causative changes in the DNA methylome, which appears relatively impervious to drought stress.

Rapid Recovery Gene Downregulation during Excess-Light Stress and Recovery in Arabidopsis

Stress recovery may prove to be a promising approach to increase plant performance and, theoretically, mRNA instability may facilitate faster recovery. Transcriptome (RNA-seq, qPCR, sRNA-seq, and PARE) and methylome profiling during repeated excess-light stress and recovery was performed at intervals as short as 3 min. We demonstrate that 87% of the stress-upregulated mRNAs analyzed exhibit very rapid recovery. For instance, HSP101 abundance declined 2-fold every 5.1 min. We term this phenomenon rapid recovery gene downregulation (RRGD), whereby mRNA abundance rapidly decreases promoting transcriptome resetting. Decay constants (k) were modeled using two strategies, linear and nonlinear least squares regressions, with the latter accounting for both transcription and degradation. This revealed extremely short half-lives ranging from 2.7 to 60.0 min for 222 genes. Ribosome footprinting using degradome data demonstrated RRGD loci undergo cotranslational decay and identified changes in the ribosome stalling index during stress and recovery. However, small RNAs and 5'-3' RNA decay were not essential for recovery of the transcripts examined, nor were any of the six excess light-associated methylome changes. We observed recovery-specific gene expression networks upon return to favorable conditions and six transcriptional memory types. In summary, rapid transcriptome resetting is reported in the context of active recovery and cellular memory.

Population scale mapping of transposable element diversity reveals links to gene regulation and epigenomic variation

Variation in the presence or absence of transposable elements (TEs) is a major source of genetic variation between individuals. Here, we identified 23,095 TE presence/absence variants between 216 Arabidopsis accessions. Most TE variants were rare, and we find these rare variants associated with local extremes of gene expression and DNA methylation levels within the population. Of the common alleles identified, two thirds were not in linkage disequilibrium with nearby SNPs, implicating these variants as a source of novel genetic diversity. Many common TE variants were associated with significantly altered expression of nearby genes, and a major fraction of inter-accession DNA methylation differences were associated with nearby TE insertions. Overall, this demonstrates that TE variants are a rich source of genetic diversity that likely plays an important role in facilitating epigenomic and transcriptional differences between individuals, and indicates a strong genetic basis for epigenetic variation.

DOI:http://dx.doi.org/10.7554/eLife.20777.001

DNA methylation profiles of diverse Brachypodium distachyon align with underlying genetic diversity

DNA methylation, a common modification of genomic DNA, is known to influence the expression of transposable elements as well as some genes. Although commonly viewed as an epigenetic mark, evidence has shown that underlying genetic variation, such as transposable element polymorphisms, often associate with differential DNA methylation states. To investigate the role of DNA methylation variation, transposable element polymorphism, and genomic diversity, whole-genome bisulfite sequencing was performed on genetically diverse lines of the model cereal Brachypodium distachyon. Although DNA methylation profiles are broadly similar, thousands of differentially methylated regions are observed between lines. An analysis of novel transposable element indel variation highlighted hundreds of new polymorphisms not seen in the reference sequence. DNA methylation and transposable element variation is correlated with the genome-wide amount of genetic variation present between samples. However, there was minimal evidence that novel transposon insertions or deletions are associated with nearby differential methylation. This study highlights unique relationships between genetic variation and DNA methylation variation within Brachypodium and provides a valuable map of DNA methylation across diverse resequenced accessions of this model cereal species.

Reconsidering plant memory: Intersections between stress recovery, RNA turnover, and epigenetics

Plants grow in dynamic environments where they can be exposed to a multitude of stressful factors, all of which affect their development, yield, and, ultimately, reproductive success. Plants are adept at rapidly acclimating to stressful conditions and are able to further fortify their defenses by retaining memories of stress to enable stronger or more rapid responses should an environmental perturbation recur. Indeed, one mechanism that is often evoked regarding environmental memories is epigenetics. Yet, there are relatively few examples of such memories; neither is there a clear understanding of their duration, considering the plethora of stresses in nature. We propose that this field would benefit from investigations into the processes and mechanisms enabling recovery from stress. An understanding of stress recovery could provide fresh insights into when, how, and why environmental memories are created and regulated. Stress memories may be maladaptive, hindering recovery and affecting development and potential yield. In some circumstances, it may be advantageous for plants to learn to forget. Accordingly, the recovery process entails a balancing act between resetting and memory formation. During recovery, RNA metabolism, posttranscriptional gene silencing, and RNA-directed DNA methylation have the potential to play key roles in resetting the epigenome and transcriptome and in altering memory. Exploration of this emerging area of research is becoming ever more tractable with advances in genomics, phenomics, and high-throughput sequencing methodology that will enable unprecedented profiling of high-resolution stress recovery time series experiments and sampling of large natural populations.

Examining the Causes and Consequences of Context-Specific Differential DNA Methylation in Maize

DNA methylation is a stable modification of chromatin that can contribute to epigenetic variation through the regulation of genes or transposons. Profiling of DNA methylation in five maize (Zea mays) inbred lines found that while DNA methylation levels for more than 99% of the analyzed genomic regions are similar, there are still 5,000 to 20,000 context-specific differentially methylated regions (DMRs) between any two genotypes. The analysis of identical-by-state genomic regions that have limited genetic variation provided evidence that DMRs can occur without local sequence variation, but they are less common than in regions with genetic variation. Characterization of the sequence specificity of DMRs, location of DMRs relative to genes and transposons, and patterns of DNA methylation in regions flanking DMRs reveals a distinct subset of DMRs. Transcriptome profiling of the same tissue revealed that only approximately 20% of genes with qualitative (on-off) differences in gene expression are associated with DMRs, and there is little evidence for association of DMRs with genes that show quantitative differences in gene expression. We also identify a set of genes that may represent cryptic information that is silenced by DNA methylation in the reference B73 genome. Many of these genes exhibit natural variation in other genotypes, suggesting the potential for selection to act upon existing epigenetic natural variation. This study provides insights into the origin and influences of DMRs in a crop species with a complex genome organization.

Post-conversion targeted capture of modified cytosines in mammalian and plant genomes

We present a capture-based approach for bisulfite-converted DNA that allows interrogation of pre-defined genomic locations, allowing quantitative and qualitative assessments of 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC) at CG dinucleotides and in non-CG contexts (CHG, CHH) in mammalian and plant genomes. We show the technique works robustly and reproducibly using as little as 500 ng of starting DNA, with results correlating well with whole genome bisulfite sequencing data, and demonstrate that human DNA can be tested in samples contaminated with microbial DNA. This targeting approach will allow cell type-specific designs to maximize the value of 5mC and 5hmC sequencing.

Minimal evidence for consistent changes in maize DNA methylation patterns following environmental stress

DNA methylation is a chromatin modification that is sometimes associated with epigenetic regulation of gene expression. As DNA methylation can be reversible at some loci, it is possible that methylation patterns may change within an organism that is subjected to environmental stress. In order to assess the effects of abiotic stress on DNA methylation patterns in maize (Zea mays), seeding plants were subjected to heat, cold, and UV stress treatments. Tissue was later collected from individual adult plants that had been subjected to stress or control treatments and used to perform DNA methylation profiling to determine whether there were consistent changes in DNA methylation triggered by specific stress treatments. DNA methylation profiling was performed by immunoprecipitation of methylated DNA followed by microarray hybridization to allow for quantitative estimates of DNA methylation abundance throughout the low-copy portion of the maize genome. By comparing the DNA methylation profiles of each individual plant to the average of the control plants it was possible to identify regions of the genome with variable DNA methylation. However, we did not find evidence of consistent DNA methylation changes resulting from the stress treatments used in this study. Instead, the data suggest that there is a low-rate of stochastic variation that is present in both control and stressed plants.

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.

Genomic distribution of H3K9me2 and DNA methylation in a maize genome

DNA methylation and dimethylation of lysine 9 of histone H3 (H3K9me2) are two chromatin modifications that can be associated with gene expression or recombination rate. The maize genome provides a complex landscape of interspersed genes and transposons. The genome-wide distribution of DNA methylation and H3K9me2 were investigated in seedling tissue for the maize inbred B73 and compared to patterns of these modifications observed in Arabidopsis thaliana. Most maize transposons are highly enriched for DNA methylation in CG and CHG contexts and for H3K9me2. In contrast to findings in Arabidopsis, maize CHH levels in transposons are generally low but some sub-families of transposons are enriched for CHH methylation and these families exhibit low levels of H3K9me2. The profile of modifications over genes reveals that DNA methylation and H3K9me2 is quite low near the beginning and end of genes. Although elevated CG and CHG methylation are found within gene bodies, CHH and H3K9me2 remain low. Maize has much higher levels of CHG methylation within gene bodies than observed in Arabidopsis and this is partially attributable to the presence of transposons within introns for some maize genes. These transposons are associated with high levels of CHG methylation and H3K9me2 but do not appear to prevent transcriptional elongation. Although the general trend is for a strong depletion of H3K9me2 and CHG near the transcription start site there are some putative genes that have high levels of these chromatin modifications. This study provides a clear view of the relationship between DNA methylation and H3K9me2 in the maize genome and how the distribution of these modifications is shaped by the interplay of genes and transposons.

Epigenetics: Beyond Chromatin Modifications and Complex Genetic Regulation

Chromatin modifications and epigenetics may play important roles in many plant processes, including developmental regulation, responses to environmental stimuli, and local adaptation. Chromatin modifications describe biochemical changes to chromatin state, such as alterations in the specific type or placement of histones, modifications of DNA or histones, or changes in the specific proteins or RNAs that associate with a genomic region. The term epigenetic is often used to describe a variety of unexpected patterns of gene regulation or inheritance. Here, we specifically define epigenetics to include the key aspects of heritability (stable transmission of gene expression states through mitotic or meiotic cell divisions) and independence from DNA sequence changes. We argue against generically equating chromatin and epigenetics; although many examples of epigenetics involve chromatin changes, those chromatin changes are not always heritable or may be influenced by genetic changes. Careful use of the terms chromatin modifications and epigenetics can help separate the biochemical mechanisms of regulation from the inheritance patterns of altered chromatin states. Here, we also highlight examples in which chromatin modifications and epigenetics affect important plant processes.

Genome-wide discovery and characterization of maize long non-coding RNAs

BACKGROUND

Long non-coding RNAs (lncRNAs) are transcripts that are 200 bp or longer, do not encode proteins, and potentially play important roles in eukaryotic gene regulation. However, the number, characteristics and expression inheritance pattern of lncRNAs in maize are still largely unknown.

RESULTS

By exploiting available public EST databases, maize whole genome sequence annotation and RNA-seq datasets from 30 different experiments, we identified 20,163 putative lncRNAs. Of these lncRNAs, more than 90% are predicted to be the precursors of small RNAs, while 1,704 are considered to be high-confidence lncRNAs. High confidence lncRNAs have an average transcript length of 463 bp and genes encoding them contain fewer exons than annotated genes. By analyzing the expression pattern of these lncRNAs in 13 distinct tissues and 105 maize recombinant inbred lines, we show that more than 50% of the high confidence lncRNAs are expressed in a tissue-specific manner, a result that is supported by epigenetic marks. Intriguingly, the inheritance of lncRNA expression patterns in 105 recombinant inbred lines reveals apparent transgressive segregation, and maize lncRNAs are less affected by cis- than by trans-genetic factors.

CONCLUSIONS

We integrate all available transcriptomic datasets to identify a comprehensive set of maize lncRNAs, provide a unique annotation resource of the maize genome and a genome-wide characterization of maize lncRNAs, and explore the genetic control of their expression using expression quantitative trait locus mapping.

Comprehensive analysis of imprinted genes in maize reveals allelic variation for imprinting and limited conservation with other species

In plants, a subset of genes exhibit imprinting in endosperm tissue such that expression is primarily from the maternal or paternal allele. Imprinting may arise as a consequence of mechanisms for silencing of transposons during reproduction, and in some cases imprinted expression of particular genes may provide a selective advantage such that it is conserved across species. Separate mechanisms for the origin of imprinted expression patterns and maintenance of these patterns may result in substantial variation in the targets of imprinting in different species. Here we present deep sequencing of RNAs isolated from reciprocal crosses of four diverse maize genotypes, providing a comprehensive analysis that allows evaluation of imprinting at more than 95% of endosperm-expressed genes. We find that over 500 genes exhibit statistically significant parent-of-origin effects in maize endosperm tissue, but focused our analyses on a subset of these genes that had >90% expression from the maternal allele (69 genes) or from the paternal allele (108 genes) in at least one reciprocal cross. Over 10% of imprinted genes show evidence of allelic variation for imprinting. A comparison of imprinting in maize and rice reveals that 13% of genes with syntenic orthologs in both species exhibit conserved imprinting. Genes that exhibit conserved imprinting between maize and rice have elevated nonsynonymous to synonymous substitution ratios compared with other imprinted genes, suggesting a history of more rapid evolution. Together, these data suggest that imprinting only has functional relevance at a subset of loci that currently exhibit imprinting in maize.

Epigenetic and genetic influences on DNA methylation variation in maize populations

DNA methylation is a chromatin modification that is frequently associated with epigenetic regulation in plants and mammals. However, genetic changes such as transposon insertions can also lead to changes in DNA methylation. Genome-wide profiles of DNA methylation for 20 maize (Zea mays) inbred lines were used to discover differentially methylated regions (DMRs). The methylation level for each of these DMRs was also assayed in 31 additional maize or teosinte genotypes, resulting in the discovery of 1966 common DMRs and 1754 rare DMRs. Analysis of recombinant inbred lines provides evidence that the majority of DMRs are heritable. A local association scan found that nearly half of the DMRs with common variation are significantly associated with single nucleotide polymorphisms found within or near the DMR. Many of the DMRs that are significantly associated with local genetic variation are found near transposable elements that may contribute to the variation in DNA methylation. Analysis of gene expression in the same samples used for DNA methylation profiling identified over 300 genes with expression patterns that are significantly associated with DNA methylation variation. Collectively, our results suggest that DNA methylation variation is influenced by genetic and epigenetic changes that are often stably inherited and can influence the expression of nearby genes.

Distribution, functional impact, and origin mechanisms of copy number variation in the barley genome

Background

There is growing evidence for the prevalence of copy number variation (CNV) and its role in phenotypic variation in many eukaryotic species. Here we use array comparative genomic hybridization to explore the extent of this type of structural variation in domesticated barley cultivars and wild barleys.

Results

A collection of 14 barley genotypes including eight cultivars and six wild barleys were used for comparative genomic hybridization. CNV affects 14.9% of all the sequences that were assessed. Higher levels of CNV diversity are present in the wild accessions relative to cultivated barley. CNVs are enriched near the ends of all chromosomes except 4H, which exhibits the lowest frequency of CNVs. CNV affects 9.5% of the coding sequences represented on the array and the genes affected by CNV are enriched for sequences annotated as disease-resistance proteins and protein kinases. Sequence-based comparisons of CNV between cultivars Barke and Morex provided evidence that DNA repair mechanisms of double-strand breaks via single-stranded annealing and synthesis-dependent strand annealing play an important role in the origin of CNV in barley.

Conclusions

We present the first catalog of CNVs in a diploid Triticeae species, which opens the door for future genome diversity research in a tribe that comprises the economically important cereal species wheat, barley, and rye. Our findings constitute a valuable resource for the identification of CNV affecting genes of agronomic importance. We also identify potential mechanisms that can generate variation in copy number in plant genomes.

Genomic distribution of maize facultative heterochromatin marked by trimethylation of H3K27

Trimethylation of histone H3 Lys-27 (H3K27me3) plays a critical role in regulating gene expression during plant and animal development. We characterized the genome-wide distribution of H3K27me3 in five developmentally distinct tissues in maize (Zea mays) plants of two genetic backgrounds, B73 and Mo17. There were more substantial differences in the genome-wide profile of H3K27me3 between different tissues than between the two genotypes. The tissue-specific patterns of H3K27me3 were often associated with differences in gene expression among the tissues and most of the imprinted genes that are expressed solely from the paternal allele in endosperm are targets of H3K27me3. A comparison of the H3K27me3 targets in rice (Oryza sativa), maize, and Arabidopsis thaliana provided evidence for conservation of the H3K27me3 targets among plant species. However, there was limited evidence for conserved targeting of H3K27me3 in the two maize subgenomes derived from whole-genome duplication, suggesting the potential for subfunctionalization of chromatin regulation of paralogs. Genomic profiling of H3K27me3 in loss-of-function mutant lines for Maize Enhancer of zeste-like2 (Mez2) and Mez3, two of the three putative H3K27me3 methyltransferases present in the maize genome, suggested partial redundancy of this gene family for maintaining H3K27me3 patterns. Only a portion of the targets of H3K27me3 required Mez2 and/or Mez3, and there was limited evidence for functional consequences of H3K27me3 at these targets.

Spreading of heterochromatin is limited to specific families of maize retrotransposons

Transposable elements (TEs) have the potential to act as controlling elements to influence the expression of genes and are often subject to heterochromatic silencing. The current paradigm suggests that heterochromatic silencing can spread beyond the borders of TEs and influence the chromatin state of neighboring low-copy sequences. This would allow TEs to condition obligatory or facilitated epialleles and act as controlling elements. The maize genome contains numerous families of class I TEs (retrotransposons) that are present in moderate to high copy numbers, and many are found in regions near genes, which provides an opportunity to test whether the spreading of heterochromatin from retrotransposons is prevalent. We have investigated the extent of heterochromatin spreading into DNA flanking each family of retrotransposons by profiling DNA methylation and di-methylation of lysine 9 of histone 3 (H3K9me2) in low-copy regions of the maize genome. The effects of different retrotransposon families on local chromatin are highly variable. Some retrotransposon families exhibit enrichment of heterochromatic marks within 800–1,200 base pairs of insertion sites, while other families exhibit very little evidence for the spreading of heterochromatic marks. The analysis of chromatin state in genotypes that lack specific insertions suggests that the heterochromatin in low-copy DNA flanking retrotransposons often results from the spreading of silencing marks rather than insertion-site preferences. Genes located near TEs that exhibit spreading of heterochromatin tend to be expressed at lower levels than other genes. Our findings suggest that a subset of retrotransposon families may act as controlling elements influencing neighboring sequences, while the majority of retrotransposons have little effect on flanking sequences.

Transposable elements comprise a substantial portion of many eukaryotic genomes. These mobile fragments of DNA can directly mutate genes through insertions into coding regions but may also affect the gene regulation through nearby insertions. There is evidence that the majority of transposable elements are epigenetically silenced, and in some cases this silencing may spread to neighboring sequences. This spreading of heterochromatin could create a significant fitness tradeoff between transposon silencing and gene expression. The maize genome has a complex organization with many genes flanked by retrotransposons, providing an opportunity to study the interaction of retrotransposons and genes. To survey the prevalence of heterochromatin spreading associated with different retrotransposon families, we profiled the spread of heterochromatin into nearby low copy sequences for 150 high copy retrotransposon families. While many retrotransposons exhibit little to no spreading of heterochromatin, there are some retrotransposon families that do exhibit spreading. Genes located near retrotransposons that spread heterochromatin have lower expression levels. The families of retrotransposons that spread heterochromatin marks to nearby low-copy sequences may have increased fitness costs for the host genome due to their suppression of genes located near insertions.

B73-Mo17 near-isogenic lines demonstrate dispersed structural variation in maize

Recombinant inbred lines developed from the maize (Zea mays ssp. mays) inbreds B73 and Mo17 have been widely used to discover quantitative trait loci controlling a wide variety of phenotypic traits and as a resource to produce high-resolution genetic maps. These two parents were used to produce a set of near-isogenic lines (NILs) with small regions of introgression into both backgrounds. A novel array-based genotyping platform was used to score genotypes of over 7,000 loci in 100 NILs with B73 as the recurrent parent and 50 NILs with Mo17 as the recurrent parent. This population contains introgressions that cover the majority of the maize genome. The set of NILs displayed an excess of residual heterozygosity relative to the amount expected based on their pedigrees, and this excess residual heterozygosity is enriched in the low-recombination regions near the centromeres. The genotyping platform provided the ability to survey copy number variants that exist in more copies in Mo17 than in B73. The majority of these Mo17-specific duplications are located in unlinked positions throughout the genome. The utility of this population for the discovery and validation of quantitative trait loci was assessed through analysis of plant height variation.

Pervasive gene content variation and copy number variation in maize and its undomesticated progenitor

Individuals of the same species are generally thought to have very similar genomes. However, there is growing evidence that structural variation in the form of copy number variation (CNV) and presence-absence variation (PAV) can lead to variation in the genome content of individuals within a species. Array comparative genomic hybridization (CGH) was used to compare gene content and copy number variation among 19 diverse maize inbreds and 14 genotypes of the wild ancestor of maize, teosinte. We identified 479 genes exhibiting higher copy number in some genotypes (UpCNV) and 3410 genes that have either fewer copies or are missing in the genome of at least one genotype relative to B73 (DownCNV/PAV). Many of these DownCNV/PAV are examples of genes present in B73, but missing from other genotypes. Over 70% of the CNV/PAV examples are identified in multiple genotypes, and the majority of events are observed in both maize and teosinte, suggesting that these variants predate domestication and that there is not strong selection acting against them. Many of the genes affected by CNV/PAV are either maize specific (thus possible annotation artifacts) or members of large gene families, suggesting that the gene loss can be tolerated through buffering by redundant functions encoded elsewhere in the genome. While this structural variation may not result in major qualitative variation due to genetic buffering, it may significantly contribute to quantitative variation.