Differential DNA Methylation Patterns Are Related to Phellogen Origin and Quality of Quercus suber Cork

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

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

Epigenetic Modifications Related to Potato Skin Russeting

Potato tuber skin is a protective corky tissue consisting of suberized phellem cells. Smooth-skinned varieties are characterized by a clean, shiny appearance compared to the darker hue of russeted potatoes. The rough skin of russeted cultivars is a desired, genetically inherited characteristic; however, unwanted russeting of smooth-skinned cultivars often occurs under suboptimal growth conditions. The involvement of epigenetic modifiers in regulating the smooth skin russeting disorder was tested. We used smooth-skin commercial cultivars with and without the russeting disorder and three lines from a breeding population segregating for russeting. Anatomically, the russet skin showed similar characteristics, whether the cause was environmentally triggered or genetically determined. The old outer layers of the corky phellem remain attached to the newly formed phellem layers instead of being sloughed off. Global DNA methylation analysis indicated a significant reduction in the percentage of 5-methylcytosine in mature vs. immature skin and russet vs. smooth skin. This was true for both the smooth-skin commercial cultivars and the russeted lines. The expression level of selected DNA methyltransferases was reduced in accordance. DNA demethylase expression did not change between the skin types and age. Hence, the reduced DNA methylation in mature and russet skin is more likely to be achieved through passive DNA demethylation and loss of methyltransferase activity.

Cork Development: What Lies Within

The cork layer present in all dicotyledonous plant species with radial growth is the result of the phellogen activity, a secondary meristem that produces phellem (cork) to the outside and phelloderm inwards. These three different tissues form the periderm, an efficient protective tissue working as a barrier against external factors such as environmental aggressions and pathogen attacks. The protective function offered by cork cells is mainly due to the abundance of suberin in their cell walls. Chemically, suberin is a complex aliphatic network of long chain fatty acids and alcohols with glycerol together with aromatic units. In most woody species growing in temperate climates, the first periderm is replaced by a new functional periderm upon a few years after being formed. One exception to this bark development can be found in cork oak (Quercus suber) which display a single periderm that grows continuously. Quercus suber stands by its thick cork layer development with continuous seasonal growth. Cork raw material has been exploited by man for centuries, especially in Portugal and Spain. Nowadays, its applications have widened vastly, from the most known product, stoppers, to purses or insulating materials used in so many industries, such as construction and car production. Research on how cork develops, and the effect environmental factors on cork oak trees is extremely important to maintain production of good-quality cork, and, by maintaining cork oak stands wealthy, we are preserving a very important ecosystem both by its biodiversity and its vital social and economic role in areas already showing a population declination.

Multiomics Molecular Research into the Recalcitrant and Orphan Quercus ilex Tree Species: Why, What for, and How

The holm oak (Quercus ilex L.) is the dominant tree species of the Mediterranean forest and the Spanish agrosilvopastoral ecosystem, “dehesa.” It has been, since the prehistoric period, an important part of the Iberian population from a social, cultural, and religious point of view, providing an ample variety of goods and services, and forming the basis of the economy in rural areas. Currently, there is renewed interest in its use for dietary diversification and sustainable food production. It is part of cultural richness, both economically (tangible) and environmentally (intangible), and must be preserved for future generations. However, a worrisome degradation of the species and associated ecosystems is occurring, observed in an increase in tree decline and mortality, which requires urgent action. Breeding programs based on the selection of elite genotypes by molecular markers is the only plausible biotechnological approach. To this end, the authors’ group started, in 2004, a research line aimed at characterizing the molecular biology of Q. ilex. It has been a challenging task due to its biological characteristics (long life cycle, allogamous, high phenotypic variability) and recalcitrant nature. The biology of this species has been characterized following the central dogma of molecular biology using the omics cascade. Molecular responses to biotic and abiotic stresses, as well as seed maturation and germination, are the two main objectives of our research. The contributions of the group to the knowledge of the species at the level of DNA-based markers, genomics, epigenomics, transcriptomics, proteomics, and metabolomics are discussed here. Moreover, data are compared with those reported for Quercus spp. All omics data generated, and the genome of Q. ilex available, will be integrated with morphological and physiological data in the systems biology direction. Thus, we will propose possible molecular markers related to resilient and productive genotypes to be used in reforestation programs. In addition, possible markers related to the nutritional value of acorn and derivate products, as well as bioactive compounds (peptides and phenolics) and allergens, will be suggested. Subsequently, the selected molecular markers will be validated by both genome-wide association and functional genomic analyses.

