Genome-wide DNA methylation profiling by modified reduced representation bisulfite sequencing in Brassica rapa suggests that epigenetic modifications play a key role in polyploid genome evolution

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.

Validation of Reference Genes in a Population of Blueberry (Vaccinium corymbosum) Plants Regenerated in Colchicine

For the first time we report the validation of reference genes in plants from a population of blueberry (Vaccinium corymbosum) clones cultured in vitro on a colchicine-supplemented medium. Nodal segment explants of the cultivar Duke were regenerated by organogenesis under different periods of colchicine 1 mg/L exposure (1, 2, 3, 5, 30 days). The clones selected for the study showed variability for phenotypic traits after 2 years of adaptation to field conditions, compared to plants of the donor genotype that were regenerated on a medium without colchicine. Vaccinium myrtillus (GAPDH) and Vaccinium macrocarpon (ATP1, NADH, RPOB and COX2) were used as reference genomes for primer design. The results show that colchicine treatments can cause genomic changes in blueberry plants. At the molecular level, exposure of plants to colchicine in early periods could promote an increase in gene expression of specific genes such as ATP1, COX2, GAPDH, MATK, NADH and RPOB. However, prolonged exposure (30 days) could decrease gene expression of the genes studied. For qPCR assays, the primers designed for ATP1, COX2, GAPDH and MATK genes showed high efficiency. In addition, the GAPDH, ATP1, NADH and COX2 genes showed high stability and could be recommended as potential reference genes for gene expression assays in Vaccinium.

Investigation of Brassica and its relative genomes in the post-genomics era

The Brassicaceae family includes many economically important crop species, as well as cosmopolitan agricultural weed species. In addition, Arabidopsis thaliana, a member of this family, is used as a molecular model plant species. The genus Brassica is mesopolyploid, and the genus comprises comparatively recently originated tetrapolyploid species. With these characteristics, Brassicas have achieved the commonly accepted status of model organisms for genomic studies. This paper reviews the rapid research progress in the Brassicaceae family from diverse omics studies, including genomics, transcriptomics, epigenomics, and three-dimensional (3D) genomics, with a focus on cultivated crops. The morphological plasticity of Brassicaceae crops is largely due to their highly variable genomes. The origin of several important Brassicaceae crops has been established. Genes or loci domesticated or contributing to important traits are summarized. Epigenetic alterations and 3D structures have been found to play roles in subgenome dominance, either in tetraploid Brassica species or their diploid ancestors. Based on this progress, we propose future directions and prospects for the genomic investigation of Brassicaceae crops.

Whole-Genome DNA Methylation Analysis in Brassica rapa subsp. perviridis in Response to Albugo candida Infection

DNA methylation is an epigenetic mark associated with several mechanisms in plants including immunity mechanisms. However, little is known about the regulatory role of DNA methylation in the resistance response of Brassica species against fungal diseases. White rust, caused by the fungus Albugo candida, is one of the most widespread and destructive diseases of all the cultivated Brassica species, particularly Brassica rapa L. and Brassica juncea (L.) Czern and Coss. Here, we investigate whole-genome DNA methylation modifications of B. rapa subsp. perviridis in response to white rust. As a result, 233 and 275 differentially methylated regions (DMRs) in the susceptible cultivar "Misugi" and the resistant cultivar "Nanane" were identified, respectively. In both cultivars, more than half of the DMRs were associated with genes (DMR-genes). Gene expression analysis showed that 13 of these genes were also differentially expressed between control and infected samples. Gene ontology enrichment analysis of DMR genes revealed their involvement in various biological processes including defense mechanisms. DMRs were unevenly distributed around genes in susceptible and resistant cultivars. In "Misugi," DMRs tended to be located within genes, while in "Nanane," DMRs tended to be located up and downstream of the genes. However, CG DMRs were predominantly located within genes in both cultivars. Transposable elements also showed association with all three sequence contexts of DMRs but predominantly with CHG and CHH DMRs in both cultivars. Our findings indicate the occurrence of DNA methylation modifications in B. rapa in response to white rust infection and suggest a potential regulatory role of DNA methylation modification in defense mechanisms which could be exploited to improve disease resistance.

An improved Raphanus sativus cv. WK10039 genome localizes centromeres, uncovers variation of DNA methylation and resolves arrangement of the ancestral Brassica genome blocks in radish chromosomes

