Optimized reduced representation bisulfite sequencing reveals tissue-specific mCHH islands in maize

Establishing an optimized ATAC-seq protocol for the maize

The advent of next-generation sequencing in crop improvement offers unprecedented insights into the chromatin landscape closely linked to gene activity governing key traits in plant development and adaptation. Particularly in maize, its dynamic chromatin structure is found to collaborate with massive transcriptional variations across tissues and developmental stages, implying intricate regulatory mechanisms, which highlights the importance of integrating chromatin information into breeding strategies for precise gene controls. The depiction of maize chromatin architecture using Assay for Transposase Accessible Chromatin with high-throughput sequencing (ATAC-seq) provides great opportunities to investigate cis-regulatory elements, which is crucial for crop improvement. In this context, we developed an easy-to-implement ATAC-seq protocol for maize with fewer nuclei and simple equipment. We demonstrate a streamlined ATAC-seq protocol with four key steps for maize in which nuclei purification can be achieved without cell sorting and using only a standard bench-top centrifuge. Our protocol, coupled with the bioinformatic analysis, including validation by read length periodicity, key metrics, and correlation with transcript abundance, provides a precise and efficient assessment of the maize chromatin landscape. Beyond its application to maize, our testing design holds the potential to be applied to other crops or other tissues, especially for those with limited size and amount, establishing a robust foundation for chromatin structure studies in diverse crop species.

Reduced non-CpG methylation is a potential epigenetic target after spinal cord injury

DNA methylation is a critical epigenetic regulator in the occurrence and development of diseases and is closely related to various functional responses in relation to spinal cord injury. To investigate the role of DNA methylation in spinal cord injury, we constructed a library with reduced-representation bisulfite sequencing data obtained at various time points (day 0-42) after spinal cord injury in mice. Global DNA methylation levels, specifically non-CpG (CHG and CHH) methylation levels, decreased modestly following spinal cord injury. Stages post-spinal cord injury were classified as early (day 0-3), intermediate (day 7-14), and late (day 28-42) based on similarity and hierarchical clustering of global DNA methylation patterns. The non-CpG methylation level, which included CHG and CHH methylation levels, was markedly reduced despite accounting for a minor proportion of total methylation abundance. At multiple genomic sites, including the 5' untranslated regions, promoter, exon, intron, and 3' untranslated regions, the non-CpG methylation level was markedly decreased following spinal cord injury, whereas the CpG methylation level remained unchanged at these locations. Approximately one-half of the differentially methylated regions were located in intergenic areas; the other differentially methylated regions in both CpG and non-CpG regions were clustered in intron regions, where the DNA methylation level was highest. The function of genes associated with differentially methylated regions in promoter regions was also investigated. From Gene Ontology analysis results, DNA methylation was implicated in a number of essential functional responses to spinal cord injury, including neuronal synaptic connection creation and axon regeneration. Notably, neither CpG methylation nor non-CpG methylation was implicated in the functional response of glial or inflammatory cells. In summary, our work elucidated the dynamic pattern of DNA methylation in the spinal cord following injury and identified reduced non-CpG methylation as an epigenetic target after spinal cord injury in mice.

Exploring the crop epigenome: a comparison of DNA methylation profiling techniques

Epigenetic modifications play a vital role in the preservation of genome integrity and in the regulation of gene expression. DNA methylation, one of the key mechanisms of epigenetic control, impacts growth, development, stress response and adaptability of all organisms, including plants. The detection of DNA methylation marks is crucial for understanding the mechanisms underlying these processes and for developing strategies to improve productivity and stress resistance of crop plants. There are different methods for detecting plant DNA methylation, such as bisulfite sequencing, methylation-sensitive amplified polymorphism, genome-wide DNA methylation analysis, methylated DNA immunoprecipitation sequencing, reduced representation bisulfite sequencing, MS and immuno-based techniques. These profiling approaches vary in many aspects, including DNA input, resolution, genomic region coverage, and bioinformatics analysis. Selecting an appropriate methylation screening approach requires an understanding of all these techniques. This review provides an overview of DNA methylation profiling methods in crop plants, along with comparisons of the efficacy of these techniques between model and crop plants. The strengths and limitations of each methodological approach are outlined, and the importance of considering both technical and biological factors are highlighted. Additionally, methods for modulating DNA methylation in model and crop species are presented. Overall, this review will assist scientists in making informed decisions when selecting an appropriate DNA methylation profiling method.

