A Comparative Epigenomic Analysis of Polyploidy-Derived Genes in Soybean and Common Bean

DNA methylation dynamics in gymnosperm duplicate genes: implications for genome evolution and stress adaptation

Duplicate genes are pivotal in driving evolutionary innovation, often exhibiting expression divergence that offers a system to investigate the role of DNA methylation in transcriptional regulation. However, previous studies have predominantly focused on angiosperms, leaving the methylation patterns in major lineages of land plants still unclear. This study explores DNA methylation evolution in duplicate genes across representative gymnosperm species with large genomes, spanning over 300 million years, using genomic, transcriptomic, and high-depth DNA methylomic data. We observed variations in DNA methylation levels along gene bodies, flanking regions, and methylation statuses of coding regions across different duplication types. Biased divergences in DNA methylation and gene expression frequently occurred between duplicate copies. Specifically, methylation divergences in the 2-kb downstream regions negatively correlated with gene expression. Both CG and CHG DNA methylation in gene bodies were positively correlated with gene length, suggesting these methylation types may function as an epigenomic buffer to mitigate the adverse impact of gene length on expression. Duplicate genes exhibiting both methylation and expression divergences were notably enriched in adaptation-related biological processes, suggesting that DNA methylation may aid adaptive evolution in gymnosperms by regulating stress response genes. Changes in expression levels correlated with switches in methylation status within coding regions of transposed duplicates. Specifically, depletion for CG methylation or enrichment for non-CG methylation significantly reduced the expression of translocated copies. This correlation suggests that DNA methylation may reduce genetic redundancy by silencing translocated copies. Our study highlights the significance of DNA methylation in plant genome evolution and stress adaptation.

Epigenetics in the modern era of crop improvements

Epigenetic mechanisms are integral to plant growth, development, and adaptation to environmental stimuli. Over the past two decades, our comprehension of these complex regulatory processes has expanded remarkably, producing a substantial body of knowledge on both locus-specific mechanisms and genome-wide regulatory patterns. Studies initially grounded in the model plant Arabidopsis have been broadened to encompass a diverse array of crop species, revealing the multifaceted roles of epigenetics in physiological and agronomic traits. With recent technological advancements, epigenetic regulations at the single-cell level and at the large-scale population level are emerging as new focuses. This review offers an in-depth synthesis of the diverse epigenetic regulations, detailing the catalytic machinery and regulatory functions. It delves into the intricate interplay among various epigenetic elements and their collective influence on the modulation of crop traits. Furthermore, it examines recent breakthroughs in technologies for epigenetic modifications and their integration into strategies for crop improvement. The review underscores the transformative potential of epigenetic strategies in bolstering crop performance, advocating for the development of efficient tools to fully exploit the agricultural benefits of epigenetic insights.

Harnessing Multi-Omics Strategies and Bioinformatics Innovations for Advancing Soybean Improvement: A Comprehensive Review

Soybean improvement has entered a new era with the advent of multi-omics strategies and bioinformatics innovations, enabling more precise and efficient breeding practices. This comprehensive review examines the application of multi-omics approaches in soybean-encompassing genomics, transcriptomics, proteomics, metabolomics, epigenomics, and phenomics. We first explore pre-breeding and genomic selection as tools that have laid the groundwork for advanced trait improvement. Subsequently, we dig into the specific contributions of each -omics field, highlighting how bioinformatics tools and resources have facilitated the generation and integration of multifaceted data. The review emphasizes the power of integrating multi-omics datasets to elucidate complex traits and drive the development of superior soybean cultivars. Emerging trends, including novel computational techniques and high-throughput technologies, are discussed in the context of their potential to revolutionize soybean breeding. Finally, we address the challenges associated with multi-omics integration and propose future directions to overcome these hurdles, aiming to accelerate the pace of soybean improvement. This review serves as a crucial resource for researchers and breeders seeking to leverage multi-omics strategies for enhanced soybean productivity and resilience.

Divergent evolution of NLR genes in the genus Glycine: impacts of annuals and perennials' life history strategies

Within the family Fabaceae, the genus Glycine is composed of two subgenera annuals (2n=40) and perennials. This life strategy transition may have differentially affected the evolution of various gene families. Its cultivated species G. max has high level of susceptibility to major pathogens including viruses, bacteria and fungi. Understanding nucleotide-binding domain leucine-rich repeat (NLR) genes evolution in soybean is of paramount importance due to their central role in plant immunity and their potential in improving disease resistance in soybean cultivars. In this study, we investigated the significance of this annual-perennial transition on the macroevolution of NLR genes in the genus Glycine. Our results reveal a remarkable distinction between annual species such as Glycine max and Glycine soja, which exhibit an expanded NLRome compared to perennial species (G. cyrtoloba, G. stenophita, G. dolichocarpa, G. falcata, G. syndetika, G. latifolia and G. tomentella). Our evolutionary timescale analysis pinpoints recent accelerated gene duplication events for this expansion, which occurred between 0.1 and 0.5 million years ago, driven predominantly by lineage-specific and terminal duplications. In contrast, perennials initially experienced significant contraction during the diploidisation phase following the Glycine-specific whole-genome duplication event (~10 million years ago). Despite the reduction in the NLRome, perennial lineages exhibit a unique and highly diversified repertoire of NLR genes with limited interspecies synteny. The investigation of gene gain and loss ratios revealed that this diversification resulted from the birth of novel genes following individual speciation events. Among perennials, G. latifolia, a well-known resistance resource, has the highest ratio of these novel genes in the tertiary gene pool. Our study suggests evolutionary mechanisms, including recombination and transposition, as potential drivers for the emergence of these novel genes. This study also provides evidence for the unbalanced expansion of the NLRome in the Dt subgenome compared with the At subgenome in the young allopolyploid G. dolichocarpa. To the best of our knowledge, this is the first study to investigate the effect of annuality and perenniality life transition on the evolution of NLR genes in the genus Glycine to identify its genomics resources for improving the resistance of soybean crop with global importance on the economy and food security.

Epigenetic variation in early and late flowering plants of the rubber-producing Russian dandelion Taraxacum koksaghyz provides insights into the regulation of flowering time

The Russian dandelion (Taraxacum koksaghyz) grows in temperate zones and produces large amounts of poly(cis-1,4-isoprene) in its roots, making it an attractive alternative source of natural rubber. Most T. koksaghyz plants require vernalization to trigger flower development, whereas early flowering varieties that have lost their vernalization dependence are more suitable for breeding and domestication. To provide insight into the regulation of flowering time in T. koksaghyz, we induced epigenetic variation by in vitro cultivation and applied epigenomic and transcriptomic analysis to the resulting early flowering plants and late flowering controls, allowing us to identify differences in methylation patterns and gene expression that correlated with flowering. This led to the identification of candidate genes homologous to vernalization and photoperiodism response genes in other plants, as well as epigenetic modifications that may contribute to the control of flower development. Some of the candidate genes were homologous to known floral regulators, including those that directly or indirectly regulate the major flowering control gene FT. Our atlas of genes can be used as a starting point to investigate mechanisms that control flowering time in T. koksaghyz in greater detail and to develop new breeding varieties that are more suited to domestication.

Frequency of cytosine methylation in the adjacent regions of soybean retrotransposon SORE-1 depends on chromosomal location