Genome-wide identification, phylogeny, and gene duplication of the epigenetic regulators in Fagaceae

Epigenetic regulators are proteins involved in controlling gene expression. Information about the epigenetic regulators within the Fagaceae, a relevant family of trees and shrubs of the northern hemisphere ecosystems, is scarce. With the intent to characterize these proteins in Fagaceae, we searched for orthologs of DNA methyltransferases (DNMTs) and demethylases (DDMEs) and Histone modifiers involved in acetylation (HATs), deacetylation (HDACs), methylation (HMTs), and demethylation (HDMTs) in Fagus, Quercus, and Castanea genera. Blast searches were performed in the available genomes, and freely available RNA-seq data were used to de novo assemble transcriptomes. We identified homologs of seven DNMTs, three DDMEs, six HATs, 11 HDACs, 32 HMTs, and 21 HDMTs proteins. Protein analysis showed that most of them have the putative characteristic domains found in these protein families, which suggests their conserved function. Additionally, to elucidate the evolutionary history of these genes within Fagaceae, paralogs were identified, and phylogenetic analyses were performed with DNA and histone modifiers. We detected duplication events in all species analyzed with higher frequency in Quercus and Castanea and discuss the evidence of transposable elements adjacent to paralogs and their involvement in gene duplication. The knowledge gathered from this work is a steppingstone to upcoming studies concerning epigenetic regulation in this economically important family of Fagaceae.

Potato Periderm Development and Tuber Skin Quality

The periderm is a corky tissue that replaces the epidermis when the latter is damaged, and is critical for preventing pathogen invasion and water loss. The periderm is formed through the meristematic activity of phellogen cells (cork cambium). The potato skin (phellem cells) composes the outer layers of the tuber periderm and is a model for studying cork development. Early in tuber development and following tuber expansion, the phellogen becomes active and produces the skin. New skin layers are continuously added by division of the phellogen cells until tuber maturation. Some physiological disorders of the potato tuber are related to abnormal development of the skin, including skinning injuries and russeting of smooth-skinned potatoes. Thus, characterizing the potato periderm contributes to modeling cork development in plants and helps to resolve critical agricultural problems. Here, we summarize the data available on potato periderm formation, highlighting tissue characteristics rather than the suberization processes.

Epigenetics at the crossroads of secondary growth regulation

The development of plant tissues and organs during post-embryonic growth occurs through the activity of both primary and secondary meristems. While primary meristems (root and shoot apical meristems) promote axial plant growth, secondary meristems (vascular and cork cambium or phellogen) promote radial thickening and plant axes strengthening. The vascular cambium forms the secondary xylem and phloem, whereas the cork cambium gives rise to the periderm that envelops stems and roots. Periderm takes on an increasingly important role in plant survival under climate change scenarios, but it is also a forest product with unique features, constituting the basis of a sustainable and profitable cork industry. There is established evidence that epigenetic mechanisms involving histone post-translational modifications, DNA methylation, and small RNAs play important roles in the activity of primary meristem cells, their maintenance, and differentiation of progeny cells. Here, we review the current knowledge on the epigenetic regulation of secondary meristems, particularly focusing on the phellogen activity. We also discuss the possible involvement of DNA methylation in the regulation of periderm contrasting phenotypes, given the potential impact of translating this knowledge into innovative breeding programs.

Quercus Conservation Genetics and Genomics: Past, Present, and Future