Antioxidant Regulation and DNA Methylation Dynamics During Mikania micrantha Seed Germination Under Cold Stress

As a primary goal, adaptation to cold climate could expand an invasion range of exotic plants. Here, we aimed to explore the regulation strategy of M. micrantha seed development under cold stress through molecular physiology and multi-omics analysis. Significant increase of hydrogen peroxide, malondialdehyde, and electrolyte leakage observed under cold stress revealed that oxidative damage within M. micrantha seed cells was induced in the initial germination phase. Proteomic data underscored an activation of antioxidant activity to maintain redox homeostasis, with a cluster of antioxidant proteins identified. Genomic-wide transcriptome, in combination with time-series whole-genome bisulfite sequencing mining, elucidated that seven candidate genes, which were the target of DNA demethylation-dependent ROS scavenging, were possibly associated with an M. micrantha germ break. Progressive gain of CHH context DNA methylation identified in an early germination phrase suggested a role of a DNA methylation pathway, while an active DNA demethylation pathway was also initiated during late seed development, which was in line with the expression trend of methylation and demethylation-related genes verified through qRT-PCR. These data pointed out that cold-dependent DNA demethylation and an antioxidant regulatory were involved together in restoring seed germination. The expression level of total 441 genes presented an opposite trend to the methylation divergence, while the expression of total 395 genes was proved to be negatively associated with their methylation levels. These data provided new insights into molecular reprograming events during M. micrantha seed development.

Divergent DNA methylation contributes to duplicated gene evolution and chilling response in tea plants

The tea plant (Camellia sinensis) is a thermophilic cash crop and contains a highly duplicated and repeat-rich genome. It is still unclear how DNA methylation regulates the evolution of duplicated genes and chilling stress in tea plants. We therefore generated a single-base-resolution DNA methylation map of tea plants under chilling stress. We found that, compared with other plants, the tea plant genome is highly methylated in all three sequence contexts, including CG, CHG and CHH (where H = A, T, or C), which is further proven to be correlated with its repeat content and genome size. We show that DNA methylation in the gene body negatively regulates the gene expression of tea plants, whereas non-CG methylation in the flanking region enables a positive regulation of gene expression. We demonstrate that transposable element-mediated methylation dynamics significantly drives the expression divergence of duplicated genes in tea plants. The DNA methylation and expression divergence of duplicated genes in the tea plant increases with evolutionary age and selective pressure. Moreover, we detect thousands of differentially methylated genes, some of which are functionally associated with chilling stress. We also experimentally reveal that DNA methyltransferase genes of tea plants are significantly downregulated, whereas demethylase genes are upregulated at the initial stage of chilling stress, which is in line with the significant loss of DNA methylation of three well-known cold-responsive genes at their promoter and gene body regions. Overall, our findings underscore the importance of DNA methylation regulation and offer new insights into duplicated gene evolution and chilling tolerance in tea plants.

Whole-genome characterization of chronological age-associated changes in methylome and circular RNAs in moso bamboo (Phyllostachys edulis) from vegetative to floral growth

In mammals, DNA methylation is associated with aging. However, age-related DNA methylation changes during phase transitions largely remain unstudied in plants. Moso bamboo (Phyllostachys edulis) requires a very long time to transition from the vegetative to the floral phase. To comprehensively investigate the association of DNA methylation with aging, we present here single-base-resolution DNA methylation profiles using both high-throughput bisulfite sequencing and single-molecule nanopore-based DNA sequencing, covering the long period of vegetative growth and transition to flowering in moso bamboo. We discovered that CHH methylation gradually accumulates from vegetative to reproductive growth in a time-dependent fashion. Differentially methylated regions, correlating with chronological aging, occurred preferentially at both transcription start sites and transcription termination sites. Genes with CG methylation changes showed an enrichment of Gene Ontology (GO) categories in 'vegetative to reproductive phase transition of meristem'. Combining methylation data with mRNA sequencing revealed that DNA methylation in promoters, introns and exons may have different roles in regulating gene expression. Finally, circular RNA (circRNA) sequencing revealed that the flanking introns of circRNAs are hypermethylated and enriched in long terminal repeat (LTR) retrotransposons. Together, the observations in this study provide insights into the dynamic DNA methylation and circRNA landscapes, correlating with chronological age, which paves the way to study further the impact of epigenetic factors on flowering in moso bamboo.

