Overexpression of REDUCED WALL ACETYLATION C increases xylan acetylation and biomass recalcitrance in Populus

Plant lignocellulosic biomass, i.e. secondary cell walls of plants, is a vital alternative source for bioenergy. However, the acetylation of xylan in secondary cell walls impedes the conversion of biomass to biofuels. Previous studies have shown that REDUCED WALL ACETYLATION (RWA) proteins are directly involved in the acetylation of xylan but the regulatory mechanism of RWAs is not fully understood. In this study, we demonstrate that overexpression of a Populus trichocarpa PtRWA-C gene increases the level of xylan acetylation and increases the lignin content and S/G ratio, ultimately yielding poplar woody biomass with reduced saccharification efficiency. Furthermore, through gene coexpression network and expression quantitative trait loci (eQTL) analysis, we found that PtRWA-C was regulated not only by the secondary cell wall hierarchical regulatory network but also by an AP2 family transcription factor HARDY (HRD). Specifically, HRD activates PtRWA-C expression by directly binding to the PtRWA-C promoter, which is also the cis-eQTL for PtRWA-C. Taken together, our findings provide insights into the functional roles of PtRWA-C in xylan acetylation and consequently saccharification and shed light on synthetic biology approaches to manipulate this gene and alter cell wall properties. These findings have substantial implications for genetic engineering of woody species, which could be used as a sustainable source of biofuels, valuable biochemicals, and biomaterials.

Expression quantitative trait loci mapping identified PtrXB38 as a key hub gene in adventitious root development in Populus

Plant establishment requires the formation and development of an extensive root system with architecture modulated by complex genetic networks. Here, we report the identification of the PtrXB38 gene as an expression quantitative trait loci (eQTL) hotspot, mapped using 390 leaf and 444 xylem Populus trichocarpa transcriptomes. Among predicted targets of this trans-eQTL were genes involved in plant hormone responses and root development. Overexpression of PtrXB38 in Populus led to significant increases in callusing and formation of both stem-born roots and base-born adventitious roots. Omics studies revealed that genes and proteins controlling auxin transport and signaling were involved in PtrXB38-mediated adventitious root formation. Protein-protein interaction assays indicated that PtrXB38 interacts with components of endosomal sorting complexes required for transport machinery, implying that PtrXB38-regulated root development may be mediated by regulating endocytosis pathway. Taken together, this work identified a crucial root development regulator and sheds light on the discovery of other plant developmental regulators through combining eQTL mapping and omics approaches.

The Plasminogen-Apple-Nematode (PAN) domain suppresses JA/ET defense pathways in plants

Suppression of immune response is a phenomenon that enables biological processes such as gamete fertilization, cell growth, cell proliferation, endophyte recruitment, parasitism, and pathogenesis. Here, we show for the first time that the Plasminogen-Apple-Nematode (PAN) domain present in G-type lectin receptor-like kinases is essential for immunosuppression in plants. Defense pathways involving jasmonic acid and ethylene are critical for plant immunity against microbes, necrotrophic pathogens, parasites, and insects. Using two Salix purpurea G-type lectin receptor kinases, we demonstrated that intact PAN domains suppress jasmonic acid and ethylene signaling in Arabidopsis and tobacco. Variants of the same receptors with mutated residues in this domain could trigger induction of both defense pathways. Assessment of signaling processes revealed significant differences between receptors with intact and mutated PAN domain in MAPK phosphorylation, global transcriptional reprogramming, induction of downstream signaling components, hormone biosynthesis and resistance to Botrytis cinerea . Further, we demonstrated that the domain is required for oligomerization, ubiquitination, and proteolytic degradation of these receptors. These processes were completely disrupted when conserved residues in the domain were mutated. Additionally, we have tested the hypothesis in recently characterized Arabidopsis mutant which has predicted PAN domain and negatively regulates plant immunity against root nematodes. ern1.1 mutant complemented with mutated PAN shows triggered immune response with elevated WRKY33 expression, hyperphosphorylation of MAPK and resistant to necrotrophic fungus Botrytis cinerea . Collectively, our results suggest that ubiquitination and proteolytic degradation mediated by the PAN domain plays a role in receptor turn-over to suppress jasmonic acid and ethylene defense signaling in plants.

Nitrogen addition alters soil fungal communities, but root fungal communities are resistant to change

Plants are colonized by numerous microorganisms serving important symbiotic functions that are vital to plant growth and success. Understanding and harnessing these interactions will be useful in both managed and natural ecosystems faced with global change, but it is still unclear how variation in environmental conditions and soils influence the trajectory of these interactions. In this study, we examine how nitrogen addition alters plant-fungal interactions within two species of Populus - Populus deltoides and P. trichocarpa. In this experiment, we manipulated plant host, starting soil (native vs. away for each tree species), and nitrogen addition in a fully factorial replicated design. After ~10 weeks of growth, we destructively harvested the plants and characterized plant growth factors and the soil and root endosphere fungal communities using targeted amplicon sequencing of the ITS2 gene region. Overall, we found nitrogen addition altered plant growth factors, e.g., plant height, chlorophyll density, and plant N content. Interestingly, nitrogen addition resulted in a lower fungal alpha diversity in soils but not plant roots. Further, there was an interactive effect of tree species, soil origin, and nitrogen addition on soil fungal community composition. Starting soils collected from Oregon and West Virginia were dominated by the ectomycorrhizal fungi Inocybe (55.8% relative abundance), but interestingly when P. deltoides was grown in its native West Virginia soil, the roots selected for a high abundance of the arbuscular mycorrhizal fungi, Rhizophagus. These results highlight the importance of soil origin and plant species on establishing plant-fungal interactions.

Development of an Experimental Approach to Achieve Spatially Resolved Plant Root-Associated Metaproteomics Using an Agar-Plate System

