Learning from methylomes: epigenomic correlates of Populus balsamifera traits based on deep learning models of natural DNA methylation.

Epigenomes have remarkable potential for the estimation of plant traits. This study tested the hypothesis that natural variation in DNA methylation can be used to estimate industrially important traits in a genetically diverse population of Populus balsamifera L. (balsam poplar) trees grown at two common garden sites. Statistical learning experiments enabled by deep learning models revealed that plant traits in novel genotypes can be modelled transparently using small numbers of methylated DNA predictors. Using this approach, tissue type, a nonheritable attribute, from which DNA methylomes were derived was assigned, and provenance, a purely heritable trait and an element of population structure, was determined. Significant proportions of phenotypic variance in quantitative wood traits, including total biomass (57.5%), wood density (40.9%), soluble lignin (25.3%) and cell wall carbohydrate (mannose: 44.8%) contents, were also explained from natural variation in DNA methylation. Modelling plant traits using DNA methylation can capture tissue-specific epigenetic mechanisms underlying plant phenotypes in natural environments. DNA methylation-based models offer new insight into natural epigenetic influence on plants and can be used as a strategy to validate the identity, provenance or quality of agroforestry products.

Author Correction: Sexual homomorphism in dioecious trees: extensive tests fail to detect sexual dimorphism in Populus.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has been fixed in the paper.

Sexual homomorphism in dioecious trees: extensive tests fail to detect sexual dimorphism in Populus †.

The evolution of sexual dimorphism and expansion of sex chromosomes are both driven through sexual conflict, arising from differing fitness optima between males and females. Here, we pair work in poplar (Populus) describing one of the smallest sex-determining regions known thus far in complex eukaryotes (~100 kbp) with comprehensive tests for sexual dimorphism using >1300 individuals from two Populus species and assessing 96 non-reproductive functional traits. Against expectation, we found sexual homomorphism (no non-reproductive trait differences between the sexes), suggesting that gender is functionally neutral with respect to non-reproductive features that affect plant survival and fitness. Combined with a small sex-determining region, we infer that sexual conflict may be effectively stymied or non-existent within these taxa. Both sexual homomorphism and the small sex-determining region occur against a background of strong environmental selection and local adaptation in Populus. This presents a powerful hypothesis for the evolution of dioecious species. Here, we suggest that environmental selection may be sufficient to suppress and stymy sexual conflict if it acts orthogonal to sexual selection, thereby placing limitations on the evolution of sexual dimorphism and genomic expansion of sex chromosomes.

Gene Expression Patterns of Wood Decay Fungi Postia placenta and Phanerochaete chrysosporium Are Influenced by Wood Substrate Composition during Degradation.

UNLABELLED: Identification of the specific genes and enzymes involved in the fungal degradation of lignocellulosic biomass derived from feedstocks with various compositions is essential to the development of improved bioenergy processes. In order to elucidate the effect of substrate composition on gene expression in wood-rotting fungi, we employed microarrays based on the annotated genomes of the brown- and white-rot fungi, Rhodonia placenta (formerly Postia placenta) and Phanerochaete chrysosporium, respectively. We monitored the expression of genes involved in the enzymatic deconstruction of the cell walls of three 4-year-old Populus trichocarpa (poplar) trees of genotypes with distinct cell wall chemistries, selected from a population of several hundred trees grown in a common garden. The woody substrates were incubated with wood decay fungi for 10, 20, and 30 days. An analysis of transcript abundance in all pairwise comparisons highlighted 64 and 84 differentially expressed genes (>2-fold, P < 0.05) in P. chrysosporium and P. placenta, respectively. Cross-fungal comparisons also revealed an array of highly differentially expressed genes (>4-fold, P < 0.01) across different substrates and time points. These results clearly demonstrate that gene expression profiles of P. chrysosporium and P. placenta are influenced by wood substrate composition and the duration of incubation. Many of the significantly expressed genes encode "proteins of unknown function," and determining their role in lignocellulose degradation presents opportunities and challenges for future research. IMPORTANCE: This study describes the variation in expression patterns of two wood-degrading fungi (brown- and white-rot fungi) during colonization and incubation on three different naturally occurring poplar substrates of differing chemical compositions, over time. The results clearly show that the two fungi respond differentially to their substrates and that several known and, more interestingly, currently unknown genes are highly misregulated in response to various substrate compositions. These findings highlight the need to characterize several unknown proteins for catalytic function but also as potential candidate proteins to improve the efficiency of enzymatic cocktails to degrade lignocellulosic substrates in industrial applications, such as in a biochemically based bioenergy platform.

