Epigenetic Variation Research Articles

The evolutionary consequences of interactions between the epigenome, the genome and the environment.

The epigenome is the suite of interacting chemical marks and molecules that helps to shape patterns of development, phenotypic plasticity and gene regulation, in part due to its responsiveness to environmental stimuli. There is increasing interest in understanding the functional and evolutionary importance of this sensitivity under ecologically realistic conditions. Observations that epigenetic variation abounds in natural populations have prompted speculation that it may facilitate evolutionary responses to rapid environmental perturbations, such as those occurring under climate change. A frequent point of contention is whether epigenetic variants reflect genetic variation or are independent of it. The genome and epigenome often appear tightly linked and interdependent. While many epigenetic changes are genetically determined, the converse is also true, with DNA sequence changes influenced by the presence of epigenetic marks. Understanding how the epigenome, genome and environment interact with one another is therefore an essential step in explaining the broader evolutionary consequences of epigenomic variation. Drawing on results from experimental and comparative studies carried out in diverse plant and animal species, we synthesize our current understanding of how these factors interact to shape phenotypic variation in natural populations, with a focus on identifying similarities and differences between taxonomic groups. We describe the main components of the epigenome and how they vary within and between taxa. We review how variation in the epigenome interacts with genetic features and environmental determinants, with a focus on the role of transposable elements (TEs) in integrating the epigenome, genome and environment. And we look at recent studies investigating the functional and evolutionary consequences of these interactions. Although epigenetic differentiation in nature is likely often a result of drift or selection on stochastic epimutations, there is growing evidence that a significant fraction of it can be stably inherited and could therefore contribute to evolution independently of genetic change.

Charting epimutation dynamics in human hematopoietic differentiation.

DNA methylation plays a critical role in hematopoietic differentiation. Epimutation is a stochastic variation in DNA methylation that induces epigenetic heterogeneity. However, the effects of epimutations on normal hematopoiesis and hematopoietic diseases remain unclear. In this study, we developed a Julia package called EpiMut that enabled rapid and accurate quantification of epimutations. EpiMut was used to evaluate and provide an epimutation landscape in steady-state hematopoietic differentiation involving 13 types of blood cells ranging from hematopoietic stem/progenitor cells to mature cells. We showed that substantial genomic regions exhibited epigenetic variations rather than significant differences in DNA methylation levels between the myeloid and lymphoid lineages. Stepwise dynamics of epimutations were observed during the differentiation of each lineage. Importantly, we found that epimutation significantly enriched signals associated with lineage differentiation. Furthermore, epimutations in hematopoietic stem cells (HSCs) derived from various sources and acute myeloid leukemia were related to the function of HSCs and malignant cell disorders. Taken together, our study comprehensively documented an epimutation map and uncovered its important roles in human hematopoiesis, thereby offering insights into hematopoietic regulation.

Epigenetic diversity of genes with copy number variations among natural populations of the three-spined stickleback.

Duplicated genes provide the opportunity for evolutionary novelty and adaptive divergence. In many cases, having more gene copies increases gene expression, which might facilitate adaptation to stressful or novel environments. Conversely, overexpression or misexpression of duplicated genes can be detrimental and subject to negative selection. In this scenario, newly duplicate genes may evade purifying selection if they are epigenetically silenced, at least temporarily, leading them to persist in populations as copy number variations (CNVs). In animals and plants, younger gene duplicates tend to have higher levels of DNA methylation and lower levels of gene expression, suggesting epigenetic regulation could promote the retention of gene duplications via expression repression or silencing. Here, we test the hypothesis that DNA methylation variation coincides with young duplicate genes that are segregating as CNVs in six populations of the three-spined stickleback that span a salinity gradient from 4 to 30 PSU. Using reduced-representation bisulfite sequencing, we found DNA methylation and CNV differentiation outliers rarely overlapped. Whereas lineage-specific genes and young duplicates were found to be highly methylated, just two gene CNVs showed a significant association between promoter methylation level and copy number, suggesting that DNA methylation might not interact with CNVs in our dataset. If most new duplications are regulated for dosage by epigenetic mechanisms, our results do not support a strong contribution from DNA methylation soon after duplication. Instead, our results are consistent with a preference to duplicate genes that are already highly methylated.