Plant-microbe interactions in the rhizosphere play a vital role in plant health and productivity. The composition and function of root-associated microbiomes is strongly influenced by their surrounding environment, which is often customized by their host. How microbiomes change with respect to space and time across plant roots remains poorly understood, and methodologies that facilitate spatiotemporal metaproteomic studies of root-associated microbiomes are yet to be realized. Here, we developed a method that provides spatially resolved metaproteome measurements along plant roots embedded in agar-plate culture systems, which have long been used to study plants. Spatially defined agar "plugs" of interest were excised and subsequently processed using a novel peptide extraction method prior to metaproteomics, which was used to infer both microbial community composition and function. As a proof-of-principle, a previously studied 10-member community constructed from a Populus root system was grown in an agar plate with a 3-week-old Populus trichocarpa plant. Metaproteomics was performed across two time points (24 and 48 h) for three distinct locations (root base, root tip, and a region distant from the root). The spatial resolution of these measurements provides evidence that microbiome composition and expression changes across the plant root interface. Interrogation of the individual microbial proteomes revealed functional profiles related to their behavioral associations with the plant root, in which chemotaxis and augmented metabolism likely supported predominance of the most abundant member. This study demonstrated a novel peptide extraction method for studying plant agar-plate culture systems, which was previously unsuitable for (meta)proteomic measurements.

Biomass formation and sugar release efficiency of Populus modified by altered expression of a NAC transcription factor

Woody biomass is an important feedstock for biofuel production. Manipulation of wood properties that enable efficient conversion of biomass to biofuel reduces cost of biofuel production. Wood cell wall composition is regulated at several levels that involve expression of transcription factors such as wood-/secondary cell wall-associated NAC domains (WND or SND). In Arabidopsis thaliana, SND1 regulates cell wall composition through activation of its down-stream targets such as MYBs. The functional aspects of SND1 homologs in the woody Populus have been studied through transgenic manipulation. In this study, we investigated the role of PdWND1B, Populus SND1 sequence ortholog, in wood formation using transgenic manipulation through over-expression or silencing under the control of a vascular-specific 4-coumarate-CoA ligase (4CL) promoter. As compared with control plants, PdWND1B-RNAi plants were shorter in height, with significantly reduced stem diameter and dry biomass, whereas there were no significant differences in growth and productivity of PdWND1B over-expression plants. Conversely, PdWND1B over-expression lines showed a significant reduction in cellulose and increase in lignin content, whereas there was no significant impact on lignin content of downregulated lines. Stem carbohydrate composition analysis revealed a decrease in glucose, mannose, arabinose, and galactose, but an increase in xylose in the over-expression lines. Transcriptome analysis revealed upregulation of several downstream transcription factors and secondary cell wall related structural genes in the PdWND1B over-expression lines, partly explaining the observed phenotypic changes in cell wall chemistry. Relative to the control, glucose release efficiency and ethanol production from stem biomass was significantly reduced in over-expression lines. Our results show that PdWND1B is an important factor determining biomass productivity, cell wall chemistry and its conversion to biofuels in Populus.

Habitat-adapted microbial communities mediate Sphagnum peatmoss resilience to warming

Sphagnum peatmosses are fundamental members of peatland ecosystems, where they contribute to the uptake and long-term storage of atmospheric carbon. Warming threatens Sphagnum mosses and is known to alter the composition of their associated microbiome. Here, we use a microbiome transfer approach to test if microbiome thermal origin influences host plant thermotolerance. We leveraged an experimental whole-ecosystem warming study to collect field-grown Sphagnum, mechanically separate the associated microbiome and then transfer onto germ-free laboratory Sphagnum for temperature experiments. Host and microbiome dynamics were assessed with growth analysis, Chla fluorescence imaging, metagenomics, metatranscriptomics and 16S rDNA profiling. Microbiomes originating from warming field conditions imparted enhanced thermotolerance and growth recovery at elevated temperatures. Metagenome and metatranscriptome analyses revealed that warming altered microbial community structure in a manner that induced the plant heat shock response, especially the HSP70 family and jasmonic acid production. The heat shock response was induced even without warming treatment in the laboratory, suggesting that the warm-microbiome isolated from the field provided the host plant with thermal preconditioning. Our results demonstrate that microbes, which respond rapidly to temperature alterations, can play key roles in host plant growth response to rapidly changing environments.

Forest stand and canopy development unaltered by 12 years of CO2 enrichment

Canopy structure-the size and distribution of tree crowns and the spatial and temporal distribution of leaves within them-exerts dominant control over primary productivity, transpiration and energy exchange. Stand structure-the spatial arrangement of trees in the forest (height, basal area and spacing)-has a strong influence on forest growth, allocation and resource use. Forest response to elevated atmospheric CO2 is likely to be dependent on the canopy and stand structure. Here, we investigated elevated CO2 effects on the forest structure of a Liquidambar styraciflua L. stand in a free-air CO2 enrichment experiment, considering leaves, tree crowns, forest canopy and stand structure. During the 12-year experiment, the trees increased in height by 5 m and basal area increased by 37%. Basal area distribution among trees shifted from a relatively narrow distribution to a much broader one, but there was little evidence of a CO2 effect on height growth or basal area distribution. The differentiation into crown classes over time led to an increase in the number of unproductive intermediate and suppressed trees and to a greater concentration of stand basal area in the largest trees. A whole-tree harvest at the end of the experiment permitted detailed analysis of canopy structure. There was little effect of CO2 enrichment on the relative leaf area distribution within tree crowns and there was little change from 1998 to 2009. Leaf characteristics (leaf mass per unit area and nitrogen content) varied with crown depth; any effects of elevated CO2 were much smaller than the variation within the crown and were consistent throughout the crown. In this young, even-aged, monoculture plantation forest, there was little evidence that elevated CO2 accelerated tree and stand development, and there were remarkably small changes in canopy structure. Questions remain as to whether a more diverse, mixed species forest would respond similarly.

Temporal dynamics of protein and post-translational modification abundances in Populus leaf across a diurnal period