Histology and cell wall biochemistry of stone cells in the physical defence of conifers against insects.

Conifers possess an array of physical and chemical defences against stem-boring insects. Stone cells provide a physical defence associated with resistance against bark beetles and weevils. In Sitka spruce (Picea sitchensis), abundance of stone cells in the cortex of apical shoots is positively correlated with resistance to white pine weevil (Pissodes strobi). We identified histological, biochemical and molecular differences in the stone cell phenotype of weevil resistant (R) or susceptible (S) Sitka spruce genotypes. R trees displayed significantly higher quantities of cortical stone cells near the apical shoot node, the primary site for weevil feeding. Lignin, cellulose, xylan and mannan were the most abundant components of stone cell secondary walls, respectively. Lignin composition of stone cells isolated from R trees contained a higher percentage of G-lignin compared with S trees. Transcript profiling revealed higher transcript abundance in the R genotype of coumarate 3-hydroxylase, a key monolignol biosynthetic gene. Developing stone cells in current year apical shoots incorporated fluorescent-tagged monolignol into the secondary cell wall, while mature stone cells of previous year apical shoots did not. Stone cell development is an ephemeral process, and fortification of shoot tips in R trees is an effective strategy against insect feeding.

High-resolution genetic mapping of allelic variants associated with cell wall chemistry in Populus.

BACKGROUND: QTL cloning for the discovery of genes underlying polygenic traits has historically been cumbersome in long-lived perennial plants like Populus. Linkage disequilibrium-based association mapping has been proposed as a cloning tool, and recent advances in high-throughput genotyping and whole-genome resequencing enable marker saturation to levels sufficient for association mapping with no a priori candidate gene selection. Here, multiyear and multienvironment evaluation of cell wall phenotypes was conducted in an interspecific P. trichocarpa x P. deltoides pseudo-backcross mapping pedigree and two partially overlapping populations of unrelated P. trichocarpa genotypes using pyrolysis molecular beam mass spectrometry, saccharification, and/ or traditional wet chemistry. QTL mapping was conducted using a high-density genetic map with 3,568 SNP markers. As a fine-mapping approach, chromosome-wide association mapping targeting a QTL hot-spot on linkage group XIV was performed in the two P. trichocarpa populations. Both populations were genotyped using the 34 K Populus Infinium SNP array and whole-genome resequencing of one of the populations facilitated marker-saturation of candidate intervals for gene identification. RESULTS: Five QTLs ranging in size from 0.6 to 1.8 Mb were mapped on linkage group XIV for lignin content, syringyl to guaiacyl (S/G) ratio, 5- and 6-carbon sugars using the mapping pedigree. Six candidate loci exhibiting significant associations with phenotypes were identified within QTL intervals. These associations were reproducible across multiple environments, two independent genotyping platforms, and different plant growth stages. cDNA sequencing for allelic variants of three of the six loci identified polymorphisms leading to variable length poly glutamine (PolyQ) stretch in a transcription factor annotated as an ANGUSTIFOLIA C-terminus Binding Protein (CtBP) and premature stop codons in a KANADI transcription factor as well as a protein kinase. Results from protoplast transient expression assays suggested that each of the polymorphisms conferred allelic differences in the activation of cellulose, hemicelluloses, and lignin pathway marker genes. CONCLUSION: This study illustrates the utility of complementary QTL and association mapping as tools for gene discovery with no a priori candidate gene selection. This proof of concept in a perennial organism opens up opportunities for discovery of novel genetic determinants of economically important but complex traits in plants.

Influence of Populus genotype on gene expression by the wood decay fungus Phanerochaete chrysosporium.

We examined gene expression patterns in the lignin-degrading fungus Phanerochaete chrysosporium when it colonizes hybrid poplar (Populus alba × tremula) and syringyl (S)-rich transgenic derivatives. A combination of microarrays and liquid chromatography-tandem mass spectrometry (LC-MS/MS) allowed detection of a total of 9,959 transcripts and 793 proteins. Comparisons of P. chrysosporium transcript abundance in medium containing poplar or glucose as a sole carbon source showed 113 regulated genes, 11 of which were significantly higher (>2-fold, P < 0.05) in transgenic line 64 relative to the parental line. Possibly related to the very large amounts of syringyl (S) units in this transgenic tree (94 mol% S), several oxidoreductases were among the upregulated genes. Peptides corresponding to a total of 18 oxidoreductases were identified in medium consisting of biomass from line 64 or 82 (85 mol% S) but not in the parental clone (65 mol% S). These results demonstrate that P. chrysosporium gene expression patterns are substantially influenced by lignin composition.