H3K9 methylation regulates heterochromatin silencing through incoherent feedforward loops.

Histone H3 lysine-9 methylation (H3K9me) is a hallmark of the condensed and transcriptionally silent heterochromatin. It remains unclear how H3K9me controls transcription silencing and how cells delimit H3K9me domains to avoid silencing essential genes. Here, using Arabidopsis genetic systems that induce H3K9me2 in genes and transposons de novo, we show that H3K9me2 accumulation paradoxically also causes the deposition of the euchromatic mark H3K36me3 by a SET domain methyltransferase, ASHH3. ASHH3-induced H3K36me3 confers anti-silencing by preventing the demethylation of H3K4me1 by LDL2, which mediates transcriptional silencing downstream of H3K9me2. These results demonstrate that H3K9me2 not only facilitates but orchestrates silencing by actuating antagonistic silencing and anti-silencing pathways, providing insights into the molecular basis underlying proper partitioning of chromatin domains and the creation of metastable epigenetic variation.

Comprehensive Analysis of the Function and Prognostic Value of TAS2Rs Family-Related Genes in Colon Cancer.

In the realm of colon carcinoma, significant genetic and epigenetic diversity is observed, underscoring the necessity for tailored prognostic features that can guide personalized therapeutic strategies. In this study, we explored the association between the type 2 bitter taste receptor (TAS2Rs) family-related genes and colon cancer using RNA-sequencing and clinical datasets from The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO). Our preliminary analysis identified seven TAS2Rs genes associated with survival using univariate Cox regression analysis, all of which were observed to be overexpressed in colon cancer. Subsequently, based on these seven TAS2Rs prognostic genes, two colon cancer molecular subtypes (Cluster A and Cluster B) were defined. These subtypes exhibited distinct prognostic and immune characteristics, with Cluster A characterized by low immune cell infiltration and less favorable outcomes, while Cluster B was associated with high immune cell infiltration and better prognosis. Finally, we developed a robust scoring system using a gradient boosting machine (GBM) approach, integrated with the gene-pairing method, to predict the prognosis of colon cancer patients. This machine learning model could improve our predictive accuracy for colon cancer outcomes, underscoring its value in the precision oncology framework.

Exploring the Significance of Somaclonal Variations in Horticultural Crops.

Genetic and epigenetic variations produced via cell and tissue culture open up new sources of variability intra-species which can be used to improve crops. The use of in vitro generated somaclonal variations for selecting novel variants aids in the development of novel genotypes having desirable agronomic traits that can be released as varieties or utilized for breeding purposes. Horticultural crops give higher yield and productivity per unit area than other crops, as well as provide good economic returns which have led to an increase in their potential benefits throughout time. The last three to four decades have seen the selection and release of a number of valuable somaclonal variants, many of which possess remarkable features including disease resistance, high yield, improved nutritional quality and abiotic stress tolerance. Generating somaclonal variations has given breeders a novel alternative option for obtaining genetic diversity in horticultural crops and without advanced technologies. The variations introduced through tissue culture process, methods to determine and validate genetic changes in vitro regenerated plantlets, along with prospective application of such variations in horticultural crops' improvement are reviewed in the present work.

ATRT-16. MAPPING THE DISORDERED METHYLOME IN ATRT

