Cellulose factories: advancing bioenergy production from forest trees

Comprehensive analysis of lncRNA-mRNA regulatory network in Populus associated with xylem development

Long noncoding RNAs (lncRNAs) play essential roles in numerous biological processes in plants, such as regulating the gene expression. However, only a few studies have looked into their potential functions in xylem development. High-throughput sequencing of P. euramericana 'Zhonglin46' developing and mature xylem was performed in this study. Through sequencing analysis, 14,028 putative lncRNA transcripts were identified, including 4525 differentially expressed lncRNAs (DELs). Additional research revealed that in mature xylem, a total of 2320 DELs were upregulated and 2205 were downregulated compared to developing xylem. Meanwhile, there were a total of 8122 differentially expressed mRNAs (DEMs) that were upregulated and 16,424 that were downregulated in mature xylem compared with developing xylem. The cis- and trans-target genes of DELs were analyzed for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment, which indicated that these DELs participate in controlling the phenylpropanoid and lignin biosynthesis pathway as well as the starch and sucrose metabolism pathway. Among the cis-regulated DELs, LNC_006291, LNC_006292, and LNC_006532 all participate in regulating multiple HCT gene family membranes. As targets, POPTR_001G045900v3 (CCR2) and POPTR_018G063500v3 (SUS) both have only one cis-regulatory lncRNA, referred to as LNC_000057 and LNC_006212, respectively. Moreover, LNC_004484 and two DELs named LNC_008014 and LNC_010781 were revealed to be important nodes in the co-expression network of trans-lncRNAs and mRNAs associated to the lignin biosynthesis pathway and cellulose and xylan biosynthetic pathways, respectively. Finally, quantitative real-time PCR (qRT-PCR) was used to confirme 34 pairs of lncRNA-mRNA. Taken together, these findings may help to clarify the regulatory role that lncRNAs play in xylem development and wood formation.

Ginseng-derived patatin-related phospholipase PgpPLAIIIβ alters plant growth and lignification of xylem in hybrid poplars

Patatin-liked phospholipase A (pPLAs) are major lipid acyl hydrolases that participate in various biological functions in plant growth and development. Previously, a ginseng-derived pPLAIII homolog was reported to reduce lignin content in Arabidopsis. This led us to evaluate its possible usefulness as a biomass source in wood plant. Herein, we report that there are six members in the pPLAIII gene family in poplar. Overexpression of pPLAIIIβ derived from ginseng resulted in a reduced plant height with radially expanded stem growth in hybrid poplars. Compared with the wild type (WT), the chlorophyll content was increased in the overexpression poplar lines, whereas the leaf size was smaller. The secondary cell wall structure in overexpression lines was also altered, exhibiting reduced lignification in the xylem. Two transcription factors, MYB92 and MYB152, which control lignin biosynthesis, were downregulated in the overexpression lines. The middle xylem of the overexpression line showed heavy thickening, making it thicker than the other xylem parts and the WT xylem, which rather could have been contributed by the presence of more cellulose in the selected surface area. Taken together, the results suggest that PgpPLAIIIβ plays a role not only in cell elongation patterns, but also in determining the secondary cell wall composition.

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.

Xylan in the Middle: Understanding Xylan Biosynthesis and Its Metabolic Dependencies Toward Improving Wood Fiber for Industrial Processing

Lignocellulosic biomass, encompassing cellulose, lignin and hemicellulose in plant secondary cell walls (SCWs), is the most abundant source of renewable materials on earth. Currently, fast-growing woody dicots such as Eucalyptus and Populus trees are major lignocellulosic (wood fiber) feedstocks for bioproducts such as pulp, paper, cellulose, textiles, bioplastics and other biomaterials. Processing wood for these products entails separating the biomass into its three main components as efficiently as possible without compromising yield. Glucuronoxylan (xylan), the main hemicellulose present in the SCWs of hardwood trees carries chemical modifications that are associated with SCW composition and ultrastructure, and affect the recalcitrance of woody biomass to industrial processing. In this review we highlight the importance of xylan properties for industrial wood fiber processing and how gaining a greater understanding of xylan biosynthesis, specifically xylan modification, could yield novel biotechnology approaches to reduce recalcitrance or introduce novel processing traits. Altering xylan modification patterns has recently become a focus of plant SCW studies due to early findings that altered modification patterns can yield beneficial biomass processing traits. Additionally, it has been noted that plants with altered xylan composition display metabolic differences linked to changes in precursor usage. We explore the possibility of using systems biology and systems genetics approaches to gain insight into the coordination of SCW formation with other interdependent biological processes. Acetyl-CoA, s-adenosylmethionine and nucleotide sugars are precursors needed for xylan modification, however, the pathways which produce metabolic pools during different stages of fiber cell wall formation still have to be identified and their co-regulation during SCW formation elucidated. The crucial dependence on precursor metabolism provides an opportunity to alter xylan modification patterns through metabolic engineering of one or more of these interdependent pathways. The complexity of xylan biosynthesis and modification is currently a stumbling point, but it may provide new avenues for woody biomass engineering that are not possible for other biopolymers.