Extensive functional pleiotropy of REVOLUTA substantiated through forward genetics.

In plants, genes may sustain extensive pleiotropic functional properties by individually affecting multiple, distinct traits. We discuss results from three genome-wide association studies of approximately 400 natural poplar (Populus trichocarpa) accessions phenotyped for 60 ecological/biomass, wood quality, and rust fungus resistance traits. Single-nucleotide polymorphisms (SNPs) in the poplar ortholog of the class III homeodomain-leucine zipper transcription factor gene REVOLUTA (PtREV) were significantly associated with three specific traits. Based on SNP associations with fungal resistance, leaf drop, and cellulose content, the PtREV gene contains three potential regulatory sites within noncoding regions at the gene's 3' end, where alternative splicing and messenger RNA processing actively occur. The polymorphisms in this region associated with leaf abscission and cellulose content are suggested to represent more recent variants, whereas the SNP associated with leaf rust resistance may be more ancient, consistent with REV's primary role in auxin signaling and its functional evolution in supporting fundamental processes of vascular plant development.

Network analysis reveals the relationship among wood properties, gene expression levels and genotypes of natural Populus trichocarpa accessions.

High-throughput approaches have been widely applied to elucidate the genetic underpinnings of industrially important wood properties. Wood traits are polygenic in nature, but gene hierarchies can be assessed to identify the most important gene variants controlling specific traits within complex networks defining the overall wood phenotype. We tested a large set of genetic, genomic, and phenotypic information in an integrative approach to predict wood properties in Populus trichocarpa. Nine-yr-old natural P. trichocarpa trees including accessions with high contrasts in six traits related to wood chemistry and ultrastructure were profiled for gene expression on 49k Nimblegen (Roche NimbleGen Inc., Madison, WI, USA) array elements and for 28,831 polymorphic single nucleotide polymorphisms (SNPs). Pre-selected transcripts and SNPs with high statistical dependence on phenotypic traits were used in Bayesian network learning procedures with a stepwise K2 algorithm to infer phenotype-centric networks. Transcripts were pre-selected at a much lower logarithm of Bayes factor (logBF) threshold than SNPs and were not accommodated in the networks. Using persistent variables, we constructed cross-validated networks for variability in wood attributes, which contained four to six variables with 94-100% predictive accuracy. Accommodated gene variants revealed the hierarchy in the genetic architecture that underpins substantial phenotypic variability, and represent new tools to support the maximization of response to selection.

Genome-wide association mapping for wood characteristics in Populus identifies an array of candidate single nucleotide polymorphisms.

Establishing links between phenotypes and molecular variants is of central importance to accelerate genetic improvement of economically important plant species. Our work represents the first genome-wide association study to the inherently complex and currently poorly understood genetic architecture of industrially relevant wood traits. Here, we employed an Illumina Infinium 34K single nucleotide polymorphism (SNP) genotyping array that generated 29,233 high-quality SNPs in c. 3500 broad-based candidate genes within a population of 334 unrelated Populus trichocarpa individuals to establish genome-wide associations. The analysis revealed 141 significant SNPs (α ≤ 0.05) associated with 16 wood chemistry/ultrastructure traits, individually explaining 3-7% of the phenotypic variance. A large set of associations (41% of all hits) occurred in candidate genes preselected for their suggested a priori involvement with secondary growth. For example, an allelic variant in the FRA8 ortholog explained 21% of the total genetic variance in fiber length, when the trait's heritability estimate was considered. The remaining associations identified SNPs in genes not previously implicated in wood or secondary wall formation. Our findings provide unique insights into wood trait architecture and support efforts for population improvement based on desirable allelic variants.

Syringyl-rich lignin renders poplars more resistant to degradation by wood decay fungi.

In order to elucidate the effects of lignin composition on the resistance of wood to degradation by decay fungi, wood specimens from two transgenic poplar lines expressing an Arabidopsis gene encoding ferulate 5-hydroxylase (F5H) driven by the cinnimate-4-hydroxylase promoter (C4H::F5H) that increased syringyl/guaiacyl (S/G) monolignol ratios relative to those in the untransformed control wood were incubated with six different wood decay fungi. Alterations in wood weight and chemical composition were monitored over the incubation period. The results showed that transgenic poplar lines extremely rich in syringyl lignin exhibited a drastically improved resistance to degradation by all decay fungi evaluated. Lignin monomer composition and its distribution among cell types and within different cell layers were the sole wood chemistry parameters determining wood durability. Since transgenic poplars with exceedingly high syringyl contents were recalcitrant to degradation, where wood durability is a critical factor, these genotypes may offer improved performance.