Abstract BACKGROUND Atypical Teratoid Rhabdoid Tumor (ATRT) is an aggressive pediatric brain tumor driven by inactivation of the chromatin regulator SMARCB1, with remarkably few other mutations. SMARCB1 loss leads to epigenetic dysregulation which plays a pivotal role in ATRT development and progression, and DNA methylation profiling has distinguished molecular subgroups of this disease. However, studies that provide a static view of the epigenome may not capture epigenetic variability that drives cellular plasticity and heterogeneity. Mapping epigenetic variability is crucial to understanding the phenotypic plasticity of this tumor and its ability to evade therapeutic interventions in the clinic. METHODS We performed whole genome bisulfite sequencing (WGBS) on 22 ATRT samples, 4 normal fetal brain controls, and 1 hindbrain neural stem cell control, and applied a novel analytic pipeline (InformME) to generate DNA methylation potential energy landscapes encapsulating epigenetic stochasticity. This characterizes mean methylation as well as methylation entropy. We identified key genes and regulatory regions with the most stochastic methylation signatures. RESULTS We find that ATRT exhibits a globally destabilized epigenome, with markedly elevated methylation entropy (a measure of variability/disorder). Methylation discordance mapped to key regulatory regions including bivalent promoters, quiescent chromatin regions, and enhancers, and specific genes, with an enrichment for targets of the Polycomb repressive complex 2 and regulators of pluripotency and differentiation. CONCLUSIONS The disordered methylome of ATRT could play a critical role in the tumor’s ability to evade therapies. Mapping the targets of epigenetic instability can identify targets for intervention. Moreover, strategies to constrain epigenetic variability may have therapeutic benefit by limiting the number of epigenetic states available to cancer cells.

A strategy for studying epigenetic diversity in natural populations: proof of concept in poplar and oak.

In the last 20 years, several techniques have been developed for quantifying DNA methylation, the most studied epigenetic marks in eukaryotes, including the gold standard method, whole-genome bisulfite sequencing (WGBS). WGBS quantifies genome-wide DNA methylation but has several inconveniences rendering it less suitable for population-scale epigenetic studies. The high cost of deep sequencing and the large amounts of data generated prompted us to seek an alternative approach. Restricting studies to parts of the genome would be a satisfactory alternative had there not been a major limitation: the need to select upstream targets corresponding to differentially methylated regions as targets. Given the need to study large numbers of samples, we propose a strategy for investigating DNA methylation variation in natural populations, taking into account the structural complexity of genomes, their size, and their content in unique coding regions versus repeated regions as transposable elements. We first identified regions of highly variable DNA methylation in a subset of genotypes representative of the biological diversity in the population by WGBS. We then analysed the variations of DNA methylation in these targeted regions at the population level by sequencing capture bisulfite (SeqCapBis). The entire strategy was then validated by applying it to another species. Our strategy was developed as a proof of concept on natural populations of two forest species: Populus nigra and Quercus petraea.

DIPG-58. TARGETING DISORDERED DNA METHYLATION IN DIPG TO CONSTRAIN VARIABILITY AND INDUCE IMMUNE SIGNALING

Abstract BACKGROUND Diffuse Intrinsic Pontine Glioma (DIPG) is a universally fatal pediatric brain cancer characterized by the histone H3 K27M mutation. The known downstream consequences of this mutation include loss of repressive chromatin marks, global DNA hypomethylation, and altered gene expression. METHODS We sought to investigate the role of epigenetic variability as a basis of DIPG cellular heterogeneity and plasticity by using Whole Genome Bisulfite Sequencing (WGBS). Existing methods of analysis of WGBS are limited in that they average/smooth methylation data and fail to capture intrinsic variability. We performed (WGBS) on 23 primary patient samples of DIPG and applied an analysis that captures the variability of DNA methylation reads to characterize methylation stochasticity. We then sought to modulate the epigenome pharmacologically using a DNA methyltransferase inhibitor and assessed changes through RNA-Seq and immunopeptidomics. RESULTS We find that DIPG has a marked increase in methylation stochasticity (quantified as methylation entropy). Additionally, while DIPG is globally hypomethylated, it retains DNA methylation at specific regulatory regions, such as bivalent promoters, and key genes, such as FOXG1, CXCR4, and KLF4. We sought to reverse these changes using the DNA methyltransferase inhibitor, decitabine, in 4 DIPG neurosphere cell lines. Treatment with decitabine reversed the entropy seen in DIPG and induced innate immune signaling and interferon signaling in DIPG cells. Using RNA-seq, we found that decitabine treatment dramatically alters gene expression, including activation of endogenous retroviral elements, activation of interferon signaling, de-repression of hypermethylated targets, and expression of putative neoantigens, such as PRAME and DAZL. Immunopeptidomics profiling of MHC Class I bound canonical and non-canonical peptides of DIPG cell lines after decitabine treatment also revealed decitabine-induced differences. CONCLUSIONS Taken together, our study finds that the methylome of DIPG is highly disordered but is also highly responsive to pharmacologic modulation by hypomethylating agents. Furthermore, hypomethylating drugs can increase the immunogenicity of DIPG, thus offering the potential for future combination with immunotherapy.

