JGI Plant Gene Atlas: an updateable transcriptome resource to improve functional gene descriptions across the plant kingdom.

Gene functional descriptions offer a crucial line of evidence for candidate genes underlying trait variation. Conversely, plant responses to environmental cues represent important resources to decipher gene function and subsequently provide molecular targets for plant improvement through gene editing. However, biological roles of large proportions of genes across the plant phylogeny are poorly annotated. Here we describe the Joint Genome Institute (JGI) Plant Gene Atlas, an updateable data resource consisting of transcript abundance assays spanning 18 diverse species. To integrate across these diverse genotypes, we analyzed expression profiles, built gene clusters that exhibited tissue/condition specific expression, and tested for transcriptional response to environmental queues. We discovered extensive phylogenetically constrained and condition-specific expression profiles for genes without any previously documented functional annotation. Such conserved expression patterns and tightly co-expressed gene clusters let us assign expression derived additional biological information to 64 495 genes with otherwise unknown functions. The ever-expanding Gene Atlas resource is available at JGI Plant Gene Atlas (https://plantgeneatlas.jgi.doe.gov) and Phytozome (https://phytozome.jgi.doe.gov/), providing bulk access to data and user-specified queries of gene sets. Combined, these web interfaces let users access differentially expressed genes, track orthologs across the Gene Atlas plants, graphically represent co-expressed genes, and visualize gene ontology and pathway enrichments.

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.

Agave REVEILLE1 regulates the onset and release of seasonal dormancy in Populus.

Deciduous woody plants like poplar (Populus spp.) have seasonal bud dormancy. It has been challenging to simultaneously delay the onset of bud dormancy in the fall and advance bud break in the spring, as bud dormancy, and bud break were thought to be controlled by different genetic factors. Here, we demonstrate that heterologous expression of the REVEILLE1 gene (named AaRVE1) from Agave (Agave americana) not only delays the onset of bud dormancy but also accelerates bud break in poplar in field trials. AaRVE1 heterologous expression increases poplar biomass yield by 166% in the greenhouse. Furthermore, we reveal that heterologous expression of AaRVE1 increases cytokinin contents, represses multiple dormancy-related genes, and up-regulates bud break-related genes, and that AaRVE1 functions as a transcriptional repressor and regulates the activity of the DORMANCY-ASSOCIATED PROTEIN 1 (DRM1) promoter. Our findings demonstrate that AaRVE1 appears to function as a regulator of bud dormancy and bud break, which has important implications for extending the growing season of deciduous trees in frost-free temperate and subtropical regions to increase crop yield.

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.

De Novo Assembly and Annotation of 11 Diverse Shrub Willow (Salix) Genomes Reveals Novel Gene Organization in Sex-Linked Regions.

Poplar and willow species in the Salicaceae are dioecious, yet have been shown to use different sex determination systems located on different chromosomes. Willows in the subgenus Vetrix are interesting for comparative studies of sex determination systems, yet genomic resources for these species are still quite limited. Only a few annotated reference genome assemblies are available, despite many species in use in breeding programs. Here we present de novo assemblies and annotations of 11 shrub willow genomes from six species. Copy number variation of candidate sex determination genes within each genome was characterized and revealed remarkable differences in putative master regulator gene duplication and deletion. We also analyzed copy number and expression of candidate genes involved in floral secondary metabolism, and identified substantial variation across genotypes, which can be used for parental selection in breeding programs. Lastly, we report on a genotype that produces only female descendants and identified gene presence/absence variation in the mitochondrial genome that may be responsible for this unusual inheritance.

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.

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.

The Ectomycorrhizal Fungus Laccaria bicolor Produces Lipochitooligosaccharides and Uses the Common Symbiosis Pathway to Colonize Populus Roots.

Mycorrhizal fungi form mutualistic associations with the roots of most land plants and provide them with mineral nutrients from the soil in exchange for fixed carbon derived from photosynthesis. The common symbiosis pathway (CSP) is a conserved molecular signaling pathway in all plants capable of associating with arbuscular mycorrhizal fungi. It is required not only for arbuscular mycorrhizal symbiosis but also for rhizobia-legume and actinorhizal symbioses. Given its role in such diverse symbiotic associations, we hypothesized that the CSP also plays a role in ectomycorrhizal associations. We showed that the ectomycorrhizal fungus Laccaria bicolor produces an array of lipochitooligosaccharides (LCOs) that can trigger both root hair branching in legumes and, most importantly, calcium spiking in the host plant Populus in a CASTOR/POLLUX-dependent manner. Nonsulfated LCOs enhanced lateral root development in Populus in a calcium/calmodulin-dependent protein kinase (CCaMK)-dependent manner, and sulfated LCOs enhanced the colonization of Populus by L. bicolor Compared with the wild-type Populus, the colonization of CASTOR/POLLUX and CCaMK RNA interference lines by L. bicolor was reduced. Our work demonstrates that similar to other root symbioses, L. bicolor uses the CSP for the full establishment of its mutualistic association with Populus.

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.

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.

Association mapping, transcriptomics, and transient expression identify candidate genes mediating plant-pathogen interactions in a tree.

Invasive microbes causing diseases such as sudden oak death negatively affect ecosystems and economies around the world. The deployment of resistant genotypes for combating introduced diseases typically relies on breeding programs that can take decades to complete. To demonstrate how this process can be accelerated, we employed a genome-wide association mapping of ca 1,000 resequenced Populus trichocarpa trees individually challenged with Sphaerulina musiva, an invasive fungal pathogen. Among significant associations, three loci associated with resistance were identified and predicted to encode one putative membrane-bound L-type receptor-like kinase and two receptor-like proteins. A susceptibility-associated locus was predicted to encode a putative G-type D-mannose-binding receptor-like kinase. Multiple lines of evidence, including allele analysis, transcriptomics, binding assays, and overexpression, support the hypothesized function of these candidate genes in the P. trichocarpa response to S. musiva.

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.

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

Following publication of the original article.

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.

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.

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.

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 CO2 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 C3, C4 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 C3, C4 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 C3 ancestors in Agave.

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.

Genome-wide analysis of lectin receptor-like kinases in Populus.

BACKGROUND: Receptor-like kinases (RLKs) belong to a large protein family with over 600 members in Arabidopsis and over 1000 in rice. Among RLKs, the lectin receptor-like kinases (LecRLKs) possess a characteristic extracellular carbohydrate-binding lectin domain and play important roles in plant development and innate immunity. There are 75 and 173 LecRLKs in Arabidopsis and rice, respectively. However, little is known about LecRLKs in perennial woody plants. RESULTS: Here we report the genome-wide analysis of classification, domain architecture and expression of LecRLKs in the perennial woody model plant Populus. We found that the LecRLK family has expanded in Populus to a total of 231, including 180 G-type, 50 L-type and 1 C-type LecRLKs. Expansion of the Populus LecRLKs (PtLecRLKs) occurred partially through tandem duplication. Based on domain architecture and orientation features, we classified PtLecRLKs into eight different classes. RNA-seq-based transcriptomics analysis revealed diverse expression patterns of PtLecRLK genes among leaves, stems, roots, buds and reproductive tissues and organs. CONCLUSIONS: This study offers a comprehensive view of LecRLKs in the perennial woody model plant Populus and provides a foundation for functional characterization of this important family of receptor-like kinases.

Knockdown of a laccase in Populus deltoides confers altered cell wall chemistry and increased sugar release.

Plant laccases are thought to function in the oxidation of monolignols which leads to higher order lignin formation. Only a hand-full of laccases in plants have been functionally evaluated, and as such little is known about the breadth of their impact on cell wall chemistry or structure. Here, we describe a previously uncharacterized laccase from Populus, encoded by locus Potri.008G064000, whose reduced expression resulted in transgenic Populus trees with changes in syringyl/guaiacyl ratios as well as altered sugar release phenotypes. These phenotypes are consistent with plant biomass exhibiting reduced recalcitrance. Interestingly, the transgene effect on recalcitrance is dependent on a mild pretreatment prior to chemical extraction of sugars. Metabolite profiling suggests the transgene modulates phenolics that are associated with the cell wall structure. We propose that this particular laccase has a range of functions related to oxidation of phenolics and conjugation of flavonoids that interact with lignin in the cell wall.