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.

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.

Design and Validation of a High-Throughput Reductive Catalytic Fractionation Method.

Reductive catalytic fractionation (RCF) is a promising method to extract and depolymerize lignin from biomass, and bench-scale studies have enabled considerable progress in the past decade. RCF experiments are typically conducted in pressurized batch reactors with volumes ranging between 50 and 1000 mL, limiting the throughput of these experiments to one to six reactions per day for an individual researcher. Here, we report a high-throughput RCF (HTP-RCF) method in which batch RCF reactions are conducted in 1 mL wells machined directly into Hastelloy reactor plates. The plate reactors can seal high pressures produced by organic solvents by vertically stacking multiple reactor plates, leading to a compact and modular system capable of performing 240 reactions per experiment. Using this setup, we screened solvent mixtures and catalyst loadings for hydrogen-free RCF using 50 mg poplar and 0.5 mL reaction solvent. The system of 1:1 isopropanol/methanol showed optimal monomer yields and selectivity to 4-propyl substituted monomers, and validation reactions using 75 mL batch reactors produced identical monomer yields. To accommodate the low material loadings, we then developed a workup procedure for parallel filtration, washing, and drying of samples and a 1H nuclear magnetic resonance spectroscopy method to measure the RCF oil yield without performing liquid-liquid extraction. As a demonstration of this experimental pipeline, 50 unique switchgrass samples were screened in RCF reactions in the HTP-RCF system, revealing a wide range of monomer yields (21-36%), S/G ratios (0.41-0.93), and oil yields (40-75%). These results were successfully validated by repeating RCF reactions in 75 mL batch reactors for a subset of samples. We anticipate that this approach can be used to rapidly screen substrates, catalysts, and reaction conditions in high-pressure batch reactions with higher throughput than standard batch reactors.

Economic and Sustainability Impacts of Yield and Composition Variation in Bioenergy Crops: Switchgrass (Panicum virgatum L.).

Economically viable production of biobased products and fuels requires high-yielding, high-quality, sustainable process-advantaged crops, developed using bioengineering or advanced breeding approaches. Understanding which crop phenotypic traits have the largest impact on biofuel economics and sustainability outcomes is important for the targeted feedstock crop development. Here, we evaluated biomass yield and cell-wall composition traits across a large natural variant population of switchgrass (Panicum virgatum L.) grown across three common garden sites. Samples from 331 switchgrass genotypes were collected and analyzed for carbohydrate and lignin components. Considering plant survival and biomass after multiple years of growth, we found that 84 of the genotypes analyzed may be suited for commercial production in the southeastern U.S. These genotypes show a range of growth and compositional traits across the population that are apparently independent of each other. We used these data to conduct techno-economic analyses and life cycle assessments evaluating the performance of each switchgrass genotype under a standard cellulosic ethanol process model with pretreatment, added enzymes, and fermentation. We find that switchgrass yield per area is the largest economic driver of the minimum fuel selling price (MSFP), ethanol yield per hectare, global warming potential (GWP), and cumulative energy demand (CED). At any yield, the carbohydrate content is significant but of secondary importance. Water use follows similar trends but has more variability due to an increased dependence on the biorefinery model. Analyses presented here highlight the primary importance of plant yield and the secondary importance of carbohydrate content when selecting a feedstock that is both economical and sustainable.

Rapid determination of sugar content in biomass hydrolysates using nuclear magnetic resonance spectroscopy.

Large populations of potential cellulosic biomass feedstocks are currently being screened for fuel and chemical applications. The monomeric sugar content, released through hydrolysis, is of particular importance and is currently measured with time-consuming HPLC methods. A method for sugar detection is presented here that employs (1)H NMR spectra regressed against primary HPLC sugar concentration data to build partial least squares (PLS) models. The PLS2 model is able to predict concentrations of both major sugar components, like glucose and xylose, and minor sugars, such as arabinose and mannose, in biomass hydrolysates. The model was built with 65 samples from a variety of different biomass species and covers a wide range of sugar concentrations. Model predictions were validated with a set of 15 samples which were all within error of both HPLC and NMR integration measurements. The data collection time for these NMR measurements is less than 20 min, offering a significant improvement to the 1 h acquisition time that is required for HPLC.

Abundance of Major Cell Wall Components in Natural Variants and Pedigrees of Populus trichocarpa.

