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.

Eucalyptus grandis AUX/INDOLE-3-ACETIC ACID 13 (EgrIAA13) is a novel transcriptional regulator of xylogenesis.

KEY MESSAGE: Our Induced Somatic Sector Analysis and protein-protein interaction experiments demonstrate that Eucalyptus grandis IAA13 regulates xylem fibre and vessel development, potentially via EgrIAA13 modules involving ARF2, ARF5, ARF6 and ARF19. Auxin is a crucial phytohormone regulating multiple aspects of plant growth and differentiation, including regulation of vascular cambium activity, xylogenesis and its responsiveness towards gravitropic stress. Although the regulation of these biological processes greatly depends on auxin and regulators of the auxin signalling pathway, many of their specific functions remain unclear. Therefore, the present study aims to functionally characterise Eucalyptus grandis AUX/INDOLE-3-ACETIC ACID 13 (EgrIAA13), a member of the auxin signalling pathway. In Eucalyptus and Populus, EgrIAA13 and its orthologs are preferentially expressed in the xylogenic tissues and downregulated in tension wood. Therefore, to further investigate EgrIAA13 and its function during xylogenesis, we conducted subcellular localisation and Induced Somatic Sector Analysis experiments using overexpression and RNAi knockdown constructs of EgrIAA13 to create transgenic tissue sectors on growing stems of Eucalyptus and Populus. Since Aux/IAAs interact with Auxin Responsive Factors (ARFs), in silico predictions of IAA13-ARF interactions were explored and experimentally validated via yeast-2-hybrid experiments. Our results demonstrate that EgrIAA13 localises to the nucleus and that downregulation of EgrIAA13 impedes Eucalyptus xylem fibre and vessel development. We also observed that EgrIAA13 interacts with Eucalyptus ARF2, ARF5, ARF6 and ARF19A. Based on these results, we conclude that EgrIAA13 is a regulator of Eucalyptus xylogenesis and postulate that the observed phenotypes are likely to result from alterations in the auxin-responsive transcriptome via IAA13-ARF modules such as EgrIAA13-EgrARF5. Our results provide the first insights into the regulatory role of EgrIAA13 during xylogenesis.

Overexpression of EgrIAA20 from Eucalyptus grandis, a Non-Canonical Aux/IAA Gene, Specifically Decouples Lignification of the Different Cell-Types in Arabidopsis Secondary Xylem.

Wood (secondary xylem) formation is regulated by auxin, which plays a pivotal role as an integrator of developmental and environmental cues. However, our current knowledge of auxin-signaling during wood formation is incomplete. Our previous genome-wide analysis of Aux/IAAs in Eucalyptus grandis showed the presence of the non-canonical paralog member EgrIAA20 that is preferentially expressed in cambium. We analyzed its cellular localization using a GFP fusion protein and its transcriptional activity using transactivation assays, and demonstrated its nuclear localization and strong auxin response repressor activity. In addition, we functionally tested the role of EgrIAA20 by constitutive overexpression in Arabidopsis to investigate for phenotypic changes in secondary xylem formation. Transgenic Arabidopsis plants overexpressing EgrIAA20 were smaller and displayed impaired development of secondary fibers, but not of other wood cell types. The inhibition in fiber development specifically affected their cell wall lignification. We performed yeast-two-hybrid assays to identify EgrIAA20 protein partners during wood formation in Eucalyptus, and identified EgrIAA9A, whose ortholog PtoIAA9 in poplar is also known to be involved in wood formation. Altogether, we showed that EgrIAA20 is an important auxin signaling component specifically involved in controlling the lignification of wood fibers.

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.

Implementing the CRISPR/Cas9 Technology in Eucalyptus Hairy Roots Using Wood-Related Genes.