This study presents an improved genome of Raphanus sativus cv. WK10039 uncovering centromeres and differentially methylated regions of radish chromosomes. Comprehensive genome comparison of radish and diploid Brassica species of U's triangle reveals that R. sativus arose from the Brassica B genome lineage and is a sibling species of B. nigra. Radish (Raphanus sativus L.) is a key root vegetable crop closely related to the Brassica crop species of the family Brassicaceae. We reported a draft genome of R. sativus cv. WK10039 (Rs1.0), which had 54.6 Mb gaps. To study the radish genome and explore previously unknown regions, we generated an improved genome assembly (Rs2.0) by long-read sequencing and high-resolution genome-wide mapping of chromatin interactions. Rs2.0 was 434.9 Mb in size with 0.27 Mb gaps, and the N50 scaffold length was 37.3 Mb (40-fold larger assembly compared to Rs1.0). Approximately 38% of Rs2.0 was comprised of repetitive sequences, and 52,768 protein-coding genes and 4845 non-protein-coding genes were predicted and annotated. The improved contiguity and coverage of Rs2.0, along with the detection of highly methylated regions, enabled localization of centromeres where R. sativus-specific centromere-associated repeats, full-length OTA and CRM LTR-Gypsy retrotransposons, hAT-Ac, CMC-EnSpm and Helitron DNA transposons, and sequences highly homologous to B. nigra centromere-specific CENH3-associated CL sequences were enriched. Whole-genome bisulfite sequencing combined with mRNA sequencing identified differential epigenetic marks in the radish genome related to tissue development. Synteny comparison and genomic distance analysis of radish and three diploid Brassica species of U's triangle suggested that the radish genome arose from the Brassica B genome lineage through unique rearrangement of the triplicated ancestral Brassica genome after splitting of the Brassica A/C and B genomes.

Reprogramming of the Hevea brasiliensis Epigenome and Transcriptome in Response to Cold Stress

Low temperature is a key factor limiting the rubber plantation extending to high latitude area. Previous work has shown that cold-induced DNA demethylation was coordinated with the expression of cold-responsive (COR) genes in Hevea brasiliensis. In this work, reduced representation bisulphite sequencing analysis of H. brasiliensis showed that cold treatment induced global genomic DNA demethylation and altered the sequence contexts of methylated cytosines, but the levels of mCG methylation in transposable elements were slightly enhanced by cold treatment. Integrated analysis of the DNA methylome and transcriptome revealed 400 genes whose expression correlated with altered DNA methylation. DNA demethylation in the upstream region of gene seems to correlate with higher gene expression, whereas demethylation in the gene body has less association. Our results suggest that cold treatment globally change the genomic DNA methylation status of the rubber tree, which might coordinate reprogramming of the transcriptome.

Distribution, expression and methylation analysis of positively selected genes provides insights into the evolution in Brassica rapa

It is believed that positive selection is one of the major evolutionary forces underlying organism phenotypic diversification. Nevertheless, the characteristics of positively selected genes (PSGs), have not been well investigated. In this study, we performed a genome-wide analysis of orthologous genes between Brassica rapa (B. rapa) and Brassica oleracea (B. oleracea), and identified 468 putative PSGs. Our data show that, (1) PSGs are enriched in plant hormone signal transduction pathway and the transcription factor family; (2) PSGs are significantly lower expressed than randomly selected non-PSGs; (3) PSGs with tissue specificity are significantly higher expressed in the callus and reproductive tissues (flower and silique) than in vegetable tissues (root, stem and leaf); (4) the proportion of PSGs is positively correlated with the number of retained triplication gene copies, but the expression level of PSGs decay with the increasing of triplication gene copies; (5) the CG and CHG methylation levels of PSGs are significantly higher in introns and UTRs than in the promoter and exon regions; (6) the percent of transposable element is in proportion to the methylation level, and DNA methylation (especially in the CG content) has the tendency to reduce the expression of PSGs. This study provides insights into the characteristics, evolution, function, expression and methylation of PSGs in B. rapa.

Epigenetic approaches to crop breeding: current status and perspectives

In order to tackle the cumulative adverse effects of global climate change, reduced farmland, and heightened needs of an ever-increasing world population, modern agriculture is in urgent search of solutions that can ensure world food security and sustainable development. Classical crop breeding is still a powerful method to obtain crops with valued agronomical traits, but its potential is gradually being compromised by the menacing decline of genetic variation. Resorting to the epigenome as a source of variation could serve as a promising alternative. Here, we discuss current status of epigenetics-mediated crop breeding (epibreeding), highlight its advances and limitations, outline currently available methodologies, and propose novel RNA-based strategies to modify the epigenome in a gene-specific and transgene-free manner.