Populus spp. are dedicated woody biomass feedstocks for advanced biofuels and bioproducts. Proper growth and fitness of poplar as a sustainable feedstock depends on timely perception and response to environmental signals (e.g., light, temperature, water). Poplar leaves, like other C3 photosynthesis plants, have evolved oscillating or circadian rhythms that play important roles in synchronizing biological processes with external cues. To characterize this phenomenon at a molecular level, we employed bottom-up proteomics using high-resolution mass spectrometry and de novo-assisted database searching to identify abundance changes in proteins and post-translational modifications in poplar leaf tissue sampled across a 12/12-hour light/dark diurnal period.

The Ancient Salicoid Genome Duplication Event: A Platform for Reconstruction of De Novo Gene Evolution in Populus trichocarpa

Orphan genes are characteristic genomic features that have no detectable homology to genes in any other species and represent an important attribute of genome evolution as sources of novel genetic functions. Here, we identified 445 genes specific to Populus trichocarpa. Of these, we performed deeper reconstruction of 13 orphan genes to provide evidence of de novo gene evolution. Populus and its sister genera Salix are particularly well suited for the study of orphan gene evolution because of the Salicoid whole-genome duplication event which resulted in highly syntenic sister chromosomal segments across the Salicaceae. We leveraged this genomic feature to reconstruct de novo gene evolution from intergenera, interspecies, and intragenomic perspectives by comparing the syntenic regions within the P. trichocarpa reference, then P. deltoides, and finally Salix purpurea. Furthermore, we demonstrated that 86.5% of the putative orphan genes had evidence of transcription. Additionally, we also utilized the Populus genome-wide association mapping panel, a collection of 1,084 undomesticated P. trichocarpa genotypes to further determine putative regulatory networks of orphan genes using expression quantitative trait loci (eQTL) mapping. Functional enrichment of these eQTL subnetworks identified common biological themes associated with orphan genes such as response to stress and defense response. We also identify a putative cis-element for a de novo gene and leverage conserved synteny to describe evolution of a putative transcription factor binding site. Overall, 45% of orphan genes were captured in trans-eQTL networks.

Overexpression of a Prefoldin β subunit gene reduces biomass recalcitrance in the bioenergy crop Populus

Prefoldin (PFD) is a group II chaperonin that is ubiquitously present in the eukaryotic kingdom. Six subunits (PFD1-6) form a jellyfish-like heterohexameric PFD complex and function in protein folding and cytoskeleton organization. However, little is known about its function in plant cell wall-related processes. Here, we report the functional characterization of a PFD gene from Populus deltoides, designated as PdPFD2.2. There are two copies of PFD2 in Populus, and PdPFD2.2 was ubiquitously expressed with high transcript abundance in the cambial region. PdPFD2.2 can physically interact with DELLA protein RGA1_8g, and its subcellular localization is affected by the interaction. In P. deltoides transgenic plants overexpressing PdPFD2.2, the lignin syringyl/guaiacyl ratio was increased, but cellulose content and crystallinity index were unchanged. In addition, the total released sugar (glucose and xylose) amounts were increased by 7.6% and 6.1%, respectively, in two transgenic lines. Transcriptomic and metabolomic analyses revealed that secondary metabolic pathways, including lignin and flavonoid biosynthesis, were affected by overexpressing PdPFD2.2. A total of eight hub transcription factors (TFs) were identified based on TF binding sites of differentially expressed genes in Populus transgenic plants overexpressing PdPFD2.2. In addition, several known cell wall-related TFs, such as MYB3, MYB4, MYB7, TT8 and XND1, were affected by overexpression of PdPFD2.2. These results suggest that overexpression of PdPFD2.2 can reduce biomass recalcitrance and PdPFD2.2 is a promising target for genetic engineering to improve feedstock characteristics to enhance biofuel conversion and reduce the cost of lignocellulosic biofuel production.

Identification of functional single nucleotide polymorphism of Populus trichocarpa PtrEPSP-TF and determination of its transcriptional effect

In plants, the phenylpropanoid pathway is responsible for the synthesis of a diverse array of secondary metabolites that include lignin monomers, flavonoids, and coumarins, many of which are essential for plant structure, biomass recalcitrance, stress defense, and nutritional quality. Our previous studies have demonstrated that Populus trichocarpa PtrEPSP-TF, an isoform of 5-enolpyruvylshikimate 3-phosphate (EPSP) synthase, has transcriptional activity and regulates phenylpropanoid biosynthesis in Populus. In this study, we report the identification of single nucleotide polymorphism (SNP) of PtrEPSP-TF that defines its functionality. Populus natural variants carrying this SNP were shown to have reduced lignin content. Here, we demonstrated that the SNP-induced substitution of 142nd amino acid (PtrEPSP-TFD142E) dramatically impairs the DNA-binding and transcriptional activity of PtrEPSP-TF. When introduced to a monocot species rice (Oryza sativa) in which an EPSP synthase isoform with the DNA-binding helix-turn-helix (HTH) motif is absent, the PtrEPSP-TF, but not PtrEPSP-TFD142E, activated genes in the phenylpropanoid pathway. More importantly, heterologous expression of PtrEPSP-TF uncovered five new transcriptional regulators of phenylpropanoid biosynthesis in rice. Collectively, this study identifies the key amino acid required for PtrEPSP-TF functionality and provides a strategy to uncover new transcriptional regulators in phenylpropanoid biosynthesis.

PdWND3A, a wood-associated NAC domain-containing protein, affects lignin biosynthesis and composition in Populus

BACKGROUND

Plant secondary cell wall is a renewable feedstock for biofuels and biomaterials production. Arabidopsis VASCULAR-RELATED NAC DOMAIN (VND) has been demonstrated to be a key transcription factor regulating secondary cell wall biosynthesis. However, less is known about its role in the woody species.

RESULTS