The rapid analysis of biopolymers including lignin and sugars in lignocellulosic biomass cell walls is essential for the analysis of the large sample populations needed for identifying heritable genetic variation in biomass feedstocks for biofuels and bioproducts. In this study, we reported the analysis of cell wall lignin content, syringyl/guaiacyl (S/G) ratio, as well as glucose and xylose content by high-throughput pyrolysis-molecular beam mass spectrometry (py-MBMS) for >3,600 samples derived from hundreds of accessions of Populus trichocarpa from natural populations, as well as pedigrees constructed from 14 parents (7 × 7). Partial Least Squares (PLS) regression models were built from the samples of known sugar composition previously determined by hydrolysis followed by nuclear magnetic resonance (NMR) analysis. Key spectral features positively correlated with glucose content consisted of m/z 126, 98, and 69, among others, deriving from pyrolyzates such as hydroxymethylfurfural, maltol, and other sugar-derived species. Xylose content positively correlated primarily with many lignin-derived ions and to a lesser degree with m/z 114, deriving from a lactone produced from xylose pyrolysis. Models were capable of predicting glucose and xylose contents with an average error of less than 4%, and accuracy was significantly improved over previously used methods. The differences in the models constructed from the two sample sets varied in training sample number, but the genetic and compositional uniformity of the pedigree set could be a potential driver in the slightly better performance of that model in comparison with the natural variants. Broad-sense heritability of glucose and xylose composition using these data was 0.32 and 0.34, respectively. In summary, we have demonstrated the use of a single high-throughput method to predict sugar and lignin composition in thousands of poplar samples to estimate the heritability and phenotypic plasticity of traits necessary to develop optimized feedstocks for bioenergy applications.

Comparison of methodologies used to determine aromatic lignin unit ratios in lignocellulosic biomass.

BACKGROUND: Multiple analytical methods have been developed to determine the ratios of aromatic lignin units, particularly the syringyl/guaiacyl (S/G) ratio, of lignin biopolymers in plant cell walls. Chemical degradation methods such as thioacidolysis produce aromatic lignin units that are released from certain linkages and may induce chemical changes rendering it difficult to distinguish and determine the source of specific aromatic lignin units released, as is the case with nitrobenzene oxidation methodology. NMR methods provide powerful tools used to analyze cell walls for lignin composition and linkage information. Pyrolysis-mass spectrometry methods are also widely used, particularly as high-throughput methodologies. However, the different techniques used to analyze aromatic lignin unit ratios frequently yield different results within and across particular studies, making it difficult to interpret and compare results. This also makes it difficult to obtain meaningful insights relating these measurements to other characteristics of plant cell walls that may impact biomass sustainability and conversion metrics for the production of bio-derived fuels and chemicals. RESULTS: The authors compared the S/G lignin unit ratios obtained from thioacidolysis, pyrolysis-molecular beam mass spectrometry (py-MBMS), HSQC liquid-state NMR and solid-state (ss) NMR methodologies of pine, several genotypes of poplar, and corn stover biomass. An underutilized approach to deconvolute ssNMR spectra was implemented to derive S/G ratios. The S/G ratios obtained for the samples did not agree across the different methods, but trends were similar with the most agreement among the py-MBMS, HSQC NMR and deconvoluted ssNMR methods. The relationship between S/G, thioacidolysis yields, and linkage analysis determined by HSQC is also addressed. CONCLUSIONS: This work demonstrates that different methods using chemical, thermal, and non-destructive NMR techniques to determine native lignin S/G ratios in plant cell walls may yield different results depending on species and linkage abundances. Spectral deconvolution can be applied to many hardwoods with lignin dominated by S and G units, but the results may not be reliable for some woody and grassy species of more diverse lignin composition. HSQC may be a better method for analyzing lignin in those species given the wealth of information provided on additional aromatic moieties and bond linkages. Additionally, trends or correlations in lignin characteristics such as S/G ratios and lignin linkages within the same species such as poplar may not necessarily exhibit the same trends or correlations made across different biomass types. Careful consideration is required when choosing a method to measure S/G ratios and the benefits and shortcomings of each method discussed here are summarized.

Selective One-Dimensional 13C-13C Spin-Diffusion Solid-State Nuclear Magnetic Resonance Methods to Probe Spatial Arrangements in Biopolymers Including Plant Cell Walls, Peptides, and Spider Silk.

