Accurate Determination of Genotypic Variance of Cell Wall Characteristics of a Populus trichocarpa Pedigree Using High-Throughput Pyrolysis-Molecular Beam Mass Spectrometry

doi:10.21203/rs.3.rs-23478/v2

Download PDF

Research

Accurate Determination of Genotypic Variance of Cell Wall Characteristics of a Populus trichocarpa Pedigree Using High-Throughput Pyrolysis-Molecular Beam Mass Spectrometry

https://doi.org/10.21203/rs.3.rs-23478/v2

This work is licensed under a CC BY 4.0 License

Journal Publication

published 06 Mar, 2021

Read the published version in Biotechnology for Biofuels and Bioproducts →

You are reading this latest preprint version

Background Pyrolysis-molecular beam mass spectrometry (py-MBMS) analysis of a pedigree of Populus trichocarpa was performed to study the phenotypic plasticity and heritability of lignin content and lignin monomer composition. Instrumental and microspatial environmental variability were observed in the spectral features and corrected to reveal underlying genetic variance of biomass composition.

Results Lignin-derived ions (including m/z 124, 154, 168, 194, 210 and others) were highly impacted by microspatial environmental variation which demonstrates phenotypic plasticity of lignin composition in Populus trichocarpa biomass. Broad-sense heritability of lignin composition after correcting for microspatial and instrumental variation was determined to be H² = 0.56 based on py-MBMS based on ions known to derive from lignin. Heritability of lignin monomeric syringyl/guaiacyl ratio (S/G) was H² = 0.81. Broad-sense heritability was also high (up to H² = 0.79) for ions derived from other components of the biomass including phenolics (e.g., salicylates) and C5 sugars (e.g., xylose). Lignin and phenolic ion abundances were primarily driven by maternal effects, and paternal effects were either similar or stronger for the most heritable carbohydrate-derived ions.

Conclusions We have shown that many biopolymer-derived ions from py-MBMS show substantial phenotypic plasticity in response to microenvironmental variation in plantations. Nevertheless, broad-sense heritability for biomass composition can be quite high after correcting for spatial environmental variation. This work outlines the importance in accounting for instrumental and microspatial environmental variation in biomass composition data for applications in heritability measurements and genomic selection for breeding poplar for renewable fuels and materials.

Biotechnology and Bioengineering

biomass composition

poplar

heritability

pyrolysis-molecular beam mass spectrometry

phenotypic plasticity

Biomass cell wall composition plays an important role in the potential use of lignocellulosic feedstocks for renewable fuels and materials. In particular, total lignin content and monomer composition can impact the technical and economic feasibility of using lignocellulosic biomass such as wood as a feedstock for biofuels or other goods [1-3]. Modifying lignin content or composition could be accomplished in several ways, including genetic engineering [4-8] or environmental priming [9, 10]. Lignin content can also be modified through breeding approaches, taking advantage of natural variation to obtain significant trait gains, using well-established protocols and experimental designs [11, 12]. Moreover, recent sophisticated techniques such as genomic selection allow for accelerated breeding with considerable reduction in costs [7, 13, 14].

Trait heritability is the cornerstone for breeding, such that a trait with no heritability is by definition unresponsive to selection, and therefore, not amenable for breeding [15]. On the contrary, traits with high heritability display large selection gains with reduced effort. However, heritablity is a not an absolute parameter, since it is codependent on the population and the experiment where it is measured [16, 17]. First, the degree of genetic control over a given trait in the population under study will determine the theoretical upper limit of heritability. Second, traits heavily affected by environmental conditions (phenotypic plasticity) will display lower heritability when measured in trials with poorly controlled environmental variation (microspatial variation) or when estimated across several trials with a range of environmental conditions (geographical variation) [17].

Accurate and affordable measurement of a trait is crucial for breeding in general and for heritability estimation in particular. Although precise methods to estimate lignin content and composition using wet chemistry are available, they are time-consuming and not practical for large sample sets [18]. High- and moderate-thoughput methods such as near-infrared spectroscopy (NIR), pyrolysis-gas chromatography (py-GC/MS), and pyrolysis-molecular beam mass spectrometry (py-MBMS) have demonstrtated the ability to rapidly measure biomass composition for a variety of studies and applications [9, 19, 20].

In particular, py-MBMS has been used as a high-throughput technique to estimate relative lignin content, lignin monomeric syringyl (S) to guaiacyl (G), (S/G) ratios, sugar composition and terpenoid content of lignocellulosic biomass [9, 10, 21-25]. The source of ions present in mass spectra derived from the pyrolysis of lignocellulosic biomass has been the focus of many investigations and has been studied on many types of instruments and scales [26-29]. Typically, py-MBMS spectra are sum (total ion chromatogram, TIC) or mean-normalized and are used to screen or elucidate cell wall composition on the basis of a small fraction of the ions present in the spectra. For example, lignin content can be estimated by regression or relative to a standard (response factor method using a single standard) by summation of ion intensities, particularly m/z 120 (H), 124 (G), 137 (G), 138 (G), 150 (H, G), 152, 154 (S), 164 (G), 167 (S), 168 (S), 178 (G), 180, 181, 182 (S), 194 (S), 208 (S) and 210 (S), where H denotes coumaryl (H) monomer derived species, G refers to species derived from coniferyl monomers and S derives from sinapyl monomers [1]. Relative S/G ratios can be estimated based on the abundance of ions known to originate from corresponding monomers. Large datasets of hundreds of samples have been analyzed in an effort to elucidate the underlying spectral data structure and its relationship to biomass composition using basic data analytics for quality control, dimension reduction and data projection, modeling and prediction. For example, principle component analysis (PCA) or discriminant analysis (DA) can be used to differentiate or group samples on the basis of py-MBMS spectra [30, 31]. Partial least squares (PLS) models have been used to predict composition based on data collected from wet chemistry and other spectroscopic methods [22, 24, 32]. One study used features in py-MBMS spectra as phenotypes to identify genes related to lignin biosynthesis in poplar using a novel multi-omic "lines of evidence" approach [33]. Ions and spectral features in the context of cell wall chemistry have also been used as traits in conjunction with genetic mapping approaches and heritability studies of poplar and other hardwoods originating from various genetic backgrounds and growing conditions [1, 3, 34]. However, thorough data analytics approaches to comprehensively analyze py-MBMS spectra from large lignocellulosic sample sets with a detailed analysis of the sources of spectral trends is lacking. As interpretation of py-MBMS spectral data in a biologically meaningful way is not trivial, particularly with respect to genetic and environmental effects, a comprehensive spectral analysis method needs to be employed for proper data interpretation.

Typically, focus has been placed on the relative abundance of a small fraction of py-MBMS ions with many assumptions in place and lack of detail outlining the variation observed in the rest of the spectra and corresponding variance attributed to instrumental drift and from metadata associated with field trials. With proper experimental design, large data sets of py-MBMS spectra within a single species consisting of various genetic or environmentally-impacted traits could potentially be efficiently analyzed to incorporate quality control, predictions of composition and properties but to also enable the generation of new hypotheses regarding the structure of the spectra as it relates to biomass properties and phenotypic plasticity.