Mobilization of transposable elements (TEs) is suppressed by epigenetic mechanisms involving cytosine methylation. However, few studies have focused on clarifying relationships between epigenetic influences of TEs on the adjacent DNA regions and time after insertion of TEs into the genome and/or their chromosomal location. Here we addressed these issues using soybean retrotransposon SORE-1. We analyzed SORE-1, inserted in exon 1 of the GmphyA2 gene, one of the newest insertions in this family so far identified. Cytosine methylation was detected in this element but was barely present in the adjacent regions. These results were correlated, respectively, with the presence and absence of the production of short interfering RNAs. Cytosine methylation profiles of 74 SORE-1 elements in the Williams 82 reference genome indicated that methylation frequency in the adjacent regions of SORE-1 was profoundly higher in pericentromeric regions than in euchromatic chromosome arms and was only weakly correlated with the length of time after insertion into the genome. Notably, the higher level of methylation in the 5' adjacent regions of SORE-1 coincided with the presence of repetitive elements in pericentromeric regions. Together, these results suggest that epigenetic influence of SORE-1 on the adjacent regions is influenced by its location on the chromosome.

Evolutionary insights into 3D genome organization and epigenetic landscape of Vigna mungo

Eukaryotic genomes show an intricate three-dimensional (3D) organization within the nucleus that regulates multiple biological processes including gene expression. Contrary to animals, understanding of 3D genome organization in plants remains at a nascent stage. Here, we investigate the evolution of 3D chromatin architecture in legumes. By using cutting-edge PacBio, Illumina, and Hi-C contact reads, we report a gap-free, chromosome-scale reference genome assembly of Vigna mungo, an important minor legume cultivated in Southeast Asia. We spatially resolved V. mungo chromosomes into euchromatic, transcriptionally active A compartment and heterochromatic, transcriptionally-dormant B compartment. We report the presence of TAD-like-regions throughout the diagonal of the HiC matrix that resembled transcriptional quiescent centers based on their genomic and epigenomic features. We observed high syntenic breakpoints but also high coverage of syntenic sequences and conserved blocks in boundary regions than in the TAD-like region domains. Our findings present unprecedented evolutionary insights into spatial 3D genome organization and epigenetic patterns and their interaction within the V. mungo genome. This will aid future genomics and epigenomics research and breeding programs of V. mungo.

Role of a LORELEI- like gene from Phaseolus vulgaris during a mutualistic interaction with Rhizobium tropici

Reactive oxygen species (ROS), produced by NADPH oxidases known as RBOHs in plants, play a key role in plant development, biotic and abiotic stress responses, hormone signaling, and reproduction. Among the subfamily of receptor-like kinases referred to as CrRLK, there is FERONIA (FER), a regulator of RBOHs, and FER requires a GPI-modified membrane protein produced by LORELEI (LRE) or LORELEI-like proteins (LLG) to reach the plasma membrane and generate ROS. In Arabidopsis, AtLLG1 is involved in interactions with microbes as AtLLG1 interacts with the flagellin receptor (FLS2) to trigger the innate immune response, but the role of LLGs in mutualistic interactions has not been examined. In this study, two Phaseolus vulgaris LLG genes were identified, PvLLG2 that was expressed in floral tissue and PvLLG1 that was expressed in vegetative tissue. Transcripts of PvLLG1 increased during rhizobial nodule formation peaking during the early period of well-developed nodules. Also, P. vulgaris roots expressing pPvLLG1:GFP-GUS showed that this promoter was highly active during rhizobium infections, and very similar to the subcellular localization using a construct pLLG1::PvLLG1-Neon. Compared to control plants, PvLLG1 silenced plants had less superoxide (O2-) at the root tip and elongation zone, spotty hydrogen peroxide (H2O2) in the elongation root zone, and significantly reduced root hair length, nodule number and nitrogen fixation. Unlike control plants, PvLLG1 overexpressing plants showed superoxide beyond the nodule meristem, and significantly increased nodule number and nodule diameter. PvLLG1 appears to play a key role during this mutualistic interaction, possibly due to the regulation of the production and distribution of ROS in roots.

Epigenomic divergence correlates with sequence polymorphism in Arabidopsis paralogs

Processes affecting rates of sequence polymorphism are fundamental to the evolution of gene duplicates. The relationship between gene activity and sequence polymorphism can influence the likelihood that functionally redundant gene copies are co-maintained in stable evolutionary equilibria vs other outcomes such as neofunctionalization. Here, we investigate genic variation in epigenome-associated polymorphism rates in Arabidopsis thaliana and consider whether these affect the evolution of gene duplicates. We compared the frequency of sequence polymorphism and patterns of genetic differentiation between genes classified by exon methylation patterns: unmethylated (unM), gene-body methylated (gbM), and transposon-like methylated (teM) states, which reflect divergence in gene expression. We found that the frequency of polymorphism was higher in teM (transcriptionally repressed, tissue-specific) genes and lower in gbM (active, constitutively expressed) genes. Comparisons of gene duplicates were largely consistent with genome-wide patterns - gene copies that exhibit teM accumulate more variation, evolve faster, and are in chromatin states associated with reduced DNA repair. This relationship between expression, the epigenome, and polymorphism may lead to the breakdown of equilibrium states that would otherwise maintain genetic redundancies. Epigenome-mediated polymorphism rate variation may facilitate the evolution of novel gene functions in duplicate paralogs maintained over evolutionary time.

Large-Scale Identification and Characterization Analysis of VQ Family Genes in Plants, Especially Gymnosperms

VQ motif-containing (VQ) proteins are a class of transcription regulatory cofactors widely present in plants, playing crucial roles in growth and development, stress response, and defense. Although there have been some reports on the member identification and functional research of VQ genes in some plants, there is still a lack of large-scale identification and clear graphical presentation of their basic characterization information to help us to better understand this family. Especially in gymnosperms, the VQ family genes and their evolutionary relationships have not yet been reported. In this study, we systematically identified 2469 VQ genes from 56 plant species, including bryophytes, gymnosperms, and angiosperms, and analyzed their molecular and evolutionary features. We found that amino acids are only highly conserved in the VQ domain, while other positions are relatively variable; most VQ genes encode relatively small proteins and do not have introns. The GC content in Poaceae plants is the highest (up to 70%); these VQ proteins can be divided into nine subgroups. In particular, we analyzed the molecular characteristics, chromosome distribution, duplication events, and expression levels of VQ genes in three gymnosperms: Ginkgo biloba, Taxus chinensis, and Pinus tabuliformis. In gymnosperms, VQ genes are classified into 11 groups, with highly similar motifs in each group; most VQ proteins have less than 300 amino acids and are predicted to be located in nucleus. Tandem duplication is an important driving force for the expansion of the VQ gene family, and the evolutionary processes of most VQ genes and duplication events are relatively independent; some candidate VQ genes are preliminarily screened, and they are likely to be involved in plant growth and stress and defense responses. These results provide detailed information and powerful references for further understanding and utilizing the VQ family genes in various plants.

Colchicine-Induced Polyploidy in Leguminous Crops Enhances Morpho-Physiological Characteristics for Drought Stress Tolerance

