Proteomic and metabolic disturbances in lignin-modified Brachypodium distachyon.

Lignin biosynthesis begins with the deamination of phenylalanine and tyrosine (Tyr) as a key branch point between primary and secondary metabolism in land plants. Here, we used a systems biology approach to investigate the global metabolic responses to lignin pathway perturbations in the model grass Brachypodium distachyon. We identified the lignin biosynthetic protein families and found that ammonia-lyases (ALs) are among the most abundant proteins in lignifying tissues in grasses. Integrated metabolomic and proteomic data support a link between lignin biosynthesis and primary metabolism mediated by the ammonia released from ALs that is recycled for the synthesis of amino acids via glutamine. RNA interference knockdown of lignin genes confirmed that the route of the canonical pathway using shikimate ester intermediates is not essential for lignin formation in Brachypodium, and there is an alternative pathway from Tyr via sinapic acid for the synthesis of syringyl lignin involving yet uncharacterized enzymatic steps. Our findings support a model in which plant ALs play a central role in coordinating the allocation of carbon for lignin synthesis and the nitrogen available for plant growth. Collectively, these data also emphasize the value of integrative multiomic analyses to advance our understanding of plant metabolism.

Ectomycorrhizal symbiosis prepares its host locally and systemically for abiotic cue signaling.

Tree growth and survival are dependent on their ability to perceive signals, integrate them, and trigger timely and fitted molecular and growth responses. While ectomycorrhizal symbiosis is a predominant tree-microbe interaction in forest ecosystems, little is known about how and to what extent it helps trees cope with environmental changes. We hypothesized that the presence of Laccaria bicolor influences abiotic cue perception by Populus trichocarpa and the ensuing signaling cascade. We submitted ectomycorrhizal or non-ectomycorrhizal P. trichocarpa cuttings to short-term cessation of watering or ozone fumigation to focus on signaling networks before the onset of any physiological damage. Poplar gene expression, metabolite levels, and hormone levels were measured in several organs (roots, leaves, mycorrhizas) and integrated into networks. We discriminated the signal responses modified or maintained by ectomycorrhization. Ectomycorrhizas buffered hormonal changes in response to short-term environmental variations systemically prepared the root system for further fungal colonization and alleviated part of the root abscisic acid (ABA) signaling. The presence of ectomycorrhizas in the roots also modified the leaf multi-omics landscape and ozone responses, most likely through rewiring of the molecular drivers of photosynthesis and the calcium signaling pathway. In conclusion, P. trichocarpa-L. bicolor symbiosis results in a systemic remodeling of the host's signaling networks in response to abiotic changes. In addition, ectomycorrhizal, hormonal, metabolic, and transcriptomic blueprints are maintained in response to abiotic cues, suggesting that ectomycorrhizas are less responsive than non-mycorrhizal roots to abiotic challenges.

Genetic underpinnings of arthropod community distributions in Populus trichocarpa.

Community genetics seeks to understand the mechanisms by which natural genetic variation in heritable host phenotypes can encompass assemblages of organisms such as bacteria, fungi, and many animals including arthropods. Prior studies that focused on plant genotypes have been unable to identify genes controlling community composition, a necessary step to predict ecosystem structure and function as underlying genes shift within plant populations. We surveyed arthropods within an association population of Populus trichocarpa in three common gardens to discover plant genes that contributed to arthropod community composition. We analyzed our surveys with traditional single-trait genome-wide association analysis (GWAS), multitrait GWAS, and functional networks built from a diverse set of plant phenotypes. Plant genotype was influential in structuring arthropod community composition among several garden sites. Candidate genes important for higher level organization of arthropod communities had broadly applicable functions, such as terpenoid biosynthesis and production of dsRNA binding proteins and protein kinases, which may be capable of targeting multiple arthropod species. We have demonstrated the ability to detect, in an uncontrolled environment, individual genes that are associated with the community assemblage of arthropods on a host plant, further enhancing our understanding of genetic mechanisms that impact ecosystem structure.

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.

Multiplex knockout of trichome-regulating MYB duplicates in hybrid poplar using a single gRNA.

