Genome-wide DNA methylation dynamics following recent polyploidy in the allotetraploid Tragopogon miscellus (Asteraceae).

Polyploidy is an important evolutionary force, yet epigenetic mechanisms, such as DNA methylation, that regulate genome-wide expression of duplicated genes remain largely unknown. Here, we use Tragopogon (Asteraceae) as a model system to discover patterns and temporal dynamics of DNA methylation in recently formed polyploids. The naturally occurring allotetraploid Tragopogon miscellus formed in the last 95-100 yr from parental diploids Tragopogon dubius and T. pratensis. We profiled the DNA methylomes of these three species using whole-genome bisulfite sequencing. Genome-wide methylation levels in T. miscellus were intermediate between its diploid parents. However, nonadditive CG and CHG methylation occurred in transposable elements (TEs), with variation among TE types. Most differentially methylated regions (DMRs) showed parental legacy, but some novel DMRs were detected in the polyploid. Differentially methylated genes (DMGs) were also identified and characterized. This study provides the first assessment of both overall and locus-specific patterns of DNA methylation in a recent natural allopolyploid and shows that novel methylation variants can be generated rapidly after polyploid formation. Together, these results demonstrate that mechanisms to regulate duplicate gene expression may arise soon after allopolyploid formation and that these mechanisms vary among genes.

Nociceptor neurons affect cancer immunosurveillance.

Solid tumours are innervated by nerve fibres that arise from the autonomic and sensory peripheral nervous systems1-5. Whether the neo-innervation of tumours by pain-initiating sensory neurons affects cancer immunosurveillance remains unclear. Here we show that melanoma cells interact with nociceptor neurons, leading to increases in their neurite outgrowth, responsiveness to noxious ligands and neuropeptide release. Calcitonin gene-related peptide (CGRP)-one such nociceptor-produced neuropeptide-directly increases the exhaustion of cytotoxic CD8+ T cells, which limits their capacity to eliminate melanoma. Genetic ablation of the TRPV1 lineage, local pharmacological silencing of nociceptors and antagonism of the CGRP receptor RAMP1 all reduced the exhaustion of tumour-infiltrating leukocytes and decreased the growth of tumours, nearly tripling the survival rate of mice that were inoculated with B16F10 melanoma cells. Conversely, CD8+ T cell exhaustion was rescued in sensory-neuron-depleted mice that were treated with local recombinant CGRP. As compared with wild-type CD8+ T cells, Ramp1-/- CD8+ T cells were protected against exhaustion when co-transplanted into tumour-bearing Rag1-deficient mice. Single-cell RNA sequencing of biopsies from patients with melanoma revealed that intratumoral RAMP1-expressing CD8+ T cells were more exhausted than their RAMP1-negative counterparts, whereas overexpression of RAMP1 correlated with a poorer clinical prognosis. Overall, our results suggest that reducing the release of CGRP from tumour-innervating nociceptors could be a strategy to improve anti-tumour immunity by eliminating the immunomodulatory effects of CGRP on cytotoxic CD8+ T cells.

Transcriptional neural-like signaling contributes to an immune-suppressive tumor microenvironment.

Tumor innervation has recently been documented and characterized in various settings and tumor types. However, the role that nerves innervating tumors play in the pathogenesis of cancer has not been clarified. In this study, we searched for neural signaling from bulk RNA sequencing from The Cancer Genome Atlas (TCGA) dataset and looked for patterns of interactions between different cell types within the tumor environment. Using a presynapse signature (PSS) as a probe, we showed that multiple stromal cell types crosstalk and/or contribute to neural signals. Based on the correlation and linear regression, we hypothesized that neural signals contribute to an immune-suppressive tumor microenvironment (TME). To test this hypothesis, we performed in vitro dorsal root ganglion (DRG)/macrophage coculture experiments. Compared to the M2 macrophage monoculture, the DRG/M2 macrophage coculture prevented anti-inflammatory M2 to pro-inflammatory M1 polarization by LPS stimulation. Finally, a survey of different TCGA tumor types indicated that higher RNA neural signature is predictive of poor patient outcomes in multiple tumor types.

High quality haplotype-resolved genome assemblies of Populus tomentosa Carr., a stabilized interspecific hybrid species widespread in Asia.

