Overexpression of REDUCED WALL ACETYLATION C increases xylan acetylation and biomass recalcitrance in Populus

Plant lignocellulosic biomass, i.e. secondary cell walls of plants, is a vital alternative source for bioenergy. However, the acetylation of xylan in secondary cell walls impedes the conversion of biomass to biofuels. Previous studies have shown that REDUCED WALL ACETYLATION (RWA) proteins are directly involved in the acetylation of xylan but the regulatory mechanism of RWAs is not fully understood. In this study, we demonstrate that overexpression of a Populus trichocarpa PtRWA-C gene increases the level of xylan acetylation and increases the lignin content and S/G ratio, ultimately yielding poplar woody biomass with reduced saccharification efficiency. Furthermore, through gene coexpression network and expression quantitative trait loci (eQTL) analysis, we found that PtRWA-C was regulated not only by the secondary cell wall hierarchical regulatory network but also by an AP2 family transcription factor HARDY (HRD). Specifically, HRD activates PtRWA-C expression by directly binding to the PtRWA-C promoter, which is also the cis-eQTL for PtRWA-C. Taken together, our findings provide insights into the functional roles of PtRWA-C in xylan acetylation and consequently saccharification and shed light on synthetic biology approaches to manipulate this gene and alter cell wall properties. These findings have substantial implications for genetic engineering of woody species, which could be used as a sustainable source of biofuels, valuable biochemicals, and biomaterials.

Overexpression of a Prefoldin β subunit gene reduces biomass recalcitrance in the bioenergy crop Populus

Prefoldin (PFD) is a group II chaperonin that is ubiquitously present in the eukaryotic kingdom. Six subunits (PFD1-6) form a jellyfish-like heterohexameric PFD complex and function in protein folding and cytoskeleton organization. However, little is known about its function in plant cell wall-related processes. Here, we report the functional characterization of a PFD gene from Populus deltoides, designated as PdPFD2.2. There are two copies of PFD2 in Populus, and PdPFD2.2 was ubiquitously expressed with high transcript abundance in the cambial region. PdPFD2.2 can physically interact with DELLA protein RGA1_8g, and its subcellular localization is affected by the interaction. In P. deltoides transgenic plants overexpressing PdPFD2.2, the lignin syringyl/guaiacyl ratio was increased, but cellulose content and crystallinity index were unchanged. In addition, the total released sugar (glucose and xylose) amounts were increased by 7.6% and 6.1%, respectively, in two transgenic lines. Transcriptomic and metabolomic analyses revealed that secondary metabolic pathways, including lignin and flavonoid biosynthesis, were affected by overexpressing PdPFD2.2. A total of eight hub transcription factors (TFs) were identified based on TF binding sites of differentially expressed genes in Populus transgenic plants overexpressing PdPFD2.2. In addition, several known cell wall-related TFs, such as MYB3, MYB4, MYB7, TT8 and XND1, were affected by overexpression of PdPFD2.2. These results suggest that overexpression of PdPFD2.2 can reduce biomass recalcitrance and PdPFD2.2 is a promising target for genetic engineering to improve feedstock characteristics to enhance biofuel conversion and reduce the cost of lignocellulosic biofuel production.

PdWND3A, a wood-associated NAC domain-containing protein, affects lignin biosynthesis and composition in Populus

BACKGROUND

Plant secondary cell wall is a renewable feedstock for biofuels and biomaterials production. Arabidopsis VASCULAR-RELATED NAC DOMAIN (VND) has been demonstrated to be a key transcription factor regulating secondary cell wall biosynthesis. However, less is known about its role in the woody species.

RESULTS

Here we report the functional characterization of Populus deltoides WOOD-ASSOCIATED NAC DOMAIN protein 3 (PdWND3A), a sequence homolog of Arabidopsis VND4 and VND5 that are members of transcription factor networks regulating secondary cell wall biosynthesis. PdWND3A was expressed at higher level in the xylem than in other tissues. The stem tissues of transgenic P. deltoides overexpressing PdWND3A (OXPdWND3A) contained more vessel cells than that of wild-type plants. Furthermore, lignin content and lignin monomer syringyl and guaiacyl (S/G) ratio were higher in OXPdWND3A transgenic plants than in wild-type plants. Consistent with these observations, the expression of FERULATE 5-HYDROXYLASE1 (F5H1), encoding an enzyme involved in the biosynthesis of sinapyl alcohol (S unit monolignol), was elevated in OXPdWND3A transgenic plants. Saccharification analysis indicated that the rate of sugar release was reduced in the transgenic plants. In addition, OXPdWND3A transgenic plants produced lower amounts of biomass than wild-type plants.

CONCLUSIONS

PdWND3A affects lignin biosynthesis and composition and negatively impacts sugar release and biomass production.

The nature of the progression of drought stress drives differential metabolomic responses in Populus deltoides

BACKGROUND AND AIMS

The use of woody crops for Quad-level (approx. 1 × 1018 J) energy production will require marginal agricultural lands that experience recurrent periods of water stress. Populus species have the capacity to increase dehydration tolerance by lowering osmotic potential via osmotic adjustment. The aim of this study was to investigate how the inherent genetic potential of a Populus clone to respond to drought interacts with the nature of the drought to determine the degree of biochemical response.

METHODS

A greenhouse drought stress study was conducted on Populus deltoides 'WV94' and the resulting metabolite profiles of leaves were determined by gas chromatography-mass spectrometry following trimethylsilylation for plants subjected to cyclic mild (-0.5 MPa pre-dawn leaf water potential) drought vs. cyclic severe (-1.26 MPa) drought in contrast to well-watered controls (-0.1 MPa) after two or four drought cycles, and in contrast to plants subjected to acute drought, where plants were desiccated for up to 8 d.

KEY RESULTS