Genome Wide MeDIP-Seq Profiling of Wild and Cultivated Olives Trees Suggests DNA Methylation Fingerprint on the Sensory Quality of Olive Oil

Secondary metabolites are particularly important to humans due to their pharmaceutical properties. Moreover, secondary metabolites are key compounds in climate change adaptation in long-living trees. Recently, it has been described that the domestication of Olea subspecies had no major selection signature on coding variants and was mainly related to changes in gene expression. In addition, the phenotypic plasticity in Olea subspecies was linked to the activation of transposable elements in the genes neighboring. Here, we investigated the imprint of DNA methylation in the unassigned fraction of the phenotypic plasticity of the Olea subspecies, using methylated DNA immuno-precipitation sequencing (MeDIP-seq) for a high-resolution genome-wide DNA methylation profiling of leaves and fruits during fruit development in wild and cultivated olives from Turkey. Notably, the methylation profiling showed a differential DNA methylation in secondary metabolism responsible for the sensory quality of olive oil. Here, we highlight for the first time the imprint of DNA methylation in modulating the activity of the Linoleate 9S lipoxygenase in the biosynthesis of volatile aromatic compounds. Unprecedently, the current study reveals the methylation status of the olive genome during fruit ripening.

In Response to Abiotic Stress, DNA Methylation Confers EpiGenetic Changes in Plants

Epigenetics involves the heritable changes in patterns of gene expression determined by developmental and abiotic stresses, i.e., drought, cold, salinity, trace metals, and heat. Gene expression is driven by changes in DNA bases, histone proteins, the biogenesis of ncRNA, and changes in the nucleotide sequence. To cope with abiotic stresses, plants adopt certain changes driven by a sophisticated biological system. DNA methylation is a primary mechanism for epigenetic variation, which can induce phenotypic alterations in plants under stress. Some of the stress-driven changes in plants are temporary, while some modifications may be stable and inheritable to the next generations to allow them to cope with such extreme stress challenges in the future. In this review, we discuss the pivotal role of epigenetically developed phenotypic characteristics in plants as an evolutionary process participating in adaptation and tolerance responses to abiotic and biotic stresses that alter their growth and development. We emphasize the molecular process underlying changes in DNA methylation, differential variation for different species, the roles of non-coding RNAs in epigenetic modification, techniques for studying DNA methylation, and its role in crop improvement in tolerance to abiotic stress (drought, salinity, and heat). We summarize DNA methylation as a significant future research priority for tailoring crops according to various challenging environmental issues.

Genomics Armed With Diversity Leads the Way in Brassica Improvement in a Changing Global Environment

Meeting the needs of a growing world population in the face of imminent climate change is a challenge; breeding of vegetable and oilseed Brassica crops is part of the race in meeting these demands. Available genetic diversity constituting the foundation of breeding is essential in plant improvement. Elite varieties, land races, and crop wild species are important resources of useful variation and are available from existing genepools or genebanks. Conservation of diversity in genepools, genebanks, and even the wild is crucial in preventing the loss of variation for future breeding efforts. In addition, the identification of suitable parental lines and alleles is critical in ensuring the development of resilient Brassica crops. During the past two decades, an increasing number of high-quality nuclear and organellar Brassica genomes have been assembled. Whole-genome re-sequencing and the development of pan-genomes are overcoming the limitations of the single reference genome and provide the basis for further exploration. Genomic and complementary omic tools such as microarrays, transcriptomics, epigenetics, and reverse genetics facilitate the study of crop evolution, breeding histories, and the discovery of loci associated with highly sought-after agronomic traits. Furthermore, in genomic selection, predicted breeding values based on phenotype and genome-wide marker scores allow the preselection of promising genotypes, enhancing genetic gains and substantially quickening the breeding cycle. It is clear that genomics, armed with diversity, is set to lead the way in Brassica improvement; however, a multidisciplinary plant breeding approach that includes phenotype = genotype × environment × management interaction will ultimately ensure the selection of resilient Brassica varieties ready for climate change.

Characterization of Histone H3 Lysine 4 and 36 Tri-methylation in Brassica rapa L