Two-dimensional (2D) and 3D through-space 13C-13C homonuclear spin-diffusion techniques are powerful solid-state nuclear magnetic resonance (NMR) tools for extracting structural information from 13C-enriched biomolecules, but necessarily long acquisition times restrict their applications. In this work, we explore the broad utility and underutilized power of a chemical shift-selective one-dimensional (1D) version of a 2D 13C-13C spin-diffusion solid-state NMR technique. The method, which is called 1D dipolar-assisted rotational resonance (DARR) difference, is applied to a variety of biomaterials including lignocellulosic plant cell walls, microcrystalline peptide fMLF, and black widow dragline spider silk. 1D 13C-13C spin-diffusion methods described here apply in select cases in which the 1D 13C solid-state NMR spectrum displays chemical shift-resolved moieties. This is analogous to the selective 1D nuclear Overhauser effect spectroscopy (NOESY) experiment utilized in liquid-state NMR as a faster (1D instead of 2D) and often less ambiguous (direct sampling of the time domain data, coupled with increased signal averaging) alternative to 2D NOESY. Selective 1D 13C-13C spin-diffusion methods are more time-efficient than their 2D counterparts such as proton-driven spin diffusion (PDSD) and dipolar-assisted rotational resonance. The additional time gained enables measurements of 13C-13C spin-diffusion buildup curves and extraction of spin-diffusion time constants TSD, yielding detailed structural information. Specifically, selective 1D DARR difference buildup curves applied to 13C-enriched hybrid poplar woody stems confirm strong spatial interaction between lignin and acetylated xylan polymers within poplar plant secondary cell walls, and an interpolymer distance of ∼0.45-0.5 nm was estimated. Additionally, Tyr/Gly long-range correlations were observed on isotopically enriched black widow spider dragline silks.

Electrocatalytic CO2 Reduction over Cu3P Nanoparticles Generated via a Molecular Precursor Route.

The design of nanoparticles (NPs) with tailored morphologies and finely tuned electronic and physical properties has become a key strategy for controlling selectivity and improving conversion efficiency in a variety of important electrocatalytic transformations. Transition metal phosphide NPs, in particular, have emerged as a versatile class of catalytic materials due to their multifunctional active sites and composition- and phase-dependent properties. Access to targeted transition metal phosphide NPs with controlled features is necessary to tune the catalytic activity. To this end, we have established a solution-synthesis route utilizing a molecular precursor containing M-P bonds to generate solid metal phosphide NPs with controlled stoichiometry and morphology. We expand here the application of molecular precursors in metal phosphide NP synthesis to include the preparation of phase-pure Cu3P NPs from the thermal decomposition of [Cu(H)(PPh3)]6. The mechanism of [Cu(H)(PPh3)]6 decomposition and subsequent formation of Cu3P was investigated through modification of the reaction parameters. Identification and optimization of the critical reaction parameters (i.e., time, temperature, and oleylamine concentration) enabled the synthesis of phase-pure 9-11 nm Cu3P NPs. To probe the multifunctionality of this materials system, Cu3P NPs were investigated as an electrocatalyst for CO2 reduction. At low overpotential (-0.30 V versus RHE) in 0.1 M KHCO3 electrolyte, Cu3P-modified carbon paper electrodes produced formate (HCOO-) at a maximum Faradaic efficiency of 8%.

Dimerization in the absence of higher-order oligomerization of the G protein-coupled secretin receptor.

Oligomerization of G protein-coupled receptors has been proposed to affect receptor function and regulation; however, little is known about the molecular nature of such complexes. We previously utilized bioluminescence resonance energy transfer (BRET) to demonstrate that the prototypic Family B secretin receptor can form oligomers. We now explore the order of oligomerization present utilizing unique bimolecular fluorescence complementation and energy transfer techniques. The non-fluorescent carboxyl-terminal and amino-terminal halves of yellow fluorescent protein (YFP) were fused to the carboxyl terminus of the secretin receptor. These constructs bound secretin normally and signaled in response to secretin like wild type receptor. When co-expressed on COS cells, these constructs physically interacted to yield typical YFP fluorescence in biosynthetic compartments and at the plasma membrane, reflecting receptor homo-dimerization. However, the addition of another potential partner in form of Rlu- or CFP-tagged secretin receptor yielded no significant BRET or FRET signal, respectively, under conditions in which intact YFP-tagged secretin receptor yielded such a signal. Absence of higher-order receptor oligomers was further confirmed using saturation BRET techniques. Absence of significant resonance transfer to the secretin receptor homo-dimer was true for carboxyl-terminally-tagged secretin receptor, as well as for receptor incorporating the transfer partner into each of the three distinct intracellular loop domains. These results suggest that the secretin receptor can exist only as a structurally-specific homo-dimer, without being present as higher-order oligomers.