Eucalypts are the most planted hardwoods worldwide. The availability of the Eucalyptus grandis genome highlighted many genes awaiting functional characterization, lagging behind because of the lack of efficient genetic transformation protocols. In order to efficiently generate knock-out mutants to study the function of eucalypts genes, we implemented the powerful CRISPR/Cas9 gene editing technology with the hairy roots transformation system. As proofs-of-concept, we targeted two wood-related genes: Cinnamoyl-CoA Reductase1 (CCR1), a key lignin biosynthetic gene and IAA9A an auxin dependent transcription factor of Aux/IAA family. Almost all transgenic hairy roots were edited but the allele-editing rates and spectra varied greatly depending on the gene targeted. Most edition events generated truncated proteins, the prevalent edition types were small deletions but large deletions were also quite frequent. By using a combination of FT-IR spectroscopy and multivariate analysis (partial least square analysis (PLS-DA)), we showed that the CCR1-edited lines, which were clearly separated from the controls. The most discriminant wave-numbers were attributed to lignin. Histochemical analyses further confirmed the decreased lignification and the presence of collapsed vessels in CCR1-edited lines, which are characteristics of CCR1 deficiency. Although the efficiency of editing could be improved, the method described here is already a powerful tool to functionally characterize eucalypts genes for both basic research and industry purposes.

A systems biology view of wood formation in Eucalyptus grandis trees submitted to different potassium and water regimes.

Wood production in fast-growing Eucalyptus grandis trees is highly dependent on both potassium (K) fertilization and water availability but the molecular processes underlying wood formation in response to the combined effects of these two limiting factors remain unknown. E. grandis trees were submitted to four combinations of K-fertilization and water supply. Weighted gene co-expression network analysis and MixOmics-based co-regulation networks were used to integrate xylem transcriptome, metabolome and complex wood traits. Functional characterization of a candidate gene was performed in transgenic E. grandis hairy roots. This integrated network-based approach enabled us to identify meaningful biological processes and regulators impacted by K-fertilization and/or water limitation. It revealed that modules of co-regulated genes and metabolites strongly correlated to wood complex traits are in the heart of a complex trade-off between biomass production and stress responses. Nested in these modules, potential new cell-wall regulators were identified, as further confirmed by the functional characterization of EgMYB137. These findings provide new insights into the regulatory mechanisms of wood formation under stressful conditions, pointing out both known and new regulators co-opted by K-fertilization and/or water limitation that may potentially promote adaptive wood traits.

The control of endopolygalacturonase expression by the sugarcane RAV transcription factor during aerenchyma formation.

The development of lysigenous aerenchyma starts with cell expansion and degradation of pectin from the middle lamella, leading to cell wall modification, and culminating with cell separation. Here we report that nutritional starvation of sugarcane induced gene expression along sections of the first 5 cm of the root and between treatments. We selected two candidate genes: a RAV transcription factor, from the ethylene response factors superfamily, and an endopolygalacturonase (EPG), a glycosyl hydrolase related to homogalacturonan hydrolysis from the middle lamella. epg1 and rav1 transcriptional patterns suggest they are essential genes at the initial steps of pectin degradation during aerenchyma development in sugarcane. Due to the high complexity of the sugarcane genome, rav1 and epg1 were sequenced from 17 bacterial artificial chromosome clones containing hom(e)ologous genomic regions, and the sequences were compared with those of Sorghum bicolor. We used one hom(e)olog sequence from each gene for transactivation assays in tobacco. rav1 was shown to bind to the epg1 promoter, repressing ß-glucuronidase activity. RAV repression upon epg1 transcription is the first reported link between ethylene regulation and pectin hydrolysis during aerenchyma formation. Our findings may help to elucidate cell wall degradation in sugarcane and therefore contribute to second-generation bioethanol production.

The Eucalyptus linker histone variant EgH1.3 cooperates with the transcription factor EgMYB1 to control lignin biosynthesis during wood formation.

Wood, also called secondary xylem, is a specialized vascular tissue constituted by different cell types that undergo a differentiation process involving deposition of thick, lignified secondary cell walls. The mechanisms needed to control the extent of lignin deposition depending on the cell type and the differentiation stage are far from being fully understood. We found that the Eucalyptus transcription factor EgMYB1, which is known to repress lignin biosynthesis, interacts specifically with a linker histone variant, EgH1.3. This interaction enhances the repression of EgMYB1's target genes, strongly limiting the amount of lignin deposited in xylem cell walls. The expression profiles of EgMYB1 and EgH1.3 overlap in xylem cells at early stages of their differentiation as well as in mature parenchymatous xylem cells, which have no or only thin lignified secondary cell walls. This suggests that a complex between EgMYB1 and EgH1.3 integrates developmental signals to prevent premature or inappropriate lignification of secondary cell walls, providing a mechanism to fine-tune the differentiation of xylem cells in time and space. We also demonstrate a role for a linker histone variant in the regulation of a specific developmental process through interaction with a transcription factor, illustrating that plant linker histones have other functions beyond chromatin organization.

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.