A New Calmodulin-Binding Protein Expresses in the Context of Secondary Cell Wall Biosynthesis and Impacts Biomass Properties in Populus

A greater understanding of biosynthesis, signaling and regulatory pathways involved in determining stem growth and secondary cell wall chemistry is important for enabling pathway engineering and genetic optimization of biomass properties. The present study describes a new functional role of PdIQD10, a Populus gene belonging to the IQ67-Domain1 family of IQD genes, in impacting biomass formation and chemistry. Expression studies showed that PdIQD10 has enhanced expression in developing xylem and tension-stressed tissues in Populus deltoides. Molecular dynamics simulation and yeast two-hybrid interaction experiments suggest interactions with two calmodulin proteins, CaM247 and CaM014, supporting the sequence-predicted functional role of the PdIQD10 as a calmodulin-binding protein. PdIQD10 was found to interact with specific Populus isoforms of the Kinesin Light Chain protein family, shown previously to function as microtubule-guided, cargo binding and delivery proteins in Arabidopsis. Subcellular localization studies showed that PdIQD10 localizes in the nucleus and plasma membrane regions. Promoter-binding assays suggest that a known master transcriptional regulator of secondary cell wall biosynthesis (PdWND1B) may be upstream of an HD-ZIP III gene that is in turn upstream of PdIQD10 gene in the transcriptional network. RNAi-mediated downregulation of PdIQD10 expression resulted in plants with altered biomass properties including higher cellulose, wall glucose content and greater biomass quantity. These results present evidence in support of a new functional role for an IQD gene family member, PdIQD10, in secondary cell wall biosynthesis and biomass formation in Populus.

Canopy and understory nitrogen addition increase the xylem tracheid size of dominant broadleaf species in a subtropical forest of China

Tree xylem anatomy is associated with carbon accumulation and wood quality. Increasing nitrogen (N) deposition can cause a significant effect on xylem anatomy, but related information is limited for subtropical broadleaf tree species. A 3-year field N addition experiment, with different N addition approaches (canopy and understory) and addition rates (0, 25, and 50 kg N ha^-1 yr^-1), was performed beginning in 2013 in a subtropical forest of China. N addition effects on xylem tracheid (wall and lumen) size, vessel, and growth of dominant broadleaf species (Schima superba Gardn. et Champ. and Castanopsis chinensis (Sprengel) Hance) were investigated. The results showed that The effect of N addition on tracheid size was dependent on the tree species and addition approaches. Canopy N addition did not affect the tracheid size of C. chinensis, while both addition approaches increased the tracheid size of S. superba. The vessel size of both species was not affected by N addition. There was no difference in radial growth or other growth-related variables between the control and N-treated trees. These findings indicated that short-term N addition could significantly affect xylem anatomy, but might not influence tree growth. Meanwhile, understory N addition may pose challenges for mechanistic understanding and forest dynamics projection.

Functional Analysis of Cellulose Synthase CesA4 and CesA6 Genes in Switchgrass (Panicum virgatum) by Overexpression and RNAi-Mediated Gene Silencing

Switchgrass (Panicum virgatum L.) is a leading lignocellulosic bioenergy feedstock. Cellulose is a major component of the plant cell walls and the primary substrate for saccharification. Accessibility of cellulose to enzymatic breakdown into fermentable sugars is limited by the presence of lignin in the plant cell wall. In this study, putatively novel switchgrass secondary cell wall cellulose synthase PvCesA4 and primary cell wall PvCesA6 genes were identified and their functional role in cellulose synthesis and cell wall composition was examined by overexpression and knockdown of the individual genes in switchgrass. The endogenous expression of PvCesA4 and PvCesA6 genes varied among including roots, leaves, stem, and reproductive tissues. Increasing or decreasing PvCesA4 and PvCesA6 expression to extreme levels in the transgenic lines resulted in decreased biomass production. PvCesA6-overexpressing lines had reduced lignin content and syringyl/guaiacyl lignin monomer ratio accompanied by increased sugar release efficiency, suggesting an impact of PvCesA6 expression levels on lignin biosynthesis. Cellulose content and cellulose crystallinity were decreased, while xylan content was increased in PvCesA4 and PvCesA6 overexpression or knockdown lines. The increase in xylan content suggests that the amount of non-cellulosic cell wall polysaccharide was modified in these plants. Taken together, the results show that the manipulation of the cellulose synthase genes alters the cell wall composition and availability of cellulose as a bioprocessing substrate.

