Functional investigation of five R2R3-MYB transcription factors associated with wood development in Eucalyptus using DAP-seq-ML

A multi-tiered transcriptional network regulates xylem differentiation and secondary cell wall (SCW) formation in plants, with evidence of both conserved and lineage-specific SCW network architecture. We aimed to elucidate the roles of selected R2R3-MYB transcription factors (TFs) linked to Eucalyptus wood formation by identifying genome-wide TF binding sites and direct target genes through an improved DAP-seq protocol combined with machine learning for target gene assignment (DAP-seq-ML). We applied this to five TFs including a well-studied SCW master regulator (EgrMYB2; homolog of AtMYB83), a repressor of lignification (EgrMYB1; homolog of AtMYB4), a TF affecting SCW thickness and vessel density (EgrMYB137; homolog of PtrMYB074) and two TFs with unclear roles in SCW regulation (EgrMYB135 and EgrMYB122). Each DAP-seq TF peak set (average 12,613 peaks) was enriched for canonical R2R3-MYB binding motifs. To improve the reliability of target gene assignment to peaks, a random forest classifier was developed from Arabidopsis DAP-seq, RNA-seq, chromatin, and conserved noncoding sequence data which demonstrated significantly higher precision and recall to the baseline method of assigning genes to proximal peaks. EgrMYB1, EgrMYB2 and EgrMYB137 predicted targets showed clear enrichment for SCW-related biological processes. As validation, EgrMYB137 overexpression in transgenic Eucalyptus hairy roots increased xylem lignification, while its dominant repression in transgenic Arabidopsis and Populus reduced xylem lignification, stunted growth, and caused downregulation of SCW genes. EgrMYB137 targets overlapped significantly with those of EgrMYB2, suggesting partial functional redundancy. Our results show that DAP-seq-ML identified biologically relevant R2R3-MYB targets supported by the finding that EgrMYB137 promotes SCW lignification in planta.

Wood Architecture and Composition Are Deeply Remodeled in Frost Sensitive Eucalyptus Overexpressing CBF/DREB1 Transcription Factors

Eucalypts are the most planted trees worldwide, but most of them are frost sensitive. Overexpressing transcription factors for CRT-repeat binding factors (CBFs) in transgenic Eucalyptus confer cold resistance both in leaves and stems. While wood plays crucial roles in trees and is affected by environmental cues, its potential role in adaptation to cold stress has been neglected. Here, we addressed this question by investigating the changes occurring in wood in response to the overexpression of two CBFs, taking advantage of available transgenic Eucalyptus lines. We performed histological, biochemical, and transcriptomic analyses on xylem samples. CBF ectopic expression led to a reduction of both primary and secondary growth, and triggered changes in xylem architecture with smaller and more frequent vessels and fibers exhibiting reduced lumens. In addition, lignin content and syringyl/guaiacyl (S/G) ratio increased. Consistently, many genes of the phenylpropanoid and lignin branch pathway were upregulated. Most of the features of xylem remodeling induced by CBF overexpression are reminiscent of those observed after long exposure of Eucalyptus trees to chilling temperatures. Altogether, these results suggest that CBF plays a central role in the cross-talk between response to cold and wood formation and that the remodeling of wood is part of the adaptive strategies to face cold stress.

Long cold exposure induces transcriptional and biochemical remodelling of xylem secondary cell wall in Eucalyptus

Although eucalypts are the most planted hardwood trees worldwide, the majority of them are frost sensitive. The recent creation of frost-tolerant hybrids such as Eucalyptus gundal plants (E. gunnii × E. dalrympleana hybrids), now enables the development of industrial plantations in northern countries. Our objective was to evaluate the impact of cold on the wood structure and composition of these hybrids, and on the biosynthetic and regulatory processes controlling their secondary cell-wall (SCW) formation. We used an integrated approach combining histology, biochemical characterization and transcriptomic profiling as well as gene co-expression analyses to investigate xylem tissues from Eucalyptus hybrids exposed to cold conditions. Chilling temperatures triggered the deposition of thicker and more lignified xylem cell walls as well as regulation at the transcriptional level of SCW genes. Most genes involved in lignin biosynthesis, except those specifically dedicated to syringyl unit biosynthesis, were up-regulated. The construction of a co-expression network enabled the identification of both known and potential new SCW transcription factors, induced by cold stress. These regulators at the crossroads between cold signalling and SCW formation are promising candidates for functional studies since they may contribute to the tolerance of E. gunnii × E. dalrympleana hybrids to cold.