Populus has a wide ecogeographical range spanning the Northern Hemisphere, and interspecific hybrids are common. Populus tomentosa Carr. is widely distributed and cultivated in the eastern region of Asia, where it plays multiple important roles in forestry, agriculture, conservation, and urban horticulture. Reference genomes are available for several Populus species, however, our goals were to produce a very high quality de novo chromosome-level genome assembly in P. tomentosa genome that could serve as a reference for evolutionary and ecological studies of hybrid speciation throughout the genus. Here, combining long-read sequencing and Hi-C scaffolding, we present a high-quality, haplotype-resolved genome assembly. The genome size was 740.2 Mb, with a contig N50 size of 5.47 Mb and a scaffold N50 size of 46.68 Mb, consisting of 38 chromosomes, as expected with the known diploid chromosome number (2n = 2x = 38). A total of 59,124 protein-coding genes were identified. Phylogenomic analyses revealed that P. tomentosa is comprised of two distinct subgenomes, which we deomonstrate is likely to have resulted from hybridization between Populus adenopoda as the female parent and Populus alba var. pyramidalis as the male parent, with an origin of approximately 3.93 Ma. Although highly colinear, significant structural variation was found between the two subgenomes. Our study provides a valuable resource for ecological genetics and forest biotechnology.

Natural variation in DNA methylation homeostasis and the emergence of epialleles.

In plants and mammals, DNA methylation plays a critical role in transcriptional silencing by delineating heterochromatin from transcriptionally active euchromatin. A homeostatic balance between heterochromatin and euchromatin is essential to genomic stability. This is evident in many diseases and mutants for heterochromatin maintenance, which are characterized by global losses of DNA methylation coupled with localized ectopic gains of DNA methylation that alter transcription. Furthermore, we have shown that genome-wide methylation patterns in Arabidopsis thaliana are highly stable over generations, with the exception of rare epialleles. However, the extent to which natural variation in the robustness of targeting DNA methylation to heterochromatin exists, and the phenotypic consequences of such variation, remain to be fully explored. Here we describe the finding that heterochromatin and genic DNA methylation are highly variable among 725 A. thaliana accessions. We found that genic DNA methylation is inversely correlated with that in heterochromatin, suggesting that certain methylation pathway(s) may be redirected to genes upon the loss of heterochromatin. This redistribution likely involves a feedback loop involving the DNA methyltransferase, CHROMOMETHYLASE 3 (CMT3), H3K9me2, and histone turnover, as highly expressed, long genes with a high density of CMT3-preferred CWG sites are more likely to be methylated. Importantly, although the presence of CG methylation in genes alone may not affect transcription, genes containing CG methylation are more likely to become methylated at non-CG sites and silenced. These findings are consistent with the hypothesis that natural variation in DNA methylation homeostasis may underlie the evolution of epialleles that alter phenotypes.

Author Correction: Widespread long-range cis-regulatory elements in the maize genome.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Widespread long-range cis-regulatory elements in the maize genome.

Genetic mapping studies on crops suggest that agronomic traits can be controlled by gene-distal intergenic loci. Despite the biological importance and the potential agronomic utility of these loci, they remain virtually uncharacterized in all crop species to date. Here, we provide genetic, epigenomic and functional molecular evidence to support the widespread existence of gene-distal (hereafter, distal) loci that act as long-range transcriptional cis-regulatory elements (CREs) in the maize genome. Such loci are enriched for euchromatic features that suggest their regulatory functions. Chromatin loops link together putative CREs with genes and recapitulate genetic interactions. Putative CREs also display elevated transcriptional enhancer activities, as measured by self-transcribing active regulatory region sequencing. These results provide functional support for the widespread existence of CREs that act over large genomic distances to control gene expression.

Multiple genes recruited from hormone pathways partition maize diterpenoid defences.

Duplication and divergence of primary pathway genes underlie the evolution of plant specialized metabolism; however, mechanisms partitioning parallel hormone and defence pathways are often speculative. For example, the primary pathway intermediate ent-kaurene is essential for gibberellin biosynthesis and is also a proposed precursor for maize antibiotics. By integrating transcriptional coregulation patterns, genome-wide association studies, combinatorial enzyme assays, proteomics and targeted mutant analyses, we show that maize kauralexin biosynthesis proceeds via the positional isomer ent-isokaurene formed by a diterpene synthase pair recruited from gibberellin metabolism. The oxygenation and subsequent desaturation of ent-isokaurene by three promiscuous cytochrome P450s and a new steroid 5α reductase indirectly yields predominant ent-kaurene-associated antibiotics required for Fusarium stalk rot resistance. The divergence and differential expression of pathway branches derived from multiple duplicated hormone-metabolic genes minimizes dysregulation of primary metabolism via the circuitous biosynthesis of ent-kaurene-related antibiotics without the production of growth hormone precursors during defence.

Monitoring the interplay between transposable element families and DNA methylation in maize.

