Lexiang Ji

Widespread natural variation of DNA methylation within angiosperms

BackgroundDNA methylation is an important feature of plant epigenomes, involved in the formation of heterochromatin and affecting gene expression. Extensive variation of DNA methylation patterns within a species has been uncovered from studies of natural variation. However, the extent to which DNA methylation varies between flowering plant species is still unclear. To understand the variation in genomic patterning of DNA methylation across flowering plant species, we compared single base resolution DNA methylomes of 34 diverse angiosperm species.ResultsBy analyzing whole-genome bisulfite sequencing data in a phylogenetic context, it becomes clear that there is extensive variation throughout angiosperms in gene body DNA methylation, euchromatic silencing of transposons and repeats, as well as silencing of heterochromatic transposons. The Brassicaceae have reduced CHG methylation levels and also reduced or loss of CG gene body methylation. The Poaceae are characterized by a lack or reduction of heterochromatic CHH methylation and enrichment of CHH methylation in genic regions. Furthermore, low levels of CHH methylation are observed in a number of species, especially in clonally propagated species.ConclusionsThese results reveal the extent of variation in DNA methylation in angiosperms and show that DNA methylation patterns are broadly a reflection of the evolutionary and life histories of plant species.

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.

On the origin and evolutionary consequences of gene body DNA methylation

Significance DNA methylation in plants is found at CG, CHG, and CHH sequence contexts. In plants, CG DNA methylation is enriched in the transcribed regions of many constitutively expressed genes (gene body methylation; gbM) and shows correlations with several chromatin modifications. Contrary to other types of DNA methylation, the evolution and function of gbM are largely unknown. Here we show two independent concomitant losses of the DNA methyltransferase CHROMOMETHYLASE 3 (CMT3) and gbM without the predicted disruption of transcription and of modifications to chromatin. This result suggests that CMT3 is required for the establishment of gbM in actively transcribed genes, and that gbM is dispensable for normal transcription as well as for the composition and modification of plant chromatin. In plants, CG DNA methylation is prevalent in the transcribed regions of many constitutively expressed genes (gene body methylation; gbM), but the origin and function of gbM remain unknown. Here we report the discovery that Eutrema salsugineum has lost gbM from its genome, to our knowledge the first instance for an angiosperm. Of all known DNA methyltransferases, only CHROMOMETHYLASE 3 (CMT3) is missing from E. salsugineum. Identification of an additional angiosperm, Conringia planisiliqua, which independently lost CMT3 and gbM, supports that CMT3 is required for the establishment of gbM. Detailed analyses of gene expression, the histone variant H2A.Z, and various histone modifications in E. salsugineum and in Arabidopsis thaliana epigenetic recombinant inbred lines found no evidence in support of any role for gbM in regulating transcription or affecting the composition and modification of chromatin over evolutionary timescales.

Maize death acids, 9-lipoxygenase–derived cyclopente(a)nones, display activity as cytotoxic phytoalexins and transcriptional mediators

Significance In plants, 12-oxo-phytodienoic acid (12-OPDA) and jasmonic acid are key 13-lipoxygenase-derived linolenate oxidation products termed jasmonates that regulate diverse processes in development and innate immunity. A less-studied metabolic pathway branch is generated by 9-lipoxygenase activity on linoleic acid, enabling the production of 10-oxo-11-phytoenoic acid (10-OPEA). In maize, fungal infection by southern leaf blight (Cochliobolus heterostrophus) results in the localized production of 10-OPEA, and a series of related 12- and 14-carbon cyclopente(a)nones, collectively termed “death acids” (DAs). DAs far exceed jasmonates in abundance within infected tissues, display direct phytoalexin activity against biotic agents, mediate defense gene expression, and can promote cytotoxicity resulting in cell death. Collectively DA activities are consistent with specialized local roles in plant defense. Plant damage promotes the interaction of lipoxygenases (LOXs) with fatty acids yielding 9-hydroperoxides, 13-hydroperoxides, and complex arrays of oxylipins. The action of 13-LOX on linolenic acid enables production of 12-oxo-phytodienoic acid (12-OPDA) and its downstream products, termed “jasmonates.” As signals, jasmonates have related yet distinct roles in the regulation of plant resistance against insect and pathogen attack. A similar pathway involving 9-LOX activity on linolenic and linoleic acid leads to the 12-OPDA positional isomer, 10-oxo-11-phytodienoic acid (10-OPDA) and 10-oxo-11-phytoenoic acid (10-OPEA), respectively; however, physiological roles for 9-LOX cyclopentenones have remained unclear. In developing maize (Zea mays) leaves, southern leaf blight (Cochliobolus heterostrophus) infection results in dying necrotic tissue and the localized accumulation of 10-OPEA, 10-OPDA, and a series of related 14- and 12-carbon metabolites, collectively termed “death acids.” 10-OPEA accumulation becomes wound inducible within fungal-infected tissues and at physiologically relevant concentrations acts as a phytoalexin by suppressing the growth of fungi and herbivores including Aspergillus flavus, Fusarium verticillioides, and Helicoverpa zea. Unlike previously established maize phytoalexins, 10-OPEA and 10-OPDA display significant phytotoxicity. Both 12-OPDA and 10-OPEA promote the transcription of defense genes encoding glutathione S transferases, cytochrome P450s, and pathogenesis-related proteins. In contrast, 10-OPEA only weakly promotes the accumulation of multiple protease inhibitor transcripts. Consistent with a role in dying tissue, 10-OPEA application promotes cysteine protease activation and cell death, which is inhibited by overexpression of the cysteine protease inhibitor maize cystatin-9. Unlike jasmonates, functions for 10-OPEA and associated death acids are consistent with specialized roles in local defense reactions.