The nature of drought (cyclic vs. acute), frequency of drought (number of cycles) and the severity of drought (mild vs. severe) all dictated the degree of osmotic adjustment and the nature of the organic solutes that accumulated. Whereas cyclic drought induced the largest responses in primary metabolism (soluble sugars, organic acids and amino acids), acute onset of prolonged drought induced the greatest osmotic adjustment and largest responses in secondary metabolism, especially populosides (hydroxycinnamic acid conjugates of salicin).

CONCLUSIONS

The differential adaptive metabolite responses in cyclic vs. acute drought suggest that stress acclimation occurs via primary metabolism in response to cyclic drought, whereas expanded metabolic plasticity occurs via secondary metabolism following severe, acute drought. The shift in carbon partitioning to aromatic metabolism with the production of a diverse suite of higher order salicylates lowers osmotic potential and increases the probability of post-stress recovery.

Mediation of plant-mycorrhizal interaction by a lectin receptor-like kinase

The molecular mechanisms underlying mycorrhizal symbioses, the most ubiquitous and impactful mutualistic plant-microbial interaction in nature, are largely unknown. Through genetic mapping, resequencing and molecular validation, we demonstrate that a G-type lectin receptor-like kinase (lecRLK) mediates the symbiotic interaction between Populus and the ectomycorrhizal fungus Laccaria bicolor. This finding uncovers an important molecular step in the establishment of symbiotic plant-fungal associations and provides a molecular target for engineering beneficial mycorrhizal relationships.

A New Calmodulin-Binding Protein Expresses in the Context of Secondary Cell Wall Biosynthesis and Impacts Biomass Properties in Populus

A greater understanding of biosynthesis, signaling and regulatory pathways involved in determining stem growth and secondary cell wall chemistry is important for enabling pathway engineering and genetic optimization of biomass properties. The present study describes a new functional role of PdIQD10, a Populus gene belonging to the IQ67-Domain1 family of IQD genes, in impacting biomass formation and chemistry. Expression studies showed that PdIQD10 has enhanced expression in developing xylem and tension-stressed tissues in Populus deltoides. Molecular dynamics simulation and yeast two-hybrid interaction experiments suggest interactions with two calmodulin proteins, CaM247 and CaM014, supporting the sequence-predicted functional role of the PdIQD10 as a calmodulin-binding protein. PdIQD10 was found to interact with specific Populus isoforms of the Kinesin Light Chain protein family, shown previously to function as microtubule-guided, cargo binding and delivery proteins in Arabidopsis. Subcellular localization studies showed that PdIQD10 localizes in the nucleus and plasma membrane regions. Promoter-binding assays suggest that a known master transcriptional regulator of secondary cell wall biosynthesis (PdWND1B) may be upstream of an HD-ZIP III gene that is in turn upstream of PdIQD10 gene in the transcriptional network. RNAi-mediated downregulation of PdIQD10 expression resulted in plants with altered biomass properties including higher cellulose, wall glucose content and greater biomass quantity. These results present evidence in support of a new functional role for an IQD gene family member, PdIQD10, in secondary cell wall biosynthesis and biomass formation in Populus.

Genome-wide association studies and expression-based quantitative trait loci analyses reveal roles of HCT2 in caffeoylquinic acid biosynthesis and its regulation by defense-responsive transcription factors in Populus

3-O-caffeoylquinic acid, also known as chlorogenic acid (CGA), functions as an intermediate in lignin biosynthesis in the phenylpropanoid pathway. It is widely distributed among numerous plant species and acts as an antioxidant in both plants and animals. Using GC-MS, we discovered consistent and extreme variation in CGA content across a population of 739 4-yr-old Populus trichocarpa accessions. We performed genome-wide association studies (GWAS) from 917 P. trichocarpa accessions and expression-based quantitative trait loci (eQTL) analyses to identify key regulators. The GWAS and eQTL analyses resolved an overlapped interval encompassing a hydroxycinnamoyl-CoA:shikimate hydroxycinnamoyl transferase 2 (PtHCT2) that was significantly associated with CGA and partially characterized metabolite abundances. PtHCT2 leaf expression was significantly correlated with CGA abundance and it was regulated by cis-eQTLs containing W-box for WRKY binding. Among all nine PtHCT homologs, PtHCT2 is the only one that responds to infection by the fungal pathogen Sphaerulina musiva (a Populus pathogen). Validation using protoplast-based transient expression system suggests that PtHCT2 is regulated by the defense-responsive WRKY. These results are consistent with reports of CGA functioning as an antioxidant in response to biotic stress. This study provides insights into data-driven and omics-based inference of gene function in woody species.

Defining the genetic components of callus formation: A GWAS approach

A characteristic feature of plant cells is the ability to form callus from parenchyma cells in response to biotic and abiotic stimuli. Tissue culture propagation of recalcitrant plant species and genetic engineering for desired phenotypes typically depends on efficient in vitro callus generation. Callus formation is under genetic regulation, and consequently, a molecular understanding of this process underlies successful generation for propagation materials and/or introduction of genetic elements in experimental or industrial applications. Herein, we identified 11 genetic loci significantly associated with callus formation in Populus trichocarpa using a genome-wide association study (GWAS) approach. Eight of the 11 significant gene associations were consistent across biological replications, exceeding a chromosome-wide-log10 (p) = 4.46 [p = 3.47E-05] Bonferroni-adjusted significance threshold. These eight genes were used as hub genes in a high-resolution co-expression network analysis to gain insight into the genome-wide basis of callus formation. A network of positively and negatively co-expressed genes, including several transcription factors, was identified. As proof-of-principle, a transient protoplast assay confirmed the negative regulation of a Chloroplast Nucleoid DNA-binding-related gene (Potri.018G014800) by the LEC2 transcription factor. Many of the candidate genes and co-expressed genes were 1) linked to cell division and cell cycling in plants and 2) showed homology to tumor and cancer-related genes in humans. The GWAS approach based on a high-resolution marker set, and the ability to manipulate targets genes in vitro, provided a catalog of high-confidence genes linked to callus formation that can serve as an important resource for successful manipulation of model and non-model plant species, and likewise, suggests a robust method of discovering common homologous functions across organisms.