DNA methylation and epigenetic silencing play important roles in the regulation of transposable elements (TEs) in many eukaryotic genomes. A majority of the maize genome is derived from TEs that can be classified into different orders and families based on their mechanism of transposition and sequence similarity, respectively. TEs themselves are highly methylated and it can be tempting to view them as a single uniform group. However, the analysis of DNA methylation profiles in flanking regions provides evidence for distinct groups of chromatin properties at different TE families. These differences among TE families are reproducible in different tissues and different inbred lines. TE families with varying levels of DNA methylation in flanking regions also show distinct patterns of chromatin accessibility and modifications within the TEs. The differences in the patterns of DNA methylation flanking TE families arise from a combination of non-random insertion preferences of TE families, changes in DNA methylation triggered by the insertion of the TE and subsequent selection pressure. A set of nearly 70,000 TE polymorphisms among four assembled maize genomes were used to monitor the level of DNA methylation at haplotypes with and without the TE insertions. In many cases, TE families with high levels of DNA methylation in flanking sequence are enriched for insertions into highly methylated regions. The majority of the >2,500 TE insertions into unmethylated regions result in changes in DNA methylation in haplotypes with the TE, suggesting the widespread potential for TE insertions to condition altered methylation in conserved regions of the genome. This study highlights the interplay between TEs and the methylome of a major crop species.

Genome-Wide Reinforcement of DNA Methylation Occurs during Somatic Embryogenesis in Soybean.

Somatic embryogenesis is an important tissue culture technique that sometimes leads to phenotypic variation via genetic and/or epigenetic changes. To understand the genomic and epigenomic impacts of somatic embryogenesis, we characterized soybean (Glycine max) epigenomes sampled from embryos at 10 different stages ranging from 6 weeks to 13 years of continuous culture. We identified genome-wide increases in DNA methylation from cultured samples, especially at CHH sites. The hypermethylation almost exclusively occurred in regions previously possessing non-CG methylation and was accompanied by increases in the expression of genes encoding the RNA-directed DNA methylation (RdDM) machinery. The epigenomic changes were similar between somatic and zygotic embryogenesis. Following the initial global wave of hypermethylation, rare decay events of maintenance methylation were observed, and the extent of the decay increased with time in culture. These losses in DNA methylation were accompanied by downregulation of genes encoding the RdDM machinery and transcriptome reprogramming reminiscent of transcriptomes during late-stage seed development. These results reveal a process for reinforcing already silenced regions to maintain genome integrity during somatic embryogenesis over the short term, which eventually decays at certain loci over longer time scales.

Rational development of transformation in Clostridium thermocellum ATCC 27405 via complete methylome analysis and evasion of native restriction-modification systems.

A major barrier to both metabolic engineering and fundamental biological studies is the lack of genetic tools in most microorganisms. One example is Clostridium thermocellum ATCC 27405T, where genetic tools are not available to help validate decades of hypotheses. A significant barrier to DNA transformation is restriction-modification systems, which defend against foreign DNA methylated differently than the host. To determine the active restriction-modification systems in this strain, we performed complete methylome analysis via single-molecule, real-time sequencing to detect 6-methyladenine and 4-methylcytosine and the rarely used whole-genome bisulfite sequencing to detect 5-methylcytosine. Multiple active systems were identified, and corresponding DNA methyltransferases were expressed from the Escherichia coli chromosome to mimic the C. thermocellum methylome. Plasmid methylation was experimentally validated and successfully electroporated into C. thermocellum ATCC 27405. This combined approach enabled genetic modification of the C. thermocellum-type strain and acts as a blueprint for transformation of other non-model microorganisms.

Epimutations are associated with CHROMOMETHYLASE 3-induced de novo DNA methylation.

In many plant species, a subset of transcribed genes are characterized by strictly CG-context DNA methylation, referred to as gene body methylation (gbM). The mechanisms that establish gbM are unclear, yet flowering plant species naturally without gbM lack the DNA methyltransferase, CMT3, which maintains CHG (H = A, C, or T) and not CG methylation at constitutive heterochromatin. Here, we identify the mechanistic basis for gbM establishment by expressing CMT3 in a species naturally lacking CMT3. CMT3 expression reconstituted gbM through a progression of de novo CHG methylation on expressed genes, followed by the accumulation of CG methylation that could be inherited even following loss of the CMT3 transgene. Thus, gbM likely originates from the simultaneous targeting of loci by pathways that promote euchromatin and heterochromatin, which primes genes for the formation of stably inherited epimutations in the form of CG DNA methylation.

Changes of gene expression but not cytosine methylation are associated with male parental care reflecting behavioural state, social context and individual flexibility.