Covalent modifications of histone proteins act as epigenetic regulators of gene expression. We report the distribution of two active histone marks (H3K4me3 and H3K36me3) in 14-day leaves in two lines of Brassica rapa L. by chromatin immunoprecipitation sequencing. Both lines were enriched with H3K4me3 and H3K36me3 marks at the transcription start site, and the transcription level of a gene was associated with the level of H3K4me3 and H3K36me3. H3K4me3- and H3K36me3-marked genes showed low tissue-specific gene expression, and genes with both H3K4me3 and H3K36me3 had a high level of expression and were constitutively expressed. Bivalent active and repressive histone modifications such as H3K4me3 and H3K27me3 marks or antagonistic coexistence of H3K36me3 and H3K27me3 marks were observed in some genes. Expression may be susceptible to changes by abiotic and biotic stresses in genes having both H3K4me3 and H3K27me3 marks. We showed that the presence of H3K36me3 marks was associated with different gene expression levels or tissue specificity between paralogous paired genes, suggesting that H3K36me3 might be involved in subfunctionalization of the subgenomes.

The effects of hybridization and genome doubling in plant evolution via allopolyploidy

Polyploidy is a pervasive and recurring phenomenon across the tree of life, which occurred at variable time scales, ecological amplitudes and cell types, and is especially prominent in the evolutionary histories of plants. Importantly, many of the world's most important crops and noxious invasive weeds are recent polyploids. Polyploidy includes two major types, autopolyploidy, referring to doubling of a single species genome, and allopolyploidy referring to doubling of two or more merged genomes via biological hybridization of distinct but related species. The prevalence of both types of polyploidy implies that both genome doubling alone and doubling coupled with hybridization confer selective advantages over their diploid progenitors under specific circumstances. In cases of allopolyploidy, the two events, genome doubling and hybridization, have both advantages and disadvantages. Accumulated studies have established that, in allopolyploidy, some advantage(s) of doubling may compensate for the disadvantage(s) of hybridity and vice versa, although further study is required to validate generality of this trend. Some studies have also revealed a variety of non-Mendelian genetic and genomic consequences induced by doubling and hybridization separately or concertedly in nascent allopolyploidy; however, the significance of which to the immediate establishment and longer-term evolutionary success of allopolyploid species remain to be empirically demonstrated and ecologically investigated. This review aims to summarize recent advances in our understanding of the roles of hybridization and genome doubling, in separation and combination, in the evolution of allopolyploid genomes, as well as fruitful future research directions that are emerging from these studies.

Single-base cytosine methylation analysis in fruits of three Capsicum species

Single-base cytosine methylation analysis across fruits of Capsicum annuum, C. chinense and C. frutescens showed global average methylation ranging from 82.8-89.1%, 77.6-83.9%, and 22.4-25% at CG, CHG and CHH contexts, respectively. High gene-body methylation at CG and CHG was observed across Capsicum species. The C. annuum showed the highest proportion (>80%) of mCs at different genomic regions compared to C. chinense and C. frutescens. Cytosine methylation for transposable-elements were lower in C. frutescens compared to C. annuum and C. chinense. A total of 510,165 CG, 583112 CHG and 277,897 CHH DMRs were identified across three Capsicum species. The differentially methylated regions (DMRs) distribution analysis revealed C. frutescens as more hypo-methylated compared to C. annuum and C. chinense, and also the presence of more intergenic DMRs in Capsicum genome. At CG and CHG context, gene expression and promoter methylation showed inverse correlations. Furthermore, the observed correlation between methylation and expression of genes suggested the potential role of methylation in Capsicum fruit development/ripening.

The genome of the Cauliflower mosaic virus, a plant pararetrovirus, is highly methylated in the nucleus

Cytosine methylation is an important defense against invasive DNAs. Here, cytosine methylation profiles of a plant pararetrovirus, Cauliflower mosaic virus (CaMV), were investigated. Nuclear CaMV DNA is highly methylated throughout the genome including at transcription regulatory regions, but the virion DNA is unmethylated. In vitro CG methylation of the viral 35S promoter reduces transcription from the downstream gene. Although nuclear CaMV DNA is highly methylated, its transcripts are accumulated in the nucleus. The data suggest that a small population of unmethylated viral genomes produced through reverse transcription are constantly delivered back to the nucleus. Small RNA profiles suggest that methylation of the CaMV DNA may be due to de novo methylation through 21-, 22-, and 24-nt small RNAs with adenines at their 5' terminus.

Changes in DNA methylation pattern of apple long-term in vitro shoot culture and acclimatized plants