Here we report the functional characterization of Populus deltoides WOOD-ASSOCIATED NAC DOMAIN protein 3 (PdWND3A), a sequence homolog of Arabidopsis VND4 and VND5 that are members of transcription factor networks regulating secondary cell wall biosynthesis. PdWND3A was expressed at higher level in the xylem than in other tissues. The stem tissues of transgenic P. deltoides overexpressing PdWND3A (OXPdWND3A) contained more vessel cells than that of wild-type plants. Furthermore, lignin content and lignin monomer syringyl and guaiacyl (S/G) ratio were higher in OXPdWND3A transgenic plants than in wild-type plants. Consistent with these observations, the expression of FERULATE 5-HYDROXYLASE1 (F5H1), encoding an enzyme involved in the biosynthesis of sinapyl alcohol (S unit monolignol), was elevated in OXPdWND3A transgenic plants. Saccharification analysis indicated that the rate of sugar release was reduced in the transgenic plants. In addition, OXPdWND3A transgenic plants produced lower amounts of biomass than wild-type plants.

CONCLUSIONS

PdWND3A affects lignin biosynthesis and composition and negatively impacts sugar release and biomass production.

Relatively rare root endophytic bacteria drive plant resource allocation patterns and tissue nutrient concentration in unpredictable ways

PREMISE

Plant endophytic bacterial strains can influence plant traits such as leaf area and root length. Yet, the influence of more complex bacterial communities in regulating overall plant phenotype is less explored. Here, in two complementary experiments, we tested whether we can predict plant phenotype response to changes in microbial community composition.

METHODS

In the first study, we inoculated a single genotype of Populus deltoides with individual root endophytic bacteria and measured plant phenotype. Next, data from this single inoculation were used to predict phenotypic traits after mixed three-strain community inoculations, which we tested in the second experiment.

RESULTS

By itself, each bacterial endophyte significantly but weakly altered plant phenotype relative to noninoculated plants. In a mixture, bacterial strain Burkholderia BT03, constituted at least 98% of community relative abundance. Yet, plant resource allocation and tissue nutrient concentrations were disproportionately influenced by Pseudomonas sp. GM17, GM30, and GM41. We found a 10% increase in leaf mass fraction and an 11% decrease in root mass fraction when replacing Pseudomonas GM17 with GM41 in communities containing both Pseudomonas GM30 and Burkholderia BT03.

CONCLUSIONS

Our results indicate that interactions among endophytic bacteria may drive plant phenotype over the contribution of each strain individually. Additionally, we have shown that low-abundance strains contribute to plant phenotype challenging the assumption that the dominant strains will drive plant function.

The nature of the progression of drought stress drives differential metabolomic responses in Populus deltoides

BACKGROUND AND AIMS

The use of woody crops for Quad-level (approx. 1 × 1018 J) energy production will require marginal agricultural lands that experience recurrent periods of water stress. Populus species have the capacity to increase dehydration tolerance by lowering osmotic potential via osmotic adjustment. The aim of this study was to investigate how the inherent genetic potential of a Populus clone to respond to drought interacts with the nature of the drought to determine the degree of biochemical response.

METHODS

A greenhouse drought stress study was conducted on Populus deltoides 'WV94' and the resulting metabolite profiles of leaves were determined by gas chromatography-mass spectrometry following trimethylsilylation for plants subjected to cyclic mild (-0.5 MPa pre-dawn leaf water potential) drought vs. cyclic severe (-1.26 MPa) drought in contrast to well-watered controls (-0.1 MPa) after two or four drought cycles, and in contrast to plants subjected to acute drought, where plants were desiccated for up to 8 d.

KEY RESULTS

The nature of drought (cyclic vs. acute), frequency of drought (number of cycles) and the severity of drought (mild vs. severe) all dictated the degree of osmotic adjustment and the nature of the organic solutes that accumulated. Whereas cyclic drought induced the largest responses in primary metabolism (soluble sugars, organic acids and amino acids), acute onset of prolonged drought induced the greatest osmotic adjustment and largest responses in secondary metabolism, especially populosides (hydroxycinnamic acid conjugates of salicin).

CONCLUSIONS

The differential adaptive metabolite responses in cyclic vs. acute drought suggest that stress acclimation occurs via primary metabolism in response to cyclic drought, whereas expanded metabolic plasticity occurs via secondary metabolism following severe, acute drought. The shift in carbon partitioning to aromatic metabolism with the production of a diverse suite of higher order salicylates lowers osmotic potential and increases the probability of post-stress recovery.

Mediation of plant-mycorrhizal interaction by a lectin receptor-like kinase

The molecular mechanisms underlying mycorrhizal symbioses, the most ubiquitous and impactful mutualistic plant-microbial interaction in nature, are largely unknown. Through genetic mapping, resequencing and molecular validation, we demonstrate that a G-type lectin receptor-like kinase (lecRLK) mediates the symbiotic interaction between Populus and the ectomycorrhizal fungus Laccaria bicolor. This finding uncovers an important molecular step in the establishment of symbiotic plant-fungal associations and provides a molecular target for engineering beneficial mycorrhizal relationships.

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.

Genome-wide association studies and expression-based quantitative trait loci analyses reveal roles of HCT2 in caffeoylquinic acid biosynthesis and its regulation by defense-responsive transcription factors in Populus

3-O-caffeoylquinic acid, also known as chlorogenic acid (CGA), functions as an intermediate in lignin biosynthesis in the phenylpropanoid pathway. It is widely distributed among numerous plant species and acts as an antioxidant in both plants and animals. Using GC-MS, we discovered consistent and extreme variation in CGA content across a population of 739 4-yr-old Populus trichocarpa accessions. We performed genome-wide association studies (GWAS) from 917 P. trichocarpa accessions and expression-based quantitative trait loci (eQTL) analyses to identify key regulators. The GWAS and eQTL analyses resolved an overlapped interval encompassing a hydroxycinnamoyl-CoA:shikimate hydroxycinnamoyl transferase 2 (PtHCT2) that was significantly associated with CGA and partially characterized metabolite abundances. PtHCT2 leaf expression was significantly correlated with CGA abundance and it was regulated by cis-eQTLs containing W-box for WRKY binding. Among all nine PtHCT homologs, PtHCT2 is the only one that responds to infection by the fungal pathogen Sphaerulina musiva (a Populus pathogen). Validation using protoplast-based transient expression system suggests that PtHCT2 is regulated by the defense-responsive WRKY. These results are consistent with reports of CGA functioning as an antioxidant in response to biotic stress. This study provides insights into data-driven and omics-based inference of gene function in woody species.