Special trends in CBF and DREB2 groups in Eucalyptus gunnii vs Eucalyptus grandis suggest that CBF are master players in the trade-off between growth and stress resistance

Annotation of the Eucalyptus grandis genome showed a large amplification of the dehydration-responsive element binding 1/C-repeat binding factor (DREB1/CBF) group without recent DREB2 gene duplication compared with other plant species. The present annotation of the CBF and DREB2 genes from a draft of the Eucalyptus gunnii genome sequence reveals at least one additional CBF copy in the E. gunnii genome compared with E. grandis, suggesting that this group is still evolving, unlike the DREB2 group. This study aims to investigate the redundancy/neo- or sub-functionalization of the duplicates and the relative involvement of the two groups in abiotic stress responses in both E. grandis and E. gunnii (lower growth but higher cold resistance). A comprehensive transcriptional analysis using high-throughput quantitative real-time polymerase chain reaction (qRT-PCR) was performed on leaves, stems and roots from the two Eucalyptus species after cold, heat or drought treatment. A large CBF cluster accounted for most of the cold response in all the organs, whereas heat and drought responses mainly involved a small CBF cluster and the DREB2 genes. In addition, CBF putative target genes, known to be involved in plant tolerance and development, were found to be cold-regulated. The higher transcript amounts of both the CBF and target genes in the cold tolerant E. gunnii contrasted with the higher CBF induction rates in the fast growing E. grandis. Altogether, the present results, in agreement with previous data about Eucalyptus transgenic lines over-expressing CBF, suggest that these factors, which promote both stress protection and growth limitation, participate in the trade-off between growth and resistance in this woody species.

Eucalyptus hairy roots, a fast, efficient and versatile tool to explore function and expression of genes involved in wood formation

Eucalyptus are of tremendous economic importance being the most planted hardwoods worldwide for pulp and paper, timber and bioenergy. The recent release of the Eucalyptus grandis genome sequence pointed out many new candidate genes potentially involved in secondary growth, wood formation or lineage-specific biosynthetic pathways. Their functional characterization is, however, hindered by the tedious, time-consuming and inefficient transformation systems available hitherto for eucalypts. To overcome this limitation, we developed a fast, reliable and efficient protocol to obtain and easily detect co-transformed E. grandis hairy roots using fluorescent markers, with an average efficiency of 62%. We set up conditions both to cultivate excised roots in vitro and to harden composite plants and verified that hairy root morphology and vascular system anatomy were similar to wild-type ones. We further demonstrated that co-transformed hairy roots are suitable for medium-throughput functional studies enabling, for instance, protein subcellular localization, gene expression patterns through RT-qPCR and promoter expression, as well as the modulation of endogenous gene expression. Down-regulation of the Eucalyptus cinnamoyl-CoA reductase1 (EgCCR1) gene, encoding a key enzyme in lignin biosynthesis, led to transgenic roots with reduced lignin levels and thinner cell walls. This gene was used as a proof of concept to demonstrate that the function of genes involved in secondary cell wall biosynthesis and wood formation can be elucidated in transgenic hairy roots using histochemical, transcriptomic and biochemical approaches. The method described here is timely because it will accelerate gene mining of the genome for both basic research and industry purposes.

The Woody-Preferential Gene EgMYB88 Regulates the Biosynthesis of Phenylpropanoid-Derived Compounds in Wood