In this study, we have analyzed the wood cell wall composition of a large full factorial 7x7 pedigree of Populus trichocarpa with a high level of technical and biological replication (n = 2721 ramets). The goal of the study was to partition the variance of the py-MBMS output given the instrumental error, the microspatial environmental heterogeneity thoughout a field trial and the genotype and familial identities of the samples. We address both the ion intensities produced by py-MBMS analysis and the compound estimations of lignin content and S/G ratio derived from appropriate ions in the spectra. This work also elucidates common trends among the ions due to microspatial environments and genotypic identity.

Assessment of instrumental error and correction of py-MBMS data

Due to the unprecedented size of the of the analysis set, quality control (QC) assessment was needed during analysis, particularly since the condition of the instrument changed between replicate analyses. Analysis of spectra from 6 types of standards monitored individually indicated no particular time-dependent trend, indicating differences were likely based on changes in instrument noise attributed to fluxuations in ion energy and the conditions of the path of pyrolysis vapors (Supplementary Table 1; Supplementary Figures 1-3). The estimated uncorrected lignin content was also reproducible between measurement replicates (Pearson correlation=0.87; R²=0.75; Supplementary Table 2). The standard deviation of lignin content between replicates ranged from 0.0–1.7 wt% lignin, all being less than 10% of the mean determined value. However, since there was minor spectral drift over the replicates of the population that was consistent among standards based on PCA and variance analysis (Supplementary Fig. 4), ions were subsequently corrected for “tray effect.” Most ions with high variance attributed to tray effects were not used in calculations for lignin composition or S/G ratio (the exceptions being m/z 167 and 181) and otherwise attributed to “noise” as outlined in the Discussion section.

Effect of microspatial environmental variation on py-MBMS spectra

We used a Thin Plate Spline (TPS) procedure to model spatial variation in py-MBMS ions in the field trial. Because of the randomized block design, genotypic effects should be randomly distributed throughout the field, so fitted values from this analysis represent environmental variation, while residuals represent genotypic effects plus error. The fitted values of the TPS models for the 421 ions displayed two distinctive patterns: simple (Fig. 1a) vs. complex surfaces (Fig. 1b). The Suface Complexity (SC) parameter was able to discriminate between these patterns, with values for simple surfaces close to 0 and values above 1 for complex surfaces. Among the 421 ions, 198 ions with null SC were free of microspatial influence (only corrected for “tray effects”), while the rest (n=223 ions) were impacted to varying degrees (Fig. 1c). The correlation between the Total Ion Chromatogram (TIC)-normalized ion intensities and the TPS residuals serves as another indication of the degree to which an ion is affected by microspatial variation (Fig. 1d). Fifteen out of the 17 ions used to quantify lignin in the spectra and all ions deriving from cell wall sugars and free phenolics had SC values in excess of 1, indicating that these cell wall components were affected by microenvironmental variation.

PCA of the ions based on their predicted TPS surface (Fig. 2), used here as a proxy for fine-scale environmental effects, yielded a PC-1 explaining 95% of the variation and PC-2 explaining 1%. The loadings for the first principal component were generally negatively correlated with lignin-derived ions with the exception of m/z 168 (primarily deriving from 4-methylsyringol) and m/z 194. When ions of TPS predicted surfaces were clustered in seven groups (Supplementary Fig. 5), the largest group (containing m/z 97) was related to phenolics and lignin-derived species, and the second largest (m/z 82) mostly consisted of sugar-derived ions and lignin dimers, again showing that these cell wall components vary spatially with the microenvironment. The rest of the clusters were small and were mostly composed of irrelevant or otherwise unannotated peaks. Peaks that are termed here as irrelevant may include noise, fragments associated with more abundant species (i.e., loss of a proton) or ions that may have many or unknown sources.

Inter-ramet variation captured in py-MBMS spectra

After TIC normalization and controlling for instrument and environmental variation, the peaks derived from cell wall components had high loadings values in PCA (Fig. 3) and were also among the most abundant and had high variance relative to the mean intensities measured across the population as shown in Figure 4a-b (Supplementary Fig. 6 shows PCA scores that are color coded corresponding to different field locations comparing before and after instrumental and environmental corrections). Approximately 120 ions were annotated based on comparisons with standards, libraries and literature where unannotated ions are either representative of an unknown component or have many sources such that their presence or source is not included or discussed (see Supplementary Table 3). The variance was highest for ions m/z 60 (C), 73 (C), 114 (C), 124 (G), 137 (G), 154 (S), 167 (S), 180 (L), 182 (S), 194 (S), and 210 (S) (where C denotes carbohydrate sugars, L for lignin, P for phenolics, G for G-lignin, S for S-lignin). These ions were also generally abundant in the spectra. However, other abundant ions such as m/z 57 (C), and 85 (C) did not exhibit as high variances relative to the mean intensity values as the former. Conversely, some ions were not particularly abundant but had high variance such as m/z 66 (P, L), 94 (P, L), 121 (P, L), and 138 (P, L).

PCA of the spatially-corrected ions also revealed negative correlations in lignin-derived ions (e.g., m/z 124, 137, 154, 210) and carbohydrate-derived ions (e.g., m/z 73, 85, 114, 126) (Fig. 3). PC-1 accounted for 41% of the spectral variation, where carbohydrate-derived ions generally were negatively correlated with lignin-derived ions. PC-2 accounted for 24% of the spectral variation, with carbohydrate and syringyl (derived from sinapyl monomers, S) ions were negatively correlated with guaiacyl (derived from coniferyl monomers, G)-derived ions (Fig. 3). Additionally, m/z 66, 94, 121 and 138 were negatively correlated with other lignin-derived species, likely indicating these ions were primarily derived from phenolics (such as phenols occurring as secondary metabolites in the case of salicylates or other lignin-like, but not true-lignin phenolics such as ferulate, coumarate, etc.) as opposed to the fragmentation of lignin-derived pyrolysates (although a positive contribution from lignin-derived analytes cannot be ruled out).

Heritability of py-MBMS spectral features

Gains in broad-sense heritability of the ions due to tray correction were marginal in most cases, though heritability of a few ions did improve noticeably with the correction (Supplementary Fig. 7). Values of broad-sense heritability for the TPS-corrected ion intensities ranged from 0 to 0.79, with annotated ions of highest heritability and noteworthiness summarized in Table 1. Permutation tests displayed thresholds of significance in heritabilities ranging from 0.028 to 0.037 for the combined tray-corrected and the TPS-corrected datasets. Although the ions with higher heritabilities were usually associated with complex surfaces for the TPS-fitted values, some ions with high heritability had simple TPS surfaces and SC values near 0 (e.g., m/z 55, 95, 167, 179, 181, 193, 195, 272, 312, and 302; Supplementary Table 4).

The most heritable ions (Table 1) also exhibited high variance in the population. However, several ions, including m/z 126 (C6 carbohydrates), 150 (G and H-lignin), 164 (G), 168 (S) 109 (P, L), 286 (G dimer) and 98 (C6 carbohydrates) were amongst the most heritable but exhibited relatively low variance. Maternal influence was almost always stronger for the most heritable ions, particularly lignin and phenolic-derived species. However, paternal effects were either more dominant or similarly influential as maternal effects for ions derived from carbohydrate sugars such as m/z 73, 97, 114 (see Supplementary Table 4 for full comparison of paternal and maternal variance associated with each ion and Supplementary Figure 8 for the % variance explained by mother vs father annotated by ion origin in biomass).