Diel rewiring and positive selection of ancient plant proteins enabled evolution of CAM photosynthesis in Agave

Background

Crassulacean acid metabolism (CAM) enhances plant water-use efficiency through an inverse day/night pattern of stomatal closure/opening that facilitates nocturnal CO₂ uptake. CAM has evolved independently in over 35 plant lineages, accounting for ~ 6% of all higher plants. Agave species are highly heat- and drought-tolerant, and have been domesticated as model CAM crops for beverage, fiber, and biofuel production in semi-arid and arid regions. However, the genomic basis of evolutionary innovation of CAM in genus Agave is largely unknown.

Results

Using an approach that integrated genomics, gene co-expression networks, comparative genomics and protein structure analyses, we investigated the molecular evolution of CAM as exemplified in Agave. Comparative genomics analyses among C₃, C₄ and CAM species revealed that core metabolic components required for CAM have ancient genomic origins traceable to non-vascular plants while regulatory proteins required for diel re-programming of metabolism have a more recent origin shared among C₃, C₄ and CAM species. We showed that accelerated evolution of key functional domains in proteins responsible for primary metabolism and signaling, together with a diel re-programming of the transcription of genes involved in carbon fixation, carbohydrate processing, redox homeostasis, and circadian control is required for the evolution of CAM in Agave. Furthermore, we highlighted the potential candidates contributing to the adaptation of CAM functional modules.

Conclusions

This work provides evidence of adaptive evolution of CAM related pathways. We showed that the core metabolic components required for CAM are shared by non-vascular plants, but regulatory proteins involved in re-reprogramming of carbon fixation and metabolite transportation appeared more recently. We propose that the accelerated evolution of key proteins together with a diel re-programming of gene expression were required for CAM evolution from C₃ ancestors in Agave.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4964-7) contains supplementary material, which is available to authorized users.

A 5-Enolpyruvylshikimate 3-Phosphate Synthase Functions as a Transcriptional Repressor in Populus

Long-lived perennial plants, with distinctive habits of inter-annual growth, defense, and physiology, are of great economic and ecological importance. However, some biological mechanisms resulting from genome duplication and functional divergence of genes in these systems remain poorly studied. Here, we discovered an association between a poplar (Populus trichocarpa) 5-enolpyruvylshikimate 3-phosphate synthase gene (PtrEPSP) and lignin biosynthesis. Functional characterization of PtrEPSP revealed that this isoform possesses a helix-turn-helix motif in the N terminus and can function as a transcriptional repressor that regulates expression of genes in the phenylpropanoid pathway in addition to performing its canonical biosynthesis function in the shikimate pathway. We demonstrated that this isoform can localize in the nucleus and specifically binds to the promoter and represses the expression of a SLEEPER-like transcriptional regulator, which itself specifically binds to the promoter and represses the expression of PtrMYB021 (known as MYB46 in Arabidopsis thaliana), a master regulator of the phenylpropanoid pathway and lignin biosynthesis. Analyses of overexpression and RNAi lines targeting PtrEPSP confirmed the predicted changes in PtrMYB021 expression patterns. These results demonstrate that PtrEPSP in its regulatory form and PtrhAT form a transcriptional hierarchy regulating phenylpropanoid pathway and lignin biosynthesis in Populus.

Quantitative proteome profile of water deficit stress responses in eastern cottonwood (Populus deltoides) leaves

Drought stress is a recurring feature of world climate and the single most important factor influencing agricultural yield worldwide. Plants display highly variable, species-specific responses to drought and these responses are multifaceted, requiring physiological and morphological changes influenced by genetic and molecular mechanisms. Moreover, the reproducibility of water deficit studies is very cumbersome, which significantly impedes research on drought tolerance, because how a plant responds is highly influenced by the timing, duration, and intensity of the water deficit. Despite progress in the identification of drought-related mechanisms in many plants, the molecular basis of drought resistance remains to be fully understood in trees, particularly in poplar species because their wide geographic distribution results in varying tolerances to drought. Herein, we aimed to better understand this complex phenomenon in eastern cottonwood (Populus deltoides) by performing a detailed contrast of the proteome changes between two different water deficit experiments to identify functional intersections and divergences in proteome responses. We investigated plants subjected to cyclic water deficit and compared these responses to plants subjected to prolonged acute water deficit. In total, we identified 108,012 peptide sequences across both experiments that provided insight into the quantitative state of 22,737 Populus gene models and 8,199 functional protein groups in response to drought. Together, these datasets provide the most comprehensive insight into proteome drought responses in poplar to date and a direct proteome comparison between short period dehydration shock and cyclic, post-drought re-watering. Overall, this investigation provides novel insights into drought avoidance mechanisms that are distinct from progressive drought stress. Additionally, we identified proteins that have been associated as drought-relevant in previous studies. Importantly, we highlight the RD26 transcription factor as a gene regulated at both the transcript and protein level, regardless of species and drought condition, and, thus, represents a key, universal drought marker for Populus species.

Overexpression of a Domain of Unknown Function 266-containing protein results in high cellulose content, reduced recalcitrance, and enhanced plant growth in the bioenergy crop Populus

BACKGROUND

Domain of Unknown Function 266 (DUF266) is a plant-specific domain. DUF266-containing proteins (DUF266 proteins) have been categorized as 'not classified glycosyltransferases (GTnc)' due to amino acid similarity with GTs. However, little is known about the function of DUF266 proteins.

RESULTS