Quercus species (oaks) have been an integral part of the landscape in the northern hemisphere for millions of years. Their ability to adapt and spread across different environments and their contributions to many ecosystem services is well documented. Human activity has placed many oak species in peril by eliminating or adversely modifying habitats through exploitative land usage and by practices that have exacerbated climate change. The goal of this review is to compile a list of oak species of conservation concern, evaluate the genetic data that is available for these species, and to highlight the gaps that exist. We compiled a list of 124 Oaks of Concern based on the Red List of Oaks 2020 and the Conservation Gap Analysis for Native U.S. Oaks and their evaluations of each species. Of these, 57% have been the subject of some genetic analysis, but for most threatened species (72%), the only genetic analysis was done as part of a phylogenetic study. While nearly half (49%) of published genetic studies involved population genetic analysis, only 16 species of concern (13%) have been the subject of these studies. This is a critical gap considering that analysis of intraspecific genetic variability and genetic structure are essential for designing conservation management strategies. We review the published population genetic studies to highlight their application to conservation. Finally, we discuss future directions in Quercus conservation genetics and genomics.

Molecular Research on Stress Responses in Quercus spp.: From Classical Biochemistry to Systems Biology through Omics Analysis

The genus Quercus (oak), family Fagaceae, comprises around 500 species, being one of the most important and dominant woody angiosperms in the Northern Hemisphere. Nowadays, it is threatened by environmental cues, which are either of biotic or abiotic origin. This causes tree decline, dieback, and deforestation, which can worsen in a climate change scenario. In the 21st century, biotechnology should take a pivotal role in facing this problem and proposing sustainable management and conservation strategies for forests. As a non-domesticated, long-lived species, the only plausible approach for tree breeding is exploiting the natural diversity present in this species and the selection of elite, more resilient genotypes, based on molecular markers. In this direction, it is important to investigate the molecular mechanisms of the tolerance or resistance to stresses, and the identification of genes, gene products, and metabolites related to this phenotype. This research is being performed by using classical biochemistry or the most recent omics (genomics, epigenomics, transcriptomics, proteomics, and metabolomics) approaches, which should be integrated with other physiological and morphological techniques in the Systems Biology direction. This review is focused on the current state-of-the-art of such approaches for describing and integrating the latest knowledge on biotic and abiotic stress responses in Quercus spp., with special reference to Quercus ilex, the system on which the authors have been working for the last 15 years. While biotic stress factors mainly include fungi and insects such as Phytophthora cinnamomi, Cerambyx welensii, and Operophtera brumata, abiotic stress factors include salinity, drought, waterlogging, soil pollutants, cold, heat, carbon dioxide, ozone, and ultraviolet radiation. The review is structured following the Central Dogma of Molecular Biology and the omic cascade, from DNA (genomics, epigenomics, and DNA-based markers) to metabolites (metabolomics), through mRNA (transcriptomics) and proteins (proteomics). An integrated view of the different approaches, challenges, and future directions is critically discussed.

Enabling reusability of plant phenomic datasets with MIAPPE 1.1

Enabling data reuse and knowledge discovery is increasingly critical in modern science, and requires an effort towards standardising data publication practices. This is particularly challenging in the plant phenotyping domain, due to its complexity and heterogeneity. We have produced the MIAPPE 1.1 release, which enhances the existing MIAPPE standard in coverage, to support perennial plants, in structure, through an explicit data model, and in clarity, through definitions and examples. We evaluated MIAPPE 1.1 by using it to express several heterogeneous phenotyping experiments in a range of different formats, to demonstrate its applicability and the interoperability between the various implementations. Furthermore, the extended coverage is demonstrated by the fact that one of the datasets could not have been described under MIAPPE 1.0. MIAPPE 1.1 marks a major step towards enabling plant phenotyping data reusability, thanks to its extended coverage, and especially the formalisation of its data model, which facilitates its implementation in different formats. Community feedback has been critical to this development, and will be a key part of ensuring adoption of the standard.

Enabling reusability of plant phenomic datasets with MIAPPE 1.1