Legumes play a significant role in the alleviation of food insecurity, maintaining soil fertility, and achieving sustainable crop production under adverse environmental conditions. The increased demand in legume production contemplates that attention on the genetic improvement of these crops through various means such as genetic engineering and mutation breeding should take a centre stage in global agriculture. Therefore, this paper provides a succinct analysis of the currently available literature on morphological and physiological traits in polyploidised leguminous plants to counter the adverse effects of drought stress. The effects of colchicine on various morphological and physiological traits of polyploidised legumes compared to their diploid counterparts were examined. Numerous reports revealed variations in these traits, such as improved root and shoot growth, plant biomass, chloroplastidic content, protein, RNA, and DNA. The differences observed were also associated with the strong relationship between plant ploidy induction and colchicine application. Furthermore, the analysis indicated that polyploidisation remains dose-dependent and may be achievable within a shorter space of time as this antimitotic chemical interferes with chromosome separations in somatic plant cells. The efficiency of this process also depends on the advancement of treatment conditions (in vitro, in vivo, or ex vitro) and the successful regeneration of polyploidised plants for adaptation under drought stress conditions. As such, the improvement in metabolite profile and other essential growth characteristics serves as a clear indication that induced polyploidy needs to be further explored to confer resilience to environmental stress and improve crop yield under drought stress conditions in leguminous plants.

Polyploidy and diploidization in soybean

Polyploidy is widespread and particularly common in angiosperms. The prevalence of polyploidy in the plant suggests it as a crucial driver of diversification and speciation. The paleopolyploid soybean (Glycine max) is one of the most important crops of plant protein and oil for humans and livestock. Soybean experienced two rounds of whole genome duplication around 13 and 59 million years ago. Due to the relatively slow process of post-polyploid diploidization, most genes are present in multiple copies across the soybean genome. Growing evidence suggests that polyploidization and diploidization could cause rapid and dramatic changes in genomic structure and epigenetic modifications, including gene loss, transposon amplification, and reorganization of chromatin architecture. This review is focused on recent progresses about genetic and epigenetic changes during polyploidization and diploidization of soybean and represents the challenges and potentials for application of polyploidy in soybean breeding.

6mA DNA Methylation on Genes in Plants Is Associated with Gene Complexity, Expression and Duplication

N6-methyladenine (6mA) DNA methylation has emerged as an important epigenetic modification in eukaryotes. Nevertheless, the evolution of the 6mA methylation of homologous genes after species and after gene duplications remains unclear in plants. To understand the evolution of 6mA methylation, we detected the genome-wide 6mA methylation patterns of four lotus plants (Nelumbo nucifera) from different geographic origins by nanopore sequencing and compared them to patterns in Arabidopsis and rice. Within lotus, the genomic distributions of 6mA sites are different from the widely studied 5mC methylation sites. Consistently, in lotus, Arabidopsis and rice, 6mA sites are enriched around transcriptional start sites, positively correlated with gene expression levels, and preferentially retained in highly and broadly expressed orthologs with longer gene lengths and more exons. Among different duplicate genes, 6mA methylation is significantly more enriched and conserved in whole-genome duplicates than in local duplicates. Overall, our study reveals the convergent patterns of 6mA methylation evolution based on both lineage and duplicate gene divergence, which underpin their potential role in gene regulatory evolution in plants.

DNA methylation signatures of duplicate gene evolution in angiosperms

Gene duplication is a source of evolutionary novelty. DNA methylation may play a role in the evolution of duplicate genes (paralogs) through its association with gene expression. While this relationship has been examined to varying extents in a few individual species, the generalizability of these results at either a broad phylogenetic scale with species of differing duplication histories or across a population remains unknown. We applied a comparative epigenomics approach to 43 angiosperm species across the phylogeny and a population of 928 Arabidopsis (Arabidopsis thaliana) accessions, examining the association of DNA methylation with paralog evolution. Genic DNA methylation was differentially associated with duplication type, the age of duplication, sequence evolution, and gene expression. Whole genome duplicates were typically enriched for CG-only gene-body methylated or unmethylated genes, while single-gene duplications were typically enriched for non-CG methylated or unmethylated genes. Non-CG methylation, in particular, was characteristic of more recent single-gene duplicates. Core angiosperm gene families differentiated into those which preferentially retain paralogs and 'duplication-resistant' families, which convergently reverted to singletons following duplication. Duplication-resistant families that still have paralogous copies were, uncharacteristically for core angiosperm genes, enriched for non-CG methylation. Non-CG methylated paralogs had higher rates of sequence evolution, higher frequency of presence-absence variation, and more limited expression. This suggests that silencing by non-CG methylation may be important to maintaining dosage following duplication and be a precursor to fractionation. Our results indicate that genic methylation marks differing evolutionary trajectories and fates between paralogous genes and have a role in maintaining dosage following duplication.

Impact of Polyploidy Induction for Salinity Stress Mitigation in Soybean (Glycine max L. Merrill)

Polyploidy induction is recognized as one of the major evolutionary processes leading to remarkable morphological, physiological, and genetic variations in plants. Soybean (Glycine max L.), also known as soja bean or soya bean, is an annual leguminous crop of the pea family (Fabaceae) that shares a paleopolypoidy history, dating back to approximately 56.5 million years ago with other leguminous crops such as cowpea and other Glycine specific polyploids. This crop has been documented as one of the polyploid complex species among legumes whose gene evolution and resultant adaptive growth characteristics following induced polyploidization has not been fully explored. Furthermore, no successfully established in vivo or in vitro based polyploidy induction protocols have been reported to date, particularly, with the intention to develop mutant plants showing strong resistance to abiotic salinity stress. This review, therefore, describes the role of synthetic polyploid plant production in soybean for the mitigation of high soil salt stress levels and how this evolving approach could be used to further enhance the nutritional, pharmaceutical and economic industrial value of soybeans. This review also addresses the challenges involved during the polyploidization process.

Evolutionary fates of gene-body methylation and its divergent association with gene expression in pigeonpea

Pigeonpea (Cajanus cajan L. Huth) is an agronomically important legume cultivated worldwide. In this study, we extensively analyzed gene-body methylation (GbM) patterns in pigeonpea. We found a bimodal distribution of CG and CHG methylation patterns. GbM features- slow evolution rate and increased length remained conserved. Genes with moderate CG body methylation showed highest expression where as highly-methylated genes showed lowest expression. Transposable element (TE)-related genes were methylated in multiple contexts and hence classified as C-methylated genes. A low expression among C-methylated genes was associated with transposons insertion in gene-body and upstream regulatory regions. The CG methylation patterns were found to be conserved in orthologs compared with non-CG methylation. By comparing methylation patterns between differentially methylated regions (DMRs) of the three genotypes, we found that variably methylated marks are less likely to target evolutionary conserved sequences. Finally, our analysis showed enrichment of nitrogen-related genes in GbM orthologs of legumes, which could be promising candidates for generating epialleles for crop improvement.

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.

A Comprehensive Study of the WRKY Transcription Factor Family in Strawberry

WRKY transcription factors play critical roles in plant growth and development or stress responses. Using up-to-date genomic data, a total of 64 and 257 WRKY genes have been identified in the diploid woodland strawberry, Fragaria vesca, and the more complex allo-octoploid commercial strawberry, Fragaria × ananassa cv. Camarosa, respectively. The completeness of the new genomes and annotations has enabled us to perform a more detailed evolutionary and functional study of the strawberry WRKY family members, particularly in the case of the cultivated hybrid, in which homoeologous and paralogous FaWRKY genes have been characterized. Analysis of the available expression profiles has revealed that many strawberry WRKY genes show preferential or tissue-specific expression. Furthermore, significant differential expression of several FaWRKY genes has been clearly detected in fruit receptacles and achenes during the ripening process and pathogen challenged, supporting a precise functional role of these strawberry genes in such processes. Further, an extensive analysis of predicted development, stress and hormone-responsive cis-acting elements in the strawberry WRKY family is shown. Our results provide a deeper and more comprehensive knowledge of the WRKY gene family in strawberry.