Comprehensive genome-wide analysis of the Aux/IAA gene family in Eucalyptus: evidence for the role of EgrIAA4 in wood formation.

Auxin plays a pivotal role in various plant growth and development processes, including vascular differentiation. The modulation of auxin responsiveness through the auxin perception and signaling machinery is believed to be a major regulatory mechanism controlling cambium activity and wood formation. To gain more insights into the roles of key Aux/IAA gene regulators of the auxin response in these processes, we identified and characterized members of the Aux/IAA family in the genome of Eucalyptus grandis, a tree of worldwide economic importance. We found that the gene family in Eucalyptus is slightly smaller than that in Populus and Arabidopsis, but all phylogenetic groups are represented. High-throughput expression profiling of different organs and tissues highlighted several Aux/IAA genes expressed in vascular cambium and/or developing xylem, some showing differential expression in response to developmental (juvenile vs. mature) and/or to environmental (tension stress) cues. Based on the expression profiles, we selected a promising candidate gene, EgrIAA4, for functional characterization. We showed that EgrIAA4 protein is localized in the nucleus and functions as an auxin-responsive repressor. Overexpressing a stabilized version of EgrIAA4 in Arabidopsis dramatically impeded plant growth and fertility and induced auxin-insensitive phenotypes such as inhibition of primary root elongation, lateral root emergence and agravitropism. Interestingly, the lignified secondary walls of the interfascicular fibers appeared very late, whereas those of the xylary fibers were virtually undetectable, suggesting that EgrIAA4 may play crucial roles in fiber development and secondary cell wall deposition.

Structural, evolutionary and functional analysis of the NAC domain protein family in Eucalyptus.

NAC domain transcription factors regulate many developmental processes and stress responses in plants and vary widely in number and family structure. We analysed the characteristics and evolution of the NAC gene family of Eucalyptus grandis, a fast-growing forest tree in the rosid order Myrtales. NAC domain genes identified in the E. grandis genome were subjected to amino acid sequence, phylogenetic and motif analyses. Transcript abundance in developing tissues and abiotic stress conditions in E. grandis and E. globulus was quantified using RNA-seq and reverse transcription quantitative PCR (RT-qPCR). One hundred and eighty-nine E. grandis NAC (EgrNAC) proteins, arranged into 22 subfamilies, are extensively duplicated in subfamilies associated with stress response. Most EgrNAC genes form tandem duplicate arrays that frequently carry signatures of purifying selection. Sixteen amino acid motifs were identified in EgrNAC proteins, eight of which are enriched in, or unique to, Eucalyptus. New candidates for the regulation of normal and tension wood development and cold responses were identified. This first description of a Myrtales NAC domain family reveals an unique history of tandem duplication in stress-related subfamilies that has likely contributed to the adaptation of eucalypts to the challenging Australian environment. Several new candidates for the regulation of stress, wood formation and tree-specific development are reported.

Genome-wide analysis of the lignin toolbox of Eucalyptus grandis.

Lignin, a major component of secondary cell walls, hinders the optimal processing of wood for industrial uses. The recent availability of the Eucalyptus grandis genome sequence allows comprehensive analysis of the genes encoding the 11 protein families specific to the lignin branch of the phenylpropanoid pathway and identification of those mainly involved in xylem developmental lignification. We performed genome-wide identification of putative members of the lignin gene families, followed by comparative phylogenetic studies focusing on bona fide clades inferred from genes functionally characterized in other species. RNA-seq and microfluid real-time quantitative PCR (RT-qPCR) expression data were used to investigate the developmental and environmental responsive expression patterns of the genes. The phylogenetic analysis revealed that 38 E. grandis genes are located in bona fide lignification clades. Four multigene families (shikimate O-hydroxycinnamoyltransferase (HCT), p-coumarate 3-hydroxylase (C3H), caffeate/5-hydroxyferulate O-methyltransferase (COMT) and phenylalanine ammonia-lyase (PAL)) are expanded by tandem gene duplication compared with other plant species. Seventeen of the 38 genes exhibited strong, preferential expression in highly lignified tissues, probably representing the E. grandis core lignification toolbox. The identification of major genes involved in lignin biosynthesis in E. grandis, the most widely planted hardwood crop world-wide, provides the foundation for the development of biotechnology approaches to develop tree varieties with enhanced processing qualities.