As the focus for CRISPR/Cas-edited plants moves from proof-of-concept to real-world applications, precise gene manipulation will increasingly require concurrent multiplex editing for polygenic traits. A common approach for editing across multiple sites is to design one guide RNA (gRNA) per target; however, this complicates construct assembly and increases the possibility of off-target mutations. In this study, we utilized one gRNA to target MYB186, a known positive trichome regulator, as well as its paralogs MYB138 and MYB38 at a consensus site for mutagenesis in hybrid poplar (Populus tremula × P. alba INRA 717-1B4). Unexpected duplications of MYB186 and MYB138 resulted in eight alleles for the three targeted genes in the hybrid poplar. Deep sequencing and polymerase chain reaction analyses confirmed editing across all eight targets in nearly all of the resultant glabrous mutants, ranging from small indels to large genomic dropouts, with no off-target activity detected at four potential sites. This highlights the effectiveness of a single gRNA targeting conserved exonic regions for multiplex editing. Additionally, cuticular wax and whole-leaf analyses showed a complete absence of triterpenes in the trichomeless mutants, hinting at a previously undescribed role for the nonglandular trichomes of poplar.

Towards engineering ectomycorrhization into switchgrass bioenergy crops via a lectin receptor-like kinase.

Soil-borne microbes can establish compatible relationships with host plants, providing a large variety of nutritive and protective compounds in exchange for photosynthesized sugars. However, the molecular mechanisms mediating the establishment of these beneficial relationships remain unclear. Our previous genetic mapping and whole-genome resequencing studies identified a gene deletion event of a Populus trichocarpa lectin receptor-like kinase gene PtLecRLK1 in Populus deltoides that was associated with poor-root colonization by the ectomycorrhizal fungus Laccaria bicolor. By introducing PtLecRLK1 into a perennial grass known to be a non-host of L. bicolor, switchgrass (Panicum virgatum L.), we found that L. bicolor colonizes ZmUbipro-PtLecRLK1 transgenic switchgrass roots, which illustrates that the introduction of PtLecRLK1 has the potential to convert a non-host to a host of L. bicolor. Furthermore, transcriptomic and proteomic analyses on inoculated-transgenic switchgrass roots revealed genes/proteins overrepresented in the compatible interaction and underrepresented in the pathogenic defence pathway, consistent with the view that pathogenic defence response is down-regulated during compatible interaction. Metabolomic profiling revealed that root colonization in the transgenic switchgrass was associated with an increase in N-containing metabolites and a decrease in organic acids, sugars, and aromatic hydroxycinnamate conjugates, which are often seen in the early steps of establishing compatible interactions. These studies illustrate that PtLecRLK1 is able to render a plant susceptible to colonization by the ectomycorrhizal fungus L. bicolor and shed light on engineering mycorrhizal symbiosis into a non-host to enhance plant productivity and fitness on marginal lands.

ALD1 accumulation in Arabidopsis epidermal plastids confers local and non-autonomous disease resistance.

The Arabidopsis plastid-localized ALD1 protein acts in the lysine catabolic pathway that produces infection-induced pipecolic acid (Pip), Pip derivatives, and basal non-Pip metabolite(s). ALD1 is indispensable for disease resistance associated with Pseudomonas syringae infections of naïve plants as well as those previously immunized by a local infection, a phenomenon called systemic acquired resistance (SAR). Pseudomonas syringae is known to associate with mesophyll as well as epidermal cells. To probe the importance of epidermal cells in conferring bacterial disease resistance, we studied plants in which ALD1 was only detectable in the epidermal cells of specific leaves. Local disease resistance and many features of SAR were restored when ALD1 preferentially accumulated in the epidermal plastids at immunization sites. Interestingly, SAR restoration occurred without appreciable accumulation of Pip or known Pip derivatives in secondary distal leaves. Our findings establish that ALD1 has a non-autonomous effect on pathogen growth and defense activation. We propose that ALD1 is sufficient in the epidermis of the immunized leaves to activate SAR, but basal ALD1 and possibly a non-Pip metabolite(s) are also needed at all infection sites to fully suppress bacterial growth. Thus, epidermal plastids that contain ALD1 play a key role in local and whole-plant immune signaling.

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.

Role of Reactive Oxygen Species and Hormones in Plant Responses to Temperature Changes.