DNA methylation is a process of epigenetic modification that can alter the functionality of a genome. Using whole-genome bisulfite sequencing, this study quantify the level of DNA methylation in the epigenomes of two diploid apple (Malus x domestica) scion cultivars ('McIntosh' and 'Húsvéti rozmaring') derived from three environmental conditions: in vivo mother plants in an orchard, in vitro culture, and acclimatized in vitro plants. The global DNA methylation levels were not dependent on the source of plant material, and the average level of DNA methylation was 49.77%, 34.65% and 8.77% in CpG, CHG and CHH contexts, respectively. Significant differences in DNA methylation were identified in 586 (specifically 334, 201 and 131 in CpG, CHG and CHH contexts, respectively) out of 45,116 genes, including promoter and coding sequences. These were classified as differentially methylated genes (DMGs). This is a 1.3% difference in the level of DNA methylation of genes in response to a change in the environment. Differential methylation was visualised by MA plots and functional genomic maps were established for biological processes, molecular functions and cellular components. When the DMGs were considered, in vitro tissue culture resulted in the highest level of methylation, but it was lower in acclimatized in vitro plants which was similar to that in the mother tree. Methylation patterns of the two scions differed, indicating cultivar-specific epigenetic regulation of gene expression during adaptation to various environments. After selecting genes that displayed differences larger than ±10% in CpG and CHG contexts, or larger than ±1.35% in the CHH context from among the DMGs, they were annotated in Blast2 GO v5.1.12 for Gene Ontology. DMGs identified as MD07G1113000 (protein transport), MD08G1041600 (extracellular space), MD09G1054800 (phosphatidic acid binding), and MD10G1265800 (not annotated) were methylated in all three contexts in in vitro shoots. These DNA methylation results suggest that epigenetic changes may contribute to the adaptation of apple to environmental changes by modifying the epigenome and thereby gene expression.

Involvement of noncoding RNAs in epigenetic modifications of esophageal cancer

Esophageal cancer (EC) is a serious digestive malignancy and is a leading cause of cancer-related mortality. Apart from genetic mutations, many epigenetic alterations including DNA methylation and histone modifications associated with chromatin remodeling have been identified in the regulation of gene expression in EC. Recently, noncoding RNAs, and mainly lncRNAs and miRNAs, have been revealed to be involved in the epigenetic regulation of EC. In this review, we focus on describing new insights on epigenetic processes associated with noncoding RNAs, which have been characterized to be responsible for the development and progression of EC.

Identification of DNA-Methylated CpG Islands Associated With Gene Silencing in the Adult Body Tissues of the Ogye Chicken Using RNA-Seq and Reduced Representation Bisulfite Sequencing

DNA methylation is an epigenetic mark that plays an essential role in regulating gene expression. CpG islands are DNA methylations regions in promoters known to regulate gene expression through transcriptional silencing of the corresponding gene. DNA methylation at CpG islands is crucial for gene expression and tissue-specific processes. At the current time, a limited number of studies have reported on gene expression associated with DNA methylation in diverse adult tissues at the genome-wide level. Expression levels are rarely affected by DNA methylation in normal adult tissues; however, statistical differences in gene expression level correlated with DNA methylation have recently been revealed. In this study, we examined 20 pairs of DNA methylomes and transcriptomes from RNA-seq and reduced representation bisulfite sequencing (RRBS) data using adult Ogye chicken tissues. A total of 3,133 CpG islands were identified from 20 tissue data in a single chicken sample which could affect downstream genes. Analyzing these CpG island and gene pairs, 121 significant units were statistically correlated. Among them, six genes (CLDN3, DECR2, EVA1B, NME4, NTSR1, and XPNPEP2) were highly significantly changed by altered DNA methylation. Finally, our data demonstrated how DNA methylation correlated to gene expression in normal adult tissues. Our source codes can be found at https://github.com/wjlim/correlation-between-rna-seq-and-RRBS.

Genome-wide characterization of DNA methylation, small RNA expression, and histone H3 lysine nine di-methylation in Brassica rapa L

Epigenetic gene regulation is crucial to plant life and can involve dynamic interactions between various histone modifications, DNA methylation, and small RNAs. Detailed analysis of epigenome information is anticipated to reveal how the DNA sequence of the genome is translated into the plant’s phenotype. The aim of this study was to map the DNA methylation state at the whole genome level and to clarify the relationship between DNA methylation and transcription, small RNA expression, and histone H3 lysine 9 di-methylation (H3K9me2) in Brassica rapa. We performed whole genome bisulfite sequencing, small RNA sequencing, and chromatin immunoprecipitation sequencing using H3K9me2 antibody in a Chinese cabbage inbred line, RJKB-T24, and examined the impact of epigenetic states on transcription. Cytosine methylation in DNA was analysed in different sequence contexts (CG, CHG, and CHH) (where H could be A, C, or T) and position (promoter, exon, intron, terminator, interspersed repeat regions), and the H3K9me2 and 24 nucleotide small interfering RNAs (24 nt-siRNA) were overlaid onto the B. rapa reference genome. The epigenome was compared with that of Arabidopsis thaliana and the relationship between the position of DNA methylation and gene expression, and the involvement of 24 nt siRNAs and H3K9me2 are discussed.