Characterization of the cork formation and production transcriptome in Quercus cerris × suber hybrids

Cork oak is the main cork-producing species worldwide, and plays a significant economic, ecological and social role in the Mediterranean countries, in particular in Portugal and Spain. The ability to produce cork is limited to a few species, hence it must involve specific regulation mechanisms that are unique to these species. However, to date, these mechanisms remain largely understudied, especially with approaches involving the use of high-throughput sequencing technology. In this study, the transcriptome of cork-producing and non-cork-producing Quercus cerris × suber hybrids was analyzed in order to elucidate the differences between the two groups of trees displaying contrasting phenotypes for cork production. The results revealed the presence of a significant number of genes exclusively associated with cork production, in the trees that developed cork. Moreover, several gene ontology subcategories, such as cell wall biogenesis, lipid metabolic processes, metal ion binding and apoplast/cell wall, were only detected in the trees with cork production. These results indicate the existence, at the transcriptome level, of mechanisms that seem to be unique and necessary for cork production, which is an advancement in our knowledge regarding the genetic regulation behind cork formation and production.

Microbial production of rhamnolipids using sugars as carbon sources

Rhamnolipids are a class of biosurfactants with effective surface-active properties. The high cost of microbial production of rhamnolipids largely affects their commercial applications. To reduce the production post, research has been carried out in screening more powerful strains, engineering microbes with higher biosurfactant yields and exploring cheaper substrates to reduce the production cost. Extensive refining is required for biosurfactant production using oils and oil-containing wastes, necessitating the use of complex and expensive biosurfactant recovery methods such as extraction with solvents or acid precipitation. As raw materials normally can account for 10-30% of the overall production cost, sugars have been proven to be an alternative carbon source for microbial production of rhamnolipids due to its lower costs and straightforward processing techniques. Studies have thus been focused on using tropical agroindustrial crop residues as renewable substrates. Herein, we reviewed studies that are using sugar-containing substrates as carbon sources for producing rhamnolipids. We speculate that sugars derived from agricultural wastes rich in cellulose and sugar-containing wastes are potential carbon sources in fermentation while challenges still remain in large scales.

Mechanocatalytic Solvent-Free Esterification of Sugarcane Bagasse

Esterification is a versatile way to produce the derivatives of lignocellulose with developed properties. However, the traditional heterogeneous esterification of lignocellulose suffered from the drawbacks of low efficiency, additional reaction medium and heating. In the present study, an efficient method was developed to produce the functionalized sugarcane bagasse (SCB) by ball milling without any additional solvents and heating. The effects of pulverization time, rotation speed, the kind of linear chain anhydrides, the ratio of anhydrides to SCB, with or without pyridine catalyst and the dosage of catalyst were investigated on weight percent gain (WPG) of SCB esters. The results indicated that the high efficiency of this mechanocatalystic esterification was probably due to the destroyed crystalline structure and the promoted penetration of the esterifying reagent onto SCB bulk caused by ball milling. The maximum WPG of SCB acetate, propionate and butyrate reached 33.3%, 33.6% and 32.4%, respectively. The physicochemical structure of the esterified SCB was characterized with Fourier transform infrared spectra (FT-IR), solid state cross-polarized magic angle spinning 13C nuclear magnetic resonance (CP/MAS 13C-NMR), X-ray diffraction (XRD), scanning electron microscopy (SEM) and thermogravimetric analysis (TGA). The direct evidence of the esterification occurrence was provided with FT-IR and solid-state CP/MAS 13C-NMR. The thermal stability of SCB increased upon the mechanocatalytic esterification. The results implied that the relatively homogeneous modification was achieved with this semi-homogeneous esterification method by ball milling.

Working towards recalcitrance mechanisms: increased xylan and homogalacturonan production by overexpression of GAlactUronosylTransferase12 (GAUT12) causes increased recalcitrance and decreased growth in Populus

BACKGROUND

The development of fast-growing hardwood trees as a source of lignocellulosic biomass for biofuel and biomaterial production requires a thorough understanding of the plant cell wall structure and function that underlie the inherent recalcitrance properties of woody biomass. Downregulation of GAUT12.1 in Populus deltoides was recently reported to result in improved biomass saccharification, plant growth, and biomass yield. To further understand GAUT12.1 function in biomass recalcitrance and plant growth, here we report the effects of P. trichocarpa GAUT12.1 overexpression in P. deltoides.

RESULTS

