Summary: 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.

The presence of lignin in secondary cell walls (SCW) is a major factor preventing hydrolytic enzymes from gaining access to cellulose, thereby limiting the saccharification potential of plant biomass. To understand how lignification is regulated is a prerequisite for selecting plant biomass better adapted to bioethanol production. Because transcriptional regulation is a major mechanism controlling the expression of genes involved in lignin biosynthesis, our aim was to identify novel transcription factors (TFs) dictating lignin profiles in the model plant Arabidopsis. To this end, we have developed a post-genomic approach by combining four independent in-house SCW-related transcriptome datasets obtained from (1) the fiber cell wall-deficient wat1 Arabidopsis mutant, (2) Arabidopsis lines over-expressing either the master regulatory activator EgMYB2 or (3) the repressor EgMYB1 and finally (4) Arabidopsis orthologs of Eucalyptus xylem-expressed genes. This allowed us to identify 502 up- or down-regulated TFs. We preferentially selected those present in more than one dataset and further analyzed their in silico expression patterns as an additional selection criteria. This selection process led to 80 candidates. Notably, 16 of them were already proven to regulate SCW formation, thereby validating the overall strategy. Then, we phenotyped 43 corresponding mutant lines focusing on histological observations of xylem and interfascicular fibers. This phenotypic screen revealed six mutant lines exhibiting altered lignification patterns. Two of them [Bel-like HomeoBox6 (blh6) and a zinc finger TF] presented hypolignified SCW. Three others (myb52, myb-like TF, hb5) showed hyperlignified SCW whereas the last one (hb15) showed ectopic lignification. In addition, our meta-analyses highlighted a reservoir of new potential regulators adding to the gene network regulating SCW but also opening new avenues to ultimately improve SCW composition for biofuel production.

Summary: 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.

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.

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. 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. 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.