Hierarchical clustering (HC) using Spearman’s rank correlation distance metric with the complete linkage criterion was used to analyze the clustering of ions in combined tray-TPS corrected spectra to elucidate spectral associations based mostly on genetic information. Eight groups were elucidated in the spectra based on K-means clustering (Supplementary Fig. 9), summarized in Table 2 (full spectral groups outlined in Supplementary Table 5). Groups separated based on biocomponent sources similarily as those in the only TPS-corrected ions (Supplementary Fig. 10), indicating the majority of ions impacted by microspatial environment, and not ions highly impacted by instrumental variation, also were impacted by genetic variation of the population. Ions in the complete tray-TPS corrected spectra generally clustered according to biopolymer source although unannotated and noise ions appeared in all clusters to some degree. Interestingly, the most heritable ions (m/z 66, 94, 121, 138), which are produced from phenolics (possibly including salicylate-like metabolites known to occur in Populus [35-39]), were clustered together in cluster EK0 along with some lignin-derived species, including lignin dimer moieties (m/z 272, 286). The rest of the most heritable ions clustered according to their biocomponent source in clusters EK4 (G-lignin), EK5 (carbohydrate sugars), and EK6 (S-lignin) (Supplementary Table 5).

Familial patterns of ramets

Clustering of the samples based on the genotypic predicted values for py-MBMS spectra revealed some of the underlying family structure present in this population. PCA shows some differentiation of maternal half-sib families (Fig. 5a). The half-sib family from female/maternal ID 1950 (See Supplementary Table 6 for additional identifier information for each parent) in the lower right quadrant of the PCA scores plot had lower S/G and lower lignin composition in comparison to the half-sib family from female 4593 in the upper left quadrant of the plot (also see Supplementary Fig. 11). Clustering by Ward’s method using Euclidean distance revealed 7 clusters (Supplementary Fig. 12), where samples were previously classified into these groups based on K-means clustering meant to elucidate at least 7 different families (Supplementary Fig. 13). These groups largely corresponded to the maternal half-sib families (clusters colored in Fig. 5) as opposed to paternal half-sib families. Interestingly, one group of siblings (from female/maternal ID 1950) also produced the highest abundance of ion m/z 94 (PCA color coded in Supplementary Fig. 14), which can come from lignin but is otherwise attributed to the presence of phenolics such as salicylates (and in this case is otherwise not correlated positively with other lignin ion abundances such as m/z 210 as described previously (Supplementary Fig. 15 for example)). The clustering of the samples in conjunction with PCA of the spectra show that compositional relationships of several families can be elucidated, particularly based on the abundance of lignin and phenolic species. In some cases, several families generally produced similar or overlapping spectral features and hence appeared to have similar biomass composition and hence differences were not noted. Additionally, clustering methods differentiating PCA projections or scores based on MBMS spectra may identify spectral groups that consist of many family types and members and don’t necessarily separate families based on spectra. Either way, the use of clustering methods in combination with PCA of the py-MBMS spectra enable the visualization and validation of the compositional relationships within or across families that are able to be captured in the spectra.

Analysis and heritability of cell wall traits from corrected py-MBMS spectra

The average lignin content of the entire population (n = 2721), after correcting for microspatial variation of the genotypes, was 25.5% (after taking replicate averages for each sample into account), ranging from 20.9-27.9% (Table 3; Fig. 6a). The S/G ratio ranged from 1.56 to 2.77, with an average of 2.10 (Table 3; Fig. 6b). These lignin metrics are typical for variants of P. trichocarpa as previously determined in Muchero et al. [1]. The broad-sense heritability of lignin composition based on TPS-corrected values was 0.56 and the heritability of S/G was 0.81.

Instrumental variance of py-MBMS spectral features

Variance associated with sample heterogeneity and changes within the analytical equipment (instrumental drift) are assumed to be the source of the variance that was removed by correcting for “tray effects.” Ions that were highly affected by the “tray variable” were mostly fragment ions typically related to a more abundant parent pyrolysate and therefore suspected to be more affected by the conditions of the instrumentation. Tray corrected ions consisted primarily of ions > 200 m/z which are generally unanottated, low abundance ions with many typically considered “noise.” Also, a significant amount of the ions, particularly those < 200 m/z, with high instrumental variance were also odd-numbered (as opposed to even), indicating their primary source as a fragment ion. For example, m/z 194 is known to originate from 2,6-dimethoxy-4-(2-propenyl) phenol, an abundant pyrolysate derived from a sinapyl moiety in the lignin polymer [9, 27] did not need correction based on “tray effects” but had significant effects from microspatial environment. Ion 193, being related to 194 after loss of a proton, had large tray effects and low associations with the microspatial environment.

Genetic and environmental influences on py-MBMS spectral features

Microspatial and genetic variance component analyses on the MBMS data separated the ions into several distinct groups, of which the following were noted: i) ions not affected by microspatial environmental variation with no significant heritability, ii) ions not affected by microspatial environmental variation with high heritability, iii) ions slightly affected by microspatial environmental variation with high heritability and iv) ions moderately affected by microspatial environmental variation with moderate heritability, where most of the ions present in groups i-iii were tray corrected as described in the previous section. Ions present in group i could mostly be attributed to “noise” and are otherwise not important for trait analyses and breeding approaches. The ions with high heritability but not affected by fine-scale environment (group ii) consisted mostly of ions that fragment from other highly heritable ions, including ions 181 (being related to 180 and 182) and 167 (an ion produced from a wide variety of S-lignin-derived species [9, 27] that are otherwise also heritable). The heritable and environmentally stable ions still need to be corrected for instrumental drift as they may be particularly sensitive to the conditions of the pyrolyzer impacting fragmentation (e.g., the “parent” ions such as 182 and 194 would have environmental correction and subsequent fragments 181 and 167 (respectively) would require instrumental correction). Otherwise, some peaks that were not impacted by microspatial environment but were impacted by instrumental variance and were highly heritable are noteworthy based on their roles in lignin structure. For example, m/z 272 and 302, believed to originate from G-G and S-G dimers in lignin [27], were heritable after tray correction but did not have high impact from environmental conditions, indicating low phenotypic plasticity for these traits. The 223 ions that were affected by microspatial environment (groups iii and iv), and subsequently TPS-corrected, included most of the ions known to originate from cell wall polymers including cellulose (e.g. m/z 126) , hemicelluloses (e.g. m/z 114) and lignin (such as m/z 124, 154, 168, 208) all together primarily occurring between m/z 50 and 210 (see Supplementary Table 3 for annotations). Therefore, ions impacted by the environment originating from the cell wall components were necessarily corrected for accurate heritability measurements and exhibited high phenotypic plasticity. After the microspatial environmental and instrumental corrections for the ions in the spectra were made, the genotypic variation associated with the fragmented ions is more likely the result of variation in the structure of the cell wall biopolymers (i.e., lignin linkages) and presence of free phenolic metabolites. However, more highly resolved information would be necessary to derive the specific biopolymer differences in heritable spectral features due to the complicated processes associated with thermal decomposition and fragmentation of bio-derived species during pyrolysis and electron ionization.