The DNA Methylome and Association of Differentially Methylated Regions with Differential Gene Expression during Heat Stress in Brassica rapa

Cytosine DNA methylation is a critical epigenetic mechanism in the silencing of transposable elements, imprinting and regulating gene expression. However, little is known about the potential role of mC in response to heat stress. To determine and explore the functions of the dynamic DNA methylome during heat stress, we characterized single-base resolution methylome maps of Brassica rapa and assessed the dynamic changes of mC under heat stress using whole genome bisulfite sequencing. On average, the DNA methylation levels of CG, CHG and CHH are 39.3%, 15.38% and 5.24% in non-heading Chinese cabbage (NHCC), respectively. We found that the patterns of methylation are similar to other eudicot plants, but with higher CHH methylation levels. Further comparative analysis revealed varying patterns for three sequence contexts (mCG, mCHG and mCHH) under heat stress indicating context- and position-dependent methylation regulation. DNA methylation near the TSS and TES may be closely associated with methylation-dependent transcriptional silencing. Association analysis of differential methylation and differential gene expression revealed a different set of methDEGs involved at early and late stages under heat stress. The systemic characterization of the dynamic DNA methylome during heat stress will improve our understanding of the mechanism of epigenetic regulation under heat stress.

Functional Divergence between Subgenomes and Gene Pairs after Whole Genome Duplications

Gene loss following whole genome duplication (WGD) is often biased, with one subgenome retaining more ancestral genes and the other sustaining more gene deletions. While bias toward the greater expression of gene copies on one subgenome can explain bias in gene loss, this raises the question to what drives differences in gene expression levels between subgenomes. Differences in chromatin modifications and epigenetic markers between subgenomes in several model species are now being identified, providing an explanation for bias in gene expression between subgenomes. WGDs can be classified into duplications with higher, biased gene loss and bias in gene expression between subgenomes versus those with lower, unbiased rates of gene loss and an absence of detectable bias between subgenomes; however, the originally proposed link between these two classes and whether WGD results from an allo- or autopolyploid event is inconsistent with recent data from the allopolyploid Capsella bursa-pastoris. The gene balance hypothesis can explain bias in the functional categories of genes retained following WGD, the difference in gene loss rates between unbiased and biased WGDs, and how plant genomes have avoided being overrun with genes encoding dose-sensitive subunits of multiprotein complexes. Comparisons of gene expression patterns between retained transcription factor pairs in maize suggest the high degree of retention for WGD-derived pairs of transcription factors may instead be explained by the older duplication-degeneration-complementation model.

Comparative epigenomics reveals evolution of duplicated genes in potato and tomato

The evolution of duplicated genes after polyploidization has been the subject of many evolutionary biology studies. Potato (Solanum tuberosum) and tomato (Solanum lycopersicum) are the first two sequenced genomes of asterids, and share a common polyploidization event. However, the epigenetic role of DNA methylation on the evolution of duplicated genes derived from polyploidization is not fully understood. Here, we explore the role of the DNA methylation in the evolution of duplicated genes in potato and tomato. The overall levels of DNA methylation are different, although patterns of DNA methylation are similar in potato and tomato. Different types of duplicated genes can display different methylation patterns in potato and tomato. In addition, we found that differences in the methylation levels between duplicated genes were associated with gene expression divergence. In particular, for the majority of duplicated gene pairs, one copy is always hyper- or hypo-methylated compared with the other copy across different tomato fruit ripening stages, and these genes are enriched for specific function related to transcription factor (TF) activity. Furthermore, transcription of hundreds of duplicated TFs was shown to be regulated by DNA methylation during fruit ripening stages in tomato, some of which are well-known fruit ripening TFs. Taken together, our results support the notion that DNA methylation may facilitate divergent evolution of duplicated genes and play roles in important biological processes such as tomato fruit ripening.

Chromosome Evolution in Connection with Repetitive Sequences and Epigenetics in Plants