Phylogenetic analysis revealed that DUF266 proteins are only present in the land plants including moss and lycophyte. We report the functional characterization of one member of DUF266 proteins in Populus, PdDUF266A. PdDUF266A was ubiquitously expressed with high abundance in the xylem. In Populus transgenic plants overexpressing PdDUF266A (OXPdDUF266A), the glucose and cellulose contents were significantly higher, while the lignin content was lower than that in the wild type. Degree of polymerization of cellulose in OXPdDUF266A transgenic plants was also higher, whereas cellulose crystallinity index remained unchanged. Gene expression analysis indicated that cellulose biosynthesis-related genes such as CESA and SUSY were upregulated in mature leaf and xylem of OXPdDUF266A transgenic plants. Moreover, PdDUF266A overexpression resulted in an increase of biomass production. Their glucose contents and biomass phenotypes were further validated via heterologous expression of PdDUF266A in Arabidopsis. Results from saccharification treatment demonstrated that the rate of sugar release was increased by approximately 38% in the OXPdDUF266A transgenic plants.

CONCLUSIONS

These results suggest that the overexpression of PdDUF266A can increase cellulose content, reduce recalcitrance, and enhance biomass production, and that PdDUF266A is a promising target for genetic manipulation for biofuel production.

Overexpression of a Domain of Unknown Function 231-containing protein increases O-xylan acetylation and cellulose biosynthesis in Populus

BACKGROUND

Domain of Unknown Function 231-containing proteins (DUF231) are plant specific and their function is largely unknown. Studies in the model plants Arabidopsis and rice suggested that some DUF231 proteins act in the process of O-acetyl substitution of hemicellulose and esterification of pectin. However, little is known about the function of DUF231 proteins in woody plant species.

RESULTS

This study provides evidence supporting that one member of DUF231 family proteins in the woody perennial plant Populus deltoides (genotype WV94), PdDUF231A, has a role in the acetylation of xylan and affects cellulose biosynthesis. A total of 52 DUF231-containing proteins were identified in the Populus genome. In P. deltoides transgenic lines overexpressing PdDUF231A (OXPdDUF231A), glucose and cellulose contents were increased. Consistent with these results, the transcript levels of cellulose biosynthesis-related genes were increased in the OXPdDUF231A transgenic lines. Furthermore, the relative content of total acetylated xylan was increased in the OXPdDUF231A transgenic lines. Enzymatic saccharification assays revealed that the rate of glucose release increased in OXPdDUF231A transgenic lines. Plant biomass productivity was also increased in OXPdDUF231A transgenic lines.

CONCLUSIONS

These results suggest that PdDUF231A affects cellulose biosynthesis and plays a role in the acetylation of xylan. PdDUF231A is a promising target for genetic modification for biofuel production because biomass productivity and compositional quality can be simultaneously improved through overexpression.

Root bacterial endophytes alter plant phenotype, but not physiology

Plant traits, such as root and leaf area, influence how plants interact with their environment and the diverse microbiota living within plants can influence plant morphology and physiology. Here, we explored how three bacterial strains isolated from the Populus root microbiome, influenced plant phenotype. We chose three bacterial strains that differed in predicted metabolic capabilities, plant hormone production and metabolism, and secondary metabolite synthesis. We inoculated each bacterial strain on a single genotype of Populus trichocarpa and measured the response of plant growth related traits (root:shoot, biomass production, root and leaf growth rates) and physiological traits (chlorophyll content, net photosynthesis, net photosynthesis at saturating light-Asat, and saturating CO2-Amax). Overall, we found that bacterial root endophyte infection increased root growth rate up to 184% and leaf growth rate up to 137% relative to non-inoculated control plants, evidence that plants respond to bacteria by modifying morphology. However, endophyte inoculation had no influence on total plant biomass and photosynthetic traits (net photosynthesis, chlorophyll content). In sum, bacterial inoculation did not significantly increase plant carbon fixation and biomass, but their presence altered where and how carbon was being allocated in the plant host.

Down-Regulation of KORRIGAN-Like Endo-β-1,4-Glucanase Genes Impacts Carbon Partitioning, Mycorrhizal Colonization and Biomass Production in Populus

A greater understanding of the genetic regulation of plant cell wall remodeling and the impact of modified cell walls on plant performance is important for the development of sustainable biofuel crops. Here, we studied the impact of down-regulating KORRIGAN-like cell wall biosynthesis genes, belonging to the endo-β-1,4-glucanase gene family, on Populus growth, metabolism and the ability to interact with symbiotic microbes. The reductions in cellulose content and lignin syringyl-to-guaiacyl unit ratio, and increase in cellulose crystallinity of cell walls of PdKOR RNAi plants corroborated the functional role of PdKOR in cell wall biosynthesis. Altered metabolism and reduced growth characteristics of RNAi plants revealed new implications on carbon allocation and partitioning. The distinctive metabolome phenotype comprised of a higher phenolic and salicylic acid content, and reduced lignin, shikimic acid and maleic acid content relative to control. Plant sustainability implications of modified cell walls on beneficial plant-microbe interactions were explored via co-culture with an ectomycorrhizal fungus, Laccaria bicolor. A significant increase in the mycorrhization rate was observed in transgenic plants, leading to measurable beneficial growth effects. These findings present new evidence for functional interconnectedness of cellulose biosynthesis pathway, metabolism and mycorrhizal association in plants, and further emphasize the consideration of the sustainability implications of plant trait improvement efforts.

A Carotenoid-Deficient Mutant in Pantoea sp. YR343, a Bacteria Isolated from the Rhizosphere of Populus deltoides, Is Defective in Root Colonization