Increasing GAUT12.1 transcript expression by 7-49% in P. deltoides PtGAUT12.1-overexpression (OE) lines resulted in a nearly complete opposite biomass saccharification and plant growth phenotype to that observed previously in PdGAUT12.1-knockdown (KD) lines. This included significantly reduced glucose, xylose, and total sugar release (12-13%), plant height (6-54%), stem diameter (8-40%), and overall total aerial biomass yield (48-61%) in 3-month-old, greenhouse-grown PtGAUT12.1-OE lines compared to controls. Total lignin content was unaffected by the gene overexpression. Importantly, selected PtGAUT12.1-OE lines retained the recalcitrance and growth phenotypes upon growth for 9 months in the greenhouse and 2.8 years in the field. PtGAUT12.1-OE plants had significantly smaller leaves with lower relative water content, and significantly reduced stem wood xylem cell numbers and size. At the cell wall level, xylose and galacturonic acid contents increased markedly in total cell walls as well as in soluble and insoluble cell wall extracts, consistent with increased amounts of xylan and homogalacturonan in the PtGAUT12.1-OE lines. This led to increased cell wall recalcitrance, as manifested by the 9-15% reduced amounts of recovered extractable wall materials and 8-15% greater amounts of final insoluble pellet in the PtGAUT12.1-OE lines compared to controls.

CONCLUSIONS

The combined phenotype and chemotype data from P. deltoides PtGAUT12.1-OE and PdGAUT12.1-KD transgenics clearly establish GAUT12.1 as a recalcitrance- and growth-associated gene in poplar. Overall, the data support the hypothesis that GAUT12.1 synthesizes either an HG-containing primer for xylan synthesis or an HG glycan required for proper xylan deposition, anchoring, and/or architecture in the wall, and the possibility of HG and xylan glycans being connected to each other by a base-sensitive covalent linkage.

Network-based integration of systems genetics data reveals pathways associated with lignocellulosic biomass accumulation and processing

As a consequence of their remarkable adaptability, fast growth, and superior wood properties, eucalypt tree plantations have emerged as key renewable feedstocks (over 20 million ha globally) for the production of pulp, paper, bioenergy, and other lignocellulosic products. However, most biomass properties such as growth, wood density, and wood chemistry are complex traits that are hard to improve in long-lived perennials. Systems genetics, a process of harnessing multiple levels of component trait information (e.g., transcript, protein, and metabolite variation) in populations that vary in complex traits, has proven effective for dissecting the genetics and biology of such traits. We have applied a network-based data integration (NBDI) method for a systems-level analysis of genes, processes and pathways underlying biomass and bioenergy-related traits using a segregating Eucalyptus hybrid population. We show that the integrative approach can link biologically meaningful sets of genes to complex traits and at the same time reveal the molecular basis of trait variation. Gene sets identified for related woody biomass traits were found to share regulatory loci, cluster in network neighborhoods, and exhibit enrichment for molecular functions such as xylan metabolism and cell wall development. These findings offer a framework for identifying the molecular underpinnings of complex biomass and bioprocessing-related traits. A more thorough understanding of the molecular basis of plant biomass traits should provide additional opportunities for the establishment of a sustainable bio-based economy.

Evaluation of Relationships between Growth Rate, Tree Size, Lignocellulose Composition, and Enzymatic Saccharification in Interspecific Corymbia Hybrids and Parental Taxa

In order for a lignocellulosic bioenergy feedstock to be considered sustainable, it must possess a high rate of growth to supply biomass for conversion. Despite the desirability of a fast growth rate for industrial application, it is unclear what effect growth rate has on biomass composition or saccharification. We characterized Klason lignin, glucan, and xylan content with response to growth in Corymbia interspecific F1 hybrid families (HF) and parental species Corymbia torelliana and C. citriodora subspecies variegata and measured the effects on enzymatic hydrolysis from hydrothermally pretreated biomass. Analysis of biomass composition within Corymbia populations found similar amounts of Klason lignin content (19.7-21.3%) among parental and hybrid populations, whereas glucan content was clearly distinguished within C. citriodora subspecies variegata (52%) and HF148 (60%) as compared to other populations (28-38%). Multiple linear regression indicates that biomass composition is significantly impacted by tree size measured at the same age, with Klason lignin content increasing with diameter breast height (DBH) (+0.12% per cm DBH increase), and glucan and xylan typically decreasing per DBH cm increase (-0.7 and -0.3%, respectively). Polysaccharide content within C. citriodora subspecies variegata and HF-148 were not significantly affected by tree size. High-throughput enzymatic saccharification of hydrothermally pretreated biomass found significant differences among Corymbia populations for total glucose production from biomass, with parental Corymbia torelliana and hybrids HF-148 and HF-51 generating the highest amounts of glucose (~180 mg/g biomass, respectively), with HF-51 undergoing the most efficient glucan-to-glucose conversion (74%). Based on growth rate, biomass composition, and further optimization of enzymatic saccharification yield, high production Corymbia hybrid trees are potentially suitable for fast-rotation bioenergy or biomaterial production.