A 5-Enolpyruvylshikimate 3-Phosphate Synthase Functions as a Transcriptional Repressor in Populus

Long-lived perennial plants, with distinctive habits of inter-annual growth, defense, and physiology, are of great economic and ecological importance. However, some biological mechanisms resulting from genome duplication and functional divergence of genes in these systems remain poorly studied. Here, we discovered an association between a poplar (Populus trichocarpa) 5-enolpyruvylshikimate 3-phosphate synthase gene (PtrEPSP) and lignin biosynthesis. Functional characterization of PtrEPSP revealed that this isoform possesses a helix-turn-helix motif in the N terminus and can function as a transcriptional repressor that regulates expression of genes in the phenylpropanoid pathway in addition to performing its canonical biosynthesis function in the shikimate pathway. We demonstrated that this isoform can localize in the nucleus and specifically binds to the promoter and represses the expression of a SLEEPER-like transcriptional regulator, which itself specifically binds to the promoter and represses the expression of PtrMYB021 (known as MYB46 in Arabidopsis thaliana), a master regulator of the phenylpropanoid pathway and lignin biosynthesis. Analyses of overexpression and RNAi lines targeting PtrEPSP confirmed the predicted changes in PtrMYB021 expression patterns. These results demonstrate that PtrEPSP in its regulatory form and PtrhAT form a transcriptional hierarchy regulating phenylpropanoid pathway and lignin biosynthesis in Populus.

Quantitative proteome profile of water deficit stress responses in eastern cottonwood (Populus deltoides) leaves

Drought stress is a recurring feature of world climate and the single most important factor influencing agricultural yield worldwide. Plants display highly variable, species-specific responses to drought and these responses are multifaceted, requiring physiological and morphological changes influenced by genetic and molecular mechanisms. Moreover, the reproducibility of water deficit studies is very cumbersome, which significantly impedes research on drought tolerance, because how a plant responds is highly influenced by the timing, duration, and intensity of the water deficit. Despite progress in the identification of drought-related mechanisms in many plants, the molecular basis of drought resistance remains to be fully understood in trees, particularly in poplar species because their wide geographic distribution results in varying tolerances to drought. Herein, we aimed to better understand this complex phenomenon in eastern cottonwood (Populus deltoides) by performing a detailed contrast of the proteome changes between two different water deficit experiments to identify functional intersections and divergences in proteome responses. We investigated plants subjected to cyclic water deficit and compared these responses to plants subjected to prolonged acute water deficit. In total, we identified 108,012 peptide sequences across both experiments that provided insight into the quantitative state of 22,737 Populus gene models and 8,199 functional protein groups in response to drought. Together, these datasets provide the most comprehensive insight into proteome drought responses in poplar to date and a direct proteome comparison between short period dehydration shock and cyclic, post-drought re-watering. Overall, this investigation provides novel insights into drought avoidance mechanisms that are distinct from progressive drought stress. Additionally, we identified proteins that have been associated as drought-relevant in previous studies. Importantly, we highlight the RD26 transcription factor as a gene regulated at both the transcript and protein level, regardless of species and drought condition, and, thus, represents a key, universal drought marker for Populus species.

The Physcomitrella patens chromosome-scale assembly reveals moss genome structure and evolution

The draft genome of the moss model, Physcomitrella patens, comprised approximately 2000 unordered scaffolds. In order to enable analyses of genome structure and evolution we generated a chromosome-scale genome assembly using genetic linkage as well as (end) sequencing of long DNA fragments. We find that 57% of the genome comprises transposable elements (TEs), some of which may be actively transposing during the life cycle. Unlike in flowering plant genomes, gene- and TE-rich regions show an overall even distribution along the chromosomes. However, the chromosomes are mono-centric with peaks of a class of Copia elements potentially coinciding with centromeres. Gene body methylation is evident in 5.7% of the protein-coding genes, typically coinciding with low GC and low expression. Some giant virus insertions are transcriptionally active and might protect gametes from viral infection via siRNA mediated silencing. Structure-based detection methods show that the genome evolved via two rounds of whole genome duplications (WGDs), apparently common in mosses but not in liverworts and hornworts. Several hundred genes are present in colinear regions conserved since the last common ancestor of plants. These syntenic regions are enriched for functions related to plant-specific cell growth and tissue organization. The P. patens genome lacks the TE-rich pericentromeric and gene-rich distal regions typical for most flowering plant genomes. More non-seed plant genomes are needed to unravel how plant genomes evolve, and to understand whether the P. patens genome structure is typical for mosses or bryophytes.

Overexpression of a Domain of Unknown Function 266-containing protein results in high cellulose content, reduced recalcitrance, and enhanced plant growth in the bioenergy crop Populus

BACKGROUND

Domain of Unknown Function 266 (DUF266) is a plant-specific domain. DUF266-containing proteins (DUF266 proteins) have been categorized as 'not classified glycosyltransferases (GTnc)' due to amino acid similarity with GTs. However, little is known about the function of DUF266 proteins.

RESULTS

Phylogenetic analysis revealed that DUF266 proteins are only present in the land plants including moss and lycophyte. We report the functional characterization of one member of DUF266 proteins in Populus, PdDUF266A. PdDUF266A was ubiquitously expressed with high abundance in the xylem. In Populus transgenic plants overexpressing PdDUF266A (OXPdDUF266A), the glucose and cellulose contents were significantly higher, while the lignin content was lower than that in the wild type. Degree of polymerization of cellulose in OXPdDUF266A transgenic plants was also higher, whereas cellulose crystallinity index remained unchanged. Gene expression analysis indicated that cellulose biosynthesis-related genes such as CESA and SUSY were upregulated in mature leaf and xylem of OXPdDUF266A transgenic plants. Moreover, PdDUF266A overexpression resulted in an increase of biomass production. Their glucose contents and biomass phenotypes were further validated via heterologous expression of PdDUF266A in Arabidopsis. Results from saccharification treatment demonstrated that the rate of sugar release was increased by approximately 38% in the OXPdDUF266A transgenic plants.