Enabling data reuse and knowledge discovery is increasingly critical in modern science, and requires an effort towards standardising data publication practices. This is particularly challenging in the plant phenotyping domain, due to its complexity and heterogeneity. We have produced the MIAPPE 1.1 release, which enhances the existing MIAPPE standard in coverage, to support perennial plants, in structure, through an explicit data model, and in clarity, through definitions and examples. We evaluated MIAPPE 1.1 by using it to express several heterogeneous phenotyping experiments in a range of different formats, to demonstrate its applicability and the interoperability between the various implementations. Furthermore, the extended coverage is demonstrated by the fact that one of the datasets could not have been described under MIAPPE 1.0. MIAPPE 1.1 marks a major step towards enabling plant phenotyping data reusability, thanks to its extended coverage, and especially the formalisation of its data model, which facilitates its implementation in different formats. Community feedback has been critical to this development, and will be a key part of ensuring adoption of the standard.

Enabling reusability of plant phenomic datasets with MIAPPE 1.1

Enabling data reuse and knowledge discovery is increasingly critical in modern science, and requires an effort towards standardising data publication practices. This is particularly challenging in the plant phenotyping domain, due to its complexity and heterogeneity. We have produced the MIAPPE 1.1 release, which enhances the existing MIAPPE standard in coverage, to support perennial plants, in structure, through an explicit data model, and in clarity, through definitions and examples. We evaluated MIAPPE 1.1 by using it to express several heterogeneous phenotyping experiments in a range of different formats, to demonstrate its applicability and the interoperability between the various implementations. Furthermore, the extended coverage is demonstrated by the fact that one of the datasets could not have been described under MIAPPE 1.0. MIAPPE 1.1 marks a major step towards enabling plant phenotyping data reusability, thanks to its extended coverage, and especially the formalisation of its data model, which facilitates its implementation in different formats. Community feedback has been critical to this development, and will be a key part of ensuring adoption of the standard.

Genome-Wide Identification of Epigenetic Regulators in Quercus suber L

Modifications of DNA and histones, including methylation and acetylation, are critical for the epigenetic regulation of gene expression during plant development, particularly during environmental adaptation processes. However, information on the enzymes catalyzing all these modifications in trees, such as Quercus suber L., is still not available. In this study, eight DNA methyltransferases (DNA Mtases) and three DNA demethylases (DDMEs) were identified in Q. suber. Histone modifiers involved in methylation (35), demethylation (26), acetylation (8), and deacetylation (22) were also identified in Q. suber. In silico analysis showed that some Q. suber DNA Mtases, DDMEs and histone modifiers have the typical domains found in the plant model Arabidopsis, which might suggest a conserved functional role. Additional phylogenetic analyses of the DNA and histone modifier proteins were performed using several plant species homologs, enabling the classification of the Q. suber proteins. A link between the expression levels of each gene in different Q. suber tissues (buds, flowers, acorns, embryos, cork, and roots) with the functions already known for their closest homologs in other species was also established. Therefore, the data generated here will be important for future studies exploring the role of epigenetic regulators in this economically important species.

Phellem versus xylem: genome-wide transcriptomic analysis reveals novel regulators of cork formation in cork oak