Temperature stress is one of the major abiotic stresses that adversely affect agricultural productivity worldwide. Temperatures beyond a plant's physiological optimum can trigger significant physiological and biochemical perturbations, reducing plant growth and tolerance to stress. Improving a plant's tolerance to these temperature fluctuations requires a deep understanding of its responses to environmental change. To adapt to temperature fluctuations, plants tailor their acclimatory signal transduction events, and specifically, cellular redox state, that are governed by plant hormones, reactive oxygen species (ROS) regulatory systems, and other molecular components. The role of ROS in plants as important signaling molecules during stress acclimation has recently been established. Here, hormone-triggered ROS produced by NADPH oxidases, feedback regulation, and integrated signaling events during temperature stress activate stress-response pathways and induce acclimation or defense mechanisms. At the other extreme, excess ROS accumulation, following temperature-induced oxidative stress, can have negative consequences on plant growth and stress acclimation. The excessive ROS is regulated by the ROS scavenging system, which subsequently promotes plant tolerance. All these signaling events, including crosstalk between hormones and ROS, modify the plant's transcriptomic, metabolomic, and biochemical states and promote plant acclimation, tolerance, and survival. Here, we provide a comprehensive review of the ROS, hormones, and their joint role in shaping a plant's responses to high and low temperatures, and we conclude by outlining hormone/ROS-regulated plant responsive strategies for developing stress-tolerant crops to combat temperature changes.

Horizontal transfer of a pathway for coumarate catabolism unexpectedly inhibits purine nucleotide biosynthesis.

A microbe's ecological niche and biotechnological utility are determined by its specific set of co-evolved metabolic pathways. The acquisition of new pathways, through horizontal gene transfer or genetic engineering, can have unpredictable consequences. Here we show that two different pathways for coumarate catabolism failed to function when initially transferred into Escherichia coli. Using laboratory evolution, we elucidated the factors limiting activity of the newly acquired pathways and the modifications required to overcome these limitations. Both pathways required host mutations to enable effective growth with coumarate, but the necessary mutations differed. In one case, a pathway intermediate inhibited purine nucleotide biosynthesis, and this inhibition was relieved by single amino acid replacements in IMP dehydrogenase. A strain that natively contains this coumarate catabolism pathway, Acinetobacter baumannii, is resistant to inhibition by the relevant intermediate, suggesting that natural pathway transfers have faced and overcome similar challenges. Molecular dynamics simulation of the wild type and a representative single-residue mutant provide insight into the structural and dynamic changes that relieve inhibition. These results demonstrate how deleterious interactions can limit pathway transfer, that these interactions can be traced to specific molecular interactions between host and pathway, and how evolution or engineering can alleviate these limitations.

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.

Development of a clostridia-based cell-free system for prototyping genetic parts and metabolic pathways.

Gas fermentation by autotrophic bacteria, such as clostridia, offers a sustainable path to numerous bioproducts from a range of local, highly abundant, waste and low-cost feedstocks, such as industrial flue gases or syngas generated from biomass or municipal waste. Unfortunately, designing and engineering clostridia remains laborious and slow. The ability to prototype individual genetic part function, gene expression patterns, and biosynthetic pathway performance in vitro before implementing designs in cells could help address these bottlenecks by speeding up design. Unfortunately, a high-yielding cell-free gene expression (CFE) system from clostridia has yet to be developed. Here, we report the development and optimization of a high-yielding (236 ± 24 µg/mL) batch CFE platform from the industrially relevant anaerobe, Clostridium autoethanogenum. A key feature of the platform is that both circular and linear DNA templates can be applied directly to the CFE reaction to program protein synthesis. We demonstrate the ability to prototype gene expression, and quantitatively map aerobic cell-free metabolism in lysates from this system. We anticipate that the C. autoethanogenum CFE platform will not only expand the protein synthesis toolkit for synthetic biology, but also serve as a platform in expediting the screening and prototyping of gene regulatory elements in non-model, industrially relevant microbes.

Arabidopsis C-terminal binding protein ANGUSTIFOLIA modulates transcriptional co-regulation of MYB46 and WRKY33.

