Providing biological context for GWAS results using eQTL regulatory and co-expression networks in Populus.

Our study utilized genome-wide association studies (GWAS) to link nucleotide variants to traits in Populus trichocarpa, a species with rapid linkage disequilibrium decay. The aim was to overcome the challenge of interpreting statistical associations at individual loci without sufficient biological context, which often leads to reliance solely on gene annotations from unrelated model organisms. We employed an integrative approach that included GWAS targeting multiple traits using three individual techniques for lignocellulose phenotyping, expression quantitative trait loci (eQTL) analysis to construct transcriptional regulatory networks around each candidate locus and co-expression analysis to provide biological context for these networks, using lignocellulose biosynthesis in Populus trichocarpa as a case study. The research identified three candidate genes potentially involved in lignocellulose formation, including one previously recognized gene (Potri.005G116800/VND1, a critical regulator of secondary cell wall formation) and two genes (Potri.012G130000/AtSAP9 and Potri.004G202900/BIC1) with newly identified putative roles in lignocellulose biosynthesis. Our integrative approach offers a framework for providing biological context to loci associated with trait variation, facilitating the discovery of new genes and regulatory networks.

Fungal elemental profiling unleashed through rapid laser-induced breakdown spectroscopy (LIBS).

Elemental profiling of fungal species as a phenotyping tool is an understudied topic and is typically performed to examine plant tissue or non-biological materials. Traditional analytical techniques such as inductively coupled plasma-optical emission spectroscopy (ICP-OES) and inductively coupled plasma-mass spectrometry (ICP-MS) have been used to identify elemental profiles of fungi; however, these techniques can be cumbersome due to the difficulty of preparing samples. Additionally, the instruments used for these techniques can be expensive to procure and operate. Laser-induced breakdown spectroscopy (LIBS) is an alternative elemental analytical technique-one that is sensitive across the periodic table, easy to use on various sample types, and is cost-effective in both procurement and operation. LIBS has not been used on axenic filamentous fungal isolates grown in substrate media. In this work, as a proof of concept, we used LIBS on two genetically distinct fungal species grown on a nutrient-rich and nutrient-poor substrate media to determine whether robust elemental profiles can be detected and whether differences between the fungal isolates can be identified. Our results demonstrate a distinct correlation between fungal species and their elemental profile, regardless of the substrate media, as the same strains shared a similar uptake of carbon, zinc, phosphorus, manganese, and magnesium, which could play a vital role in their survival and propagation. Independently, each fungal species exhibited a unique elemental profile. This work demonstrates a unique and valuable approach to rapidly phenotype fungi through optical spectroscopy, and this approach can be critical in understanding these fungi's behavior and interactions with the environment. IMPORTANCE: Historically, ionomics, the elemental profiling of an organism or materials, has been used to understand the elemental composition in waste materials to identify and recycle heavy metals or rare earth elements, identify the soil composition in space exploration on the moon or Mars, or understand human disorders or disease. To our knowledge, ionomic profiling of microbes, particularly fungi, has not been investigated to answer applied and fundamental biological questions. The reason is that current ionomic analytical techniques can be laborious in sample preparation, fail to measure all potential elements accurately, are cost-prohibitive, or provide inconsistent results across replications. In our previous efforts, we explored whether laser-induced breakdown spectroscopy (LIBS) could be used in determining the elemental profiles of poplar tissue, which was successful. In this proof-of-concept endeavor, we undertook a transdisciplinary effort between applied and fundamental mycology and elemental analytical techniques to address the biological question of how LIBS can used for fungi grown axenically in a nutrient-rich and nutrient-poor environment.

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.

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.

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.

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.

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.

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.