DNA methylation variations underlie lettuce domestication and divergence

BackgroundLettuce (Lactuca sativa L.) is an economically important vegetable crop worldwide. Lettuce is believed to be domesticated from a single wild ancestor Lactuca serriola and subsequently diverged into two major morphologically distinct vegetable types: leafy lettuce and stem lettuce. However, the role of epigenetic variation in lettuce domestication and divergence remains largely unknown.ResultsTo understand the genetic and epigenetic basis underlying lettuce domestication and divergence, we generate single-base resolution DNA methylomes from 52 Lactuca accessions, including major lettuce cultivars and wild relatives. We find a significant increase of DNA methylation during lettuce domestication and uncover abundant epigenetic variations associated with lettuce domestication and divergence. Interestingly, DNA methylation variations specifically associated with leafy and stem lettuce are related to regulation and metabolic processes, respectively, while those associated with both types are enriched in stress responses. Moreover, we reveal that domestication-induced DNA methylation changes could influence expression levels of nearby and distal genes possibly through affecting chromatin accessibility and chromatin loop.ConclusionOur study provides population epigenomic insights into crop domestication and divergence and valuable resources for further domestication for diversity and epigenetic breeding to boost crop improvement.

Engineered lignocellulosic based biochar to remove endocrine-disrupting chemicals: Assessment of binding mechanism

The safety and health of aquatic organisms and humans are threatened by the increasing presence of pollutants in the environment. Endocrine disrupting chemicals are common pollutants which affect the function of endocrine and causes adverse effects on human health. These chemicals can disrupt metabolic processes by interacting with hormone receptors upon consumptions by humans or aquatic species. Several studies have reported the presence of endocrine disrupting chemicals in waterbodies, food, air and soil. These chemicals are associated with increasing occurrence of obesity, metabolic disorders, reproductive abnormalities, autism, cancer, epigenetic variation and cardiovascular risk. Conventional treatment processes are expensive, not environment friendly and unable to achieve complete removal of these harmful chemicals. In recent years, biochar from different sources has gained a considerable interest due to their adsorption efficiency with porous structure and large surface areas. biochar derived from lignocellulosic biomass are widely used as sustainable catalysts in soil remediation, carbon sequestration, removal of organic and inorganic pollutants and wastewater treatment. This review conceptualizes the production techniques of biochar from lignocellulosic biomass and explores the functionalization and interaction of biochar with endocrine-disrupting chemicals. This review also identifies the further needs of research. Overall, the environmental and health risks of endocrine-disrupting chemicals can be dealt with by biochar produced from lignocellulosic biomass as a sustainable and prominent approach.

Epigenetic differences between wild and cultivated grapevines highlight the contribution of DNA methylation during crop domestication