Clustering of the ions on the TPS-predicted values identified groups of ions with similar responses to microspatial variation in the field trial. The two largest groups that represented the principal components of the TPS-predicted values primarily corresponded to lignin and phenolic-derived ions (FK1) and lignin dimer as well as sugar-derived ions (FK6). The patterns of spatial variation in predicted values (e.g., Fig. 1b) roughly corresponds to areas of the field that had high water content and standing water during the winter (personal observation). PC-1, accounting for 95% of the variation in the data, may therefore reflect this variation in soil water content. Therefore, most ions were affected by this single environmental factor, but in opposing directions: FK1 ions (phenolics, G and S lignin) had lower values in the areas under excess water, whereas FK6 (sugar and lignin dimers) had higher values with excess water. The other five clusters were driven by other minor PCs, and no clear pattern could be inferred. Independent clustering on TPS-corrected values resulted in groups that were largely congruent with those identified by the TPS-predicted values. However, clustering of ions based on the TPS-corrected values should be driven by intrinsic factors of each ramet, including genetics, and was hence more effective at resolving genotypes on the basis of lignin monomer composition. Thus, ions in the same cluster are likely to share some of the same underlying genetic control, and may be derived from related pathways. This may provide a valuable clue about the identities of some of the unknown ions in these clusters warranting further investigation on their sources.

Sources of ions that were highly heritable and impacted by microspatial variation

Hierarchical clustering of the ions showed groups of ions relating to certain cell wall components and Ward’s clustering of the samples was able to differentiate the samples based on familial relationships that were related by corresponding cell wall compositional differences. Lignin-derived ions were highly heritable, particularly those originating from the sinapyl (S) monomer present in the lignin polymer. For example, m/z 210, 194, 182, 167, 154 are all ions present in the spectra of sinapyl alcohol consisting of fragments corresponding to 2,6-dimethoxy-4-(2-propenyl)phenol (m/z 194), 4-ethylsyringol (m/z 167), and syringol (m/z 154), each of which are subsequently produced from various syringyl monomers present in lignin [40, 41]. The ion m/z 208 originates from sinapylaldehyde which can be generated upon pyrolysis of S-lignins but if present in native lignin, can have dramatic implications relating to genetic and recalcitrance properties [42]. Ions 182 and 181 are due to the presence of syringaldehyde, which can be produced upon pyrolysis of syringyl lignin but may also be present in native lignin structures [42]. Ions 124 and 137 derive from a great number of guaicyl moieties (from coniferyl, G, monomer) in the lignin polymer whereas 164 derives primarily from eugenol which was likely produced upon pyrolysis of labile (i.e., b-O-4) guaiacyl units in the lignin [43]. Other guaiacyl-related ions consisted of m/z 151 and 152, indicating the presence of vanillins which are produced from the pyrolysis of guaiacyl lignins but may also be present in the native lignin polymers [44]. Ferulic acid and 4-vinylguaiacol are likely the sources of m/z 150 and 135. There are many S- and G-derived pyrolysates attributing to the presence of m/z 168 and 180, hence their high abundance in the spectra and high variance. Interestingly, m/z 272 is likely derived from a guaiacyl stilbene dimer that has been observed from the pyrolysis of lignins [27, 43] but has recently been identified as a monomer in endocarp lignins [45]. Additionally, since the samples were analyzed without an extraction of low-molecular weight phenolics, the heritability of phenolic-derived ions and the relationship between phenolic and lignin derived ions was observable. One set of half-siblings in particular, exhibited low S/G and low lignin abundance with higher production of phenolic/salicylate derived ions.

The most heritable and most variable ions originating from sugars were generated from many sugar-derived species including m/z 60 (derived from hydroxyacetaldehyde and acetate) and 73 (3-hydroxypropionaldehyde), each of which are also fragments of levoglucosan produced upon the pyrolysis of sugars [27]. Otherwise, m/z 114 which is particularly produced from C5 sugars (i.e., xylose [24, 27]), was amongst the most heritable carbohydrate-derived ions. C6 sugars (i.e., glucose from cellulose [24, 27]) are known to produce m/z 126, which was also found to be heritable but to a lesser degree than the C5-derived m/z 114. These sugar-derived ions were driven similarly by both maternal and paternal effects, but were amongst those of highest variance from paternal effects.

Maternal effects

Ions originating from phenolics and salicylate-like species (and as fragments from lignin), including m/z 66 (P, L), 94 (P, L), 121 (P), and 138 (P, L), were primarily driven by maternal parent-of-origin effects and these were among the most heritable as well as the most variable ions produced from the population.

Parent-of-origin effects are when offspring phenotypes are not in line with an additive pattern of inheritance of the two progenitors where deviation is mostly biased to either the female progenitors or to the male progentors. For example, chill tolerance in cucumber (Cucumis sativus) displays strong maternal effects (when the phenotype of the mother prevails over the father [46]) and seems to be linked to chloroplast elements [47, 48]. Parent-of-origin effects have been shown to have important ecological and evolutionary implications (for example, genomic imprinting where the silencing or overexpression of alleles is based on the sex of parent of origin, review provided by [49]). Given the asymmetric dispersal distance between ovules (maternal) and pollen (paternal), maternal effects seem to play a role in local adaptation of offspring to environments that normally are more similar to that of the mother than the father [50]. Also, parent-of-origin effects have been implicated in interspecific reproductive barriers [51, 52] with a high rate of diversifying evolution among plant species through genomic imprinting [53]. However, imprinted expression is conserved in other genes across plant lineages, suggesting stabilizing selection for these loci [53, 54].