The regulatory landscape of early maize inflorescence development

BACKGROUND

The functional genome of agronomically important plant species remains largely unexplored, yet presents a virtually untapped resource for targeted crop improvement. Functional elements of regulatory DNA revealed through profiles of chromatin accessibility can be harnessed for fine-tuning gene expression to optimal phenotypes in specific environments.

RESULT

Here, we investigate the non-coding regulatory space in the maize (Zea mays) genome during early reproductive development of pollen- and grain-bearing inflorescences. Using an assay for differential sensitivity of chromatin to micrococcal nuclease (MNase) digestion, we profile accessible chromatin and nucleosome occupancy in these largely undifferentiated tissues and classify at least 1.6% of the genome as accessible, with the majority of MNase hypersensitive sites marking proximal promoters, but also 3' ends of maize genes. This approach maps regulatory elements to footprint-level resolution. Integration of complementary transcriptome profiles and transcription factor occupancy data are used to annotate regulatory factors, such as combinatorial transcription factor binding motifs and long non-coding RNAs, that potentially contribute to organogenesis, including tissue-specific regulation between male and female inflorescence structures. Finally, genome-wide association studies for inflorescence architecture traits based solely on functional regions delineated by MNase hypersensitivity reveals new SNP-trait associations in known regulators of inflorescence development as well as new candidates.

CONCLUSIONS

These analyses provide a comprehensive look into the cis-regulatory landscape during inflorescence differentiation in a major cereal crop, which ultimately shapes architecture and influences yield potential.

Epigenomic landscape and epigenetic regulation in maize

Epigenetic regulation has been implicated in the control of multiple agronomic traits in maize. Here, we review current advances in our understanding of epigenetic regulation, which has great potential for improving agronomic traits and the environmental adaptability of crops. Epigenetic regulation plays vital role in the control of complex agronomic traits. Epigenetic variation could contribute to phenotypic diversity and can be used to improve the quality and productivity of crops. Maize (Zea mays L.), one of the most widely cultivated crops for human food, animal feed, and ethanol biofuel, is a model plant for genetic studies. Recent advances in high-throughput sequencing technology have made possible the study of epigenetic regulation in maize on a genome-wide scale. In this review, we discuss recent epigenetic studies in maize many achieved by Chinese research groups. These studies have explored the roles of DNA methylation, posttranslational modifications of histones, chromatin remodeling, and noncoding RNAs in the regulation of gene expression in plant development and environment response. We also provide our future prospects for manipulating epigenetic regulation to improve crops.

Meiocyte-Specific and AtSPO11-1-Dependent Small RNAs and Their Association with Meiotic Gene Expression and Recombination

Meiotic recombination ensures accurate chromosome segregation and results in genetic diversity in sexually reproducing eukaryotes. Over the last few decades, the genetic regulation of meiotic recombination has been extensively studied in many organisms. However, the role of endogenous meiocyte-specific small RNAs (ms-sRNAs; 21-24 nucleotide [nt]) and their involvement in meiotic recombination are unclear. Here, we sequenced the total small RNA (sRNA) and messenger RNA populations from meiocytes and leaves of wild type Arabidopsis (Arabidopsis thaliana) and meiocytes of spo11-1, a mutant defective in double-strand break formation, and we discovered 2,409 ms-sRNA clusters, 1,660 of which areSPORULATION 11-1 (AtSPO11-1)-dependent. Unlike mitotic small interfering RNAs that are enriched in intergenic regions and associated with gene silencing, ms-sRNAs are significantly enriched in genic regions and exhibit a positive correlation with genes that are preferentially expressed in meiocytes (i.e. Arabidopsis SKP1-LIKE1 and RAD51), in a fashion unrelated to DNA methylation. We also found that AtSPO11-1-dependent sRNAs have distinct characteristics compared with ms-sRNAs and tend to be associated with two known types of meiotic recombination hotspot motifs (i.e. CTT-repeat and A-rich motifs). These results reveal different meiotic and mitotic sRNA landscapes and provide new insights into how sRNAs relate to gene expression in meiocytes and meiotic recombination.