Comparative phylogenetic analyses of the R2R3-MYB transcription factor family revealed that five subgroups were preferentially found in woody species and were totally absent from Brassicaceae and monocots (Soler et al., 2015). Here, we analyzed one of these subgroups (WPS-I) for which no gene had been yet characterized. Most Eucalyptus members of WPS-I are preferentially expressed in the vascular cambium, the secondary meristem responsible for tree radial growth. We focused on EgMYB88, which is the most specifically and highly expressed in vascular tissues, and showed that it behaves as a transcriptional activator in yeast. Then, we functionally characterized EgMYB88 in both transgenic Arabidopsis and poplar plants overexpressing either the native or the dominant repression form (fused to the Ethylene-responsive element binding factor-associated Amphiphilic Repression motif, EAR). The transgenic Arabidopsis lines had no phenotype whereas the poplar lines overexpressing EgMYB88 exhibited a substantial increase in the levels of the flavonoid catechin and of some salicinoid phenolic glycosides (salicortin, salireposide, and tremulacin), in agreement with the increase of the transcript levels of landmark biosynthetic genes. A change in the lignin structure (increase in the syringyl vs. guaiacyl, S/G ratio) was also observed. Poplar lines overexpressing the EgMYB88 dominant repression form did not show a strict opposite phenotype. The level of catechin was reduced, but the levels of the salicinoid phenolic glycosides and the S/G ratio remained unchanged. In addition, they showed a reduction in soluble oligolignols containing sinapyl p-hydroxybenzoate accompanied by a mild reduction of the insoluble lignin content. Altogether, these results suggest that EgMYB88, and more largely members of the WPS-I group, could control in cambium and in the first layers of differentiating xylem the biosynthesis of some phenylpropanoid-derived secondary metabolites including lignin.

Foliar treatments with Gaultheria procumbens essential oil induce defense responses and resistance against a fungal pathogen in Arabidopsis

Essential oil from Gaultheria procumbens is mainly composed of methylsalicylate (MeSA) (>96%), a compound which can be metabolized in plant tissues to salicylic acid, a phytohormone inducing plant immunity against microbial pathogens. The potential use of G. procumbens essential oil as a biocontrol agent was evaluated on the model plant Arabidopsis thaliana. Expression of a selection of defense genes was detected 1, 6, and 24 h after essential oil treatment (0.1 ml/L) using a high-throughput qPCR-based microfluidic technology. Control treatments included methyl jasmonate and a commercialized salicylic acid (SA) analog, benzo(1,2,3)-thiadiazole-7carbothiolic acid (BTH). Strong induction of defense markers known to be regulated by the SA pathway was observed after the treatment with G. procumbens essential oil. Treatment induced the accumulation of total SA in the wild-type Arabidopsis line Col-0 and analysis of the Arabidopsis line sid2, mutated in a SA biosynthetic gene, revealed that approximately 30% of MeSA sprayed on the leaves penetrated inside plant tissues and was demethylated by endogenous esterases. Induction of plant resistance by G. procumbens essential oil was tested following inoculation with a GFP-expressing strain of the Arabidopsis fungal pathogen Colletotrichum higginsianum. Fluorescence measurement of infected tissues revealed that treatments led to a strong reduction (60%) of pathogen development and that the efficacy of the G. procumbens essential oil was similar to the commercial product BION(®). Together, these results show that the G. procubens essential oil is a natural source of MeSA which can be formulated to develop new biocontrol products.

Pathogen-associated molecular pattern-triggered immunity and resistance to the root pathogen Phytophthora parasitica in Arabidopsis

The cellulose binding elicitor lectin (CBEL) of the genus Phytophthora induces necrosis and immune responses in several plant species, including Arabidopsis thaliana. However, the role of CBEL-induced responses in the outcome of the interaction is still unclear. This study shows that some of CBEL-induced defence responses, but not necrosis, required the receptor-like kinase BAK1, a general regulator of basal immunity in Arabidopsis, and the production of a reactive oxygen burst mediated by respiratory burst oxidases homologues (RBOH). Screening of a core collection of 48 Arabidopsis ecotypes using CBEL uncovered a large variability in CBEL-induced necrotic responses. Analysis of non-responsive CBEL lines Ws-4, Oy-0, and Bla-1 revealed that Ws-4 and Oy-0 were also impaired in the production of the oxidative burst and expression of defence genes, whereas Bla-1 was partially affected in these responses. Infection tests using two Phytophthora parasitica strains, Pp310 and Ppn0, virulent and avirulent, respectively, on the Col-0 line showed that BAK1 and RBOH mutants were susceptible to Ppn0, suggesting that some immune responses controlled by these genes, but not CBEL-induced cell death, are required for Phytophthora parasitica resistance. However, Ws-4, Oy-0, and Bla-1 lines were not affected in Ppn0 resistance, showing that natural variability in CBEL responsiveness is not correlated to Phytophthora susceptibility. Overall, the results uncover a BAK1- and RBOH-dependent CBEL-triggered immunity essential for Phytophthora resistance and suggest that natural quantitative variation of basal immunity triggered by conserved general elicitors such as CBEL does not correlate to Phytophthora susceptibility.