The genome and methylome of a beetle with complex social behavior, Nicrophorus vespilloides (Coleoptera: Silphidae)

Testing for conserved and novel mechanisms underlying phenotypic evolution requires a diversity of genomes available for comparison spanning multiple independent lineages. For example, complex social behavior in insects has been investigated primarily with eusocial lineages, nearly all of which are Hymenoptera. If conserved genomic influences on sociality do exist, we need data from a wider range of taxa that also vary in their levels of sociality. Here, we present the assembled and annotated genome of the subsocial beetle Nicrophorus vespilloides, a species long used to investigate evolutionary questions of complex social behavior. We used this genome to address two questions. First, do aspects of life history, such as using a carcass to breed, predict overlap in gene models more strongly than phylogeny? We found that the overlap in gene models was similar between N. vespilloides and all other insect groups regardless of life history. Second, like other insects with highly developed social behavior but unlike other beetles, does N. vespilloides have DNA methylation? We found strong evidence for an active DNA methylation system. The distribution of methylation was similar to other insects with exons having the most methylated CpGs. Methylation status appears highly conserved; 85% of the methylated genes in N. vespilloides are also methylated in the hymentopteran Nasonia vitripennis. The addition of this genome adds a coleopteran resource to answer questions about the evolution and mechanistic basis of sociality and to address questions about the potential role of methylation in social behavior.

Methylated DNA is over-represented in whole-genome bisulfite sequencing data

The development of whole-genome bisulfite sequencing (WGBS) has resulted in a number of exciting discoveries about the role of DNA methylation leading to a plethora of novel testable hypotheses. Methods for constructing sodium bisulfite-converted and amplified libraries have recently advanced to the point that the bottleneck for experiments that use WGBS has shifted to data analysis and interpretation. Here we present empirical evidence for an over-representation of reads from methylated DNA in WGBS. This enrichment for methylated DNA is exacerbated by higher cycles of PCR and is influenced by the type of uracil-insensitive DNA polymerase used for amplifying the sequencing library. Future efforts to computationally correct for this enrichment bias will be essential to increasing the accuracy of determining methylation levels for individual cytosines. It is especially critical for studies that seek to accurately quantify DNA methylation levels in populations that may segregate for allelic DNA methylation states.

The evolution of CHROMOMETHYLASES and gene body DNA methylation in plants

BackgroundThe evolution of gene body methylation (gbM), its origins, and its functional consequences are poorly understood. By pairing the largest collection of transcriptomes (>1000) and methylomes (77) across Viridiplantae, we provide novel insights into the evolution of gbM and its relationship to CHROMOMETHYLASE (CMT) proteins.ResultsCMTs are evolutionary conserved DNA methyltransferases in Viridiplantae. Duplication events gave rise to what are now referred to as CMT1, 2 and 3. Independent losses of CMT1, 2, and 3 in eudicots, CMT2 and ZMET in monocots and monocots/commelinids, variation in copy number, and non-neutral evolution suggests overlapping or fluid functional evolution of this gene family. DNA methylation within genes is widespread and is found in all major taxonomic groups of Viridiplantae investigated. Genes enriched with methylated CGs (mCG) were also identified in species sister to angiosperms. The proportion of genes and DNA methylation patterns associated with gbM are restricted to angiosperms with a functional CMT3 or ortholog. However, mCG-enriched genes in the gymnosperm Pinus taeda shared some similarities with gbM genes in Amborella trichopoda. Additionally, gymnosperms and ferns share a CMT homolog closely related to CMT2 and 3. Hence, the dependency of gbM on a CMT most likely extends to all angiosperms and possibly gymnosperms and ferns.ConclusionsThe resulting gene family phylogeny of CMT transcripts from the most diverse sampling of plants to date redefines our understanding of CMT evolution and its evolutionary consequences on DNA methylation. Future, functional tests of homologous and paralogous CMTs will uncover novel roles and consequences to the epigenome.