The apparent antagonism between salicylic acid (SA) and jasmonic acid (JA)/ethylene (ET) signalling resulting in trade-offs between defence against (hemi)biotrophic and necrotrophic pathogens has been widely described across multiple plant species. However, the underlying mechanism remains to be fully established. The molecular and cellular functions of ANGUSTIFOLIA (AN) were characterised, and its role in regulating the pathogenic response was studied in Arabidopsis. We demonstrated that AN, a plant homologue of mammalian C-TERMINAL BINDING PROTEIN (CtBP), antagonistically regulates plant resistance to the hemibiotrophic pathogen Pseudomonas syringae and the necrotrophic pathogen Botrytis cinerea. Consistent with phenotypic observations, transcription of genes involved in SA and JA/ET pathways was antagonistically regulated by AN. By interacting with another nuclear protein TYROSYL-DNA PHOSPHODIESTERASE1 (TDP1), AN imposes transcriptional repression on MYB46, encoding a transcriptional activator of PHENYLALANINE AMMONIA-LYASE (PAL) genes which are required for SA biosynthesis, while releasing TDP1-imposed transcriptional repression on WRKY33, a master regulator of the JA/ET signalling pathway. These findings demonstrate that transcriptional co-regulation of MYB46 and WRKY33 by AN mediates the coordination of SA and JA/ET pathways to optimise defences against (hemi)biotrophic and necrotrophic pathogens.

Ectopic Defense Gene Expression Is Associated with Growth Defects in Medicago truncatula Lignin Pathway Mutants.

Lignin provides essential mechanical support for plant cell walls but decreases the digestibility of forage crops and increases the recalcitrance of biofuel crops. Attempts to modify lignin content and/or composition by genetic modification often result in negative growth effects. Although several studies have attempted to address the basis for such effects in individual transgenic lines, no common mechanism linking lignin modification with perturbations in plant growth and development has yet been identified. To address whether a common mechanism exists, we have analyzed transposon insertion mutants resulting in independent loss of function of five enzymes of the monolignol pathway, as well as one double mutant, in the model legume Medicago truncatula These plants exhibit growth phenotypes from essentially wild type to severely retarded. Extensive phenotypic, transcriptomic, and metabolomics analyses, including structural characterization of differentially expressed compounds, revealed diverse phenotypic consequences of lignin pathway perturbation that were perceived early in plant development but were not predicted by lignin content or composition alone. Notable phenotypes among the mutants with severe growth impairment were increased trichome numbers, accumulation of a variety of triterpene saponins, and extensive but differential ectopic expression of defense response genes. No currently proposed model explains the observed phenotypes across all lines. We propose that reallocation of resources into defense pathways is linked to the severity of the final growth phenotype in monolignol pathway mutants of M. truncatula, although it remains unclear whether this is a cause or an effect of the growth impairment.

The genome of Eucalyptus grandis.

Eucalypts are the world's most widely planted hardwood trees. Their outstanding diversity, adaptability and growth have made them a global renewable resource of fibre and energy. We sequenced and assembled >94% of the 640-megabase genome of Eucalyptus grandis. Of 36,376 predicted protein-coding genes, 34% occur in tandem duplications, the largest proportion thus far in plant genomes. Eucalyptus also shows the highest diversity of genes for specialized metabolites such as terpenes that act as chemical defence and provide unique pharmaceutical oils. Genome sequencing of the E. grandis sister species E. globulus and a set of inbred E. grandis tree genomes reveals dynamic genome evolution and hotspots of inbreeding depression. The E. grandis genome is the first reference for the eudicot order Myrtales and is placed here sister to the eurosids. This resource expands our understanding of the unique biology of large woody perennials and provides a powerful tool to accelerate comparative biology, breeding and biotechnology.

Underground Azelaic Acid-Conferred Resistance to Pseudomonas syringae in Arabidopsis.

Local interactions between individual plant organs and diverse microorganisms can lead to whole plant immunity via the mobilization of defense signals. One such signal is the plastid lipid-derived oxylipin azelaic acid (AZA). Arabidopsis lacking AZI1 or EARLI1, related lipid transfer family proteins, exhibit reduced AZA transport among leaves and cannot mount systemic immunity. AZA has been detected in roots as well as leaves. Therefore, the present study addresses the effects on plants of AZA application to roots. AZA but not the structurally related suberic acid inhibits root growth when directly in contact with roots. Treatment of roots with AZA also induces resistance to Pseudomonas syringae in aerial tissues. These effects of AZA on root growth and disease resistance depend, at least partially, on AZI1 and EARLI1. AZI1 in roots localizes to plastids, similar to its known location in leaves. Interestingly, kinases previously shown to modify AZI1 in vitro, MPK3 and MPK6, are also needed for AZA-induced root-growth inhibition and aboveground immunity. Finally, deuterium-labeled AZA applied to the roots does not move to aerial tissues. Thus, AZA application to roots triggers systemic immunity through an AZI1/EARLI1/MPK3/MPK6-dependent pathway and AZA effects may involve one or more additional mobile signals.