Reference genes for high-throughput quantitative reverse transcription-PCR analysis of gene expression in organs and tissues of Eucalyptus grown in various environmental conditions

Interest in the genomics of Eucalyptus has skyrocketed thanks to the recent sequencing of the genome of Eucalyptus grandis and to a growing number of large-scale transcriptomic studies. Quantitative reverse transcription-PCR (RT-PCR) is the method of choice for gene expression analysis and can now also be used as a high-throughput method. The selection of appropriate internal controls is becoming of utmost importance to ensure accurate expression results in Eucalyptus. To this end, we selected 21 candidate reference genes and used high-throughput microfluidic dynamic arrays to assess their expression among a large panel of developmental and environmental conditions with a special focus on wood-forming tissues. We analyzed the expression stability of these genes by using three distinct statistical algorithms (geNorm, NormFinder and ΔCt), and used principal component analysis to compare methods and rankings. We showed that the most stable genes identified depended not only on the panel of biological samples considered but also on the statistical method used. We then developed a comprehensive integration of the rankings generated by the three methods and identified the optimal reference genes for 17 distinct experimental sets covering 13 organs and tissues, as well as various developmental and environmental conditions. The expression patterns of Eucalyptus master genes EgMYB1 and EgMYB2 experimentally validated our selection. Our findings provide an important resource for the selection of appropriate reference genes for accurate and reliable normalization of gene expression data in the organs and tissues of Eucalyptus trees grown in a range of conditions including abiotic stresses.

EgMYB1, an R2R3 MYB transcription factor from eucalyptus negatively regulates secondary cell wall formation in Arabidopsis and poplar

• The eucalyptus R2R3 transcription factor, EgMYB1 contains an active repressor motif in the regulatory domain of the predicted protein. It is preferentially expressed in differentiating xylem and is capable of repressing the transcription of two key lignin genes in vivo. • In order to investigate in planta the role of this putative transcriptional repressor of the lignin biosynthetic pathway, we overexpressed the EgMYB1 gene in Arabidopsis and poplar. • Expression of EgMYB1 produced similar phenotypes in both species, with stronger effects in transgenic Arabidopsis plants than in poplar. Vascular development was altered in overexpressors showing fewer lignified fibres (in phloem and interfascicular zones in poplar and Arabidopsis, respectively) and reduced secondary wall thickening. Klason lignin content was moderately but significantly reduced in both species. Decreased transcript accumulation was observed for genes involved in the biosynthesis of lignins, cellulose and xylan, the three main polymers of secondary cell walls. Transcriptomic profiles of transgenic poplars were reminiscent of those reported when lignin biosynthetic genes are disrupted. • Together, these results strongly suggest that EgMYB1 is a repressor of secondary wall formation and provide new opportunities to dissect the transcriptional regulation of secondary wall biosynthesis.

A new genomic resource dedicated to wood formation in Eucalyptus

BACKGROUND

Renowned for their fast growth, valuable wood properties and wide adaptability, Eucalyptus species are amongst the most planted hardwoods in the world, yet they are still at the early stages of domestication because conventional breeding is slow and costly. Thus, there is huge potential for marker-assisted breeding programs to improve traits such as wood properties. To this end, the sequencing, analysis and annotation of a large collection of expressed sequences tags (ESTs) from genes involved in wood formation in Eucalyptus would provide a valuable resource.

RESULTS

We report here the normalization and sequencing of a cDNA library from developing Eucalyptus secondary xylem, as well as the construction and sequencing of two subtractive libraries (juvenile versus mature wood and vice versa). A total of 9,222 high quality sequences were collected from about 10,000 cDNA clones. The EST assembly generated a set of 3,857 wood-related unigenes including 2,461 contigs (Cg) and 1,396 singletons (Sg) that we named 'EUCAWOOD'. About 65% of the EUCAWOOD sequences produced matches with poplar, grapevine, Arabidopsis and rice protein sequence databases. BlastX searches of the Uniref100 protein database allowed us to allocate gene ontology (GO) and protein family terms to the EUCAWOOD unigenes. This annotation of the EUCAWOOD set revealed key functional categories involved in xylogenesis. For instance, 422 sequences matched various gene families involved in biosynthesis and assembly of primary and secondary cell walls. Interestingly, 141 sequences were annotated as transcription factors, some of them being orthologs of regulators known to be involved in xylogenesis. The EUCAWOOD dataset was also mined for genomic simple sequence repeat markers, yielding a total of 639 putative microsatellites. Finally, a publicly accessible database was created, supporting multiple queries on the EUCAWOOD dataset.