Progress in soybean functional genomics over the past decade

Soybean is one of the most important oilseed and fodder crops. Benefiting from the efforts of soybean breeders and the development of breeding technology, large number of germplasm has been generated over the last 100 years. Nevertheless, soybean breeding needs to be accelerated to meet the needs of a growing world population, to promote sustainable agriculture and to address future environmental changes. The acceleration is highly reliant on the discoveries in gene functional studies. The release of the reference soybean genome in 2010 has significantly facilitated the advance in soybean functional genomics. Here, we review the research progress in soybean omics (genomics, transcriptomics, epigenomics and proteomics), germplasm development (germplasm resources and databases), gene discovery (genes that are responsible for important soybean traits including yield, flowering and maturity, seed quality, stress resistance, nodulation and domestication) and transformation technology during the past decade. At the end, we also briefly discuss current challenges and future directions.

Fab Advances in Fabaceae for Abiotic Stress Resilience: From 'Omics' to Artificial Intelligence

Legumes are a better source of proteins and are richer in diverse micronutrients over the nutritional profile of widely consumed cereals. However, when exposed to a diverse range of abiotic stresses, their overall productivity and quality are hugely impacted. Our limited understanding of genetic determinants and novel variants associated with the abiotic stress response in food legume crops restricts its amelioration. Therefore, it is imperative to understand different molecular approaches in food legume crops that can be utilized in crop improvement programs to minimize the economic loss. 'Omics'-based molecular breeding provides better opportunities over conventional breeding for diversifying the natural germplasm together with improving yield and quality parameters. Due to molecular advancements, the technique is now equipped with novel 'omics' approaches such as ionomics, epigenomics, fluxomics, RNomics, glycomics, glycoproteomics, phosphoproteomics, lipidomics, regulomics, and secretomics. Pan-omics-which utilizes the molecular bases of the stress response to identify genes (genomics), mRNAs (transcriptomics), proteins (proteomics), and biomolecules (metabolomics) associated with stress regulation-has been widely used for abiotic stress amelioration in food legume crops. Integration of pan-omics with novel omics approaches will fast-track legume breeding programs. Moreover, artificial intelligence (AI)-based algorithms can be utilized for simulating crop yield under changing environments, which can help in predicting the genetic gain beforehand. Application of machine learning (ML) in quantitative trait loci (QTL) mining will further help in determining the genetic determinants of abiotic stress tolerance in pulses.

Epimutations Define a Fast-Ticking Molecular Clock in Plants

Stochastic gains and losses of DNA methylation at CG dinucleotides are a frequent occurrence in plants. These spontaneous 'epimutations' occur at a rate that is 100 000 times higher than the genetic mutation rate, are effectively neutral at the genome-wide scale, and are stably inherited across mitotic and meiotic cell divisions. Mathematical models have been extraordinarily successful at describing how epimutations accumulate in plant genomes over time, making this process one of the most predictable epigenetic phenomena to date. Here, we propose that their high rate and effective neutrality make epimutations a powerful new molecular clock for timing evolutionary events of the recent past and for age dating of long-lived perennials such as trees.

Altered chromatin architecture and gene expression during polyploidization and domestication of soybean

Polyploidy or whole-genome duplication (WGD) is widespread in plants and is a key driver of evolution and speciation, accompanied by rapid and dynamic changes in genomic structure and gene expression. The 3D structure of the genome is intricately linked to gene expression, but its role in transcription regulation following polyploidy and domestication remains unclear. Here, we generated high-resolution (∼2 kb) Hi-C maps for cultivated soybean (Glycine max), wild soybean (Glycine soja), and common bean (Phaseolus vulgaris). We found polyploidization in soybean may induce architecture changes of topologically associating domains and subsequent diploidization led to chromatin topology alteration around chromosome-rearrangement sites. Compared with single-copy and small-scale duplicated genes, WGD genes displayed more long-range chromosomal interactions and were coupled with higher levels of gene expression and chromatin accessibilities but void of DNA methylation. Interestingly, chromatin loop reorganization was involved in expression divergence of the genes during soybean domestication. Genes with chromatin loops were under stronger artificial selection than genes without loops. These findings provide insights into the roles of dynamic chromatin structures on gene expression during polyploidization, diploidization, and domestication of soybean.

Intertwined evolution of plant epigenomes and genomes

DNA methylation is found across eukaryotes; however, plants have evolved patterns and pathways of DNA methylation that are distinct from animals and fungi. DNA methylation shapes the evolution of genomes through its direct roles in transposon silencing, gene expression, genome stability, and its impact on mutation rates. In return the diversity of DNA methylation across species is shaped by genome sequence evolution. Extensive diversification of key DNA methylation pathways has continued in plants through gene duplication and loss. Meanwhile, frequent movement of transposons has altered local DNA methylation patterns and the genes affected. Only recently has the diversity and evolutionary history of plant DNA methylation become evident with the availability of increasing genomic and epigenomic data. However, much remains unresolved regarding the evolutionary forces that have shaped the dynamics of the complex and intertwined history of plant genome and epigenome evolution.

Evolutionary Contribution of Duplicated Genes to Genome Evolution in the Ginseng Species Complex

Genes duplicated by whole genome duplication (WGD) and small-scale duplication (SSD) have played important roles in adaptive evolution of all flowering plants. However, it still remains underinvestigated how the distinct models of duplication events and their contending evolutionary patterns have shaped the genome and epigenomes of extant plant species. In this study, we investigated the contribution of the WGD- and SSD-derived duplicate genes to the genome evolution of one diploid and three closely related allotetraploid Panax species based on genome, methylome, and proteome data sets. Our genome-wide comparative analyses revealed that although the ginseng species complex was recently diverged, they have evolved distinct overall patterns of nucleotide variation, cytosine methylation, and protein-level expression. In particular, genetic and epigenetic asymmetries observed in the recent WGD-derived genes are largely consistent across the ginseng species complex. In addition, our results revealed that gene duplicates generated by ancient WGD and SSD mechanisms exhibited distinct evolutionary patterns. We found the ancient WGD-derived genes (i.e., ancient collinear gene) are genetically more conserved and hypomethylated at the cytosine sites. In contrast, some of the SSD-derived genes (i.e., dispersal duplicated gene) showed hypermethylation and high variance in nucleotide variation pattern. Functional enrichment analyses of the duplicated genes indicated that adaptation-related traits (i.e., photosynthesis) created during the distant ancient WGDs are further strengthened by both the more recent WGD and SSD. Together, our findings suggest that different types of duplicated genes may have played distinct but relaying evolutionary roles in the polyploidization and speciation processes in the ginseng species complex.

A new decade and new data at SoyBase, the USDA-ARS soybean genetics and genomics database