The domestication process in grapevines has facilitated the fixation of desired traits. Nowadays, vegetative propagation through cuttings enables easier preservation of these genotypes compared to sexual reproduction. Nonetheless, even with vegetative propagation, various phenotypes are often present within the same vineyard due to the accumulation of somatic mutations. These mutations are not the sole factors influencing phenotype. Alongside somatic variations, epigenetic variation has been proposed as a pivotal player in regulating phenotypic variability acquired during domestication. The emergence of these epialleles might have significantly influenced grapevine domestication over time. This study aims to investigate the impact of domestication on methylation patterns in cultivated grapevines. Reduced-representation bisulfite sequencing was conducted on 18 cultivated and wild accessions. Results revealed that cultivated grapevines exhibited higher methylation levels than their wild counterparts. Differential Methylation Analysis between wild and cultivated grapevines identified a total of 9955 differentially methylated cytosines, of which 78% were hypermethylated in cultivated grapevines. Functional analysis shows that core methylated genes (consistently methylated in both wild and cultivated accessions) are associated with stress response and terpenoid/isoprenoid metabolic processes. Meanwhile, genes with differential methylation are linked to protein targeting to the peroxisome, ethylene regulation, histone modifications, and defense response. Collectively, our results highlight the significant roles that epialleles may have played throughout the domestication history of grapevines.

Microbiome and epigenetic variation in wild fish with low genetic diversity

Non-genetic sources of phenotypic variation, such as the epigenome and the microbiome, could be important contributors to adaptive variation for species with low genetic diversity. However, little is known about the complex interaction between these factors and the genetic diversity of the host, particularly in wild populations. Here, we examine the skin microbiome composition of two closely-related mangrove killifish species with different mating systems (self-fertilising and outcrossing) under sympatric and allopatric conditions. This allows us to partition the influence of the genotype and the environment on their microbiome and (previously described) epigenetic profiles. We find the diversity and community composition of the skin microbiome are strongly shaped by the environment and, to a lesser extent, by species-specific influences. Heterozygosity and microbiome alpha diversity, but not epigenetic variation, are associated with the fluctuating asymmetry of traits related to performance (vision) and behaviour (aggression). Our study identifies that a proportion of the epigenetic diversity and microbiome differentiation is unrelated to genetic variation, and we find evidence for an associative relationship between microbiome and epigenetic diversity in these wild populations. This suggests that both mechanisms could potentially contribute to variation in species with low genetic diversity.

APOA1/C3/A4/A5 Gene Cluster at 11q23.3 and Lipid Metabolism Disorders: From Epigenetic Mechanisms to Clinical Practices.

The APOA1/C3/A4/A5 cluster is an essential component in regulating lipoprotein metabolism and maintaining plasma lipid homeostasis. A genome-wide association analysis and Mendelian randomization have revealed potential associations between genetic variants within this cluster and lipid metabolism disorders, including hyperlipidemia and cardiovascular events. An enhanced understanding of the complexity of gene regulation has led to growing recognition regarding the role of epigenetic variation in modulating APOA1/C3/A4/A5 gene expression. Intensive research into the epigenetic regulatory patterns of the APOA1/C3/A4/A5 cluster will help increase our understanding of the pathogenesis of lipid metabolism disorders and facilitate the development of new therapeutic approaches. This review discusses the biology of how the APOA1/C3/A4/A5 cluster affects circulating lipoproteins and the current progress in the epigenetic regulation of the APOA1/C3/A4/A5 cluster.

Recent advances in epigenomic techniques: Analysis of DNA base modifications

Epigenetics refers to the study of modifications in DNA bases, histone proteins and/or alterations in non-coding-RNA biogenesis that cause changes in gene expression without affecting the underlying nucleotide sequence. Epigenetic variation is emerging as one of the regulators of developmental processes and stress responses in animals and plants. Identification and interrogation of epigenetic changes provide unprecedented opportunities for resolving some of the enigmas that could not be solved based on genetic principles. Moreover, epigenome editing might become a preferred technique for the manipulation of gene expression in the desired organism. Although epigenetics is still in its infancy, the need of the day is to identify epigenetic marks and associate them with heritable phenotypes. Several techniques have been developed to analyze the epigenetic changes and associate them with molecular/physiological processes. It is necessary to examine them in a site- and cell-specific manner. Though considerable progress has been made in the last decade towards mapping DNA modifications as well as other epigenetic changes, there is a need to improve the specificity, sensitivity, and resolving power of the techniques to conduct epigenetic research in a cell-specific manner. The recent technological advancements are currently driving the field of epigenetics, and facilitating a better understanding of the field of functional genomics. The present review outlines the recent advances in epigenomic tools and techniques and presents the future perspectives of epigenomic studies.