CONCLUSIONS

These results suggest that the overexpression of PdDUF266A can increase cellulose content, reduce recalcitrance, and enhance biomass production, and that PdDUF266A is a promising target for genetic manipulation for biofuel production.

Overexpression of a Domain of Unknown Function 231-containing protein increases O-xylan acetylation and cellulose biosynthesis in Populus

BACKGROUND

Domain of Unknown Function 231-containing proteins (DUF231) are plant specific and their function is largely unknown. Studies in the model plants Arabidopsis and rice suggested that some DUF231 proteins act in the process of O-acetyl substitution of hemicellulose and esterification of pectin. However, little is known about the function of DUF231 proteins in woody plant species.

RESULTS

This study provides evidence supporting that one member of DUF231 family proteins in the woody perennial plant Populus deltoides (genotype WV94), PdDUF231A, has a role in the acetylation of xylan and affects cellulose biosynthesis. A total of 52 DUF231-containing proteins were identified in the Populus genome. In P. deltoides transgenic lines overexpressing PdDUF231A (OXPdDUF231A), glucose and cellulose contents were increased. Consistent with these results, the transcript levels of cellulose biosynthesis-related genes were increased in the OXPdDUF231A transgenic lines. Furthermore, the relative content of total acetylated xylan was increased in the OXPdDUF231A transgenic lines. Enzymatic saccharification assays revealed that the rate of glucose release increased in OXPdDUF231A transgenic lines. Plant biomass productivity was also increased in OXPdDUF231A transgenic lines.

CONCLUSIONS

These results suggest that PdDUF231A affects cellulose biosynthesis and plays a role in the acetylation of xylan. PdDUF231A is a promising target for genetic modification for biofuel production because biomass productivity and compositional quality can be simultaneously improved through overexpression.

Agronomic performance of Populus deltoides trees engineered for biofuel production

BACKGROUND

One of the major barriers to the development of lignocellulosic feedstocks is the recalcitrance of plant cell walls to deconstruction and saccharification. Recalcitrance can be reduced by targeting genes involved in cell wall biosynthesis, but this can have unintended consequences that compromise the agronomic performance of the trees under field conditions. Here we report the results of a field trial of fourteen distinct transgenic Populus deltoides lines that had previously demonstrated reduced recalcitrance without yield penalties under greenhouse conditions.

RESULTS

Survival and productivity of the trial were excellent in the first year, and there was little evidence for reduced performance of the transgenic lines with modified target gene expression. Surprisingly, the most striking phenotypic effects in this trial were for two empty-vector control lines that had modified bud set and bud flush. This is most likely due to somaclonal variation or insertional mutagenesis. Traits related to yield, crown architecture, herbivory, pathogen response, and frost damage showed few significant differences between target gene transgenics and empty vector controls. However, there were a few interesting exceptions. Lines overexpressing the DUF231 gene, a putative O-acetyltransferase, showed early bud flush and marginally increased height growth. Lines overexpressing the DUF266 gene, a putative glycosyltransferase, had significantly decreased stem internode length and slightly higher volume index. Finally, lines overexpressing the PFD2 gene, a putative member of the prefoldin complex, had a slightly reduced volume index.

CONCLUSIONS

This field trial demonstrates that these cell wall modifications, which decreased cell wall recalcitrance under laboratory conditions, did not seriously compromise first-year performance in the field, despite substantial challenges, including an outbreak of a stem boring insect (Gypsonoma haimbachiana), attack by a leaf rust pathogen (Melampsora spp.), and a late frost event. This bodes well for the potential utility of these lines as advanced biofuels feedstocks.

Two Poplar-Associated Bacterial Isolates Induce Additive Favorable Responses in a Constructed Plant-Microbiome System

The biological function of the plant-microbiome system is the result of contributions from the host plant and microbiome members. The Populus root microbiome is a diverse community that has high abundance of β- and γ-Proteobacteria, both classes which include multiple plant-growth promoting representatives. To understand the contribution of individual microbiome members in a community, we studied the function of a simplified community consisting of Pseudomonas and Burkholderia bacterial strains isolated from Populus hosts and inoculated on axenic Populus cutting in controlled laboratory conditions. Both strains increased lateral root formation and root hair production in Arabidopsis plate assays and are predicted to encode for different functions related to growth and plant growth promotion in Populus hosts. Inoculation individually, with either bacterial isolate, increased root growth relative to uninoculated controls, and while root area was increased in mixed inoculation, the interaction term was insignificant indicating additive effects of root phenotype. Complementary data including photosynthetic efficiency, whole-transcriptome gene expression and GC-MS metabolite expression data in individual and mixed inoculated treatments indicate that the effects of these bacterial strains are unique and additive. These results suggest that the function of a microbiome community may be predicted from the additive functions of the individual members.

Scaling nitrogen and carbon interactions: what are the consequences of biological buffering?