Systems genetics of wood formation

In woody plants, xylogenesis is an exceptionally strong carbon sink requiring robust transcriptional control and dynamic coordination of cellular and metabolic processes directing carbon allocation and partitioning into secondary cell wall biosynthesis. As a biological process, wood formation is an excellent candidate for systems modeling due to the strong correlation patterns and interconnectedness observed for transcriptional and metabolic component traits contributing to complex phenotypes such as cell wall chemistry and ultrastructure. Genetic variation in undomesticated tree populations provides abundant perturbation of systems components, adding another dimension to plant systems biology (besides spatial and temporal variation). High-throughput analysis of molecular component traits in adult trees has provided the first insights into the systems genetics of wood, an important renewable feedstock for biomaterials and bioenergy.

Genetic architecture of wood properties based on association analysis and co-expression networks in white spruce

Association studies are widely utilized to analyze complex traits but their ability to disclose genetic architectures is often limited by statistical constraints, and functional insights are usually minimal in nonmodel organisms like forest trees. We developed an approach to integrate association mapping results with co-expression networks. We tested single nucleotide polymorphisms (SNPs) in 2652 candidate genes for statistical associations with wood density, stiffness, microfibril angle and ring width in a population of 1694 white spruce trees (Picea glauca). Associations mapping identified 229-292 genes per wood trait using a statistical significance level of P < 0.05 to maximize discovery. Over-representation of genes associated for nearly all traits was found in a xylem preferential co-expression group developed in independent experiments. A xylem co-expression network was reconstructed with 180 wood associated genes and several known MYB and NAC regulators were identified as network hubs. The network revealed a link between the gene PgNAC8, wood stiffness and microfibril angle, as well as considerable within-season variation for both genetic control of wood traits and gene expression. Trait associations were distributed throughout the network suggesting complex interactions and pleiotropic effects. Our findings indicate that integration of association mapping and co-expression networks enhances our understanding of complex wood traits.

Free-Flow Electrophoresis of Plasma Membrane Vesicles Enriched by Two-Phase Partitioning Enhances the Quality of the Proteome from Arabidopsis Seedlings

The plant plasma membrane is the interface between the cell and its environment undertaking a range of important functions related to transport, signaling, cell wall biosynthesis, and secretion. Multiple proteomic studies have attempted to capture the diversity of proteins in the plasma membrane using biochemical fractionation techniques. In this study, two-phase partitioning was combined with free-flow electrophoresis to produce a population of highly purified plasma membrane vesicles that were subsequently characterized by tandem mass spectroscopy. This combined high-quality plasma membrane isolation technique produced a reproducible proteomic library of over 1000 proteins with an extended dynamic range including plasma membrane-associated proteins. The approach enabled the detection of a number of putative plasma membrane proteins not previously identified by other studies, including peripheral membrane proteins. Utilizing multiple data sources, we developed a PM-confidence score to provide a value indicating association to the plasma membrane. This study highlights over 700 proteins that, while seemingly abundant at the plasma membrane, are mostly unstudied. To validate this data set, we selected 14 candidates and transiently localized 13 to the plasma membrane using a fluorescent tag. Given the importance of the plasma membrane, this data set provides a valuable tool to further investigate important proteins. The mass spectrometry data are available via ProteomeXchange, identifier PXD001795.

Modular organization of the white spruce (Picea glauca) transcriptome reveals functional organization and evolutionary signatures

Transcript profiling has shown the molecular bases of several biological processes in plants but few studies have developed an understanding of overall transcriptome variation. We investigated transcriptome structure in white spruce (Picea glauca), aiming to delineate its modular organization and associated functional and evolutionary attributes. Microarray analyses were used to: identify and functionally characterize groups of co-expressed genes; investigate expressional and functional diversity of vascular tissue preferential genes which were conserved among Picea species, and identify expression networks underlying wood formation. We classified 22 857 genes as variable (79%; 22 coexpression groups) or invariant (21%) by profiling across several vegetative tissues. Modular organization and complex transcriptome restructuring among vascular tissue preferential genes was revealed by their assignment to coexpression groups with partially overlapping profiles and partially distinct functions. Integrated analyses of tissue-based and temporally variable profiles identified secondary xylem gene networks, showed their remodelling over a growing season and identified PgNAC-7 (no apical meristerm (NAM), Arabidopsis transcription activation factor (ATAF) and cup-shaped cotyledon (CUC) transcription factor 007 in Picea glauca) as a major hub gene specific to earlywood formation. Reference profiling identified comprehensive, statistically robust coexpressed groups, revealing that modular organization underpins the evolutionary conservation of the transcriptome structure.