SoyBase, a USDA genetic and genomics database, holds professionally curated soybean genetic and genomic data, which is integrated and made accessible to researchers and breeders. The site holds several reference genome assemblies, as well as genetic maps, thousands of mapped traits, expression and epigenetic data, pedigree information, and extensive variant and genotyping data sets. SoyBase displays include genetic, genomic, and epigenetic maps of the soybean genome. Gene expression data is presented in the genome viewer as heat maps and pictorial and tabular displays in gene report pages. Millions of sequence variants have been added, representing variations across various collections of cultivars. This variant data is explorable using new interactive tools to visualize the distribution of those variants across the genome, between selected accessions. SoyBase holds several reference-quality soybean genome assemblies, accessible via various query tools and browsers, including a new visualization system for exploring the soybean pan-genome. SoyBase also serves as a nexus of announcements pertinent to the greater soybean research community. The database also includes a soybean-specific anatomic and biochemical trait ontology. The database can be accessed at https://soybase.org.

Suppressed Methionine γ-Lyase Expression Causes Hyperaccumulation of S-Methylmethionine in Soybean Seeds

Several soybean (Glycine max) germplasms, such as Nishiyamahitashi 98-5 (NH), have an intense seaweed-like flavor after cooking because of their high seed S-methylmethionine (SMM) content. In this study, we compared the amounts of amino acids in the phloem sap, leaves, pods, and seeds between NH and the common soybean cultivar Fukuyutaka. This revealed a comparably higher SMM content alongside a higher free Met content in NH seeds, suggesting that the SMM-hyperaccumulation phenotype of NH soybean was related to Met metabolism in seeds. To investigate the molecular mechanism behind SMM hyperaccumulation, we examined the phenotype-associated gene locus in NH plants. Analyses of the quantitative trait loci in segregated offspring of the cross between NH and the common soybean cultivar Williams 82 indicated that one locus on chromosome 10 explains 71.4% of SMM hyperaccumulation. Subsequent fine-mapping revealed that a transposon insertion into the intron of a gene, Glyma.10g172700, is associated with the SMM-hyperaccumulation phenotype. The Glyma.10g172700-encoded recombinant protein showed Met-γ-lyase (MGL) activity in vitro, and the transposon-insertion mutation in NH efficiently suppressed Glyma.10g172700 expression in developing seeds. Exogenous administration of Met to sections of developing soybean seeds resulted in transient increases in Met levels, followed by continuous increases in SMM concentrations, which was likely caused by Met methyltransferase activity in the seeds. Accordingly, we propose that the SMM-hyperaccumulation phenotype is caused by suppressed MGL expression in developing soybean seeds, resulting in transient accumulation of Met, which is converted into SMM to avoid the harmful effects caused by excess free Met.

Transcriptomic Insights into Mechanisms of Early Seed Maturation in the Garden Pea (Pisum sativum L.)

The garden pea (Pisum sativum L.) is a legume crop of immense economic value. Extensive breeding has led to the emergence of numerous pea varieties, of which some are distinguished by accelerated development in various stages of ontogenesis. One such trait is rapid seed maturation, which, despite novel insights into the genetic control of seed development in legumes, remains poorly studied. This article presents an attempt to dissect mechanisms of early maturation in the pea line Sprint-2 by means of whole transcriptome RNA sequencing in two developmental stages. By using a de novo assembly approach, we have obtained a reference transcriptome of 25,756 non-redundant entries expressed in pea seeds at either 10 or 20 days after pollination. Differential expression in Sprint-2 seeds has affected 13,056 transcripts. A comparison of the two pea lines with a common maturation rate demonstrates that while at 10 days after pollination, Sprint-2 seeds show development retardation linked to intensive photosynthesis, morphogenesis, and cell division, and those at 20 days show a rapid onset of desiccation marked by the cessation of translation and cell anabolism and accumulation of dehydration-protective and -storage moieties. Further inspection of certain transcript functional categories, including the chromatin constituent, transcription regulation, protein turnover, and hormonal regulation, has revealed transcriptomic trends unique to specific stages and cultivars. Among other remarkable features, Sprint-2 demonstrated an enhanced expression of transposable element-associated open reading frames and an altered expression of major maturation regulators and DNA methyltransferase genes. To the best of our knowledge, this is the first comparative transcriptomic study in which the issue of the seed maturation rate is addressed.

DNA methylation of LDOX gene contributes to the floral colour variegation in peach

Peach is an important fruit and ornamental plant around the globe. Variegation in flowers often captures consumers' attention, and variegated plants are of high ornamental value. To determine the relationship between DNA methylation and phenotype, we obtained the first single-nucleotide resolution DNA methylation of variegation cultivars in peach through bisulfite sequencing. In this study, a similar methylation rate of 12.90 % in variegated flower buds (VF) and 11.96 % in red flower buds (RF) was determined. The methyl-CG (mCG) was the main context in both samples. We identified 503 differentially methylated regions (DMRs) in all chromosomes. These DMRs were focused on 96 genes and 156 promoters. Associated with the transcriptional and proteome analysis, 106 differently expressed genes and 52 different proteins had varying degrees of methylation. Silent genes exhibited higher methylation levels than expressed genes. The methylation state of the leucoanthocyanidin dioxygenase (LDOX) promoter in VF was higher than RF at flower stages 2 (FS2) based on bisulfite sequencing PCR (BSP) results. Moreover, further experiments showed LDOX gene expression and enzyme activity in RF was higher than VF. The DNA methylation trend consisted of the gene expression and flower colour phenotype. Several cis-acting regulatory elements on BSP sequences were involved in phytohormones, transcription factors, and light responsiveness, which could affect gene expression. The higher level of LDOX gene expression promoted synthesis of colourful anthocyanidins, which caused red spots on the petal. Together, this study identified the context and level of methylation at each site with bisulfite sequencing (BS). These results are helpful in uncovering the mechanism of variegated flower petal formation in peach.

Whole-Genome DNA Methylation Associated With Differentially Expressed Genes Regulated Anthocyanin Biosynthesis Within Flower Color Chimera of Ornamental Tree Prunus mume

DNA methylation is one of the best-studied epigenetic modifications involved in many biological processes. However, little is known about the epigenetic mechanism for flower color chimera of Prunus mume (Japanese apricot, mei). Using bisulfate sequencing and RNA sequencing, we analyzed the white (FBW) and red (FBR) petals collected from an individual tree of Japanese apricot cv. ‘Fuban Tiaozhi’ mei to reveal the different changes in methylation patterns associated with gene expression leading to significant difference in anthocyanins accumulation of FBW (0.012 ± 0.005 mg/g) and FBR (0.078 ± 0.013 mg/g). It was found that gene expression levels were positively correlated with DNA methylation levels within gene-bodies of FBW and FBR genomes; however, negative correlations between gene expression and DNA methylation levels were detected within promoter domains. In general, the methylation level within methylome of FBW was higher; and in total, 4,618 differentially methylated regions (DMRs) and 1,212 differentially expressed genes (DEGs) were detected from FBW vs. FBR. We also identified 82 DMR-associated DEGs, and 13 of them, including PmBAHD, PmCYP450, and PmABC, were playing critical roles in phenylalanine metabolism pathway, glycosyltransferase activity, and ABC transporter. The evidence exhibited DNA methylation may regulate gene expression resulting in flower color chimera of Japanese apricot.

Establishment, maintenance, and biological roles of non-CG methylation in plants

