BackgroundCork oak (Quercus suber) is one of the rare trees with the ability to produce cork, a material widely used to make wine bottle stoppers, flooring and insulation materials, among many other uses. The molecular mechanisms of cork formation are still poorly understood, in great part due to the difficulty in studying a species with a long life-cycle and for which there is scarce molecular/genomic information. Cork oak forests are of great ecological importance and represent a major economic and social resource in Southern Europe and Northern Africa. However, global warming is threatening the cork oak forests by imposing thermal, hydric and many types of novel biotic stresses. Despite the economic and social value of the Q. suber species, few genomic resources have been developed, useful for biotechnological applications and improved forest management.ResultsWe generated in excess of 7 million sequence reads, by pyrosequencing 21 normalized cDNA libraries derived from multiple Q. suber tissues and organs, developmental stages and physiological conditions. We deployed a stringent sequence processing and assembly pipeline that resulted in the identification of ~159,000 unigenes. These were annotated according to their similarity to known plant genes, to known Interpro domains, GO classes and E.C. numbers. The phylogenetic extent of this ESTs set was investigated, and we found that cork oak revealed a significant new gene space that is not covered by other model species or EST sequencing projects. The raw data, as well as the full annotated assembly, are now available to the community in a dedicated web portal at http://www.corkoakdb.org.ConclusionsThis genomic resource represents the first trancriptome study in a cork producing species. It can be explored to develop new tools and approaches to understand stress responses and developmental processes in forest trees, as well as the molecular cascades underlying cork differentiation and disease response.

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.

Cork oak (Quercus suber ) is an important Portuguese species, mainly due to the economic value of the cork it produces. Cork results from phellogen, a meristematic tissue, which can locally produce lenticels or have discontinuities, originating "defects": pores and nail inclusions that are detrimental to cork industrial use. Epigenetic processes control plant development and its deregulation can lead to altered phenotypes; therefore, the study of epigenetic players in the phellogen is important to understand the emergence of cork's defects. DNA methyltransferases (DNMTs) and one protein associated to MET1 (DMAP1) were characterized in Q. suber, and their gene expression was analyzed in phellogen and contiguous differentiating cell layers of trees producing high and low quality cork, after the evaluation of their defects by physical and image analysis methods. All classes of DNMTs (MET, DRM, and CMT) with the respective canonical motifs were identified in Q. suber. The expression analyses of these genes showed that QsDRM2 was the most active methyltransferases in the cells analyzed, and that all the genes were differentially expressed in trees with distinct cork quality, with a tendency for higher expression levels in low quality producers. Interestingly, the global methylation level was higher in cells with low expression of DNA methyltransferases. A positive and significant correlation was obtained between QsDMAP1 gene expression and the percentage of cork defects. This work provides the first evidence that cork quality in Q. suber is likely influenced by epigenetic mechanisms.

BackgroundGene activity is largely controlled by transcriptional regulation through the action of transcription factors and other regulators. QsMYB1 is a member of the R2R3-MYB transcription factor family related to secondary growth, and in particular, with the cork development process. In order to identify the putative gene targets of QsMYB1 across the cork oak genome we developed a ChIP-Seq strategy.ResultsResults provide direct evidence that QsMY1B targets genes encoding for enzymes involved in the lignin and suberin pathways as well as gene encoding for ABCG transporters and LTPs implicated in the transport of monomeric suberin units across the cellular membrane. These results highlight the role of QsMYB1 as a regulator of lignin and suberin biosynthesis, transport and assembly.ConclusionTo our knowledge, this work constitutes the first ChIP-Seq experiment performed in cork oak, a non-model plant species with a long-life cycle, and these results will contribute to deepen the knowledge about the molecular mechanisms of cork formation and differentiation.Electronic supplementary materialThe online version of this article (10.1186/s12870-018-1403-5) contains supplementary material, which is available to authorized users.

The 35S ribosomal DNA (rDNA) units, repeated in tandem at one or more chromosomal loci, are separated by an intergenic spacer (IGS) containing functional elements involved in the regulation of transcription of downstream rRNA genes. In the present work, we have compared the IGS molecular organizations in two divergent species of Fagaceae, Fagus sylvatica and Quercus suber, aiming to comprehend the evolution of the IGS sequences within the family. Self- and cross-hybridization FISH was done on representative species of the Fagaceae. The IGS length variability and the methylation level of 18 and 25S rRNA genes were assessed in representatives of three genera of this family: Fagus, Quercus and Castanea. The intergenic spacers in Beech and Cork Oak showed similar overall organizations comprising putative functional elements needed for rRNA gene activity and containing a non-transcribed spacer (NTS), a promoter region, and a 5′-external transcribed spacer. In the NTS: the sub-repeats structure in Beech is more organized than in Cork Oak, sharing some short motifs which results in the lowest sequence similarity of the entire IGS; the AT-rich region differed in both spacers by a GC-rich block inserted in Cork Oak. The 5′-ETS is the region with the higher similarity, having nonetheless different lengths. FISH with the NTS-5′-ETS revealed fainter signals in cross-hybridization in agreement with the divergence between genera. The diversity of IGS lengths revealed variants from ∼2 kb in Fagus, and Quercus up to 5.3 kb in Castanea, and a lack of correlation between the number of variants and the number of rDNA loci in several species. Methylation of 25S Bam HI site was confirmed in all species and detected for the first time in the 18S of Q. suber and Q. faginea. These results provide important clues for the evolutionary trends of the rDNA 25S-18S IGS in the Fagaceae family.