Chromosome evolution is a fundamental aspect of evolutionary biology. The evolution of chromosome size, structure and shape, number, and the change in DNA composition suggest the high plasticity of nuclear genomes at the chromosomal level. Repetitive DNA sequences, which represent a conspicuous fraction of every eukaryotic genome, particularly in plants, are found to be tightly linked with plant chromosome evolution. Different classes of repetitive sequences have distinct distribution patterns on the chromosomes. Mounting evidence shows that repetitive sequences may play multiple generative roles in shaping the chromosome karyotypes in plants. Furthermore, recent development in our understanding of the repetitive sequences and plant chromosome evolution has elucidated the involvement of a spectrum of epigenetic modification. In this review, we focused on the recent evidence relating to the distribution pattern of repetitive sequences in plant chromosomes and highlighted their potential relevance to chromosome evolution in plants. We also discussed the possible connections between evolution and epigenetic alterations in chromosome structure and repatterning, such as heterochromatin formation, centromere function, and epigenetic-associated transposable element inactivation.

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

Background

DNA methylation plays important roles in many regulatory processes in plants. It is economically infeasible to profile genome-wide DNA methylation at a single-base resolution in maize, given its genome size of ~2.5 Gb. As an alternative, we adapted region of interest (ROI)-directed reduced representation bisulfite sequencing (RRBS) to survey genome-wide methylation in maize.

Results

We developed a pipeline for selecting restriction enzymes in silico and experimentally showed that, in the maize genome, MseI- and CviQI-digested fragments are precisely enriched in promoters and gene bodies, respectively. We proceeded with comparisons of epigenomes and transcriptomes between shoots and tassels and found that the occurrences of highly methylated, tissue-specific, mCHH islands upstream of transcription start sites (TSSs) were positively correlated with differential gene expression. Furthermore, 5′ regulatory regions between TSS and mCHH islands often contain putative binding sites of known transcription factors (TFs) that regulate the flowering process and the timing of the transition from the vegetative to the reproductive phase. By integrating MNase-seq and siRNA-seq data, we found that regions of mCHH islands accumulate 21nt-siRNAs in a tissue-specific manner, marking the transition to open chromatin, thereby ensuring the accessibility of TFs for tissue-specific gene regulation.

Conclusions

Our ROI-directed RRBS pipeline is eminently applicable to DNA methylation profiling of large genomes. Our results provide novel insights into the tissue-specific epigenomic landscapes in maize, demonstrating that DNA methylation and siRNA and chromatin accessibility constitute a critical, interdependent component that orchestrates the transition from the vegetative to the reproductive phase.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-017-0148-y) contains supplementary material, which is available to authorized users.

Plant-RRBS, a bisulfite and next-generation sequencing-based methylome profiling method enriching for coverage of cytosine positions

BACKGROUND

Cytosine methylation in plant genomes is important for the regulation of gene transcription and transposon activity. Genome-wide methylomes are studied upon mutation of the DNA methyltransferases, adaptation to environmental stresses or during development. However, from basic biology to breeding programs, there is a need to monitor multiple samples to determine transgenerational methylation inheritance or differential cytosine methylation. Methylome data obtained by sodium hydrogen sulfite (bisulfite)-conversion and next-generation sequencing (NGS) provide genome-wide information on cytosine methylation. However, a profiling method that detects cytosine methylation state dispersed over the genome would allow high-throughput analysis of multiple plant samples with distinct epigenetic signatures. We use specific restriction endonucleases to enrich for cytosine coverage in a bisulfite and NGS-based profiling method, which was compared to whole-genome bisulfite sequencing of the same plant material.

METHODS

We established an effective methylome profiling method in plants, termed plant-reduced representation bisulfite sequencing (plant-RRBS), using optimized double restriction endonuclease digestion, fragment end repair, adapter ligation, followed by bisulfite conversion, PCR amplification and NGS. We report a performant laboratory protocol and a straightforward bioinformatics data analysis pipeline for plant-RRBS, applicable for any reference-sequenced plant species.

RESULTS

As a proof of concept, methylome profiling was performed using an Oryza sativa ssp. indica pure breeding line and a derived epigenetically altered line (epiline). Plant-RRBS detects methylation levels at tens of millions of cytosine positions deduced from bisulfite conversion in multiple samples. To evaluate the method, the coverage of cytosine positions, the intra-line similarity and the differential cytosine methylation levels between the pure breeding line and the epiline were determined. Plant-RRBS reproducibly covers commonly up to one fourth of the cytosine positions in the rice genome when using MspI-DpnII within a group of five biological replicates of a line. The method predominantly detects cytosine methylation in putative promoter regions and not-annotated regions in rice.

CONCLUSIONS

Plant-RRBS offers high-throughput and broad, genome-dispersed methylation detection by effective read number generation obtained from reproducibly covered genome fractions using optimized endonuclease combinations, facilitating comparative analyses of multi-sample studies for cytosine methylation and transgenerational stability in experimental material and plant breeding populations.