Somatic epigenetic drift during shoot branching: a cell lineage-based model.

Plant architecture is shaped by the production of new organs, most of which emerge postembryonically. This process includes the formation of new lateral branches along existing shoots. Current evidence supports a detached-meristem model as the cellular basis of lateral shoot initiation. In this model, a small number of undifferentiated cells are sampled from the periphery of the shoot apical meristem (SAM) to act as precursors for axillary buds, which eventually develop into new shoots. Repeated branching thus creates cellular bottlenecks (i.e. somatic drift) that affect how de novo (epi)genetic mutations propagate through the plant body during development. Somatic drift could be particularly relevant for stochastic DNA methylation gains and losses (i.e. spontaneous epimutations), as they have been shown to arise rapidly with each cell division. Here, we formalize a special case of the detached-meristem model, where precursor cells are randomly sampled from the SAM periphery in a way that maximizes cell lineage independence. We show that somatic drift during repeated branching gives rise to a mixture of cellular phylogenies within the SAM over time. This process is dependent on the number of branch points, the strength of drift as well as the epimutation rate. Our model predicts that cell-to-cell DNA methylation heterogeneity in the SAM converges to nonzero states during development, suggesting that epigenetic variation is an inherent property of the SAM cell population. Our insights have direct implications for empirical studies of somatic (epi)genomic diversity in long-lived perennial and clonal species using bulk or single-cell sequencing approaches.

Phenotypical variability affecting the commercial value of the stalked barnacle Pollicipes pollicipes: No evidence for epigenetic variation

The stalked barnacle Pollicipes pollicipes is an important economic resource in Portugal and Spain. Two extreme phenotypes can be identified, based on their morphology. More elongated barnacles are associated with bad quality and have a lower commercial value.The fishers perception about the existence, definition and causes for this phenotypical/quality variation was evaluated through a survey performed in Portugal and Galicia, Spain. The existence of two extreme commercial qualities was validated. Good quality barnacles were mainly defined as thick and short in both countries. In Spain (Galicia), the definition of bad quality animals corresponded mainly to the terms long and thin, while in Portugal, fishers used a wider variety of terms including watery, thin and soft. The characteristics of the rock and the hydrodynamics were the causes most referred by the fishers for this variation.The morphological variation of P. pollicipes was described by the ratio between maximal rostro-carinal length (RC) and total height (TH): RC/TH values > 0.4 indicate good quality barnacles; and values < 0.4 indicate bad quality barnacles.Although morphological variation between the two extreme qualities/phenotypes was found, no genetic (amplified fragment length polymorphism - AFLP) or epigenetic (methylation sensitive amplification polymorphism - MSAP) differences were detected.

An epigenome-wide association study of waist circumference in Chinese monozygotic twins.