The effect of coumaryl alcohol incorporation on the structure and composition of lignin dehydrogenation polymers.

BACKGROUND: Lignin dehydrogenation polymers (DHPs) are polymers generated from phenolic precursors for the purpose of studying lignin structure and polymerization processes. METHODS: Here, DHPs were synthesized using a Zutropfverfahren method with horseradish peroxidase and three lignin monomers, sinapyl (S), coumaryl (H), and coniferyl (G) alcohols, in the presence of hydrogen peroxide. The H monomer was reacted with G and a 1:1 molar mixture of S:G monomers at H molar compositions of 0, 5, 10, and 20 mol% to study how the presence of the H monomer affected the structure and composition of the recovered polymers. RESULTS: At low H concentrations, solid-state NMR spectra suggest that the H and G monomers interact to form G:H polymers that have a lower average molecular weight than the solely G-based polymer or the G:H polymer produced at higher H concentrations. Solid-state NMR and pyrolysis-MBMS analyses suggest that at higher H concentrations, the H monomer primarily self-polymerizes to produce clusters of H-based polymer that are segregated from clusters of G- or S:G-based polymers. Thioacidolysis generally showed higher recoveries of thioethylated products from S:G or S:G:H polymers made with higher H content, indicating an increase in the linear ether linkages. CONCLUSIONS: Overall, the experimental results support theoretical predictions for the reactivity and structural influences of the H monomer on the formation of lignin-like polymers.

A thioacidolysis method tailored for higher-throughput quantitative analysis of lignin monomers.

Thioacidolysis is a method used to measure the relative content of lignin monomers bound by ß-O-4 linkages. Current thioacidolysis methods are low-throughput as they require tedious steps for reaction product concentration prior to analysis using standard GC methods. A quantitative thioacidolysis method that is accessible with general laboratory equipment and uses a non-chlorinated organic solvent and is tailored for higher-throughput analysis is reported. The method utilizes lignin arylglycerol monomer standards for calibration, requires 1-2 mg of biomass per assay and has been quantified using fast-GC techniques including a Low Thermal Mass Modular Accelerated Column Heater (LTM MACH). Cumbersome steps, including standard purification, sample concentrating and drying have been eliminated to help aid in consecutive day-to-day analyses needed to sustain a high sample throughput for large screening experiments without the loss of quantitation accuracy. The method reported in this manuscript has been quantitatively validated against a commonly used thioacidolysis method and across two different research sites with three common biomass varieties to represent hardwoods, softwoods, and grasses.

O-glycosylation effects on family 1 carbohydrate-binding module solution structures.

UNLABELLED: Family 1 carbohydrate-binding modules (CBMs) are ubiquitous components of multimodular fungal enzymes that degrade plant cell wall polysaccharides and bind specifically to cellulose. Native glycosylation of family 1 CBMs has been shown to substantially impact multiple physical properties, including thermal and proteolytic stability and cellulose binding affinity. To gain molecular insights into the changes in CBM properties upon glycosylation, solution structures of two glycoforms of a Trichoderma reesei family 1 CBM were studied by NMR spectroscopy: a glycosylated family 1 CBM with a mannose group attached to both Thr1 and Ser3 and a second family 1 CBM with single mannose groups attached to Thr1, Ser3 and Ser14. The structures clearly reveal that monosaccharides at both Ser3 and Ser14 on family 1 CBMs present additional cellulose binding platforms, similar to well-characterized aromatic residues at the binding interface, which align to the cellulose surface. These results are in agreement with previous experimental work demonstrating that glycans at Ser3 and Ser14 impart significant improvements in binding affinity. Additionally, detailed analysis of the NMR structures and molecular simulations indicates that the protein backbone of the CBM is not significantly altered by attachment of monosaccharides, and that the mannose attached to Ser14 may be more flexible than the mannose at Ser3. Overall, the present study reveals how family 1 CBM structures are affected by covalent attachment of monosaccharides, which are likely important post-translational modifications of these common subdomains of fungal plant cell wall degrading enzymes. DATABASE: Structural data have been deposited in the RCSB Protein Data Bank (PDB codes: 2MWJ and 2MWK) and the BioMagRes Bank (BMRB codes: 25331 and 25332) for CBM_M2 and CBM_M3, respectively.