Comparative interrogation of the developing xylem transcriptomes of two wood-forming species: Populus trichocarpa and Eucalyptus grandis

Wood formation is a complex developmental process governed by genetic and environmental stimuli. Populus and Eucalyptus are fast-growing, high-yielding tree genera that represent ecologically and economically important species suitable for generating significant lignocellulosic biomass. Comparative analysis of the developing xylem and leaf transcriptomes of Populus trichocarpa and Eucalyptus grandis together with phylogenetic analyses identified clusters of homologous genes preferentially expressed during xylem formation in both species. A conserved set of 336 single gene pairs showed highly similar xylem preferential expression patterns, as well as evidence of high functional constraint. Individual members of multi-gene orthologous clusters known to be involved in secondary cell wall biosynthesis also showed conserved xylem expression profiles. However, species-specific expression as well as opposite (xylem versus leaf) expression patterns observed for a subset of genes suggest subtle differences in the transcriptional regulation important for xylem development in each species. Using sequence similarity and gene expression status, we identified functional homologs likely to be involved in xylem developmental and biosynthetic processes in Populus and Eucalyptus. Our study suggests that, while genes involved in secondary cell wall biosynthesis show high levels of gene expression conservation, differential regulation of some xylem development genes may give rise to unique xylem properties.

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.

Plant lignin content altered by soil microbial community

Questions have been raised in various fields of research about the consequences of plants with modified lignin production. As a result of their roles in nutrient cycling and plant diversity, plant-soil interactions should be a major focus of ecological studies on lignin-modified plants. However, most studies have been decomposition studies conducted in a single soil or in sterile soil. Thus, we understand little about plant-soil interactions in living lignin-modified plants. In lignin mutants of three different barley (Hordeum vulgare) cultivars and their corresponding wild-types associated with three different soil microbial communities, we asked: do plant-soil microbiome interactions influence the lignin content of plants?; does a mutation in lignin production alter the outcome of plant-soil microbiome interactions?; does the outcome of plant-soil microbiome interactions depend on host genotype or the presence of a mutation altering lignin production? In roots, the soil community explained 6% of the variation in lignin content, but, in shoots, the soil community explained 21% of the variation in lignin content and was the only factor influencing lignin content. Neither genotype nor mutations in lignin production explained associations with fungi. Lignin content changes in response to a plant's soil microbial community, and may be a defensive response to particular components of the soil community.

Systems and synthetic biology approaches to alter plant cell walls and reduce biomass recalcitrance

Fine-tuning plant cell wall properties to render plant biomass more amenable to biofuel conversion is a colossal challenge. A deep knowledge of the biosynthesis and regulation of plant cell wall and a high-precision genome engineering toolset are the two essential pillars of efforts to alter plant cell walls and reduce biomass recalcitrance. The past decade has seen a meteoric rise in use of transcriptomics and high-resolution imaging methods resulting in fresh insights into composition, structure, formation and deconstruction of plant cell walls. Subsequent gene manipulation approaches, however, commonly include ubiquitous mis-expression of a single candidate gene in a host that carries an intact copy of the native gene. The challenges posed by pleiotropic and unintended changes resulting from such an approach are moving the field towards synthetic biology approaches. Synthetic biology builds on a systems biology knowledge base and leverages high-precision tools for high-throughput assembly of multigene constructs and pathways, precision genome editing and site-specific gene stacking, silencing and/or removal. Here, we summarize the recent breakthroughs in biosynthesis and remodelling of major secondary cell wall components, assess the impediments in obtaining a systems-level understanding and explore the potential opportunities in leveraging synthetic biology approaches to reduce biomass recalcitrance.

Flavonoid supplementation affects the expression of genes involved in cell wall formation and lignification metabolism and increases sugar content and saccharification in the fast-growing eucalyptus hybrid E. urophylla x E. grandis

BACKGROUND

Eucalyptus species are the most widely planted hardwood species in the world and are renowned for their rapid growth and adaptability. In Brazil, one of the most widely grown Eucalyptus cultivars is the fast-growing Eucalyptus urophylla x Eucalyptus grandis hybrid. In a previous study, we described a chemical characterization of these hybrids when subjected to flavonoid supplementation on 2 distinct timetables, and our results revealed marked differences between the wood composition of the treated and untreated trees.

RESULTS

In this work, we report the transcriptional responses occurring in these trees that may be related to the observed chemical differences. Gene expression was analysed through mRNA-sequencing, and notably, compared to control trees, the treated trees display differential down-regulation of cell wall formation pathways such as phenylpropanoid metabolism as well as differential expression of genes involved in sucrose, starch and minor CHO metabolism and genes that play a role in several stress and environmental responses. We also performed enzymatic hydrolysis of wood samples from the different treatments, and the results indicated higher sugar contents and glucose yields in the flavonoid-treated plants.