CONCLUSION

In this work, we have identified a large set of wood-related Eucalyptus unigenes called EUCAWOOD, thus creating a valuable resource for functional genomics studies of wood formation and molecular breeding in this economically important genus. This set of publicly available annotated sequences will be instrumental for candidate gene approaches, custom array development and marker-assisted selection programs aimed at improving and modulating wood properties.

Transcript profiling of a xylem vs phloem cDNA subtractive library identifies new genes expressed during xylogenesis in Eucalyptus

Eucalyptus is one of the world's main sources of biomass. The genus includes species representing the principle hardwood trees used for pulp and paper. Here, we aimed to identify genes specifically expressed in differentiating secondary xylem compared with phloem. We constructed a xylem vs phloem subtractive library (Xp) that generated 263 unique sequences. By transcript profiling of xylem, phloem, vascular cambium and leaves using macroarrays, we classified the 263 unigenes into distinct tissue-specific groups. Reverse transcription-polymerase chain reaction (RT-PCR) confirmed the differential expression of representative expressed sequence tags (ESTs). A total of 87 unigenes were preferentially expressed in xylem. They were involved in functional categories known to play roles in xylogenesis, such as hormone signaling and metabolism, secondary cell wall thickening and proteolysis. Some of these genes, including unknown genes, may be considered xylem-specific and they are likely to control important functions in xylogenesis. These data shed light on the cellular functions of xylem cells and, importantly, provide us with a portfolio of Eucalyptus xylem genes that may be major players in the control of wood formation and quality.

Comparative analysis of early embryonic sunflower cDNA libraries

To gain information concerning cell functions and activities during sunflower embryogenesis, an expressed sequence tag (EST) approach was used to analyse gene expression in the early stages of sunflower embryos development. Confocal microscopy observations of whole-mounted embryos allowed us to identify precisely the major steps of the zygotic embryonic development. A time-course analysis was then employed to collect the embryonic material. Three cDNA libraries were constructed from microdissected embryos, and three other cDNA libraries were created using a classical day after pollination schedule. A total of 7106 ESTs were produced and assembled. The total number of putative different genes represents about 43.1 (3064 tentative contigs and singlets) of the analysed sequences. The unigenes that showed similarity to proteins with known or predicted functions (50.3) were classified into 15 different functional categories. The functional profiles were found to be quite similar for all studied embryo stages but statistical analysis revealed that successive and coordinate sets of genes are expressed at each embryonic stage. The analysis allowed us to identify abundant and differentially expressed genes at the early stages of embryos development as well as some putatively interesting genes, showing strong similarities with genes playing key roles in plant and animal embryogenesis. The data presented in this study not only provide a first global overview of the genes expression profile during sunflower embryogenesis but also represent an original and valuable tool for developmental genomics studies on exalbuminous dicots.

Identification of genes preferentially expressed during wood formation in Eucalyptus

Wood is the most abundant biological resource on earth and it is also an important raw material for a major global industry with rapidly increasing demand. The genus Eucalyptus includes the most widely used tree species for industrial plantation, mainly for making pulp and paper. With the aim of identifying major genes involved in wood formation in Eucalyptus , we have developed a targeted approach of functional genomics based on the isolation of xylem preferentially expressed genes by subtractive PCR. Transcript profiling using cDNA arrays and analysis of variance (ANOVA) were used to identify differentially expressed ESTs between secondary xylem and leaves. Real-time RT-PCR was performed to confirm the differential expression of representative EST. Of 224 independent EST sequences obtained, 81% were preferentially expressed in xylem. One-third of the ESTs exhibiting homologies with proteins of known function fell into two main classes highlighting the importance of the auxin signalling through ubiquitin-dependent proteolysis on one hand, and of the enzymes involved in cell wall biosynthesis and remodelling, on the other. The functions of the genes represented by the remaining 61% of ESTs should be of great interest for future research. This systematic analysis of genes involved in wood formation in Eucalyptus provides valuable insights into the molecular mechanisms involved in secondary xylem differentiation as well as new candidate-genes for wood quality improvement.