Opportunities and limitations of reduced representation bisulfite sequencing in plant ecological epigenomics

Contents Summary 738 I. Introduction 738 II. RRBS loci as genome-wide epigenetic markers 739 III. Exploiting functional annotation of RRBS loci 739 IV. Limitations of RRBS methods for nonmodel species 740 V. Maximising the impact of RRBS in plants 741 VI. Conclusions 741 Acknowledgements 741 SUMMARY: Investigating the features and implications of epigenetic mechanisms across the breadth of organisms and ecosystems is important for understanding the ecological relevance of epigenetics. Several cost-effective reduced representation bisulfite sequencing approaches (RRBS) have been recently developed and applied to different organisms that lack a well annotated reference genome. These new approaches improve the assessment of epigenetic diversity in ecological settings and may provide functional insights. We assess here the opportunities and limitations of RRBS in nonmodel plant species. Well thought out experimental designs that include complementary gene expression studies, and the improvement of genomics resources for the target group, promise to maximize the effect of future RRBS studies.

DNA methylation and genetic degeneration of the Y chromosome in the dioecious plant Silene latifolia

Background

S. latifolia is a model organism for the study of sex chromosome evolution in plants. Its sex chromosomes include large regions in which recombination became gradually suppressed. The regions tend to expand over time resulting in the formation of evolutionary strata. Non-recombination and later accumulation of repetitive sequences is a putative cause of the size increase in the Y chromosome. Gene decay and accumulation of repetitive DNA are identified as key evolutionary events. Transposons in the X and Y chromosomes are distributed differently and there is a regulation of transposon insertion by DNA methylation of the target sequences, this points to an important role of DNA methylation during sex chromosome evolution in Silene latifolia. The aim of this study was to elucidate whether the reduced expression of the Y allele in S. latifolia is caused by genetic degeneration or if the cause is methylation triggered by transposons and repetitive sequences.

Results

Gene expression analysis in S. latifolia males has shown expression bias in both X and Y alleles. To determine whether these differences are caused by genetic degeneration or methylation spread by transposons and repetitive sequences, we selected several sex-linked genes with varying degrees of degeneration and from different evolutionary strata. Immunoprecipitation of methylated DNA (MeDIP) from promoter, exon and intron regions was used and validated through bisulfite sequencing. We found DNA methylation in males, and only in the promoter of genes of stratum I (older). The Y alleles in genes of stratum I were methylation enriched compared to X alleles. There was also abundant and high percentage methylation in the CHH context in most sequences, indicating de novo methylation through the RdDM pathway.

Conclusions

We speculate that TE accumulation and not gene decay is the cause of DNA methylation in the S. latifolia Y sex chromosome with influence on the process of heterochromatinization.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4936-y) contains supplementary material, which is available to authorized users.

Reduced Representation Bisulfite Sequencing in Maize

DNA methylation is an epigenetic modification that regulates plant development (Law and Jacobsen, 2010). Whole genome bisulfite sequencing (WGBS) is a state-of-the-art method for profiling genome-wide methylation patterns with single-base resolution ( Cokus et al., 2008 ). However, for an organism with a large genome, e.g., the 2.1 Gb genome of maize, WGBS may be very expensive. Reduced representation bisulfite sequencing (RRBS) has been developed in mammalian studies ( Smith et al., 2009 ). By digesting the genome with MspI with a size selection range of approximately 40-220 bp, CG-rich regions covering only ~1% of the human genome can be specifically sequenced. However, unlike mammalian genomes, plant genomes do not exhibit clear CpG islands. Therefore the original RRBS protocol is not suitable for plants. Accordingly, we developed an in silico pipeline to select specific enzymes to generate a region of interest (ROI)-enriched, e.g., promoter-enriched, reduced representation genome in plants ( Hsu et al., 2017 ). By digesting the maize genome with MseI and selecting 40-300 bp segments, we sequenced about one-fourth of the maize genome while preserving 84.3% of the promoter information. The protocol has been successfully established in maize and can be broadly used in any genome. Our in silico pipeline is combined with the RRBS library preparation protocol, allowing for the computational analysis and experimental validation.