Correction to: Diel rewiring and positive selection of ancient plant proteins enabled evolution of CAM photosynthesis in Agave.

Following publication of the original article.

Multitrait genome-wide association analysis of Populus trichocarpa identifies key polymorphisms controlling morphological and physiological traits.

Genome-wide association studies (GWAS) have great promise for identifying the loci that contribute to adaptive variation, but the complex genetic architecture of many quantitative traits presents a substantial challenge. We measured 14 morphological and physiological traits and identified single nucleotide polymorphism (SNP)-phenotype associations in a Populus trichocarpa population distributed from California, USA to British Columbia, Canada. We used whole-genome resequencing data of 882 trees with more than 6.78 million SNPs, coupled with multitrait association to detect polymorphisms with potentially pleiotropic effects. Candidate genes were validated with functional data. Broad-sense heritability (H2 ) ranged from 0.30 to 0.56 for morphological traits and 0.08 to 0.36 for physiological traits. In total, 4 and 20 gene models were detected using the single-trait and multitrait association methods, respectively. Several of these associations were corroborated by additional lines of evidence, including co-expression networks, metabolite analyses, and direct confirmation of gene function through RNAi. Multitrait association identified many more significant associations than single-trait association, potentially revealing pleiotropic effects of individual genes. This approach can be particularly useful for challenging physiological traits such as water-use efficiency or complex traits such as leaf morphology, for which we were able to identify credible candidate genes by combining multitrait association with gene co-expression and co-methylation data.

Comparative genomics can provide new insights into the evolutionary mechanisms and gene function in CAM plants.

Crassulacean acid metabolism (CAM) photosynthesis is an important biological innovation enabling plant adaptation to hot and dry environments. CAM plants feature high water-use efficiency, with potential for sustainable crop production under water-limited conditions. A deep understanding of CAM-related gene function and molecular evolution of CAM plants is critical for exploiting the potential of engineering CAM into C3 crops to enhance crop production on semi-arid or marginal agricultural lands. With the newly emerging genomics resources for multiple CAM species, progress has been made in comparative genomics studies on the molecular basis and subsequently on the evolution of CAM. Here, recent advances in CAM comparative genomics research in constitutive and facultative CAM plants are reviewed, with a focus on the analyses of DNA/protein sequences and gene expression to provide new insights into the path and driving force of CAM evolution and to identify candidate genes involved in CAM-related biological processes. Potential applications of new computational and experimental technologies (e.g. CRISPR/Cas-mediated genome-editing technology) to the comparative and evolutionary genomics research on CAM plants are offered.

Poplar PtabZIP1-like enhances lateral root formation and biomass growth under drought stress.

Developing drought-resistance varieties is a major goal for bioenergy crops, such as poplar (Populus), which will be grown on marginal lands with little or no water input. Root architecture can affect drought resistance, but few genes that affect root architecture in relation to water availability have been identified. Here, using activation tagging in the prime bioenergy crop poplar, we have identified a mutant that overcomes the block of lateral root (LR) formation under osmotic stress. Positioning of the tag, validation of the activation and recapitulation showed that the phenotype is caused by the poplar PtabZIP1-like (PtabZIP1L) gene with highest homology to bZIP1 from Arabidopsis. PtabZIP1L is predominantly expressed in roots, particularly in zones where lateral root primordia (LRP) initiate and LR differentiate and emerge. Transgenics overexpressing PtabZIP1L showed precocious LRP and LR development, while PtabZIP1L suppression significantly delayed both LRP and LR formation. Transgenic overexpression and suppression of PtabZIP1L also resulted in modulation of key metabolites like proline, asparagine, valine and several flavonoids. Consistently, expression of both of the poplar Proline Dehydrogenase orthologs and two of the Flavonol Synthases genes was also increased and decreased in overexpressed and suppressed transgenics, respectively. These findings suggest that PtabZIP1L mediates LR development and drought resistance through modulation of multiple metabolic pathways.