DNA methylation is thought to influence Quercus suber cork quality, which is the main constraint for its economic valorisation. However, a deep knowledge of the cytosine methylation patterns disclosing the epigenetic variability of trees with different cork quality types is totally missing. This study investigates the hypothesis that variations in DNA methylation contribute to differences in cork cellular characteristics directly related to original or traumatic phellogen activity. We used MSAPs (Methylation Sensitive Amplified Polymorphism) to assess DNA methylation patterns of cork and leaf tissues of Q. suber adult trees growing in three cork oak stands. The relationship between the detected polymorphisms and the diversity of cork quality traits was explored by a marker-trait analysis focusing on the most relevant quality characteristics. Populations differed widely in cork quality, but only slightly in degree of epigenetic differentiation. Four MSAP markers (1.3% of the total) were significantly associated with the most noteworthy quality traits: wood inclusions (nails) and porosity. This evidence supports the potential role of cytosine methylation in the modulation of differential phellogen activity either involved in localized cell death or in pore production, resulting in different cork qualities. Although, the underlying basis of the methylation polymorphism of loci affecting cork quality traits remain unclear, the disclosure of markers statistically associated with cork quality strengthens the potential role of DNA methylation in the regulation of these traits, namely at the phellogen level.

Oaks, chestnuts, and beeches are economically important species of the Fagaceae. To understand the relationship between these members of this family, a deep knowledge of their genome composition and organization is needed. In this work, we have isolated and characterized several AFLP fragments obtained from Quercus rotundifolia Lam. through homology searches in available databases. Genomic polymorphisms involving some of these sequences were evaluated in two species of Quercus, one of Castanea, and one of Fagus with specific primers. Comparative FISH analysis with generated sequences was performed in interphase nuclei of the four species, and the co-immunolocalization of 5-methylcytosine was also studied. Some of the sequences isolated proved to be genus-specific, while others were present in all the genera. Retroelements, either gypsy-like of the Tat/Athila clade or copia-like, are well represented, and most are dispersed in euchromatic regions of these species with no DNA methylation associated, pointing to an interspersed arrangement of these retroelements with potential gene-rich regions. A particular gypsy-sequence is dispersed in oaks and chestnut nuclei, but its confinement to chromocenters in beech evidences genome restructuring events during evolution of Fagaceae. Several sequences generated in this study proved to be good tools to comparatively study Fagaceae genome organization.

Vascular plants with secondary growth develop a periderm mostly composed of dead suberized cork cells to face environmental hostile conditions. Cork oak has a highly active and long-living phellogen forming a remarkably thick periderm that is periodically debarked for industrial purposes. This wounding originates the quick formation of a new traumatic periderm, making cork oak an exceptional model to study the first periderm differentiation during normal development in young sprigs and traumatic (wound) periderm formation after debarking. Here we studied the poorly known first periderm differentiation steps that involve cell wall suberization, polyphenolic accumulation, and programmed cell death (PCD) by combining TEM, histochemical, and molecular methods in periderms from young sprigs. These processes were further compared with traumatic periderms formed after wounding using molecular and histochemical techniques, such as the polyphenolic accumulation. In the first periderms from young sprigs, four distinct differentiation stages were defined according to the presence of PCD morphological features. First young and traumatic periderms showed an upregulation of genes related to suberin biosynthesis, proanthocyanidins biosynthesis and transport, autophagy, and PCD. Traumatic periderms revealed an overall upregulation of these genes, likely resulting from ontogeny differences and distinct phellogen origin associated with a faster metabolism, highlighting the impact of wounding on phellogen activity after debarking. First periderms from young sprigs showed gradual accumulation of proanthocyanidins in the vacuoles throughout PCD stages until total filled lumens, whereas in traumatic periderms these compounds were found cell wall linked in already empty cells. This work enabled a comprehensive overview of the cork cells differentiation processes contributing to deepening the knowledge of the fundamental ontogenic program of this protective tissue, which is also a unique forest product, constituting the basis of a sustainable and profitable industry.