The Eucalyptus grandis R2R3-MYB transcription factor family: evidence for woody growth-related evolution and function.

The R2R3-MYB family, one of the largest transcription factor families in higher plants, controls a wide variety of plant-specific processes including, notably, phenylpropanoid metabolism and secondary cell wall formation. We performed a genome-wide analysis of this superfamily in Eucalyptus, one of the most planted hardwood trees world-wide. A total of 141 predicted R2R3-MYB sequences identified in the Eucalyptus grandis genome sequence were subjected to comparative phylogenetic analyses with Arabidopsis thaliana, Oryza sativa, Populus trichocarpa and Vitis vinifera. We analysed features such as gene structure, conserved motifs and genome location. Transcript abundance patterns were assessed by RNAseq and validated by high-throughput quantitative PCR. We found some R2R3-MYB subgroups with expanded membership in E. grandis, V. vinifera and P. trichocarpa, and others preferentially found in woody species, suggesting diversification of specific functions in woody plants. By contrast, subgroups containing key genes regulating lignin biosynthesis and secondary cell wall formation are more conserved across all of the species analysed. In Eucalyptus, R2R3-MYB tandem gene duplications seem to disproportionately affect woody-preferential and woody-expanded subgroups. Interestingly, some of the genes belonging to woody-preferential subgroups show higher expression in the cambial region, suggesting a putative role in the regulation of secondary growth.

Contrasting nitrogen fertilization treatments impact xylem gene expression and secondary cell wall lignification in Eucalyptus.

BACKGROUND: Nitrogen (N) is a main nutrient required for tree growth and biomass accumulation. In this study, we analyzed the effects of contrasting nitrogen fertilization treatments on the phenotypes of fast growing Eucalyptus hybrids (E. urophylla x E. grandis) with a special focus on xylem secondary cell walls and global gene expression patterns. RESULTS: Histological observations of the xylem secondary cell walls further confirmed by chemical analyses showed that lignin was reduced by luxuriant fertilization, whereas a consistent lignin deposition was observed in trees grown in N-limiting conditions. Also, the syringyl/guaiacyl (S/G) ratio was significantly lower in luxuriant nitrogen samples. Deep sequencing RNAseq analyses allowed us to identify a high number of differentially expressed genes (1,469) between contrasting N treatments. This number is dramatically higher than those obtained in similar studies performed in poplar but using microarrays. Remarkably, all the genes involved the general phenylpropanoid metabolism and lignin pathway were found to be down-regulated in response to high N availability. These findings further confirmed by RT-qPCR are in agreement with the reduced amount of lignin in xylem secondary cell walls of these plants. CONCLUSIONS: This work enabled us to identify, at the whole genome level, xylem genes differentially regulated by N availability, some of which are involved in the environmental control of xylogenesis. It further illustrates that N fertilization can be used to alter the quantity and quality of lignocellulosic biomass in Eucalyptus, offering exciting prospects for the pulp and paper industry and for the use of short coppices plantations to produce second generation biofuels.

Targeted linkage map densification to improve cell wall related QTL detection and interpretation in maize.