Central obesity poses significant health risks because it increases susceptibility to multiple chronic diseases. Epigenetic features such as DNA methylation may be associated with specific obesity traits, which could help us understand how genetic and environmental factors interact to influence the development of obesity. This study aims to identify DNA methylation sites associated with the waist circumference (WC) in Northern Han Chinese population, and to elucidate potential causal relationships. A total of 59 pairs of WC discordant monozygotic twins (ΔWC >0) were selected from the Qingdao Twin Registry in China. Generalized estimated equation model was employed to estimate the methylation levels of CpG sites on WC. Causal relationships between methylation and WC were assessed through the examination of family confounding factors using FAmiliaL CONfounding (ICE FALCON). Additionally, the findings of the epigenome-wide analysis were corroborated in the validation stage. We identified 26 CpG sites with differential methylation reached false discovery rate (FDR) < 0.05 and 22 differentially methylated regions (slk-corrected p < 0.05) strongly linked to WC. These findings provided annotations for 26 genes, with notable emphasis on MMP17, ITGA11, COL23A1, TFPI, A2ML1-AS1, MRGPRE, C2orf82, and NINJ2. ICE FALCON analysis indicated the DNA methylation of ITGA11 and TFPI had a causal effect on WC and vice versa (p < 0.05). Subsequent validation analysis successfully replicated 10 (p < 0.05) out of the 26 identified sites. Our research has ascertained an association between specific epigenetic variations and WC in the Northern Han Chinese population. These DNA methylation features can offer fresh insights into the epigenetic regulation of obesity and WC as well as hints to plausible biological mechanisms.

Genome wide inherited modifications of the tomato epigenome by trans-activated bacterial CG methyltransferase

BackgroundEpigenetic variation is mediated by epigenetic marks such as DNA methylation occurring in all cytosine contexts in plants. CG methylation plays a critical role in silencing transposable elements and regulating gene expression. The establishment of CG methylation occurs via the RNA-directed DNA methylation pathway and CG methylation maintenance relies on METHYLTRANSFERASE1, the homologue of the mammalian DNMT1.PurposeHere, we examined the capacity to stably alter the tomato genome methylome by a bacterial CG-specific M.SssI methyltransferase expressed through the LhG4/pOP transactivation system.ResultsMethylome analysis of M.SssI expressing plants revealed that their euchromatic genome regions are specifically hypermethylated in the CG context, and so are most of their genes. However, changes in gene expression were observed only with a set of genes exhibiting a greater susceptibility to CG hypermethylation near their transcription start site. Unlike gene rich genomic regions, our analysis revealed that heterochromatic regions are slightly hypomethylated at CGs only. Notably, some M.SssI-induced hypermethylation persisted even without the methylase or transgenes, indicating inheritable epigenetic modification.ConclusionCollectively our findings suggest that heterologous expression of M.SssI can create new inherited epigenetic variations and changes in the methylation profiles on a genome wide scale. This open avenues for the conception of epigenetic recombinant inbred line populations with the potential to unveil agriculturally valuable tomato epialleles.

Sequence and epigenetic variations of R2R3-MYB transcription factors determine the diversity of taproot skin and flesh colors in different cultivated types of radish (Raphanus sativus L.).

This study found that three paralogous R2R3-MYB transcription factors exhibit functional divergence among different subspecies and cultivated types in radish. Cultivated radish taproots exhibit a wide range of color variations due to unique anthocyanin accumulation patterns in various tissues. This study investigated the universal principles of taproot color regulation that developed during domestication of different subspecies and cultivated types. The key candidate genes RsMYB1 and RsMYB2, which control anthocyanin accumulation in radish taproots, were identified using bulked segregant analysis in two genetic populations. We introduced the RsMYB1-RsF3'H-RsMYB1Met genetic model to elucidate the complex and unstable genetic regulation of taproot flesh color in Xinlimei radish. Furthermore, we analyzed the expression patterns of three R2R3-MYB transcription factors in lines with different taproot colors and investigated the relationship between RsMYB haplotypes and anthocyanin accumulation in a natural population of 56 germplasms. The results revealed that three paralogous RsMYBs underwent functional divergence during radish domestication, with RsMYB1 regulating the red flesh of Xinlimei radish, and RsMYB2 and RsMYB3 regulating the red skin of East Asian big long radish (R. sativus var. hortensis) and European small radish (R. sativus var. sativus), respectively. Moreover, RsMYB1-H1, RsMYB2-H10, and RsMYB3-H6 were identified as the primary haplotypes exerting regulatory functions on anthocyanin synthesis. These findings provide an understanding of the genetic mechanisms regulating anthocyanin synthesis in radish and offer a potential strategy for early prediction of color variations in breeding programs.