CONCLUSIONS

Our results further illustrate the potential use of flavonoids as a nutritional complement for modifying Eucalyptus wood, since, supplementation with flavonoids alters its chemical composition, gene expression and increases saccharification probably as part of a stress response.

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.

Conservation and functional influence of alternative splicing in wood formation of Populus and Eucalyptus

BACKGROUND

Wood formation in tree species is regulated by multiple factors at various layers. Alternative splicing (AS) occurs within a large number of genes in wood formation. However, the functional implications and conservation of the AS occurrence are not well understood.

RESULTS

In this study, we profiled AS events in wood-forming tissues of Populus and Eucalyptus, and analyzed their functional implications as well as inter-species conservation. 28.3% and 20.7% of highly expressed transcripts in the developing xylem of Populus and Eucalyptus respectively were affected by AS events. Around 42% of the AS events resulted in changes to the original reading frame. 25.0% (in Populus) and 26.8% (in Eucalyptus) of the AS events may cause protein domain modification. In the process of wood formation, about 28% of AS-occurring genes were putative orthologs and 71 conserved AS events were identified in the two species.

CONCLUSION

Through analysis of AS events in developing xylem of two tree species, this study reveals an array of new information regarding AS occurrence and function in tree development.

Uncovering the defence responses of Eucalyptus to pests and pathogens in the genomics age

Long-lived tree species are subject to attack by various pests and pathogens during their lifetime. This problem is exacerbated by climate change, which may increase the host range for pathogens and extend the period of infestation by pests. Plant defences may involve preformed barriers or induced resistance mechanisms based on recognition of the invader, complex signalling cascades, hormone signalling, activation of transcription factors and production of pathogenesis-related (PR) proteins with direct antimicrobial or anti-insect activity. Trees have evolved some unique defence mechanisms compared with well-studied model plants, which are mostly herbaceous annuals. The genome sequence of Eucalyptus grandis W. Hill ex Maiden has recently become available and provides a resource to extend our understanding of defence in large woody perennials. This review synthesizes existing knowledge of defence mechanisms in model plants and tree species and features mechanisms that may be important for defence in Eucalyptus, such as anatomical variants and the role of chemicals and proteins. Based on the E. grandis genome sequence, we have identified putative PR proteins based on sequence identity to the previously described plant PR proteins. Putative orthologues for PR-1, PR-2, PR-4, PR-5, PR-6, PR-7, PR-8, PR-9, PR-10, PR-12, PR-14, PR-15 and PR-17 have been identified and compared with their orthologues in Populus trichocarpa Torr. & A. Gray ex Hook and Arabidopsis thaliana (L.) Heynh. The survey of PR genes in Eucalyptus provides a first step in identifying defence gene targets that may be employed for protection of the species in future. Genomic resources available for Eucalyptus are discussed and approaches for improving resistance in these hardwood trees, earmarked as a bioenergy source in future, are considered.

Cell wall-related proteins of unknown function: missing links in plant cell wall development

Lignocellulosic biomass is an important feedstock for the pulp and paper industry as well as emerging biofuel and biomaterial industries. However, the recalcitrance of the secondary cell wall to chemical or enzymatic degradation remains a major hurdle for efficient extraction of economically important biopolymers such as cellulose. It has been estimated that approximately 10-15% of about 27,000 protein-coding genes in the Arabidopsis genome are dedicated to cell wall development; however, only about 130 Arabidopsis genes thus far have experimental evidence validating cell wall function. While many genes have been implicated through co-expression analysis with known genes, a large number are broadly classified as proteins of unknown function (PUFs). Recently the functionality of some of these unknown proteins in cell wall development has been revealed using reverse genetic approaches. Given the large number of cell wall-related PUFs, how do we approach and subsequently prioritize the investigation of such unknown genes that may be essential to or influence plant cell wall development and structure? Here, we address the aforementioned question in two parts; we first identify the different kinds of PUFs based on known and predicted features such as protein domains. Knowledge of inherent features of PUFs may allow for functional inference and a concomitant link to biological context. Secondly, we discuss omics-based technologies and approaches that are helping identify and prioritize cell wall-related PUFs by functional association. In this way, hypothesis-driven experiments can be designed for functional elucidation of many proteins that remain missing links in our understanding of plant cell wall biosynthesis.

Navigating the transcriptional roadmap regulating plant secondary cell wall deposition