Cork cambium (or phellogen) is a secondary meristem responsible for the formation of phelloderm and phellem/cork, which together compose the periderm. In Quercus suber L., the phellogen is active throughout the entire life of the tree, producing a continuous and renewable outer bark of cork. To identify specific candidate genes associated with cork cambium activity and phellem differentiation, we performed a comparative transcriptomic study of Q. suber secondary growth tissues (xylem and phellogen/phellem) using RNA-seq. The present work provides a high-resolution map of all the transcripts identified in the phellogen/phellem tissues. A total of 6013 differentially expressed genes were identified, with 2875 of the transcripts being specifically enriched during the cork formation process versus secondary xylem formation. Furthermore, cork samples originating from the original phellogen (`virgin' cork) and from a traumatic phellogen (`amadia' cork) were also compared. Our results point to a shortlist of potentially relevant candidate genes regulating phellogen activity and phellem differentiation, including novel genes involved in the suberization process, as well as genes associated to ethylene and jasmonate signaling and to meristem function. The future functional characterization of some of the identified candidate genes will help to elucidate the molecular mechanisms underlying cork cambium activity and phellem differentiation.

Phellem Cell-Wall Components Are Discriminants of Cork Quality in Quercus suber

Cork is a renewable, non-wood high valued forest product, with relevant ecological and economic impact in the Mediterranean-type ecosystems. Currently, cork is ranked according to its commercial quality. The most valuable planks are chosen for cork stoppers production. Cork planks with adequate thickness and porosity are classified as stoppable quality cork (SQC). The chemical composition of cork is known, but the regulation of metabolic pathways responsible of cork production and composition, hence of cork quality, is largely unknown. Here, we tested the hypothesis that post-genomic events may be responsible for the development of SQC and N-SQC (non-stoppable quality cork). Here, we show that combined proteomics and targeted metabolomics (namely soluble and cell wall bound phenolics) analyzed on recently formed phellem allows discriminate cork planks of different quality. Phellem cells of SQC and N-SQC displayed different reducing capacity, with consequential impact on both enzymatic pathways (e.g., glycolysis) and other cellular functions, including cell wall assembly and suberization. Glycolysis and respiration related proteins were abundant in both cork quality groups, whereas the level of several proteins associated to mitochondrial metabolism was higher in N-SQC. The soluble and cell wall-bound phenolics in recently formed phellem clearly discriminated SQC from N-SCQ. In our study, SQC was characterized by a high incorporation of aromatic components of the phenylpropanoid pathway in the cell wall, together with a lower content of hydrolysable tannins. Here, we propose that the level of hydrolysable tannins may represent a valuable diagnostic tool for screening recently formed phellem, and used as a proxy for the quality grade of cork plank produced by each tree.

Cork Oak Young and Traumatic Periderms Show PCD Typical Chromatin Patterns but Different Chromatin-Modifying Genes Expression

Plants are subjected to adverse conditions being outer protective tissues fundamental to their survival. Tree stems are enveloped by a periderm made of cork cells, resulting from the activity of the meristem phellogen. DNA methylation and histone modifications have important roles in the regulation of plant cell differentiation. However, studies on its involvement in cork differentiation are scarce despite periderm importance. Cork oak periderm development was used as a model to study the formation and differentiation of secondary protective tissues, and their behavior after traumatic wounding (traumatic periderm). Nuclei structural changes, dynamics of DNA methylation, and posttranslational histone modifications were assessed in young and traumatic periderms, after cork harvesting. Lenticular phellogen producing atypical non-suberized cells that disaggregate and form pores was also studied, due to high impact for cork industrial uses. Immunolocalization of active and repressive marks, transcription analysis of the corresponding genes, and correlations between gene expression and cork porosity were investigated. During young periderm development, a reduction in nuclei area along with high levels of DNA methylation occurred throughout epidermis disruption. As cork cells became more differentiated, whole nuclei progressive chromatin condensation with accumulation in the nuclear periphery and increasing DNA methylation was observed. Lenticular cells nuclei were highly fragmented with faint 5-mC labeling. Phellogen nuclei were less methylated than in cork cells, and in lenticular phellogen were even lower. No significant differences were detected in H3K4me3 and H3K18ac signals between cork cells layers, although an increase in H3K4me3 signals was found from the phellogen to cork cells. Distinct gene expression patterns in young and traumatic periderms suggest that cork differentiation might be under specific silencing regulatory pathways. Significant correlations were found between QsMET1, QsMET2, and QsSUVH4 gene expression and cork porosity. This work evidences that DNA methylation and histone modifications play a role in cork differentiation and epidermis induced tension-stress. It also provides the first insights into chromatin dynamics during cork and lenticular cells differentiation pointing to a distinct type of remodeling associated with cell death.