Understanding the consequences of elevated CO₂ (eCO₂; 800 ppm) on terrestrial ecosystems is a central theme in global change biology, but relatively little is known about how altered plant C and N metabolism influences higher levels of biological organization. Here, we investigate the consequences of C and N interactions by genetically modifying the N-assimilation pathway in Arabidopsis and initiating growth chamber and mesocosm competition studies at current CO₂ (cCO₂; 400 ppm) and eCO₂ over multiple generations. Using a suite of ecological, physiological, and molecular genomic tools, we show that a single-gene mutant of a key enzyme (nia2) elicited a highly orchestrated buffering response starting with a fivefold increase in the expression of a gene paralog (nia1) and a 63% increase in the expression of gene network module enriched for N-assimilation genes. The genetic perturbation reduced amino acids, protein, and TCA-cycle intermediate concentrations in the nia2 mutant compared to the wild-type, while eCO₂ mainly increased carbohydrate concentrations. The mutant had reduced net photosynthetic rates due to a 27% decrease in carboxylation capacity and an 18% decrease in electron transport rates. The expression of these buffering mechanisms resulted in a penalty that negatively correlated with fitness and population dynamics yet showed only minor alterations in our estimates of population function, including total per unit area biomass, ground cover, and leaf area index. This study provides insight into the consequences of buffering mechanisms that occur post-genetic perturbations in the N pathway and the associated outcomes these buffering systems have on plant populations relative to eCO₂.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1215-z) contains supplementary material, which is available to authorized users.

A 34K SNP genotyping array for Populus trichocarpa: design, application to the study of natural populations and transferability to other Populus species

Genetic mapping of quantitative traits requires genotypic data for large numbers of markers in many individuals. For such studies, the use of large single nucleotide polymorphism (SNP) genotyping arrays still offers the most cost-effective solution. Herein we report on the design and performance of a SNP genotyping array for Populus trichocarpa (black cottonwood). This genotyping array was designed with SNPs pre-ascertained in 34 wild accessions covering most of the species latitudinal range. We adopted a candidate gene approach to the array design that resulted in the selection of 34 131 SNPs, the majority of which are located in, or within 2 kb of, 3543 candidate genes. A subset of the SNPs on the array (539) was selected based on patterns of variation among the SNP discovery accessions. We show that more than 95% of the loci produce high quality genotypes and that the genotyping error rate for these is likely below 2%. We demonstrate that even among small numbers of samples (n = 10) from local populations over 84% of loci are polymorphic. We also tested the applicability of the array to other species in the genus and found that the number of polymorphic loci decreases rapidly with genetic distance, with the largest numbers detected in other species in section Tacamahaca. Finally, we provide evidence for the utility of the array to address evolutionary questions such as intraspecific studies of genetic differentiation, species assignment and the detection of natural hybrids.

Initial characterization of shade avoidance response suggests functional diversity between Populus phytochrome B genes

Shade avoidance signaling involves perception of incident red/far-red (R/FR) light by phytochromes (PHYs) and modulation of downstream transcriptional networks. Although these responses are well studied in Arabidopsis, little is known about the role of PHYs and the transcriptional responses to shade in the woody perennial Populus. Tissue expression and subcellular localization of Populus PHYs was studied by quantitative real-time PCR (qRT-PCR) and protoplast transient assay. Transgenic lines with altered PHYB1 and/or PHYB2 expression were used in phenotypic assays and transcript profiling with qRT-PCR. RNA-Seq was used to identify transcriptional responses to enriched FR light. All three PHYs were differentially expressed among tissue types and PHYBs were targeted to the nucleus under white light. Populus PHYB1 rescued Arabidopsis phyB mutant phenotypes. Phenotypes of Populus transgenic lines and the expression of candidate shade response genes suggested that PHYB1 and PHYB2 have distinct yet overlapping functions. RNA-Seq analysis indicated that genes associated with cell wall modification and brassinosteroid signaling were induced under enriched FR light in Populus. This study is an initial attempt at deciphering the role of Populus PHYs by evaluating transcriptional reprogramming to enriched FR and demonstrates functional diversity and overlap of the Populus PHYB1 and PHYB2 in regulating shade responses.

Genome resequencing reveals multiscale geographic structure and extensive linkage disequilibrium in the forest tree Populus trichocarpa

• Plant population genomics informs evolutionary biology, breeding, conservation and bioenergy feedstock development. For example, the detection of reliable phenotype-genotype associations and molecular signatures of selection requires a detailed knowledge about genome-wide patterns of allele frequency variation, linkage disequilibrium and recombination. • We resequenced 16 genomes of the model tree Populus trichocarpa and genotyped 120 trees from 10 subpopulations using 29,213 single-nucleotide polymorphisms. • Significant geographic differentiation was present at multiple spatial scales, and range-wide latitudinal allele frequency gradients were strikingly common across the genome. The decay of linkage disequilibrium with physical distance was slower than expected from previous studies in Populus, with r(2) dropping below 0.2 within 3-6 kb. Consistent with this, estimates of recent effective population size from linkage disequilibrium (N(e) ≈ 4000-6000) were remarkably low relative to the large census sizes of P. trichocarpa stands. Fine-scale rates of recombination varied widely across the genome, but were largely predictable on the basis of DNA sequence and methylation features. • Our results suggest that genetic drift has played a significant role in the recent evolutionary history of P. trichocarpa. Most importantly, the extensive linkage disequilibrium detected suggests that genome-wide association studies and genomic selection in undomesticated populations may be more feasible in Populus than previously assumed.

Comparative physiology and transcriptional networks underlying the heat shock response in Populus trichocarpa, Arabidopsis thaliana and Glycine max

The heat shock response continues to be layered with additional complexity as interactions and crosstalk among heat shock proteins (HSPs), the reactive oxygen network and hormonal signalling are discovered. However, comparative analyses exploring variation in each of these processes among species remain relatively unexplored. In controlled environment experiments, photosynthetic response curves were conducted from 22 to 42 °C and indicated that temperature optimum of light-saturated photosynthesis was greater for Glycine max relative to Arabidopsis thaliana or Populus trichocarpa. Transcript profiles were taken at defined states along the temperature response curves, and inferred pathway analysis revealed species-specific variation in the abiotic stress and the minor carbohydrate raffinose/galactinol pathways. A weighted gene co-expression network approach was used to group individual genes into network modules linking biochemical measures of the antioxidant system to leaf-level photosynthesis among P. trichocarpa, G. max and A. thaliana. Network-enabled results revealed an expansion in the G. max HSP17 protein family and divergence in the regulation of the antioxidant and heat shock modules relative to P. trichocarpa and A. thaliana. These results indicate that although the heat shock response is highly conserved, there is considerable species-specific variation in its regulation.