Several QTLs for cell wall degradability and lignin content were previously detected in the F288 × F271 maize RIL progeny, including a set of major QTLs located in bin 6.06. Unexpectedly, allelic sequencing of genes located around the bin 6.06 QTL positions revealed a monomorphous region, suggesting that these QTLs were likely "ghost" QTLs. Refining the positions of all QTLs detected in this population was thus considered, based on a linkage map densification in most important QTL regions, and in several large still unmarked regions. Re-analysis of data with an improved genetic map (173 markers instead of 108) showed that ghost QTLs located in bin 6.06 were then fractionated over two QTL positions located upstream and downstream of the monomorphic region. The area located upstream of bin 6.06 position carried the major QTLs, which explained from 37 to 59 % of the phenotypic variation for per se values and extended on only 6 cM, corresponding to a physical distance of 2.2 Mbp. Among the 92 genes present in the corresponding area of the B73 maize reference genome, nine could putatively be considered as involved in the formation of the secondary cell wall [bHLH, FKBP, laccase, fasciclin, zinc finger C2H2-type and C3HC4-type (two genes), NF-YB, and WRKY]. In addition, based on the currently improved genetic map, eight QTLs were detected in bin 4.09, while only one QTL was highlighted in the initial investigation. Moreover, significant epistatic interaction effects were shown for all traits between these QTLs located in bin 4.09 and the major QTLs located in bin 6.05. Three genes related to secondary cell wall assembly (ZmMYB42, COV1-like, PAL-like) underlay QTL support intervals in this newly identified bin 4.09 region. The current investigations, even if they were based only on one RIL progeny, illustrated the interest of a targeted marker mapping on a genetic map to improve QTL position.

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.

Comprehensive genetic dissection of wood properties in a widely-grown tropical tree: Eucalyptus.

BACKGROUND: Eucalyptus is an important genus in industrial plantations throughout the world and is grown for use as timber, pulp, paper and charcoal. Several breeding programmes have been launched worldwide to concomitantly improve growth performance and wood properties (WPs). In this study, an interspecific cross between Eucalyptus urophylla and E. grandis was used to identify major genomic regions (Quantitative Trait Loci, QTL) controlling the variability of WPs. RESULTS: Linkage maps were generated for both parent species. A total of 117 QTLs were detected for a series of wood and end-use related traits, including chemical, technological, physical, mechanical and anatomical properties. The QTLs were mainly clustered into five linkage groups. In terms of distribution of QTL effects, our result agrees with the typical L-shape reported in most QTL studies, i.e. most WP QTLs had limited effects and only a few (13) had major effects (phenotypic variance explained > 15%). The co-locations of QTLs for different WPs as well as QTLs and candidate genes are discussed in terms of phenotypic correlations between traits, and of the function of the candidate genes. The major wood property QTL harbours a gene encoding a Cinnamoyl CoA reductase (CCR), a structural enzyme of the monolignol-specific biosynthesis pathway. CONCLUSIONS: Given the number of traits analysed, this study provides a comprehensive understanding of the genetic architecture of wood properties in this Eucalyptus full-sib pedigree. At the dawn of Eucalyptus genome sequence, it will provide a framework to identify the nature of genes underlying these important quantitative traits.

Advancing Eucalyptus genomics: identification and sequencing of lignin biosynthesis genes from deep-coverage BAC libraries.

BACKGROUND: Eucalyptus species are among the most planted hardwoods in the world because of their rapid growth, adaptability and valuable wood properties. The development and integration of genomic resources into breeding practice will be increasingly important in the decades to come. Bacterial artificial chromosome (BAC) libraries are key genomic tools that enable positional cloning of important traits, synteny evaluation, and the development of genome framework physical maps for genetic linkage and genome sequencing. RESULTS: We describe the construction and characterization of two deep-coverage BAC libraries EG_Ba and EG_Bb obtained from nuclear DNA fragments of E. grandis (clone BRASUZ1) digested with HindIII and BstYI, respectively. Genome coverages of 17 and 15 haploid genome equivalents were estimated for EG_Ba and EG_Bb, respectively. Both libraries contained large inserts, with average sizes ranging from 135 Kb (Eg_Bb) to 157 Kb (Eg_Ba), very low extra-nuclear genome contamination providing a probability of finding a single copy gene ≥ 99.99%. Libraries were screened for the presence of several genes of interest via hybridizations to high-density BAC filters followed by PCR validation. Five selected BAC clones were sequenced and assembled using the Roche GS FLX technology providing the whole sequence of the E. grandis chloroplast genome, and complete genomic sequences of important lignin biosynthesis genes. CONCLUSIONS: The two E. grandis BAC libraries described in this study represent an important milestone for the advancement of Eucalyptus genomics and forest tree research. These BAC resources have a highly redundant genome coverage (> 15×), contain large average inserts and have a very low percentage of clones with organellar DNA or empty vectors. These publicly available BAC libraries are thus suitable for a broad range of applications in genetic and genomic research in Eucalyptus and possibly in related species of Myrtaceae, including genome sequencing, gene isolation, functional and comparative genomics. Because they have been constructed using the same tree (E. grandis BRASUZ1) whose full genome is being sequenced, they should prove instrumental for assembly and gap filling of the upcoming Eucalyptus reference genome sequence.