Behaviour is often a front line response to changing environments. Recent studies show behavioural changes are associated with changes of gene expression; however, these studies have primarily focused on discrete behavioural states. We build on these studies by addressing additional contexts that produce qualitatively similar behavioural changes. We measured levels of gene expression and cytosine methylation, which is hypothesized to regulate the transcriptional architecture of behavioural transitions, within the brain during male parental care of the burying beetle Nicrophorus vespilloides in a factorial design. Male parenting is a suitably plastic behaviour because although male N. vespilloides typically do not provide direct care (i.e. feed offspring) when females are present, levels of feeding by a male equivalent to the female can be induced by removing the female. We examined three different factors: behavioural state (caring versus non-caring), social context (with or without a female mate) and individual flexibility (if a male switched to direct care after his mate was removed). The greatest number of differentially expressed genes were associated with behavioural state, followed by social context and individual flexibility. Cytosine methylation was not associated with changes of gene expression in any of the factors. Our results suggest a hierarchical association between gene expression and the different factors, but that this process is not controlled by cytosine methylation. Our results further suggest that the extent a behaviour is transient plays an underappreciated role in determining its underpinning molecular mechanisms.

TET-mediated epimutagenesis of the Arabidopsis thaliana methylome.

DNA methylation in the promoters of plant genes sometimes leads to transcriptional repression, and the loss of DNA methylation in methyltransferase mutants results in altered gene expression and severe developmental defects. However, many cases of naturally occurring DNA methylation variations have been reported, whereby altered expression of differentially methylated genes is responsible for agronomically important traits. The ability to manipulate plant methylomes to generate epigenetically distinct individuals could be invaluable for breeding and research purposes. Here, we describe "epimutagenesis," a method to rapidly generate DNA methylation variation through random demethylation of the Arabidopsis thaliana genome. This method involves the expression of a human ten-eleven translocation (TET) enzyme, and results in widespread hypomethylation that can be inherited to subsequent generations, mimicking mutants in the maintenance of DNA methyltransferase met1. Application of epimutagenesis to agriculturally significant plants may result in differential expression of alleles normally silenced by DNA methylation, uncovering previously hidden phenotypic variations.

A novel requirement for DROSHA in maintenance of mammalian CG methylation.

In mammals, faithful inheritance of genomic methylation patterns ensures proper gene regulation and cell behaviour, impacting normal development and fertility. Following establishment, genomic methylation patterns are transmitted through S-phase by the maintenance methyltransferase Dnmt1. Using a protein interaction screen, we identify Microprocessor component DROSHA as a novel DNMT1-interactor. Drosha-deficient embryonic stem (ES) cells display genomic hypomethylation that is not accounted for by changes in the levels of DNMT proteins. DNMT1-mediated methyltransferase activity is also reduced in these cells. We identify two transcripts that are specifically upregulated in Drosha- but not Dicer-deficient ES cells. Regions within these transcripts predicted to form stem-loop structures are processed by Microprocessor and can inhibit DNMT1-mediated methylation in vitro. Our results highlight DROSHA as a novel regulator of mammalian DNA methylation and we propose that DROSHA-mediated processing of RNA is necessary to ensure full DNMT1 activity. This adds to the DROSHA repertoire of non-miRNA dependent functions as well as implicating RNA in regulating DNMT1 activity and correct levels of genomic methylation.

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

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.

The evolution of CHROMOMETHYLASES and gene body DNA methylation in plants.

BACKGROUND: The 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. RESULTS: CMTs 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. CONCLUSIONS: The 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.

Dynamic changes in the transcriptome of Populus hopeiensis in response to abscisic acid.

Abscisic acid (ABA) plays a fundamental role in plant response and adaptation to abiotic stresses, such as drought, high salinity and low temperature. Populus hopeiensis exhibits exceptional tolerance to water-deficit environments and is therefore an excellent choice for studying drought tolerance in trees. This study provides a global view of transcriptome dynamics in P. hopeiensis in response to exogenous ABA using Illumina RNA-sequencing. Endogenous ABA content increased and reached a peak at 8 h after ABA treatment and then significantly decreased at latter time points. Differential expression analysis and Gene ontology enrichment revealed that the number of transcripts exhibited significant increase during the first 8 hours after ABA treatment, which then significantly decreased at 12 and 24 h. Transcription factors (TFs) analysis showed that six different patterns were observed based on the expression of the six TFs families (AP2/ERF, NAC, MYB, MYB-related, bZIP and WRKY) and the majority of differentially expressed TFs increased rapidly after ABA treatment. This study provides a robust resource for investigating the functions of genes induced by ABA and will help to develop a better understanding of the molecular regulatory mechanism in response to drought in poplar.

Widespread natural variation of DNA methylation within angiosperms.

BACKGROUND: DNA 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. RESULTS: By 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. CONCLUSIONS: These 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.