GWAS supported by computer vision identifies large numbers of candidate regulators of in planta regeneration in Populus trichocarpa.

Plant regeneration is an important dimension of plant propagation and a key step in the production of transgenic plants. However, regeneration capacity varies widely among genotypes and species, the molecular basis of which is largely unknown. Association mapping methods such as genome-wide association studies (GWAS) have long demonstrated abilities to help uncover the genetic basis of trait variation in plants; however, the performance of these methods depends on the accuracy and scale of phenotyping. To enable a large-scale GWAS of in planta callus and shoot regeneration in the model tree Populus, we developed a phenomics workflow involving semantic segmentation to quantify regenerating plant tissues over time. We found that the resulting statistics were of highly non-normal distributions, and thus employed transformations or permutations to avoid violating assumptions of linear models used in GWAS. We report over 200 statistically supported quantitative trait loci (QTLs), with genes encompassing or near to top QTLs including regulators of cell adhesion, stress signaling, and hormone signaling pathways, as well as other diverse functions. Our results encourage models of hormonal signaling during plant regeneration to consider keystone roles of stress-related signaling (e.g. involving jasmonates and salicylic acid), in addition to the auxin and cytokinin pathways commonly considered. The putative regulatory genes and biological processes we identified provide new insights into the biological complexity of plant regeneration, and may serve as new reagents for improving regeneration and transformation of recalcitrant genotypes and species.

Interactions with microbial consortia have variable effects in organic carbon and production of exometabolites among genotypes of Populus trichocarpa.

Poplar is a short-rotation woody crop frequently studied for its significance as a sustainable bioenergy source. The successful establishment of a poplar plantation partially depends on its rhizosphere-a dynamic zone governed by complex interactions between plant roots and a plethora of commensal, mutualistic, symbiotic, or pathogenic microbes that shape plant fitness. In an exploratory endeavor, we investigated the effects of a consortium consisting of ectomycorrhizal fungi and a beneficial Pseudomonas sp. strain GM41 on plant growth (including height, stem girth, leaf, and root growth) and as well as growth rate over time, across four Populus trichocarpa genotypes. Additionally, we compared the level of total organic carbon and plant exometabolite profiles across different poplar genotypes in the presence of the microbial consortium. These data revealed no significant difference in plant growth parameters between the treatments and the control across four different poplar genotypes at 7 weeks post-inoculation. However, total organic carbon and exometabolite profiles were significantly different between the genotypes and the treatments. These findings suggest that this microbial consortium has the potential to trigger early signaling responses in poplar, influencing its metabolism in ways crucial for later developmental processes and stress tolerance.

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.

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.

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.

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.

Newly identified sex chromosomes in the Sphagnum (peat moss) genome alter carbon sequestration and ecosystem dynamics.

Peatlands are crucial sinks for atmospheric carbon but are critically threatened due to warming climates. Sphagnum (peat moss) species are keystone members of peatland communities where they actively engineer hyperacidic conditions, which improves their competitive advantage and accelerates ecosystem-level carbon sequestration. To dissect the molecular and physiological sources of this unique biology, we generated chromosome-scale genomes of two Sphagnum species: S. divinum and S. angustifolium. Sphagnum genomes show no gene colinearity with any other reference genome to date, demonstrating that Sphagnum represents an unsampled lineage of land plant evolution. The genomes also revealed an average recombination rate an order of magnitude higher than vascular land plants and short putative U/V sex chromosomes. These newly described sex chromosomes interact with autosomal loci that significantly impact growth across diverse pH conditions. This discovery demonstrates that the ability of Sphagnum to sequester carbon in acidic peat bogs is mediated by interactions between sex, autosomes and environment.

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.

Draft Metagenome Sequences of the Sphagnum (Peat Moss) Microbiome from Ambient and Warmed Environments across Europe.

We present 49 metagenome assemblies of the microbiome associated with Sphagnum (peat moss) collected from ambient, artificially warmed, and geothermally warmed conditions across Europe. These data will enable further research regarding the impact of climate change on plant-microbe symbiosis, ecology, and ecosystem functioning of northern peatland ecosystems.

Comparative genomics analysis of drought response between obligate CAM and C3 photosynthesis plants.

Crassulacean acid metabolism (CAM) plants exhibit elevated drought and heat tolerance compared to C3 and C4 plants through an inverted pattern of day/night stomatal closure and opening for CO2 assimilation. However, the molecular responses to water-deficit conditions remain unclear in obligate CAM species. In this study, we presented genome-wide transcription sequencing analysis using leaf samples of an obligate CAM species Kalanchoë fedtschenkoi under moderate and severe drought treatments at two-time points of dawn (2-h before the start of light period) and dusk (2-h before the dark period). Differentially expressed genes were identified in response to environmental drought stress and a whole genome wide co-expression network was created as well. We found that the expression of CAM-related genes was not regulated by drought stimuli in K. fedtschenkoi. Our comparative analysis revealed that CAM species (K. fedtschenkoi) and C3 species (Arabidopsis thaliana, Populus deltoides 'WV94') share some common transcriptional changes in genes involved in multiple biological processes in response to drought stress, including ABA signaling and biosynthesis of secondary metabolites.

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.

Reporter genes confer new-to-nature ornamental traits in plants.

This manuscript has been authored by UT-Battelle, LLC under Contract No. DE-AC05-00OR22725 with the U.S. Department of Energy. The United States Government retains and the publisher, by accepting the article for publication, acknowledges that the United States Government retains a non-exclusive, paid-up, irrevocable, worldwide license to publish or reproduce the published form of this manuscript, or allow others to do so, for United States Government purposes. The Department of Energy will provide public access to these results of federally sponsored research in accordance with the DOE Public Access Plan (http://energy.gov/downloads/doe-public-access-plan).

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.

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.

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.

Novel metabolic interactions and environmental conditions mediate the boreal peatmoss-cyanobacteria mutualism.

Interactions between Sphagnum (peat moss) and cyanobacteria play critical roles in terrestrial carbon and nitrogen cycling processes. Knowledge of the metabolites exchanged, the physiological processes involved, and the environmental conditions allowing the formation of symbiosis is important for a better understanding of the mechanisms underlying these interactions. In this study, we used a cross-feeding approach with spatially resolved metabolite profiling and metatranscriptomics to characterize the symbiosis between Sphagnum and Nostoc cyanobacteria. A pH gradient study revealed that the Sphagnum-Nostoc symbiosis was driven by pH, with mutualism occurring only at low pH. Metabolic cross-feeding studies along with spatially resolved matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI) identified trehalose as the main carbohydrate source released by Sphagnum, which were depleted by Nostoc along with sulfur-containing choline-O-sulfate, taurine and sulfoacetate. In exchange, Nostoc increased exudation of purines and amino acids. Metatranscriptome analysis indicated that Sphagnum host defense was downregulated when in direct contact with the Nostoc symbiont, but not as a result of chemical contact alone. The observations in this study elucidated environmental, metabolic, and physiological underpinnings of the widespread plant-cyanobacterial symbioses with important implications for predicting carbon and nitrogen cycling in peatland ecosystems as well as the basis of general host-microbe interactions.

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.

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.