The complex interactions between plants and their microbiome can have a profound effect on the health and productivity of the plant host. A better understanding of the microbial mechanisms that promote plant health and stress tolerance will enable strategies for improving the productivity of economically important plants. Pantoea sp. YR343 is a motile, rod-shaped bacterium isolated from the roots of Populus deltoides that possesses the ability to solubilize phosphate and produce the phytohormone indole-3-acetic acid (IAA). Pantoea sp. YR343 readily colonizes plant roots and does not appear to be pathogenic when applied to the leaves or roots of selected plant hosts. To better understand the molecular mechanisms involved in plant association and rhizosphere survival by Pantoea sp. YR343, we constructed a mutant in which the crtB gene encoding phytoene synthase was deleted. Phytoene synthase is responsible for converting geranylgeranyl pyrophosphate to phytoene, an important precursor to the production of carotenoids. As predicted, the ΔcrtB mutant is defective in carotenoid production, and shows increased sensitivity to oxidative stress. Moreover, we find that the ΔcrtB mutant is impaired in biofilm formation and production of IAA. Finally we demonstrate that the ΔcrtB mutant shows reduced colonization of plant roots. Taken together, these data suggest that carotenoids are important for plant association and/or rhizosphere survival in Pantoea sp. YR343.

Two Poplar-Associated Bacterial Isolates Induce Additive Favorable Responses in a Constructed Plant-Microbiome System

The biological function of the plant-microbiome system is the result of contributions from the host plant and microbiome members. The Populus root microbiome is a diverse community that has high abundance of β- and γ-Proteobacteria, both classes which include multiple plant-growth promoting representatives. To understand the contribution of individual microbiome members in a community, we studied the function of a simplified community consisting of Pseudomonas and Burkholderia bacterial strains isolated from Populus hosts and inoculated on axenic Populus cutting in controlled laboratory conditions. Both strains increased lateral root formation and root hair production in Arabidopsis plate assays and are predicted to encode for different functions related to growth and plant growth promotion in Populus hosts. Inoculation individually, with either bacterial isolate, increased root growth relative to uninoculated controls, and while root area was increased in mixed inoculation, the interaction term was insignificant indicating additive effects of root phenotype. Complementary data including photosynthetic efficiency, whole-transcriptome gene expression and GC-MS metabolite expression data in individual and mixed inoculated treatments indicate that the effects of these bacterial strains are unique and additive. These results suggest that the function of a microbiome community may be predicted from the additive functions of the individual members.

Simultaneous knockdown of six non-family genes using a single synthetic RNAi fragment in Arabidopsis thaliana

BACKGROUND

Genetic engineering of plants that results in successful establishment of new biochemical or regulatory pathways requires stable introduction of one or more genes into the plant genome. It might also be necessary to down-regulate or turn off expression of endogenous genes in order to reduce activity of competing pathways. An established way to knockdown gene expression in plants is expressing a hairpin-RNAi construct, eventually leading to degradation of a specifically targeted mRNA. Knockdown of multiple genes that do not share homologous sequences is still challenging and involves either sophisticated cloning strategies to create vectors with different serial expression constructs or multiple transformation events that is often restricted by a lack of available transformation markers.

RESULTS

Synthetic RNAi fragments were assembled in yeast carrying homologous sequences to six or seven non-family genes and introduced into pAGRIKOLA. Transformation of Arabidopsis thaliana and subsequent expression analysis of targeted genes proved efficient knockdown of all target genes.

CONCLUSIONS

We present a simple and cost-effective method to create constructs to simultaneously knockdown multiple non-family genes or genes that do not share sequence homology. The presented method can be applied in plant and animal synthetic biology as well as traditional plant and animal genetic engineering.

Metabolic profiling reveals altered sugar and secondary metabolism in response to UGPase overexpression in Populus

BACKGROUND

UDP-glucose pyrophosphorylase (UGPase) is a sugar-metabolizing enzyme (E.C. 2.7.7.9) that catalyzes a reversible reaction of UDP-glucose and pyrophosphate from glucose-1-phosphate and UTP. UDP-glucose is a key intermediate sugar that is channeled to multiple metabolic pathways. The functional role of UGPase in perennial woody plants is poorly understood.

RESULTS

We characterized the functional role of a UGPase gene in Populus deltoides, PdUGPase2. Overexpression of the native gene resulted in increased leaf area and leaf-to-shoot biomass ratio but decreased shoot and root growth. Metabolomic analyses showed that manipulation of PdUGPase2 results in perturbations in primary, as well as secondary metabolism, resulting in reduced sugar and starch levels and increased phenolics, such as caffeoyl and feruloyl conjugates. While cellulose and lignin levels in the cell walls were not significantly altered, the syringyl-to-guaiacyl ratio was significantly reduced.

CONCLUSIONS

These results demonstrate that PdUGPase2 plays a key role in the tightly coupled primary and secondary metabolic pathways and perturbation in its function results in pronounced effects on growth and metabolism beyond cell wall biosynthesis of Populus.

Initial characterization of shade avoidance response suggests functional diversity between Populus phytochrome B genes

Shade avoidance signaling involves perception of incident red/far-red (R/FR) light by phytochromes (PHYs) and modulation of downstream transcriptional networks. Although these responses are well studied in Arabidopsis, little is known about the role of PHYs and the transcriptional responses to shade in the woody perennial Populus. Tissue expression and subcellular localization of Populus PHYs was studied by quantitative real-time PCR (qRT-PCR) and protoplast transient assay. Transgenic lines with altered PHYB1 and/or PHYB2 expression were used in phenotypic assays and transcript profiling with qRT-PCR. RNA-Seq was used to identify transcriptional responses to enriched FR light. All three PHYs were differentially expressed among tissue types and PHYBs were targeted to the nucleus under white light. Populus PHYB1 rescued Arabidopsis phyB mutant phenotypes. Phenotypes of Populus transgenic lines and the expression of candidate shade response genes suggested that PHYB1 and PHYB2 have distinct yet overlapping functions. RNA-Seq analysis indicated that genes associated with cell wall modification and brassinosteroid signaling were induced under enriched FR light in Populus. This study is an initial attempt at deciphering the role of Populus PHYs by evaluating transcriptional reprogramming to enriched FR and demonstrates functional diversity and overlap of the Populus PHYB1 and PHYB2 in regulating shade responses.