Discovery and annotation of small proteins using genomics, proteomics, and computational approaches

Small proteins (10-200 amino acids [aa] in length) encoded by short open reading frames (sORF) play important regulatory roles in various biological processes, including tumor progression, stress response, flowering, and hormone signaling. However, ab initio discovery of small proteins has been relatively overlooked. Recent advances in deep transcriptome sequencing make it possible to efficiently identify sORFs at the genome level. In this study, we obtained ~2.6 million expressed sequence tag (EST) reads from Populus deltoides leaf transcriptome and reconstructed full-length transcripts from the EST sequences. We identified an initial set of 12,852 sORFs encoding proteins of 10-200 aa in length. Three computational approaches were then used to enrich for bona fide protein-coding sORFs from the initial sORF set: (1) coding-potential prediction, (2) evolutionary conservation between P. deltoides and other plant species, and (3) gene family clustering within P. deltoides. As a result, a high-confidence sORF candidate set containing 1469 genes was obtained. Analysis of the protein domains, non-protein-coding RNA motifs, sequence length distribution, and protein mass spectrometry data supported this high-confidence sORF set. In the high-confidence sORF candidate set, known protein domains were identified in 1282 genes (higher-confidence sORF candidate set), out of which 611 genes, designated as highest-confidence candidate sORF set, were supported by proteomics data. Of the 611 highest-confidence candidate sORF genes, 56 were new to the current Populus genome annotation. This study not only demonstrates that there are potential sORF candidates to be annotated in sequenced genomes, but also presents an efficient strategy for discovery of sORFs in species with no genome annotation yet available.

Microsatellite primer resource for Populus developed from the mapped sequence scaffolds of the Nisqually-1 genome

In this study, 148 428 simple sequence repeat (SSR) primer pairs were designed from the unambiguously mapped sequence scaffolds of the Nisqually-1 genome. The physical position of the priming sites were identified along each of the 19 Populus chromosomes, and it was specified whether the priming sequences belong to intronic, intergenic, exonic or UTR regions. A subset of 150 SSR loci were amplified and a high amplification success rate (72%) was obtained in P. tremuloides, which belongs to a divergent subgenus of Populus relative to Nisqually-1. PCR reactions showed that the amplification success rate of exonic primer pairs was much higher than that of the intronic/intergenic primer pairs. Applying ANOVA and regression analyses to the flanking sequences of microsatellites, the repeat lengths, the GC contents of the repeats, the repeat motif numbers, the repeat motif length and the base composition of the repeat motif, it was determined that only the base composition of the repeat motif and the repeat motif length significantly affect the microsatellite variability in P. tremuloides samples. The SSR primer resource developed in this study provides a database for selecting highly transferable SSR markers with known physical position in the Populus genome and provides a comprehensive genetic tool to extend the genome sequence of Nisqually-1 to genetic studies in different Populus species.

The F-box gene family is expanded in herbaceous annual plants relative to woody perennial plants

F-box proteins are generally responsible for substrate recognition in the Skp1-Cullin-F-box complexes that are involved in protein degradation via the ubiquitin-26S proteasome pathway. In plants, F-box genes influence a variety of biological processes, such as leaf senescence, branching, self-incompatibility, and responses to biotic and abiotic stresses. The number of F-box genes in Populus (Populus trichocarpa; approximately 320) is less than half that found in Arabidopsis (Arabidopsis thaliana; approximately 660) or Oryza (Oryza sativa; approximately 680), even though the total number of genes in Populus is equivalent to that in Oryza and 1.5 times that in Arabidopsis. We performed comparative genomics analysis between the woody perennial plant Populus and the herbaceous annual plants Arabidopsis and Oryza in order to explicate the functional implications of this large gene family. Our analyses reveal interspecific differences in genomic distribution, orthologous relationship, intron evolution, protein domain structure, and gene expression. The set of F-box genes shared by these species appear to be involved in core biological processes essential for plant growth and development; lineage-specific differences primarily occurred because of an expansion of the F-box genes via tandem duplications in Arabidopsis and Oryza. The number of F-box genes in the newly sequenced woody species Vitis (Vitis vinifera; 156) and Carica (Carica papaya; 139) is similar to that in Populus, supporting the hypothesis that the F-box gene family is expanded in herbaceous annual plants relative to woody perennial plants. This study provides insights into the relationship between the structure and composition of the F-box gene family in herbaceous and woody species and their associated developmental and physiological features.

Genome structure and emerging evidence of an incipient sex chromosome in Populus

The genus Populus consists of dioecious woody species with largely unknown genetic mechanisms for gender determination. We have discovered genetic and genomic features in the peritelomeric region of chromosome XIX that suggest this region of the Populus genome is in the process of developing characteristics of a sex chromosome. We have identified a gender-associated locus that consistently maps to this region. Furthermore, comparison of genetic maps across multiple Populus families reveals consistently distorted segregation within this region. We have intensively characterized this region using an F(1) interspecific cross involving the female genotype that was used for genome sequencing. This region shows suppressed recombination and high divergence between the alternate haplotypes, as revealed by dense map-based genome assembly using microsatellite markers. The suppressed recombination, distorted segregation, and haplotype divergence were observed only for the maternal parent in this cross. Furthermore, the progeny of this cross showed a strongly male-biased sex ratio, in agreement with Haldane's rule that postulates that the heterogametic sex is more likely to be absent, rare, or sterile in interspecific crosses. Together, these results support the role of chromosome XIX in sex determination and suggest that sex determination in Populus occurs through a ZW system in which the female is the heterogametic gender.