Crop Epigenomics: Identifying, Unlocking, and Harnessing Cryptic Variation in Crop Genomes

DNA methylation is a key chromatin modification in plant genomes that is meiotically and mitotically heritable, and at times is associated with gene expression and morphological variation. Benefiting from the increased availability of high-quality reference genome assemblies and methods to profile single-base resolution DNA methylation states, DNA methylomes for many crop species are available. These efforts are making it possible to begin answering crucial questions, including understanding the role of DNA methylation in developmental processes, its role in crop species evolution, and whether DNA methylation is dynamically altered and heritable in response to changes in the environment. These genome-wide maps provide evidence for the existence of silent epialleles in plant genomes which, once identified, can be targeted for reactivation leading to phenotypic variation.

Genome-wide redistribution of H3K27me3 is linked to genotoxic stress and defective growth

Significance Regulators of chromatin structure play critical roles in DNA-based processes. Lysine (K) Methyltransferase 1 (KMT1) homologs perform methylation of H3 lysine-9 and are best known for their essential role in heterochromatin formation and transcriptional silencing. Heterochromatin formation is also important for maintenance of genome stability, although the mechanisms are not well understood. We report that altered activity of Polycomb repressive complex-2 (PRC2), a histone lysine-27 methyltransferase complex, is responsible for genotoxic stress, poor growth, and defective development in KMT1-deficient mutants of Neurospora crassa. Mammalian KMT1 and PRC2 are required for development and are frequently mutated in cancer. This work provides information about the cellular consequences of KMT1 and PRC2 deficiency and provides insights into the regulatory and functional relationships of these conserved enzymes. H3K9 methylation directs heterochromatin formation by recruiting multiple heterochromatin protein 1 (HP1)-containing complexes that deacetylate histones and methylate cytosine bases in DNA. In Neurospora crassa, a single H3K9 methyltransferase complex, called the DIM-5,-7,-9, CUL4, DDB1 Complex (DCDC), is required for normal growth and development. DCDC-deficient mutants are hypersensitive to the genotoxic agent methyl methanesulfonate (MMS), but the molecular basis of genotoxic stress is unclear. We found that both the MMS sensitivity and growth phenotypes of DCDC-deficient strains are suppressed by mutation of embryonic ectoderm development or Su-(var)3-9; E(z); Trithorax (set)-7, encoding components of the H3K27 methyltransferase Polycomb repressive complex-2 (PRC2). Trimethylated histone H3K27 (H3K27me3) undergoes genome-wide redistribution to constitutive heterochromatin in DCDC- or HP1-deficient mutants, and introduction of an H3K27 missense mutation is sufficient to rescue phenotypes of DCDC-deficient strains. Accumulation of H3K27me3 in heterochromatin does not compensate for silencing; rather, strains deficient for both DCDC and PRC2 exhibit synthetic sensitivity to the topoisomerase I inhibitor Camptothecin and accumulate γH2A at heterochromatin. Together, these data suggest that PRC2 modulates the response to genotoxic stress.

Subgenome Dominance in an Interspecific Hybrid, Synthetic Allopolyploid, and a 140-Year-Old Naturally Established Neo-Allopolyploid Monkeyflower

Subgenome differences in levels of transposon methylation mirror increases in subgenome-specific expression bias observed over generations following hybridization and allopolyploidization in Mimulus. Recent studies have shown that one of the parental subgenomes in ancient polyploids is generally more dominant, having retained more genes and being more highly expressed, a phenomenon termed subgenome dominance. The genomic features that determine how quickly and which subgenome dominates within a newly formed polyploid remain poorly understood. To investigate the rate of emergence of subgenome dominance, we examined gene expression, gene methylation, and transposable element (TE) methylation in a natural, <140-year-old allopolyploid (Mimulus peregrinus), a resynthesized interspecies triploid hybrid (M. robertsii), a resynthesized allopolyploid (M. peregrinus), and progenitor species (M. guttatus and M. luteus). We show that subgenome expression dominance occurs instantly following the hybridization of divergent genomes and significantly increases over generations. Additionally, CHH methylation levels are reduced in regions near genes and within TEs in the first-generation hybrid, intermediate in the resynthesized allopolyploid, and are repatterned differently between the dominant and recessive subgenomes in the natural allopolyploid. Subgenome differences in levels of TE methylation mirror the increase in expression bias observed over the generations following hybridization. These findings provide important insights into genomic and epigenomic shock that occurs following hybridization and polyploid events and may also contribute to uncovering the mechanistic basis of heterosis and subgenome dominance.