Populus trichocarpa cell wall chemistry and ultrastructure trait variation, genetic control and genetic correlations.

The increasing ecological and economical importance of Populus species and hybrids has stimulated research into the investigation of the natural variation of the species and the estimation of the extent of genetic control over its wood quality traits for traditional forestry activities as well as the emerging bioenergy sector. A realized kinship matrix based on informative, high-density, biallelic single nucleotide polymorphism (SNP) genetic markers was constructed to estimate trait variance components, heritabilities, and genetic and phenotypic correlations. Seventeen traits related to wood chemistry and ultrastructure were examined in 334 9-yr-old Populus trichocarpa grown in a common-garden plot representing populations spanning the latitudinal range 44° to 58.6°. In these individuals, 9342 SNPs that conformed to Hardy-Weinberg expectations were employed to assess the genomic pair-wise kinship to estimate narrow-sense heritabilities and genetic correlations among traits. The range-wide phenotypic variation in all traits was substantial and several trait heritabilities were > 0.6. In total, 61 significant genetic and phenotypic correlations and a network of highly interrelated traits were identified. The high trait variation, the evidence for moderate to high heritabilities and the identification of advantageous trait combinations of industrially important characteristics should aid in providing the foundation for the enhancement of poplar tree breeding strategies for modern industrial use.

Significant alteration of gene expression in wood decay fungi Postia placenta and Phanerochaete chrysosporium by plant species.

Identification of specific genes and enzymes involved in conversion of lignocellulosics from an expanding number of potential feedstocks is of growing interest to bioenergy process development. The basidiomycetous wood decay fungi Phanerochaete chrysosporium and Postia placenta are promising in this regard because they are able to utilize a wide range of simple and complex carbon compounds. However, systematic comparative studies with different woody substrates have not been reported. To address this issue, we examined gene expression of these fungi colonizing aspen (Populus grandidentata) and pine (Pinus strobus). Transcript levels of genes encoding extracellular glycoside hydrolases, thought to be important for hydrolytic cleavage of hemicelluloses and cellulose, showed little difference for P. placenta colonizing pine versus aspen as the sole carbon source. However, 164 genes exhibited significant differences in transcript accumulation for these substrates. Among these, 15 cytochrome P450s were upregulated in pine relative to aspen. Of 72 P. placenta extracellular proteins identified unambiguously by mass spectrometry, 52 were detected while colonizing both substrates and 10 were identified in pine but not aspen cultures. Most of the 178 P. chrysosporium glycoside hydrolase genes showed similar transcript levels on both substrates, but 13 accumulated >2-fold higher levels on aspen than on pine. Of 118 confidently identified proteins, 31 were identified in both substrates and 57 were identified in pine but not aspen cultures. Thus, P. placenta and P. chrysosporium gene expression patterns are influenced substantially by wood species. Such adaptations to the carbon source may also reflect fundamental differences in the mechanisms by which these fungi attack plant cell walls.

Comparative transcriptome and secretome analysis of wood decay fungi Postia placenta and Phanerochaete chrysosporium.

Cellulose degradation by brown rot fungi, such as Postia placenta, is poorly understood relative to the phylogenetically related white rot basidiomycete, Phanerochaete chrysosporium. To elucidate the number, structure, and regulation of genes involved in lignocellulosic cell wall attack, secretome and transcriptome analyses were performed on both wood decay fungi cultured for 5 days in media containing ball-milled aspen or glucose as the sole carbon source. Using liquid chromatography-tandem mass spectrometry (LC-MS/MS), a total of 67 and 79 proteins were identified in the extracellular fluids of P. placenta and P. chrysosporium cultures, respectively. Viewed together with transcript profiles, P. chrysosporium employs an array of extracellular glycosyl hydrolases to simultaneously attack cellulose and hemicelluloses. In contrast, under these same conditions, P. placenta secretes an array of hemicellulases but few potential cellulases. The two species display distinct expression patterns for oxidoreductase-encoding genes. In P. placenta, these patterns are consistent with an extracellular Fenton system and include the upregulation of genes involved in iron acquisition, in the synthesis of low-molecular-weight quinones, and possibly in redox cycling reactions.