Cytosine DNA methylation is prevalent throughout eukaryotes and prokaryotes. While most commonly thought of as being localized to dinucleotide CpG sites, non-CG sites can also be modified. Such non-CG methylation is widespread in plants, occurring at trinucleotide CHG and CHH (H = A, T, or C) sequence contexts. The prevalence of non-CG methylation in plants is due to the plant-specific CHROMOMETHYLASE (CMT) and RNA-directed DNA Methylation (RdDM) pathways. These pathways have evolved through multiple rounds of gene duplication and gene loss, generating epigenomic variation both within and between species. They regulate both transposable elements and genes, ensure genome integrity, and ultimately influence development and environmental responses. In these capacities, non-CG methylation influence and shape plant genomes.

DNA Methylation: Toward Crop Disease Resistance Improvement

Crop diseases, in conjunction with climate change, are a major threat to global crop production. DNA methylation is an epigenetic mark and is involved in plants' biological processes, including development, stress adaptation, and genome evolution. By providing a new source of variation, DNA methylation introduces novel direction to both scientists and breeders with its potential in disease resistance enhancement. Here, we discuss the impact of pathogen-induced DNA methylation modifications on a host's transcriptome reprogramming and genome stability, as part of the plant's defense mechanisms. We also highlight the knowledge gaps that need to be investigated for understanding the entire role of DNA methylation in plant pathogen interactions. This will ultimately assist breeders toward improving resistance and decreasing yield losses.

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.

MethBank 3.0: a database of DNA methylomes across a variety of species