Cryptogein, a fungal elicitor, remodels the phenylpropanoid metabolism of tobacco cell suspension cultures in a calcium-dependent manner.

Plant cells use calcium-based signalling pathways to transduce biotic and/or abiotic stimuli into adaptive responses. However, little is known about the coupling between calcium signalling, transcriptional regulation and the downstream biochemical processes. To understand these relationships better, we challenged tobacco BY-2 cells with cryptogein and evaluated how calcium transients (monitored through the calcium sensor aequorin) impact (1) transcript levels of phenylpropanoid genes (assessed by RT-qPCR); and (2) derived-phenolic compounds (analysed by mass spectrometry). Most genes of the phenylpropanoid pathway were up-regulated by cryptogein and cell wall-bound phenolic compounds accumulated (mainly 5-hydroxyferulic acid). The accumulation of both transcripts and phenolics was calcium-dependent. The transcriptional regulation of phenylpropanoid genes was correlated in a non-linear manner with stimulus intensity and with components of the cryptogein-induced calcium signature. In addition, calmodulin inhibitors increased the sensitivity of cells to low concentrations of cryptogein. These results led us to propose a model of coupling between the cryptogein signal, calcium signalling and the transcriptional response, exerting control of transcription through the coordinated action of two decoding modules exerting opposite effects.

Coordinated transcriptional regulation of two key genes in the lignin branch pathway--CAD and CCR--is mediated through MYB- binding sites.

BACKGROUND: Cinnamoyl CoA reductase (CCR) and cinnamyl alcohol dehydrogenase (CAD) catalyze the final steps in the biosynthesis of monolignols, the monomeric units of the phenolic lignin polymers which confer rigidity, imperviousness and resistance to biodegradation to cell walls. We have previously shown that the Eucalyptus gunnii CCR and CAD2 promoters direct similar expression patterns in vascular tissues suggesting that monolignol production is controlled, at least in part, by the coordinated transcriptional regulation of these two genes. Although consensus motifs for MYB transcription factors occur in most gene promoters of the whole phenylpropanoid pathway, functional evidence for their contribution to promoter activity has only been demonstrated for a few of them. Here, in the lignin-specific branch, we studied the functional role of MYB elements as well as other cis-elements identified in the regulatory regions of EgCAD2 and EgCCR promoters, in the transcriptional activity of these gene promoters. RESULTS: By using promoter deletion analysis and in vivo footprinting, we identified an 80 bp regulatory region in the Eucalyptus gunnii EgCAD2 promoter that contains two MYB elements, each arranged in a distinct module with newly identified cis-elements. A directed mutagenesis approach was used to introduce block mutations in all putative cis-elements of the EgCAD2 promoter and in those of the 50 bp regulatory region previously delineated in the EgCCR promoter. We showed that the conserved MYB elements in EgCAD2 and EgCCR promoters are crucial both for the formation of DNA-protein complexes in EMSA experiments and for the transcriptional activation of EgCAD2 and EgCCR promoters in vascular tissues in planta. In addition, a new regulatory cis-element that modulates the balance between two DNA-protein complexes in vitro was found to be important for EgCAD2 expression in the cambial zone. CONCLUSIONS: Our assignment of functional roles to the identified cis-elements clearly demonstrates the importance of MYB cis-elements in the transcriptional regulation of two genes of the lignin-specific pathway and support the hypothesis that MYB elements serve as a common means for the coordinated regulation of genes in the entire lignin biosynthetic pathway.