Connecting genes, coexpression modules, and molecular signatures to environmental stress phenotypes in plants

Background

One of the eminent opportunities afforded by modern genomic technologies is the potential to provide a mechanistic understanding of the processes by which genetic change translates to phenotypic variation and the resultant appearance of distinct physiological traits. Indeed much progress has been made in this area, particularly in biomedicine where functional genomic information can be used to determine the physiological state (e.g., diagnosis) and predict phenotypic outcome (e.g., patient survival). Ecology currently lacks an analogous approach where genomic information can be used to diagnose the presence of a given physiological state (e.g., stress response) and then predict likely phenotypic outcomes (e.g., stress duration and tolerance, fitness).

Results

Here, we demonstrate that a compendium of genomic signatures can be used to classify the plant abiotic stress phenotype in Arabidopsis according to the architecture of the transcriptome, and then be linked with gene coexpression network analysis to determine the underlying genes governing the phenotypic response. Using this approach, we confirm the existence of known stress responsive pathways and marker genes, report a common abiotic stress responsive transcriptome and relate phenotypic classification to stress duration.

Conclusion

Linking genomic signatures to gene coexpression analysis provides a unique method of relating an observed plant phenotype to changes in gene expression that underlie that phenotype. Such information is critical to current and future investigations in plant biology and, in particular, to evolutionary ecology, where a mechanistic understanding of adaptive physiological responses to abiotic stress can provide researchers with a tool of great predictive value in understanding species and population level adaptation to climate change.

Development of AFLP and RAPD markers linked to a locus associated with twisted growth in corkscrew willow (Salix matsudana 'Tortuosa')

Salix matsudana Koidz. cultivar 'Tortuosa' (corkscrew willow) is characterized by extensive stem bending and curling of leaves. To investigate the genetic basis of this trait, controlled crosses were made between a corkscrew female (S. matsudana 'Tortuosa') and a straight-stemmed, wild-type male (Salix alba L. Clone 99010). Seventy-seven seedlings from this family (ID 99270) were grown in the field for phenotypic observation. Among the progeny, 39 had straight stems and leaves and 38 had bent stems and curled leaves, suggesting that a dominant allele at a single locus controls this phenotype. As a first step in characterizing the locus, we searched for amplified fragment length polymorphism (AFLP) and randomly amplified polymorphic DNA (RAPD) markers linked to the tortuosa allele using bulked segregant analysis. Samples of DNA from 10 corkscrew individuals were combined to produce a corkscrew pool, and DNA from 10 straight progeny was combined to make a wild-type pool. Sixty-four AFLP primer combinations and 640 RAPD primers were screened to identify marker bands amplified from the corkscrew parent and progeny pool, but not from the wild-type parent or progeny pool. An AFLP marker and a RAPD marker linked to and flanking the tortuosa locus were placed on a preliminary linkage map constructed based on segregation among the 77 progeny. Sectioning and analysis of shoot tips revealed that the corkscrew phenotype is associated with vascular cell collapse, smaller cell size in regions near the cambium and less developed phloem fibers than in wild-type progeny. Identification of a gene associated with this trait could lead to greater understanding of the control of normal stem development in woody plants.

Genetic and physical mapping of Melampsora rust resistance genes in Populus and characterization of linkage disequilibrium and flanking genomic sequence

•  In an attempt to elucidate the molecular mechanisms of Melampsora rust resistance in Populus trichocarpa, we have mapped two resistance loci, MXC3 and MER, and intensively characterized the flanking genomic sequence for the MXC3 locus and the level of linkage disequilibrium (LD) in natural populations. •  We used an interspecific backcross pedigree and a genetic map that was highly saturated with AFLP and SSR markers, and assembled shotgun-sequence data in the region containing markers linked to MXC3. •  The two loci were mapped to different linkage groups. Linkage disequilibrium for MXC3 was confined to two closely linked regions spanning 34 and 16 kb, respectively. The MXC3 region also contained six disease-resistance candidate genes. •  The MER and MXC3 loci are clearly distinct, and may have different mechanisms of resistance, as different classes of putative resistance genes were present near each locus. The suppressed recombination previously observed in the MXC3 region was possibly caused by extensive hemizygous rearrangements confined to the original parent tree. The relatively low observed LD may facilitate association studies using candidate genes for rust resistance, but will probably inhibit marker-aided selection.

Large-scale heterospecific segregation distortion in Populus revealed by a dense genetic map

We report the most complete genetic map to have been constructed for the genus Populus. This map includes 544 markers mapped onto 19 linkage groups, equivalent to the Populus chromosome number, with all markers displaying internally consistent linkage patterns. We estimate the genome length to be between 2,300 and 2,500 cM, based both on the observed number of crossovers in the maternal haplotypes, as well as the total observed map length. Genome coverage was estimated to be greater than 99.9% at 20 cM per marker. We did not detect obvious recombination repression in the maternal tree (a hybrid of Populus trichocarpa Hooker x P. deltoides Marsh.) compared to the paternal tree (pure P. deltoides). Finally, most markers exhibiting segregation distortion were derived from the donor parent in this backcross, and generally occurred in large contiguous blocks on two linkage groups. We hypothesize that divergent selection has occurred on chromosomal scales among the parental species used to create this pedigree, and explore the evolutionary implications of this observation. This genetic linkage map provides the most comprehensive view of the Populus genome reported to date and will prove invaluable for future inquiries into the structural and functional genomics, evolutionary biology, and genetic improvement of this ecologically important model species.

The development of two flanking SCAR markers linked to a sex determination locus in Salix viminalis L