MethBank (http://bigd.big.ac.cn/methbank) is a database that integrates high-quality DNA methylomes across a variety of species and provides an interactive browser for visualization of methylation data. Here, we present an updated implementation of MethBank (version 3.0) by incorporating more DNA methylomes from multiple species and equipping with more enhanced functionalities for data annotation and more friendly web interfaces for data presentation, search and visualization. MethBank 3.0 features large-scale integration of high-quality methylomes, involving 34 consensus reference methylomes derived from a large number of human samples, 336 single-base resolution methylomes from different developmental stages and/or tissues of five plants, and 18 single-base resolution methylomes from gametes and early embryos at multiple stages of two animals. Additionally, it is enhanced by improving the functionalities for data annotation, which accordingly enables systematic identification of methylation sites closely associated with age, sites with constant methylation levels across different ages, differentially methylated promoters, age-specific differentially methylated cytosines/regions, and methylated CpG islands. Moreover, MethBank provides tools to estimate human methylation age online and to identify differentially methylated promoters, respectively. Taken together, MethBank is upgraded with significant improvements and advances over the previous version, which is of great help for deciphering DNA methylation regulatory mechanisms for epigenetic studies.

Towards Exploitation of Adaptive Traits for Climate-Resilient Smart Pulses

Pulses are the main source of protein and minerals in the vegetarian diet. These are primarily cultivated on marginal lands with few inputs in several resource-poor countries of the world, including several in South Asia. Their cultivation in resource-scarce conditions exposes them to various abiotic and biotic stresses, leading to significant yield losses. Furthermore, climate change due to global warming has increased their vulnerability to emerging new insect pests and abiotic stresses that can become even more serious in the coming years. The changing climate scenario has made it more challenging to breed and develop climate-resilient smart pulses. Although pulses are climate smart, as they simultaneously adapt to and mitigate the effects of climate change, their narrow genetic diversity has always been a major constraint to their improvement for adaptability. However, existing genetic diversity still provides opportunities to exploit novel attributes for developing climate-resilient cultivars. The mining and exploitation of adaptive traits imparting tolerance/resistance to climate-smart pulses can be accelerated further by using cutting-edge approaches of biotechnology such as transgenics, genome editing, and epigenetics. This review discusses various classical and molecular approaches and strategies to exploit adaptive traits for breeding climate-smart pulses.

De Novo Domestication: An Alternative Route toward New Crops for the Future

Current global agricultural production must feed over 7 billion people. However, productivity varies greatly across the globe and is under threat from both increased competitions for land and climate change and associated environmental deterioration. Moreover, the increase in human population size and dietary changes are putting an ever greater burden on agriculture. The majority of this burden is met by the cultivation of a very small number of species, largely in locations that differ from their origin of domestication. Recent technological advances have raised the possibility of de novo domestication of wild plants as a viable solution for designing ideal crops while maintaining food security and a more sustainable low-input agriculture. Here we discuss how the discovery of multiple key domestication genes alongside the development of technologies for accurate manipulation of several target genes simultaneously renders de novo domestication a route toward crops for the future.

Genome-wide nucleotide patterns and potential mechanisms of genome divergence following domestication in maize and soybean

BACKGROUND

Plant domestication provides a unique model to study genome evolution. Many studies have been conducted to examine genes, genetic diversity, genome structure, and epigenome changes associated with domestication. Interestingly, domesticated accessions have significantly higher [A] and [T] values across genome-wide polymorphic sites than accessions sampled from the corresponding progenitor species. However, the relative contributions of different genomic regions to this genome divergence pattern and underlying mechanisms have not been well characterized.

RESULTS

Here, we investigate the genome-wide base-composition patterns by analyzing millions of SNPs segregating among 100 accessions from a teosinte-maize comparison set and among 302 accessions from a wild-domesticated soybean comparison set. We show that non-genic part of the genome has a greater contribution than genic SNPs to the [AT]-increase observed between wild and domesticated accessions in maize and soybean. The separation between wild and domesticated accessions in [AT] values is significantly enlarged in non-genic and pericentromeric regions. Motif frequency and sequence context analyses show the motifs (PyCG) related to solar-UV signature are enriched in these regions, particularly when they are methylated. Additional analysis using population-private SNPs also implicates the role of these motifs in relatively recent mutations. With base-composition across polymorphic sites as a genome phenotype, genome scans identify a set of putative candidate genes involved in UV damage repair pathways.

CONCLUSIONS

The [AT]-increase is more pronounced in genomic regions that are non-genic, pericentromeric, transposable elements; methylated; and with low recombination. Our findings establish important links among UV radiation, mutation, DNA repair, methylation, and genome evolution.

DNA methylome analysis provides evidence that the expansion of the tea genome is linked to TE bursts

DNA methylation is essential for gene regulation, imprinting and silencing of transposable elements (TEs). Although bursts of transposable elements are common in many plant lineages, how plant DNA methylation is related to transposon bursts remains unclear. Here we explore the landscape of DNA methylation of tea, a species thought to have experienced a recent transposon burst event. This species possesses more transposable elements than any other sequenced asterids (potato, tomato, coffee, pepper and tobacco). The overall average DNA methylation levels were found to differ among the tea, potato and tomato genomes, and methylation at CHG sequence sites was found to be significantly higher in tea than that in potato or tomato. Moreover, the abundant TEs resulting from burst events not only resulted in tea developing a very large genome size, but also affected many genes involved in importantly biological processes, including caffeine, theanine and flavonoid metabolic pathway genes. In addition, recently transposed TEs were more heavily methylated than ancient ones, implying that DNA methylation is proportionate to the degree of TE silencing, especially on recent active ones. Taken together, our results show that DNA methylation regulates transposon silencing and may play a role in genome size expansion.

Effects of ploidy variation on promoter DNA methylation and gene expression in rice (Oryza sativa L.)

BACKGROUND

Polyploidy, or whole-genome duplication (WGD) promotes genetic diversification in plants. However, whether WGD is accompanied by epigenetic regulation especially DNA methylation remains yet elusive. Methylation of different region in genomic DNA play discrete role in gene regulation and developmental processes in plants.

RESULTS

In our study, we used an apomictic rice line (SARII-628) that produces twin seedlings of different ploidy for methylated DNA immunoprecipitation sequencing (MeDIP-seq). We compared the level of methylation and mRNA expression in three different (CG, CHG, and CHH) sequence contexts of promoter region among haploid (1X), diploid (2X), and triploid (3X) seedling. We used MeDIP-Seq analysis of 14 genes to investigate whole genome DNA methylation and found that relative level of DNA methylation across different ploidy was in following order e.g. diploid > triploid > haploid. GO functional classification of differentially methylated genes into 9 comparisons group of promoter, intergenic and intragenic region discovered, these genes were mostly enriched for cellular component, molecular function, and biological process. By the comparison of methylome data, digital gene expression (DGE), mRNA expression profile, and Q-PCR findings LOC_ Os07g31450 and LOC_ Os01g59320 were analyzed for BS-Seq (Bisulphite sequencing).

CONCLUSIONS

We found that (1) The level of the promoter DNA methylation is negatively correlated with gene expression within each ploidy level. (2) Among all ploidy levels, CG sequence context had highest methylation frequency, and demonstrated that the high CG methylation did reduce gene expression change suggesting that DNA methylation exert repressive function and ensure genome stability during WGD. (3) Alteration in ploidy (from diploid to haploid, or diploid to triploid) reveals supreme changes in methylation frequency of CHH sequence context. Our finding will contribute an understanding towards lower stability of CHH sequence context and educate the effect of promoter region methylation during change in ploidy state in rice.

LncRNAs in polyploid cotton interspecific hybrids are derived from transposon neofunctionalization

BACKGROUND

Interspecific hybridization and whole genome duplication are driving forces of genomic and organism diversification. But the effect of interspecific hybridization and whole genome duplication on the non-coding portion of the genome in particular remains largely unknown. In this study, we examine the profile of long non-coding RNAs (lncRNAs), comparing them with that of coding genes in allotetraploid cotton (Gossypium hirsutum), its putative diploid ancestors (G. arboreum; G. raimondii), and an F1 hybrid (G. arboreum × G. raimondii, AD).

RESULTS

We find that most lncRNAs (80%) that were allelic expressed in the allotetraploid genome. Moreover, the genome shock of hybridization reprograms the non-coding transcriptome in the F1 hybrid. Interestingly, the activated lncRNAs are predominantly transcribed from demethylated TE regions, especially from long interspersed nuclear elements (LINEs). The DNA methylation dynamics in the interspecies hybridization are predominantly associated with the drastic expression variation of lncRNAs. Similar trends of lncRNA bursting are also observed in the progress of polyploidization. Additionally, we find that a representative novel lncRNA XLOC_409583 activated after polyploidization from a LINE in the A subgenome of allotetraploid cotton was involved in control of cotton seedling height.

CONCLUSION

Our results reveal that the processes of hybridization and polyploidization enable the neofunctionalization of lncRNA transcripts, acting as important sources of increased plasticity for plants.

Genetic and epigenetic divergence of duplicate genes in two legume species

Soybean (Glycine max) and common bean (Phaseolus vulgaris) share a polyploidy event ~59 MYA, followed by a Glycine-specific whole genome duplication (WGD) ~8-13 MYA. Duplicated genes were classified into five categories: singletons, dispersed, proximal, tandem, or WGD/segmental and found strong correlations between gene category and functional annotation. Photosynthesis and transcriptional regulation-related Gene Ontology terms were significantly over-represented in singletons and WGD genes, respectively, aligning with the gene balance hypothesis. We found that the divergence of gene expression and DNA methylation between WGD-derived paralogs increased with age and that WGD genes, initially retained via dosage constraints, subsequently underwent expression divergence, associated with other factors such as DNA methylation. Genes derived from different modes of duplication differed in breadth, level, and specificity of expression in both species. Orthologous genes and ungrouped genes (genes not in an ortholog group) differed in expression patterns. The protein divergence rates of WGD paralog pairs containing an ungrouped gene were higher than those for which both copies had orthologs. We propose that many ungrouped genes are derived from divergent and redundant gene copies, concordant with the neofunctionalization hypothesis. Tandemly duplicated genes were distinct from WGD-derived genes, indicating that mode of duplication contributes to the evolutionary fate of duplicated genes.

Substantial Epigenetic Variation Causing Flower Color Chimerism in the Ornamental Tree Prunus mume Revealed by Single Base Resolution Methylome Detection and Transcriptome Sequencing

Epigenetic changes caused by methylcytosine modification participate in gene regulation and transposable element (TE) repression, resulting in phenotypic variation. Although the effects of DNA methylation and TE repression on flower, fruit, seed coat, and leaf pigmentation have been investigated, little is known about the relationship between methylation and flower color chimerism. In this study, we used a comparative methylomic⁻transcriptomic approach to explore the molecular mechanism responsible for chimeric flowers in Prunus mume "Danban Tiaozhi". High-performance liquid chromatography-electrospray ionization mass spectrometry revealed that the variation in white (WT) and red (RT) petal tissues in this species is directly due to the accumulation of anthocyanins, i.e., cyanidin 3,5-O-diglucoside, cyanidin 3-O-glucoside, and peonidin 3-O-glucoside. We next mapped the first-ever generated methylomes of P. mume, and found that 11.29⁻14.83% of the genomic cytosine sites were methylated. We also determined that gene expression was negatively correlated with methylcytosine level in general, and uncovered significant epigenetic variation between WT and RT. Furthermore, we detected differentially methylated regions (DMRs) and DMR-related genes between WT and RT, and concluded that many of these genes, including differentially expressed genes (DEGs) and transcription factor genes, are critical participants in the anthocyanin regulatory pathway. Importantly, some of the associated DEGs harbored TE insertions that were also modified by methylcytosine. The above evidence suggest that flower color chimerism in P. mume is induced by the DNA methylation of critical genes and TEs.

Genome-wide bisulphite-sequencing reveals organ-specific methylation patterns in chickpea

DNA methylation is widely known to regulate gene expression in eukaryotes. Here, we unraveled DNA methylation patterns in cultivated chickpea to understand the regulation of gene expression in different organs. We analyzed the methylation pattern in leaf tissue of wild chickpea too, and compared it with cultivated chickpea. Our analysis indicated abundant CG methylation within gene-body and CHH methylation in intergenic regions of the chickpea genome in all the organs examined. Analysis of differentially methylated regions (DMRs) demonstrated a higher number of CG context DMRs in wild chickpea and CHH context DMRs in cultivated chickpea. We observed increased preponderance of hypermethylated DMRs in the promoter regions and hypomethylated DMRs in the genic regions in cultivated chickpea. Genomic location and context of the DMRs correlated well with expression of proximal genes. Our results put forth a positive correlation of promoter hypermethylation with increased transcript abundance via identification of DMR-associated genes involved in flower development in cultivated chickpea. The atypical correlation observed between promoter hypermethylation and increased transcript abundance might be dependent on 24-nt small RNAs and transcription factors binding to the promoter region. This study provides novel insights into DNA methylation patterns in chickpea and their role in regulation of gene expression.

Contribution of transposable elements in the plant's genome

Plants maintain extensive growth flexibility under different environmental conditions, allowing them to continuously and rapidly adapt to alterations in their environment. A large portion of many plant genomes consists of transposable elements (TEs) that create new genetic variations within plant species. Different types of mutations may be created by TEs in plants. Many TEs can avoid the host's defense mechanisms and survive alterations in transposition activity, internal sequence and target site. Thus, plant genomes are expected to utilize a variety of mechanisms to tolerate TEs that are near or within genes. TEs affect the expression of not only nearby genes but also unlinked inserted genes. TEs can create new promoters, leading to novel expression patterns or alternative coding regions to generate alternate transcripts in plant species. TEs can also provide novel cis-acting regulatory elements that act as enhancers or inserts within original enhancers that are required for transcription. Thus, the regulation of plant gene expression is strongly managed by the insertion of TEs into nearby genes. TEs can also lead to chromatin modifications and thereby affect gene expression in plants. TEs are able to generate new genes and modify existing gene structures by duplicating, mobilizing and recombining gene fragments. They can also facilitate cellular functions by sharing their transposase-coding regions. Hence, TE insertions can not only act as simple mutagens but can also alter the elementary functions of the plant genome. Here, we review recent discoveries concerning the contribution of TEs to gene expression in plant genomes and discuss the different mechanisms by which TEs can affect plant gene expression and reduce host defense mechanisms.

Epigenetic perspectives on the evolution and domestication of polyploid plant and crops

Polyploidy or whole genome duplication (WGD) is a prominent feature for genome evolution of some animals and all flowering plants, including many important crops such as wheat, cotton, and canola. In autopolyploids, genome duplication often perturbs dosage regulation on biological networks. In allopolyploids, interspecific hybridization could induce genetic and epigenetic changes, the effects of which could be amplified by genome doubling (ploidy changes). Albeit the importance of genetic changes, some epigenetic changes can be stabilized and transmitted as epialleles into the progeny, which are subject to natural selection, adaptation, and domestication. Here we review recent advances for general and specific roles of epigenetic changes in the evolution of flowering plants and domestication of agricultural crops.

Genomic and epigenomic immunity in common bean: the unusual features of NB-LRR gene family

In plants, a key class of genes comprising most of disease resistance (R) genes encodes Nucleotide-binding leucine-rich repeat (NL) proteins. Access to common bean (Phaseolus vulgaris) genome sequence provides unparalleled insight into the organization and evolution of this large gene family (∼400 NL) in this important crop. As observed in other plant species, most common bean NL are organized in cluster of genes. However, a particularity of common bean is that these clusters are often located in subtelomeric regions close to terminal knobs containing the satellite DNA khipu. Phylogenetically related NL are spread between different chromosome ends, suggesting frequent exchanges between non-homologous chromosomes. NL peculiar location, in proximity to heterochromatic regions, led us to study their DNA methylation status using a whole-genome cytosine methylation map. In common bean, NL genes displayed an unusual body methylation pattern since half of them are methylated in the three contexts, reminiscent of the DNA methylation pattern of repeated sequences. Moreover, 90 NL were also abundantly targeted by 24 nt siRNA, with 90% corresponding to methylated NL genes. This suggests the existence of a transcriptional gene silencing mechanism of NL through the RdDM (RNA-directed DNA methylation) pathway in common bean that has not been described in other plant species.

Genic C-Methylation in Soybean Is Associated with Gene Paralogs Relocated to Transposable Element-Rich Pericentromeres

Most plants are polyploid due to whole-genome duplications (WGD) and can thus have duplicated genes. Following a WGD, paralogs are often fractionated (lost) and few duplicate pairs remain. Little attention has been paid to the role of DNA methylation in the functional divergence of paralogous genes. Using high-resolution methylation maps of accessions of domesticated and wild soybean, we show that in soybean, a recent paleopolyploid with many paralogs, DNA methylation likely contributed to the elimination of genetic redundancy of polyploidy-derived gene paralogs. Transcriptionally silenced paralogs exhibit particular genomic features as they are often associated with proximal transposable elements (TEs) and are preferentially located in pericentromeres, likely due to gene movement during evolution. Additionally, we provide evidence that gene methylation associated with proximal TEs is implicated in the divergence of expression profiles between orthologous genes of wild and domesticated soybean, and within populations.

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.

DNA methylation profiles in cancer diagnosis and therapeutics

Cancer initiation and proliferation is regulated by both epigenetic and genetic events with epigenetic modifications being increasingly identified as important targets for cancer research. DNA methylation catalyzed by DNA methyltransferases (DNMTs) is one of the essential epigenetic mechanisms that control cell proliferation, apoptosis, differentiation, cell cycle, and transformation in eukaryotes. Recent progress in epigenetics revealed a deeper understanding of the mechanisms of tumorigenesis and provided biomarkers for early detection, diagnosis, and prognosis in cancer patients. Although DNA methylation biomarker possesses potential contributing to precision medicine, there are still limitations to be overcome before it reaches clinical setting. Hence, the current status of DNA methylation biomarkers was reviewed and the future use in clinic was also predicted.

Interaction and Regulation of Carbon, Nitrogen, and Phosphorus Metabolisms in Root Nodules of Legumes

Members of the plant family Leguminosae (Fabaceae) are unique in that they have evolved a symbiotic relationship with rhizobia (a group of soil bacteria that can fix atmospheric nitrogen). Rhizobia infect and form root nodules on their specific host plants before differentiating into bacteroids, the symbiotic form of rhizobia. This complex relationship involves the supply of C₄-dicarboxylate and phosphate by the host plants to the microsymbionts that utilize them in the energy-intensive process of fixing atmospheric nitrogen into ammonium, which is in turn made available to the host plants as a source of nitrogen, a macronutrient for growth. Although nitrogen-fixing bacteroids are no longer growing, they are metabolically active. The symbiotic process is complex and tightly regulated by both the host plants and the bacteroids. The metabolic pathways of carbon, nitrogen, and phosphate are heavily regulated in the host plants, as they need to strike a fine balance between satisfying their own needs as well as those of the microsymbionts. A network of transporters for the various metabolites are responsible for the trafficking of these essential molecules between the two partners through the symbiosome membrane (plant-derived membrane surrounding the bacteroid), and these are in turn regulated by various transcription factors that control their expressions under different environmental conditions. Understanding this complex process of symbiotic nitrogen fixation is vital in promoting sustainable agriculture and enhancing soil fertility.

Patterns and Consequences of Subgenome Differentiation Provide Insights into the Nature of Paleopolyploidy in Plants

Polyploidy is an important feature of plant genomes, but the nature of many polyploidization events remains to be elucidated. Here, we demonstrate that the evolutionary fates of the subgenomes in maize (Zea mays) and soybean (Glycine max) have followed different trajectories. One subgenome has been subject to relaxed selection, lower levels of gene expression, higher rates of transposable element accumulation, more small interfering RNAs and DNA methylation around genes, and higher rates of gene loss in maize, whereas none of these features were observed in soybean. Nevertheless, individual gene pairs exhibit differentiation with respect to these features in both species. In addition, we observed a higher number of chromosomal rearrangements and higher frequency of retention of duplicated genes in soybean than in maize. Furthermore, soybean "singletons" were found to be more frequently tandemly duplicated than "duplicates" in soybean, which may, to some extent, counteract the genome imbalance caused by gene loss. We propose that unlike in maize, in which two subgenomes were distinct prior to the allotetraploidization event and thus experienced global differences in selective constraints, in soybean, the two subgenomes were far less distinct prior to polyploidization, such that individual gene pairs, rather than subgenomes, experienced stochastic differences over longer periods of time, resulting in retention of the majority of duplicates.