Epigenomic diversification within the genus Lupinus

Deciphering the various chemical modifications of both DNA and the histone compound of chromatin not only leads to a better understanding of the genome-wide organisation of epigenetic landmarks and their impact on gene expression but may also provide some insights into the evolutionary processes. Although both histone modifications and DNA methylation have been widely investigated in various plant genomes, here we present the first study for the genus Lupinus. Lupins, which are members of grain legumes (pulses), are beneficial for food security, nutrition, health and the environment. In order to gain a better understanding of the epigenetic organisation of genomes in lupins we applied the immunostaining of methylated histone H3 and DNA methylation as well as whole-genome bisulfite sequencing. We revealed variations in the patterns of chromatin modifications at the chromosomal level among three crop lupins, i.e. L. angustifolius (2n = 40), L. albus (2n = 50) and L. luteus (2n = 52), and the legume model plant Medicago truncatula (2n = 16). Different chromosomal patterns were found depending on the specific modification, e.g. H3K4me2 was localised in the terminal parts of L. angustifolius and M. truncatula chromosomes, which is in agreement with the results that have been obtained for other species. Interestingly, in L. albus and L. luteus this modification was limited to one arm in the case of all of the chromosomes in the complement. Additionally, H3K9me2 was detected in all of the analysed species except L. luteus. DNA methylation sequencing (CG, CHG and CHH contexts) of aforementioned crop but also wild lupins such as L. cosentinii (2n = 32), L. digitatus (2n = 36), L. micranthus (2n = 52) and L. pilosus (2n = 42) supported the range of interspecific diversity. The examples of epigenetic modifications illustrate the diversity of lupin genomes and could be helpful for elucidating further epigenetic changes in the evolution of the lupin genome.

Facilitation of environmental adaptation and evolution by epigenetic phenotype variation: insights from clonal, invasive, polyploid, and domesticated animals

There is increasing evidence, particularly from plants, that epigenetic mechanisms can contribute to environmental adaptation and evolution. The present article provides an overview on this topic for animals and highlights the special suitability of clonal, invasive, hybrid, polyploid, and domesticated species for environmental and evolutionary epigenetics. Laboratory and field studies with asexually reproducing animals have shown that epigenetically diverse phenotypes can be produced from the same genome either by developmental stochasticity or environmental induction. The analysis of invasions revealed that epigenetic phenotype variation may help to overcome genetic barriers typically associated with invasions such as bottlenecks and inbreeding. Research with hybrids and polyploids established that epigenetic mechanisms are involved in consolidation of speciation by contributing to reproductive isolation and restructuring of the genome in the neo-species. Epigenetic mechanisms may even have the potential to trigger speciation but evidence is still meager. The comparison of domesticated animals and their wild ancestors demonstrated heritability and selectability of phenotype modulating DNA methylation patterns. Hypotheses, model predictions, and empirical results are presented to explain how epigenetic phenotype variation could facilitate adaptation and speciation. Clonal laboratory lineages, monoclonal invaders, and adaptive radiations of different evolutionary age seem particularly suitable to empirically test the proposed ideas. A respective research agenda is presented.

Profiling genome-wide DNA methylation

DNA methylation is an epigenetic modification that plays an important role in regulating gene expression and therefore a broad range of biological processes and diseases. DNA methylation is tissue-specific, dynamic, sequence-context-dependent and trans-generationally heritable, and these complex patterns of methylation highlight the significance of profiling DNA methylation to answer biological questions. In this review, we surveyed major methylation assays, along with comparisons and biological examples, to provide an overview of DNA methylation profiling techniques. The advances in microarray and sequencing technologies make genome-wide profiling possible at a single-nucleotide or even a single-cell resolution. These profiling approaches vary in many aspects, such as DNA input, resolution, genomic region coverage, and bioinformatics analysis, and selecting a feasible method requires knowledge of these methods. We first introduce the biological background of DNA methylation and its pattern in plants, animals and fungi. We present an overview of major experimental approaches to profiling genome-wide DNA methylation and hydroxymethylation and then extend to the single-cell methylome. To evaluate these methods, we outline their strengths and weaknesses and perform comparisons across the different platforms. Due to the increasing need to compute high-throughput epigenomic data, we interrogate the computational pipeline for bisulfite sequencing data and also discuss the concept of identifying differentially methylated regions (DMRs). This review summarizes the experimental and computational concepts for profiling genome-wide DNA methylation, followed by biological examples. Overall, this review provides researchers useful guidance for the selection of a profiling method suited to specific research questions.