Pseudomonas fluorescens induces strain-dependent and strain-independent host plant responses in defense networks, primary metabolism, photosynthesis, and fitness

Colonization of plants by nonpathogenic Pseudomonas fluorescens strains can confer enhanced defense capacity against a broad spectrum of pathogens. Few studies, however, have linked defense pathway regulation to primary metabolism and physiology. In this study, physiological data, metabolites, and transcript profiles are integrated to elucidate how molecular networks initiated at the root-microbe interface influence shoot metabolism and whole-plant performance. Experiments with Arabidopsis thaliana were performed using the newly identified P. fluorescens GM30 or P. fluorescens Pf-5 strains. Co-expression networks indicated that Pf-5 and GM30 induced a subnetwork specific to roots enriched for genes participating in RNA regulation, protein degradation, and hormonal metabolism. In contrast, only GM30 induced a subnetwork enriched for calcium signaling, sugar and nutrient signaling, and auxin metabolism, suggesting strain dependence in network architecture. In addition, one subnetwork present in shoots was enriched for genes in secondary metabolism, photosynthetic light reactions, and hormone metabolism. Metabolite analysis indicated that this network initiated changes in carbohydrate and amino acid metabolism. Consistent with this, we observed strain-specific responses in tryptophan and phenylalanine abundance. Both strains reduced host plant carbon gain and fitness, yet provided a clear fitness benefit when plants were challenged with the pathogen P. syringae DC3000.

Comparative physiology and transcriptional networks underlying the heat shock response in Populus trichocarpa, Arabidopsis thaliana and Glycine max

The heat shock response continues to be layered with additional complexity as interactions and crosstalk among heat shock proteins (HSPs), the reactive oxygen network and hormonal signalling are discovered. However, comparative analyses exploring variation in each of these processes among species remain relatively unexplored. In controlled environment experiments, photosynthetic response curves were conducted from 22 to 42 °C and indicated that temperature optimum of light-saturated photosynthesis was greater for Glycine max relative to Arabidopsis thaliana or Populus trichocarpa. Transcript profiles were taken at defined states along the temperature response curves, and inferred pathway analysis revealed species-specific variation in the abiotic stress and the minor carbohydrate raffinose/galactinol pathways. A weighted gene co-expression network approach was used to group individual genes into network modules linking biochemical measures of the antioxidant system to leaf-level photosynthesis among P. trichocarpa, G. max and A. thaliana. Network-enabled results revealed an expansion in the G. max HSP17 protein family and divergence in the regulation of the antioxidant and heat shock modules relative to P. trichocarpa and A. thaliana. These results indicate that although the heat shock response is highly conserved, there is considerable species-specific variation in its regulation.

Genetic and physical mapping of Melampsora rust resistance genes in Populus and characterization of linkage disequilibrium and flanking genomic sequence

•  In an attempt to elucidate the molecular mechanisms of Melampsora rust resistance in Populus trichocarpa, we have mapped two resistance loci, MXC3 and MER, and intensively characterized the flanking genomic sequence for the MXC3 locus and the level of linkage disequilibrium (LD) in natural populations. •  We used an interspecific backcross pedigree and a genetic map that was highly saturated with AFLP and SSR markers, and assembled shotgun-sequence data in the region containing markers linked to MXC3. •  The two loci were mapped to different linkage groups. Linkage disequilibrium for MXC3 was confined to two closely linked regions spanning 34 and 16 kb, respectively. The MXC3 region also contained six disease-resistance candidate genes. •  The MER and MXC3 loci are clearly distinct, and may have different mechanisms of resistance, as different classes of putative resistance genes were present near each locus. The suppressed recombination previously observed in the MXC3 region was possibly caused by extensive hemizygous rearrangements confined to the original parent tree. The relatively low observed LD may facilitate association studies using candidate genes for rust resistance, but will probably inhibit marker-aided selection.

Leaf dynamics of a deciduous forest canopy: no response to elevated CO2

Leaf area index (LAI) and its seasonal dynamics are key determinants of terrestrial productivity and, therefore, of the response of ecosystems to a rising atmospheric CO(2) concentration. Despite the central importance of LAI, there is very little evidence from which to assess how forest LAI will respond to increasing [CO(2)]. We assessed LAI and related leaf indices of a closed-canopy deciduous forest for 4 years in 25-m-diameter plots that were exposed to ambient or elevated CO(2) (542 ppm) in a free-air CO(2) enrichment (FACE) experiment. LAI of this Liquidambar styraciflua (sweetgum) stand was about 6 and was relatively constant year-to-year, including the 2 years prior to the onset of CO(2) treatment. LAI throughout the 1999-2002 growing seasons was assessed through a combination of data on photosynthetically active radiation (PAR) transmittance, mass of litter collected in traps, and leaf mass per unit area (LMA). There was no effect of [CO(2)] on any expression of leaf area, including peak LAI, average LAI, or leaf area duration. Canopy mass and LMA, however, were significantly increased by CO(2) enrichment. The hypothesized connection between light compensation point (LCP) and LAI was rejected because LCP was reduced by [CO(2)] enrichment only in leaves under full sun, but not in shaded leaves. Data on PAR interception also permitted calculation of absorbed PAR (APAR) and light use efficiency (LUE), which are key parameters connecting satellite assessments of terrestrial productivity with ecosystem models of future productivity. There was no effect of [CO(2)] on APAR, and the observed increase in net primary productivity in elevated [CO(2)] was ascribed to an increase in LUE, which ranged from 1.4 to 2.4 g MJ(-1). The current evidence seems convincing that LAI of non-expanding forest stands will not be different in a future CO(2)-enriched atmosphere and that increases in LUE and productivity in elevated [CO(2)] are driven primarily by functional responses rather than by structural changes. Ecosystem or regional models that incorporate feedbacks on resource use through LAI should not assume that LAI will increase with CO(2) enrichment of the atmosphere.