The current status of lignocellulosic biomass as an invaluable resource in industry, agriculture, and health has spurred increased interest in understanding the transcriptional regulation of secondary cell wall (SCW) biosynthesis. The last decade of research has revealed an extensive network of NAC, MYB and other families of transcription factors regulating Arabidopsis SCW biosynthesis, and numerous studies have explored SCW-related transcription factors in other dicots and monocots. Whilst the general structure of the Arabidopsis network has been a topic of several reviews, they have not comprehensively represented the detailed protein-DNA and protein-protein interactions described in the literature, and an understanding of network dynamics and functionality has not yet been achieved for SCW formation. Furthermore the methodologies employed in studies of SCW transcriptional regulation have not received much attention, especially in the case of non-model organisms. In this review, we have reconstructed the most exhaustive literature-based network representations to date of SCW transcriptional regulation in Arabidopsis. We include a manipulable Cytoscape representation of the Arabidopsis SCW transcriptional network to aid in future studies, along with a list of supporting literature for each documented interaction. Amongst other topics, we discuss the various components of the network, its evolutionary conservation in plants, putative modules and dynamic mechanisms that may influence network function, and the approaches that have been employed in network inference. Future research should aim to better understand network function and its response to dynamic perturbations, whilst the development and application of genome-wide approaches such as ChIP-seq and systems genetics are in progress for the study of SCW transcriptional regulation in non-model organisms.

Identification of molecular processes needed for vascular formation through transcriptome analysis of different vascular systems

Background

Vascular system formation has been studied through molecular and genetic approaches in Arabidopsis, a herbaceous dicot that is used as a model system. Different vascular systems have developed in other plants such as crops and trees. Uncovering shared mechanisms underlying vascular development by transcriptome analysis of different vascular systems may help to transfer knowledge acquired from Arabidopsis to other economically important species.

Results

Conserved vascular genes and biological processes fundamental to vascular development were explored across various plants. Through comparative transcriptome analysis, 226 genes from Arabidopsis, 217 genes from poplar and 281 genes from rice were identified as constituting 107 conserved vascular gene groups. These gene groups are expressed mainly in vascular tissues and form a complex coexpression network with multiple functional connections. To date, only half of the groups have been experimentally investigated. The conserved vascular gene groups were classified into 9 essential processes for vascular development. 18 groups (17%) lack of annotations were classified as having unknown functions.

Conclusion

The study provides a map of fundamental biological processes conserved across different vascular systems. It identifies gaps in the experimental investigation of pathways active in vascular formation, which if explored, could lead to a more complete understanding of vascular development.

Induced somatic sector analysis of cellulose synthase (CesA) promoter regions in woody stem tissues

The increasing focus on plantation forestry as a renewable source of cellulosic biomass has emphasized the need for tools to study the unique biology of woody genera such as Eucalyptus, Populus and Pinus. The domestication of these woody crops is hampered by long generation times, and breeders are now looking to molecular approaches such as marker-assisted breeding and genetic modification to accelerate tree improvement. Much of what is known about genes involved in the growth and development of plants has come from studies of herbaceous models such as Arabidopsis and rice. However, transferring this information to woody plants often proves difficult, especially for genes expressed in woody stems. Here we report the use of induced somatic sector analysis (ISSA) for characterization of promoter expression patterns directly in the stems of Populus and Eucalyptus trees. As a case study, we used previously characterized primary and secondary cell wall-related cellulose synthase (CesA) promoters cloned from Eucalyptus grandis. We show that ISSA can be used to elucidate the phloem and xylem expression patterns of the CesA genes in Eucalyptus and Populus stems and also show that the staining patterns differ in Eucalyptus and Populus stems. These findings show that ISSA is an efficient approach to investigate promoter function in the developmental context of woody plant tissues and raise questions about the suitability of heterologous promoters for genetic manipulation in plant species.

Identification of quantitative trait loci controlling fibre length and lignin content in Arabidopsis thaliana stems

Fibre properties and the biochemical composition of cell walls are important traits in many applications. For example, the lengths of fibres define the strength and quality of paper, and lignin content is a critical parameter for the use of biomass in biofuel production. Identifying genes controlling these traits is comparatively difficult in woody species, because of long generation times and limited amenability to high-resolution genetic mapping. To address this problem, this study mapped quantitative trait loci (QTLs) defining fibre length and lignin content in the Arabidopsis recombinant inbred line population Col-4 × Ler-0. Adapting high-throughput phenotyping techniques for both traits for measurements in Arabidopsis inflorescence stems identified significant QTLs for fibre length on chromosomes 2 and 5, as well as one significant QTL affecting lignin content on chromosome 2. For fibre length, total variation within the population was 208% higher than between parental lines and the identified QTLs explained 50.58% of the observed variation. For lignin content, the values were 261 and 26.51%, respectively. Bioinformatics analysis of the associated intervals identified a number of candidate genes for fibre length and lignin content. This study demonstrates that molecular mapping of QTLs pertaining to wood and fibre properties is possible in Arabidopsis, which substantially broadens the use of Arabidopsis as a model species for the functional characterization of plant genes.