Most studies of sex determination systems in plants involve dioecious annuals that have known sex chromosomes. Despite the absence of such structures in the majority of dioecious plants, gender seems to be under relatively strict genetic control in some species. Genetic markers linked to a female sex-determination locus in Salix viminalis L. have been discovered through bulked segregant analysis of three full-sib families using approximately 1,000 arbitrary primers. Two RAPD markers that were present in the common female parent as well as in predominantly female progeny of these families were subsequently sequenced and converted to sequence characterized amplified region (SCAR) markers. The two SCAR markers are correlated with gender in the three full-sib families and are present in 96.4% of the female progeny and 2.2% of the males, providing evidence of linkage to a putative female-specific locus associated with gender determination in S. viminalis. Estimates of recombination suggest that the two markers are flanking a putative sex determination locus, SDL-II, in certain families of S. viminalis.

Demethylation of satellite I DNA during senescence of bovine adrenocortical cells in culture

Over the finite proliferative life span of cultured bovine adrenocortical cells, satellite I DNA shows a progressive and extensive loss of methylation at CCGG sites. This was shown by Southern blotting after digestion with the methylation-sensitive enzyme HpaII alone, which provides a sensitive indicator of methylation loss, or digestion with the combination of EcoRI and HpaII, which provides a quantitative indication of loss of methylation. Bovine tissues, including adrenal cortex, all showed a much higher level of satellite methylation than cultured adrenocortical cells. After adrenocortical cells are placed in culture, some demethylation of satellite I is seen as early as 10 population doublings. By 80 population doublings, loss of satellite DNA methylation is extensive. The loss does not appear to prevent continued cell division, since an extended life span clone of bovine adrenocortical cells transfected with SV40 T antigen showed a similar pattern of extensive demethylation. Satellite demethylation has been reported in aging in vivo and the present cell culture system may provide an in vitro model for this form of genetic instability.

Translation initiation in Drosophila melanogaster is reduced by mutations upstream of the AUG initiator codon

The importance to in vivo translation of sequences immediately upstream of the Drosophila alcohol dehydrogenase (Adh) start codon was examined at two developmental stages. Mutations were introduced into the Adh gene in vitro, and the mutant gene was inserted into the genome via germ line transformation. An A-to-T substitution at the -3 position did not affect relative translation rates of the ADH protein at the second-instar larval stage but resulted in a 2.4-fold drop in translation of ADH at the adult stage. A second mutant gene, containing five mutations in the region -1 to -9, was designed to completely block translation initiation. However, transformant lines bearing these mutations still exhibit detectable ADH, albeit at substantially reduced levels. The average fold reduction at the second-instar larval stage was 5.9, while at the adult stage a 12.5-fold reduction was observed.

Characterization of cell surface polypeptides of unfertilized, fertilized, and protease-treated zona-free mouse eggs

The polypeptide composition of unfertilized, fertilized, and protease-treated zona-free mouse eggs was evaluated in this study. Zona-free eggs were radioiodinated by an Iodogen-catalyzed reaction. Light microscopic autoradiography of egg sections revealed that labeling was restricted to the cell surface. Labeled eggs were solubilized, and cell surface polypeptides were identified by one-dimensional SDS polyacrylamide gel electrophoresis and autoradiography. The unfertilized egg demonstrated 8-10 peptides that incorporated 125I, with major bands observed at approximately 145-150, 94, and 23 kilodaltons (kD). Zona-free eggs fertilized in vitro and then radiolabeled demonstrated several new bands in comparison to unfertilized eggs, with a major band appearing at approximately 36 kD. Treatment of radiolabeled unfertilized eggs with either trypsin or chymotrypsin (1 mg/ml for 5-20 min) caused enzyme-specific modifications in labeled polypeptides. Trypsin (T) treatment resulted in time-dependant modification of the three major peptides at 145-150, 94, and 23 kD. Chymotrypsin (CT) treatment, in contrast, was associated with loss or modification of the 94 kD band, with no apparent effect on either the 145-150 or 23 kD band. Taken together with previous data indicating that T or CT egg treatment interferes with sperm-egg attachment and fusion (Boldt et al.: Biol Reprod 39:19-27, 1988), these results suggest a possible role for the 94 kD protein in sperm-egg interaction.

Carbohydrate involvement in sperm-egg fusion in mice

The potential involvement of cell-surface carbohydrates in sperm-egg fusion in mice was evaluated in this study. Zona-free mouse eggs were inseminated in the presence of a variety of simple saccharides to determine if certain sugars would act as competitive inhibitors of sperm-egg fusion. Of the sugars tested, L-fucose, galactose, and N-acetylglucosamine caused the greatest inhibition of sperm penetration levels relative to controls. A number of complex saccharides or glycoproteins with differing carbohydrate structures, including fucoidan, ascophyllan, ovomucoid, ovalbumin, fetuin, asialofetuin, and chondroitin sulfate, were also tested as competitive inhibitors of fusion. Only the L-fucose containing saccharides fucoidan and ascophyllan caused significant inhibition of fusion at concentrations of 0.05-1.0 mg/ml and 0.1-5.0 mg/ml, respectively. None of the other compounds tested had any inhibitory effect on fusion when tested at concentrations up to 5 mg/ml. The effect of the inhibitory saccharides was not due to the presence of residual zona material on the surface of the zona-free eggs, since zona-free eggs did not bind an 125I-labeled antibody directed against the ZP3 protein of the mouse zona pellucida. Pretreatment of either sperm or eggs with fucoidan (1 mg/ml) for 60 min prior to insemination caused only small decreases in sperm penetration levels, indicating that fucoidan exerted its major inhibitory effect on fusion only when present during insemination. Treatment of sperm, but not zona-free eggs, with fucosidase prior to insemination caused significant reductions in sperm penetration levels. Other glycosidic enzymes, including glucosidase, galactosidase, and N-acetylglucosaminidase, had no inhibitory effect on the sperm. These data suggest that an L-fucose component of the sperm surface is involved in sperm-egg fusion in the mouse.