Despite the apparent evolutionary and ecological importance of parent-of-origin effects and their common occurrence in quantitative genetics studies of many traits of diverse plants [55], further investigation of causality is lacking. Furthermore, most in-depth studies have focused on the early stage of plant life cycle (endosperm and seed development; [46, 56, 57] and not much is known about genomic-imprinting and organelle effects in later stages of life in plants. In this study we show how parent-of-origin effects (mostly maternal) have an important role in wood chemistry profiles. This may suggest that wood composition could play important roles in environmental adaptation in plants.

Implications in breeding and genomic selection

The high level of accuracy and reproducilibity of py-MBMS instrumentation showcased by this experiment and the high heritability of many biologically relevant ions provide many possibilities for the use of this analytical technique in tree breeding. First, precise and rapid estimation of lignin content and S/G ratio allowed for large sample size (i.e., breeding population or more replication) and the reduction of error, that translate immediately to gains in selection efficiency. The broad-sense heritability values observed in this study for total lignin (0.56) and S:G (0.81) were comparable to or higher than those observed in previous studies using similar methodology (range of 0.23-0.58 and 0.42-0.81 for lignin and S/G respectively; Supplementary Table 7) [3, 58-60]. Second, this method and data treatment also permits larger training populations to create models for genomic selection with a lower level of uncertainty. Finally, recent research has pointed towards the integration of phenomics layers into the genomic selection process in order to increase predictive power of the models [61], making use of the genetic correlation between phenotypes and also being able to improve associations between the genome and the traits. This latter approach takes advantage of the normally simpler genetic architecture of some phenotypes (such as specific pyrolysates or metabolites) to detect true positive associated loci. In this sense, in order to select for a trait, other layers of phenotypes are added to help to establish true correlations between the genome and the phenotype of interest. Among the many approaches that are under investigation currently, the use of selection indices seems promising [62]. Selection index approaches could be used to select for several phenotypes simultaneously, but also to integrate other phenotypes in the selection process that are not of ultimate interest but due to high correlation to the phenotype of interest and high heritability can improve accuracy substantially.

Comprehensive analysis of a large pedigree of Populus trichocarpa was streamlined in a high-throughput py-MBMS analysis and data analytics framework that elucidated variation and heritability of biomass composition on the basis of spectral features. Ions that were sensitive to micro-spatial variation originated from phenolics and cell wall biopolymers such as lignin and hemicelluloses and cellulose whereas ions that were primarily impacted by instrumental variance were primarily attributed to noise, fragments and unannotated ions. The most variable and heritable ions originated from phenolics and S-derived lignin monomers with fewer highly heritable ions originating from G-derived monomers and carbohydrates within the cell wall. Maternal effects were generally higher for most heritable ions, particularly those derived from phenolics and lignin whereas paternal effects were similar to maternal or more impactful for carbohydrate-derived ions.

Progeny trial design and sample collection

Seven females and seven males of black cottonwood (Populus trichocarpa) were collected from naturally occurring stands in Oregon and Washington, USA and were selected based on extremes and intermediate values in lignin content and S/G ratio phenotypes. These trees were grown in a common garden in Westport, OR on the lower Columbia River. These trees were crossed to generate 49 full-sib families in a full factorial design (i.e., full factorial; 7×7), parent identifier information is provided in Supplementary Table 6. Additionally, seven open-pollinated families were created from the seven mothers. In total, 986 offspring were obtained and propagated. Three clonal replicates of these progeny plus the fourteen parents (n = 1000) were planted in April of 2015 in three complete, randomized blocks in Westport, OR. The trial consisted of 30 rows with 100 trees per row and the trial was surrounded by a double border consisting of extra ramets from the same crosses (Supplementary Figure 16). Spacing was 3 m between rows and 1 m within rows. The location of the trial (46.130742° N, 123.370276° W) falls within the core of the natural black cottonwood distribution, with favorable climatic and edaphic conditions for this species.

In January 2018, wood cores were collected at breast height (1.3 m) from the southern face of the trunks of all live ramets (n = 2721). The cores were stored in envelopes and dried at 50°C for 72 hr. They were then ground in a Wiley Mini Cutting Mill (Thomas Scientific) and filtered through a 20 mesh sieve to 1 mm for py-MBMS analysis.

Py-MBMS experiments

Py-MBMS was performed according to previously described methods [9, 18, 23, 63]. Briefly, 4 mg of debarked, dried and ground material (1 mm mesh) was pyrolyzed in a Frontier PY2020 at 500°C for 30 seconds. Pyrolysis vapors were analyzed using an Extrel Super Sonic MBMS Model Max 1000 and spectra were collected at 17 eV from m/z 30-450. MBMS data was processed using Merlin 3.0 software to produce average spectrum acquired during sample pyrolysis (a total ion chromatogram (TIC) peak over approximately 0.5 min at 0.5 s/scan rate) and spectral analysis was performed on (TIC) sum-normalized data as described in the following sections. Almost all samples of ground wood were run twice (n = 5440). Additionally, six standards were run throughout (n = 509) for quality control monitoring. Runs were carried out in a total of 125 forty-eight-vial trays. Samples were analyzed in one replicate prior to instrument cleaning and retuning and the second replicate (both sets in random order) was collected after cleaning to validate the reproducibility of the spectra. Lignin content was estimated by summing m/z 120 (H), 124 (G), 137 (G), 138 (G), 150 (H, G), 152, 154 (S), 164 (G), 167 (S), 168 (S), 178 (G), 180, 181, 182 (S), 194 (S), 208 (S) and 210 (S), where H denotes coumaryl (H) monomer derived species, G refers to species derived from coniferyl monomers and S derives from sinapyl monomers. S/G was determined by dividing sum of S-derived ions (154, 167, 168, 182, 194, 208, 210) by the sum of G-derived ions (124, 137, 138, 150, 164, 178). Coumaryl (H) content was not determined based on known overlaps with coumarate-derived ions. Ion annotations were made based on [9, 24, 27, 64], comparisons to NIST databases (online as there is not one available for MBMS software) and by comparisons to standards analyzed in-house (spectral information available upon request).

Correction of py-MBMS instrumental drift (tray effect)

This uniquely large sample set from a population developed for quantitative genetic analyses required the development of a data analysis workflow that encompassed QC and comprehensive spectral analysis. To assess and correct for the batch effect due to tray, the following linear mixed model was run for every single ion for the six replicated standards:

where, y_i,s,t is the TIC-normalized intensity of the ion i in the tray t for the standard s, μ is the overall mean, Tray_t is the tray, Standard_s is the identitity of the standard and Ɛ_i,t,s is the error. Both tray and standards were treated as random effects. The effect of each tray was used to correct the spectrum reads of each ion in the whole dataset (henceforth tray-corrected dataset). For one tray that lacked standards, a model was run on the whole dataset just having tray as random effect. These models were solved with R package lme4.

Thin plate spline model for microspatial variation analysis and correction of py-MBMS spectra

The means of each ramet of the spectral values for each ion were regressed to their spatial position in the field trial using a Thin plate spline model (TPS). This model estimates a smoothed surface, in our case matching with the two-dimensional spatial layout of the trial, but minimizing the variance of the residuals. Therefore, the predicted values of the model serve as an estimation of how the fine-scale environmental variation affects the ion intensity value, and the residuals quantify the variation that is attributed to other factors, such as genotype. Thus, we used the model to: (i) pinpoint which ions were not affected by fine-scale environment; those ions could be putative noise (i.e., no biological significance), metabolites under very strong genetic control, or stochastic or hyperresponsive metabolites, (ii) cluster ions based on their response to the environment and (iii) correct for the effects of fine-scale environmental variation. When the TPS model of an ion predicts a simple surface, it implies that the ion is not affected by fine-scale environmental variation (Fig. 1a). Complete lack of association between environment and phenotype is unlikely. Therefore, this criterion could be used as a line of evidence to identify noise. However, linkage between environment and phenotype could be obscured for metabolites that are hyper-responsive to the environment, or which vary stochastically. Finally, metabolites under strong genetic control would not be expected to vary spatially. The complexity of the TPS-fitted values was assessed with a ad-hoc parameter, dubbed as Surface complexity (SC), employing the following formula:

Predicted values in simple surfaces vary linearly within row, so the absolute value of the correlation between the predicted value and the tree position within row is almost one. The product of this correlation across all the rows should therefore be close to one, and the log of this value should approach zero. Conversely, the higher complexity the predicted surface has, the lower the correlations will be and the higher the parameter SC will be. The R package fields [65] was used to fit the TPS models.

To group the ions based on response to the environment, the predicted values from the TPS model were submitted to a Principal Component Analysis (PCA) and to hierarchical clustering (HC) as described below. Finally, the residuals of the models were used as estimators of the variation of the ions due to factors not spatially driven (e.g., genotype). The residuals for each were added to the ion mean value to project them to the original scale (henceforth TPS-corrected dataset). To study common non-enviromental trends among ions, PCA on standardized values and HC were carried out as described below.

Broad-sense heritability

Broad-sense heritabilities were calculated on the TIC-normalized, combined tray-corrected and the TPS-corrected full spectra and TPS-corrected values of lignin and S/G ratios determined prior to any spectral corrections. These variance components were estimated using a linear mixed model having ion intensities as response variable and genotype as a random effect by means of the R package lme4. Heritability was calculated from the variance explained by genotype divided by the total variance (i.e., the sum of the variance explained by genotype plus the variance due to residuals of the model). To assess significant departures from null values, we performed permutation tests on the heritability of each ion. We ran 1000 permutations for each ion setting significance threshold α<0.05. Maternal and paternal effects were estimated in a similar way, adding mother and father as random effects to the aforementioned model. For example, maternal effects were estimated by the proportion of the variance explained by the mother relative to the total variance (sum of variance explained by genotype, mother, father and residuals).

Multivariate analytical procedures

Descriptive statistical analysis (i.e., mean, variance, standard deviation), principal component analysis (PCA) and clustering were used to explore spectral data with the software R [66] and Unscrambler X V.10.5 (Camo Software). PCA was performed using 100 iterations of the NIPALS algorithm, with 6 or 7 principal components depending on convergence of data, 20 random cross validation segments, and mean-centered data. Hierarchical clustering by complete linkage using Spearman’s rank correlation distance was performed for clustering of ions in spectra and Ward’s method using Euclidean distance for clustering of biomass samples. K-means clustering of ions and samples (ramets) was performed in R using the packages cluster and factoextra.

C6 – carbohydrates with 6 carbon sugars; C5 – carbohydrates with 5 carbon sugars; HC- hierarchical clustering; PCA – principle component analysis; py-MBMS – pyrolysis-molecular beam mass spectrometry; SC – surface complexity; S/G – syringyl/guaiacyl ratio; TPS – thin plate spline

Table 1. Most heritable and informative ions with corresponding sources observed in py-MBMS spectra of 7x7 P. trichocarpa pedigree. Sources obtained based on analysis of standards and from comparison to NIST online databases and also from [9, 18, 27].

Ion (m/z)	Source	Heritability
94	Phenolics, lignin	0.79
138	Phenolics, G-lignin	0.76
66	Phenolics, lignin	0.74
121	Phenolics, lignin	0.74
167*	S-lignin	0.72
182	S-lignin	0.70
210	S-lignin	0.68
181*	S-lignin	0.65
194	S-lignin	0.64
154	S-lignin	0.59
208	S-lignin	0.58
124	G-lignin	0.56
150	G and H-lignin	0.55
164	G-lignin	0.53
137	G-lignin	0.53
168	S-lignin	0.48
180	Lignin	0.47
73	Carbohydrates	0.45
109	Phenolics, lignin	0.45
60	Acetyl	0.45
114	C5 carbohydrates	0.44
272*	G-G lignin dimer	0.41
286	Lignin dimer	0.38
151	G-lignin	0.36
57	Carbohydrates	0.37
126	C6 carbohydrates	0.32
98	C6 carbohydrates	0.32
85	C5 carbohydrates	0.31

*denotes surface complexity < 1

Table 2. Classes of clusters of ions (complete spectra of Thin Plate Spline and tray corrected ions, clustered by Spearman’s rank correlation complete linkage heiracrchal clustering) annotated for cell wall composition relevance.

Cluster class (EK#)	Number of ions	Annotations
0	79	Phenolics, lignin dimers
1	103	Lignin dimers, large mass noise
2	40	Low and moderate mass noise
3	32	Low and high mass noise
4	58	G-lignin
5	50	Carbohydrates
6	34	S-lignin
7	25	Lignin, carbohydrate, and phenolic fragments, moderate mass noise

Table 3. Summary of composition analysis for Populus trichocarpa pedigree samples based on Thin Plate Spline corrected data, parenthesis indicate standard deviation.

Populus family analysis	Population (n=2721)	Heritability (H²)
Average lignin content (%)	25.5 (± 0.9)	0.56
Average S/G ratio	2.10 (± 0.17)	0.81

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Competing interests

The authors declare that they have no competing interests.

Funding

This work was authored in part by Alliance for Sustainable Energy, LLC, the manager and operator of the National Renewable Energy Laboratory for the U.S. Department of Energy (DOE) under Contract No. DE-AC36-08GO28308. Support at NREL was provided in part by the BioEnergy Science Center (BESC) and the Center for Bioenergy Innovation (CBI) from the U.S. Department of Energy Bioenergy Research Centers (BER) supported by the Office of Biological and Environmental Research in the DOE Office of Science. Additional funding was provided by the Departments of Transportation and Department of Energy’s Western Sun Grant program (T0013A). The views expressed in the article do not necessarily represent the views of the DOE or the U.S. Government. The U.S. Government retains and the publisher, by accepting the article for publication, acknowledges that the U.S. Government retains a nonexclusive, paid-up, irrevocable, worldwide license to publish or reproduce the published form of this work, or allow others to do so, for U.S. Government purposes.

Author’s contributions

AEHW operated py-MBMS, analyzed and interpreted py-MBMS data, contributed text, figures and editing. DMS collected the cores, performed data analysis and contributed text, figures and editing. CRA collected the cores, processed samples and provided consultation. CD processed samples and operated py-MBMS. KH performed breeding crosses and managed the field trial. GT provided oversight and editing. BS designed, generated and managed the crosses. SPD designed the experiment, collected the cores, supervised wood grinding and edited the manuscript. MFD provided oversight and editing of the manuscript. All authors read and approved of the final manuscript.

Acknowledgments

The authors would like to thank Sean Cornell for his assistance in preparation of the biomass samples as well as Ambarish Nag and Alida Gerristen for helpful discussion regarding data analysis. We would also like to thank Jason Mack, Rich Shuren, Austin Himes and Carlos Gantz of Greenwood Resources for their help in establishing and maintaining the field trial.

Muchero W, Guo J, DiFazio SP, Chen J-G, Ranjan P, Slavov GT, et al. High-resolution genetic mapping of allelic variants associated with cell wall chemistry in Populus. BMC Genomics. 2015;16(1):24.
Tuskan GA, Muchero W, Tschaplinski TJ, Ragauskas AJ. Population-level approaches reveal novel aspects of lignin biosynthesis, content, composition and structure. Current Opinion in Biotechnology. 2019;56:250-7.
Ohlsson JA, Hallingbäck HR, Jebrane M, Harman-Ware AE, Shollenberger T, Decker SR, et al. Genetic variation of biomass recalcitrance in a natural Salix viminalis (L.) population. Biotechnology for Biofuels. 2019;12(1):135.
Sykes RW, Gjersing EL, Foutz K, Rottmann WH, Kuhn SA, Foster CE, et al. Down-regulation of p-coumaroyl quinate/shikimate 3′-hydroxylase (C3′H) and cinnamate 4-hydroxylase (C4H) genes in the lignin biosynthetic pathway of Eucalyptus urophylla × E. grandis leads to improved sugar release. Biotechnology for Biofuels. 2015;8(1):1-10.
Edmunds CW, Peralta P, Kelley SS, Chiang VL, Sharma-Shivappa RR, Davis MF, et al. Characterization and enzymatic hydrolysis of wood from transgenic Pinus taeda engineered with syringyl lignin or reduced lignin content. Cellulose. 2017;24(4):1901-14.
Biswal AK, Hao Z, Pattathil S, Yang X, Winkeler K, Collins C, et al. Downregulation of GAUT12 in Populus deltoides by RNA silencing results in reduced recalcitrance, increased growth and reduced xylan and pectin in a woody biofuel feedstock. Biotechnology for Biofuels. 2015;8(1):41.
Chanoca A, de Vries L, Boerjan W. Lignin Engineering in Forest Trees. Frontiers in Plant Science. 2019;10(912).
Wang JP, Matthews ML, Williams CM, Shi R, Yang C, Tunlaya-Anukit S, et al. Improving wood properties for wood utilization through multi-omics integration in lignin biosynthesis. Nature Communications. 2018;9(1):1579.
Sykes R, Yung M, Novaes E, Kirst M, Peter G, Davis M. High-Throughput Screening of Plant Cell-Wall Composition Using Pyrolysis Molecular Beam Mass Spectroscopy. In: Biofuels: Methods and Protocols. Mielenz JR, editor. Totowa, NJ: Humana Press; 2009. p. 169-83.
Mann DGJ, Labbé N, Sykes RW, Gracom K, Kline L, Swamidoss IM, et al. Rapid assessment of lignin content and structure in switchgrass (Panicum virgatum L.) grown under different environmental conditions. BioEnergy Research. 2009;2(4):246-56.
Hinchee M, Rottmann W, Mullinax L, Zhang C, Chang S, Cunningham M, et al. Short-rotation woody crops for bioenergy and biofuels applications. In vitro cellular & developmental biology Plant : journal of the Tissue Culture Association. 2009;45(6):619-29.
Dinus RJ, Payne P, Sewell MM, Chiang VL, Tuskan GA. Genetic Modification of Short Rotation Popular Wood: Properties for Ethanol Fuel and Fiber Productions. Critical Reviews in Plant Sciences. 2001;20(1):51-69.
Grattapaglia D, Resende MDV. Genomic selection in forest tree breeding. Tree Genetics & Genomes. 2011;7(2):241-55.
Neale DB, Kremer A. Forest tree genomics: growing resources and applications. Nature Reviews Genetics. 2011;12(2):111-22.
Frankham R, Falconer DS, Mackay Longman TFC. Introduction to quantitative genetics. Trends in Genetics. 1996;12(7):280.
Falconer DS. The Problem of Environment and Selection. The American Naturalist. 1952;86(830):293-8.
Visscher PM, Hill WG, Wray NR. Heritability in the genomics era — concepts and misconceptions. Nature Reviews Genetics. 2008;9(4):255-66.
Decker SR, Harman-Ware AE, Happs RM, Wolfrum EJ, Tuskan GA, Kainer D, et al. High Throughput Screening Technologies in Biomass Characterization. Frontiers in Energy Research. 2018;6.
Payne CE, Wolfrum EJ. Rapid analysis of composition and reactivity in cellulosic biomass feedstocks with near-infrared spectroscopy. Biotechnology for Biofuels. 2015;8(1):43.
Harman-Ware A, Crocker M, Pace R, Placido A, Morton S, III, DeBolt S. Characterization of Endocarp Biomass and Extracted Lignin Using Pyrolysis and Spectroscopic Methods. BioEnergy Research. 2014:1-19.
Harman-Ware AE, Davis MF, Peter GF, Wang Y, Sykes RW. Estimation of terpene content in loblolly pine biomass using a hybrid fast-GC and pyrolysis-molecular beam mass spectrometry method. Journal of Analytical and Applied Pyrolysis. 2017.
Penning BW, Sykes RW, Babcock NC, Dugard CK, Klimek JF, Gamblin D, et al. Validation of PyMBMS as a High-throughput Screen for Lignin Abundance in Lignocellulosic Biomass of Grasses. BioEnergy Research. 2014;7(3):899-908.
Sykes R, Kodrzycki B, Tuskan G, Foutz K, Davis M. Within tree variability of lignin composition in Populus. Wood science and technology. 2008;42(8):649-61.
Sykes RW, Gjersing EL, Doeppke CL, Davis MF. High-Throughput Method for Determining the Sugar Content in Biomass with Pyrolysis Molecular Beam Mass Spectrometry. BioEnergy Research. 2015;8(3):964-72.
Li M, Williams DL, Heckwolf M, de Leon N, Kaeppler S, Sykes RW, et al. Prediction of Cell Wall Properties and Response to Deconstruction Using Alkaline Pretreatment in Diverse Maize Genotypes Using Py-MBMS and NIR. BioEnergy Research. 2017;10(2):329-43.
Christensen E, Evans RJ, Carpenter D. High-resolution mass spectrometric analysis of biomass pyrolysis vapors. Journal of Analytical and Applied Pyrolysis. 2017;124:327-34.
Evans RJ, Milne TA. Molecular characterization of the pyrolysis of biomass. Energy & Fuels. 1987;1(2):123-37.
Jarvis JM, Page-Dumroese DS, Anderson NM, Corilo Y, Rodgers RP. Characterization of Fast Pyrolysis Products Generated from Several Western USA Woody Species. Energy & Fuels. 2014;28(10):6438-46.
Tessarolo NS, Silva RC, Vanini G, Pinho A, Romão W, de Castro EVR, et al. Assessing the chemical composition of bio-oils using FT-ICR mass spectrometry and comprehensive two-dimensional gas chromatography with time-of-flight mass spectrometry. Microchemical Journal. 2014;117:68-76.
Davis MF, Johnson DK, Deutch SP, Agblevor FA, Fennell J, Ashley P, editors. Variability in the Composition of Short Rotation Woody Feedstocks. Second Biomass Conference of the Americas; 1995; Portland, OR.
Johnson D, Ashley P, Deutch S, Davis M, Fennell J, Wiselogel A, editors. Compositional variability in herbaceous energy crops. Second Biomass Conference of the Americas: Energy, Environment, Agriculture, and Industry Proceedings Golden, CO: National Renewable Energy Laboratory; 1995.
Tuskan G, West D, Bradshaw HD, Neale D, Sewell M, Wheeler N, et al. Two high-throughput techniques for determining wood properties as part of a molecular genetics analysis of hybrid poplar and loblolly pine. Applied Biochemistry and Biotechnology. 1999;77(1):55-65.
Weighill D, Jones P, Shah M, Ranjan P, Muchero W, Schmutz J, et al. Pleiotropic and Epistatic Network-Based Discovery: Integrated Networks for Target Gene Discovery. Frontiers in Energy Research. 2018;6(30).
Yin T, Zhang X, Gunter L, Priya R, Sykes R, Davis M, et al. Differential Detection of Genetic Loci Underlying Stem and Root Lignin Content in Populus. PLOS ONE. 2010;5(11):e14021.
Morse AM, Tschaplinski TJ, Dervinis C, Pijut PM, Schmelz EA, Day W, et al. Salicylate and catechol levels are maintained in nahG transgenic poplar. Phytochemistry. 2007;68(15):2043-52.
Tschaplinski TJ, Abraham PE, Jawdy SS, Gunter LE, Martin MZ, Engle NL, et al. The nature of the progression of drought stress drives differential metabolomic responses in Populus deltoides. Annals of Botany. 2019.
Abreu IN, Ahnlund M, Moritz T, Albrectsen BR. UHPLC-ESI/TOFMS determination of salicylate-like phenolic gycosides in Populus tremula leaves. Journal of chemical ecology. 2011;37(8):857-70.
Tsai CJ, Guo W, Babst B, Nyamdari B, Yuan Y, Payyavula R, et al. Salicylate metabolism in Populus. BMC Proceedings. 2011;5(7):I9.
Zhao N, Guan J, Forouhar F, Tschaplinski TJ, Cheng Z-M, Tong L, et al. Two poplar methyl salicylate esterases display comparable biochemical properties but divergent expression patterns. Phytochemistry. 2009;70(1):32-9.
Harman-Ware AE, Crocker M, Kaur AP, Meier MS, Kato D, Lynn B. Pyrolysis–GC/MS of sinapyl and coniferyl alcohol. Journal of analytical and applied pyrolysis. 2013;99:161-9.
Ware AE. Application of pyrolysis-GC/MS to the study of biomass and biomass constituents. 2013.
Lapierre C, Pilate G, Pollet B, Mila I, Leplé J-C, Jouanin L, et al. Signatures of cinnamyl alcohol dehydrogenase deficiency in poplar lignins. Phytochemistry. 2004;65(3):313-21.
Kawamoto H, Horigoshi S, Saka S. Pyrolysis reactions of various lignin model dimers. Journal of Wood Science. 2007;53(2):168-74.
Kim H, Ralph J, Lu F, Ralph SA, Boudet A-M, MacKay JJ, et al. NMR analysis of lignins in CAD-deficient plants. Part 1. Incorporation of hydroxycinnamaldehydes and hydroxybenzaldehydes into lignins. Organic & biomolecular chemistry. 2003;1(2):268-81.
Carlos del Río J, Rencoret J, Gutiérrez A, Kim H, Ralph J. Hydroxystilbenes Are Monomers in Palm Fruit Endocarp Lignins. Plant Physiology. 2017;174(4):2072.
Autran D, Huanca-Mamani W, Vielle-Calzada J-P. Genomic imprinting in plants: the epigenetic version of an Oedipus complex. Current Opinion in Plant Biology. 2005;8(1):19-25.
Jia S, Rob D, Michael JH. Diallel Crossing Among Doubled Haploids of Cucumber Reveals Significant Reciprocal-cross Differences. Journal of the American Society for Horticultural Science J Amer Soc Hort Sci. 2015;140(2):178-82.
Chung S-M, Staub J, Fazio G. Inheritance of Chilling Injury: A Maternally Inherited Trait in Cucumber. J Am Soc Hortic Sci. 2003;128.
Patten MM, Ross L, Curley JP, Queller DC, Bonduriansky R, Wolf JB. The evolution of genomic imprinting: theories, predictions and empirical tests. Heredity. 2014;113(2):119-28.
Galloway LF. Maternal effects provide phenotypic adaptation to local environmental conditions. New Phytologist. 2005;166(1):93-100.
Lafon-Placette C, Hatorangan MR, Steige KA, Cornille A, Lascoux M, Slotte T, et al. Paternally expressed imprinted genes associate with hybridization barriers in Capsella. Nature Plants. 2018;4(6):352-7.
Garner AG, Kenney AM, Fishman L, Sweigart AL. Genetic loci with parent-of-origin effects cause hybrid seed lethality in crosses between Mimulus species. New Phytologist. 2016;211(1):319-31.
Hatorangan MR, Laenen B, Steige KA, Slotte T, Köhler C. Rapid Evolution of Genomic Imprinting in Two Species of the Brassicaceae. The Plant Cell. 2016;28(8):1815.
Chen C, Li T, Zhu S, Liu Z, Shi Z, Zheng X, et al. Characterization of Imprinted Genes in Rice Reveals Conservation of Regulation and Imprinting with Other Plant Species. Plant Physiology. 2018;177(4):1754.
Jumbo MB, Carena MJ. Combining ability, maternal, and reciprocal effects of elite early-maturing maize population hybrids. Euphytica. 2008;162(3):325-33.
Jiang H, Köhler C. Evolution, function, and regulation of genomic imprinting in plant seed development. Journal of Experimental Botany. 2012;63(13):4713-22.
Schmid B, Dolt C. Effects of Maternal and Paaternal Environment and Genotype on Offspring Phenotype in Solidago Altissma L. . Evolution. 1994;48(5):1525-49.
Zhang J, Novaes E, Kirst M, Peter GF. Comparison of Pyrolysis Mass Spectrometry and Near Infrared Spectroscopy for Genetic Analysis of Lignocellulose Chemical Composition in Populus. Forests. 2014;5:466-81.
Guerra FP, Richards JH, Fiehn O, Famula R, Stanton BJ, Shuren R, et al. Analysis of the genetic variation in growth, ecophysiology, and chemical and metabolomic composition of wood of Populus trichocarpa provenances. Tree Genetics & Genomes. 2016;12(1):6.
Guerra FP, Wegrzyn JL, Sykes R, Davis MF, Stanton BJ, Neale DB. Association genetics of chemical wood properties in black poplar (Populus nigra). New Phytologist. 2013;197(1):162-76.
Schrag TA, Westhues M, Schipprack W, Seifert F, Thiemann A, Scholten S, et al. Beyond Genomic Prediction: Combining Different Types of <em>omics</em> Data Can Improve Prediction of Hybrid Performance in Maize. Genetics. 2018;208(4):1373.
Davey CL, Robson P, Hawkins S, Farrar K, Clifton-Brown JC, Donnison IS, et al. Genetic relationships between spring emergence, canopy phenology, and biomass yield increase the accuracy of genomic prediction in Miscanthus. Journal of Experimental Botany. 2017;68(18):5093-102.
Hu Z, Sykes R, Davis MF, Charles Brummer E, Ragauskas AJ. Chemical profiles of switchgrass. Bioresource Technology. 2010;101(9):3253-7.
Decker SR, Harman-Ware AE, Happs RM, Wolfrum EJ, Tuskan GA, Kainer D, et al. High Throughput Screening Technologies in Biomass Characterization. Frontiers in Energy Research. 2018;6.
Furrer R, Nychka D, Sain S. Fields: Tools for spatial data. 2017.
Team RC. R: A language and environment for statistical computing. 2013.

7x7MBMSheritSupplementaryMaterialfinal.pdf
Supplementary Material
SuppTable3MZannotations.xlsx
Supplementary Table 3
SuppTable4final.xlsx
Supplementary Table 4
SuppTable5final.xlsx
Supplementary Table 5

Download PDF

Journal Publication

published 06 Mar, 2021

Read the published version in Biotechnology for Biofuels and Bioproducts →

Editorial decision: Accept
15 Feb, 2021
Review #1 received at journal
10 Jan, 2021
Reviewer #1 agreed at journal
03 Jan, 2021
Reviewers invited by journal
27 Dec, 2020
Editor assigned by journal
26 Dec, 2020
Submission checks completed at journal
26 Dec, 2020
Editor invited by journal
26 Dec, 2020

You are reading this latest preprint version

Accurate Determination of Genotypic Variance of Cell Wall Characteristics of a Populus trichocarpa Pedigree Using High-Throughput Pyrolysis-Molecular Beam Mass Spectrometry

Status:

Journal Publication

Version 2

Abstract

Figures

Background

Results

Discussion

Conclusions

Methods

List Of Abbreviations

Tables

Declarations

References

Supplementary Files

Status:

Journal Publication

Version 2